Correcting method of eddy covariance fluxes over non-flat surfaces and its application in ChinaFLUX

Although Eddy Covariance (EC) technique is one of the best methods for estimating the energy and mass exchanges between underlying surface and atmosphere in micrometeorology, errors and uncertainties still exist without necessary corrections. In this paper, we will focus on the effect of coordinate system on the eddy fluxes. Based on the data observed over four sites (one farmland site, one grassland site and two forest sites), the effects of three coordinate system transforming methods (Double Rotation-DR, Triple Rotation-TR and Planar Fit-PF)on the turbulent fluxes are analyzed. It shows that (i) the corrected fluxes are more or less than the uncorrected fluxes, which is related mainly to the sloping degree of surface, wind speed and wind direction; and (ii) pitch angle has a sinusoidal dependence on wind direction, especially in the regular sloping terrain; and (iii) PF method is something like the simplification of TR or DR,and there are not obvious distinctions in correction in sloping grassland and flat farmland, but PF method is not suitable for uneven and irregular forest sites.     

Establishment of apparent quantum yield and maximum ecosystem assimilation on Tibetan Plateau alpine meadow ecosystem

The alpine meadow is widely distributed on the Tibetan Plateau with an area of about 1.2×106kn2. Damxung County, located in the hinterland of the Tibetan Plateau, is the place covered with this typical vegetation. An open-path eddy covariance system was set up in Damxung rangeland station to measure the carbon flux of alpine meadow from July to October,2003. The continuous carbon flux data were used to analyze the relationship between net ecosystem carbon dioxide exchange (NEE) and photosynthetically active radiation (PAR), as well as the seasonal patterns of apparent quantum yield (α) and maximum ecosystem assimilation (Pmax).Results showed that the daytime NEE fitted fairly well with the PAR in a rectangular hyperbola function, with α declining in the order of peak growth period (0.0244 μmolCO2 · μmol-1pAR) ＞early growth period ＞ seed maturing period ＞ withering period (0.0098 μmolCO2 · μmol-1pAR).The Pmax did not change greatly during the first three periods, with an average of 0.433mgCO2· m-2· s-1, i.e. 9.829 μmolCO2· m-2· s-1. However, during the withering period, Pmax was only 0.35 mgCO2 · m-2 · s-1, i.e. 7.945 μmolCO2 · m-2 · s-1. Compared with other grassland ecosystems, the α of the Tibetan Plateau alpine meadow ecosystem was much lower.     

Large Eddy Simulation on Interaction Between in-Line and Cross-Flow Oscillation of A Circular Cylinder

This paper discusses numerical results from three-dimensional large eddy simulations of an oscillating cylinder under prescribed movements in uniform flow. Six cases,namely pure in-line,pure cross-flow and two groups of 'Figure of Eight' oscillation patterns are under investigation at Reynolds number Re=24000. The 'Figure of Eight' pattern in each group is with identical shape but opposite orbital directions. The numerical results on hydrodynamic forces,higher order force components,and vortex shedding mode...     

南海中尺度涡强度的季节和年际变化分析

为了研究南海中尺度涡强度的季节和年际变化规律,利用Matlab提取50 a(1958~2007年)简单海洋资料同化(Simple Ocean Data Assimilation,SODA)月平均数据集中流场和海表面高度场数据,应用一个涡旋自动探测算法对南海中尺度涡初始生成位置进行分析,并分析了海表面高度异常均方根值的季节变化和年际变化。结果表明:50 a里南海中尺度涡主要分布在吕宋岛西北海域、吕宋岛西南海域和越南以东广大海域,秋、冬季中尺度涡能量较高,春季中尺度涡最弱,中尺度涡强度高值区年际变化明显。从季节变化上看,海面高度异常均方根春、夏季最小,秋冬季最大;从年际变化上看,与同时期Nino3指数有显著负相关,周期大约为3 a。     

北极薄冰-煎饼冰区域中不同重力波驱动形成海流的敏感性分析

A coupled ocean-ice-wave model is used to study ice-edge jet and eddy genesis during surface gravity wave dissipation in a frazil-pancake ice zone. With observational data from the Beaufort Sea, possible wave dissipation processes are evaluated using sensitivity experiments. As wave energy dissipated, energy was transferred into ice floe through radiation stress. Later, energy was in turn transferred into current through ocean-ice interfacial stress. Since most of the wave energy is dissipated at the ice edge, ice-edge jets, which contained strong horizontal shear, appeared both in the ice zone and the ocean. Meanwhile, the wave propagation direction determines the velocity partition in the along-ice-edge and cross-ice-edge directions, which in turn determines the strength of the along-ice-edge jet and cross-ice-edge velocity. The momentum applied in the along-ice-edge(cross-ice-edge)direction increased(decreased) with larger incident angle, which is favorable condition for producing stronger mesoscale eddies, vice versa. The dissipation rate increases(decreases) with larger(smaller) wavenumber, which enhances(reduces) the jet strength and the strength of the mesoscale eddy. The strong along-ice-edge jet may extend to a deep layer(> 200 m). If the water depth is too shallow(e.g., 80 m), the jet may be largely dampened by bottom drag, and no visible mesoscale eddies are found. The results suggest that the bathymetry and incident wavenumber(magnitude and propagation direction) are important for wave-driven current and mesoscale eddy genesis.     

巴芬湾西部和格陵兰岛东部大气经向热量输送协同变化对夏季北极海冰的影响

利用欧洲中心气候再分析资料和美国国家冰雪数据中心北极海冰面积资料,分析了夏季北极海冰面积与前期大气经向热量输送年际变化的联系。结果表明：6月北半球中高纬大气的经向热量输送以瞬变热量形式为主,其中巴芬湾西部（B区）和格陵兰岛东部（G区）是瞬变热量向极区传输的两个通道,二者之间存在反位相的协同变化,且这种协同变化与夏季北极海冰面积变化密切相关。可能的机制为：6月,AD、AO和NAO三种北极大气环流型能够引起巴芬湾西部和格陵兰岛东部瞬变热量输送的协同变化,这种协同变化通过涡旋动力作用激发夏季极区大气表现为AD异常,同时影响途经区域的气温,从而通过热动力作用影响夏季北极海冰。将向极区输送的热量称为暖输送,从极区输出的热量为冷输送,则上述两个区域的瞬变热量协同输送可分为三种情况：B暖G冷、B冷G暖、B和G均冷,而B和G均暖的情况十分罕见。当B区向极区输入、G区输出热量时,有利于太平洋扇区和喀拉海的海冰偏少;当G区输入、B区输出热量时,利于喀拉海和拉普捷夫海海冰偏少;当B区和G区均输出热量时,利于波佛特海南部、喀拉海和拉普捷夫海海冰偏多,反之则相反。     

Effects of cold eddy on phytoplankton production and assemblages in Luzon strait bordering the South China Sea

The biochemical effects of a cold-core eddy that was shed from the Kuroshio Current at the Luzon Strait bordering the South China Sea (SCS) were studied in late spring, a relatively unproductive season in the SCS. The extent of the eddy was determined by time-series images of SeaWiFS ocean color, AVHRR sea surface temperature, and TOPEX/Jason-1 sea surface height anomaly. Nutrient budgets, nitrate-based new production, primary production, and phytoplankton assemblages were compared between the eddy and its surrounding Kuroshio and SCS waters. The enhanced productivity in the eddy was comparable to wintertime productivity in the SCS basin, which is supported by upwelled subsurface nitrate under the prevailing Northeastern Monsoon. There were more Synechococcus, pico-eucaryotes, and diatoms, but less Trichodesmium in the surface water inside the eddy than outside. Prochlorococcus and Richelia intracellularis showed no spatial differences. Water column-integrated primary production (IPP) inside the eddy was 2–3 times that outside the eddy in the SCS (1.09 vs. 0.59 g C m⁻²d⁻¹), as was nitrate-based new production (INP) (0.67 vs. 0.25 g C m⁻²d⁻¹). INP in the eddy was 6 times that in the Kuroshio (0.12 g C m⁻²d⁻¹). IPP and INP in the eddy were higher than the maximum production values ever measured in the SCS basin. Surface chlorophyll a concentration (0.40 mg m⁻³) in the eddy equaled the maximum concentration registered for the SCS basin and was higher than the wintertime average (0.29 ± 0.04 mg m⁻³). INP was 3.5 times as great and IPP was doubled in the eddy compared to the wintertime SCS basin. As cold core eddies form intermittently all year round as the Kuroshio invades the SCS, their effects on phytoplankton productivity and assemblages are likely to have important influences on the biogeochemical cycle of the region.     

Numerical modeling investigation on turbulent oscillatory flow over a plane rough bed composed by randomly arrayed particles

A three-dimensional numerical model is established to simulate the turbulent oscillatory boundary layer over a fixed and rough bed composed by randomly arrayed solid spheres based on the lattice Boltzmann method and the large eddy simulation model.The equivalent roughness height,the location of the theoretical bed and the time variation of the friction velocity are investigated using the log-fit method.The time series of turbulent intensity and Reynolds stress are also investigated.The equivalent roughness height of cases with Reynolds numbers of 1×10~4–6×10~4 is approximately 2.81 d(grain size).The time variation of the friction velocity in an oscillatory cycle exhibits sinusoidal-like behavior.The friction factor depends on the relative roughness in the rough turbulent regime,and the pattern of solid particles arrayed as the rough bed in the numerical simulations has no obvious effect on the friction factor.     

Reconstruction of eddies by assimilating satellite altimeter data into Princeton Ocean Model

An optimal interpolation assimilation model for satellite altimetry data is developed based on Princeton Ocean Model (POM), which is applied in a quasi-global domain, by the method of isotropic correlation between sea level anomaly (SLA) and sea temperature anomaly. The perfor- mance of this assimilation model is validated by the modeled results of SLA and the current patterns. Comparisons between modeling and satellite data show that both the magnitudes and distribution patterns of the simulated SLA are improved by assimilation. The most signiˉcant improvement is that meso-scale systems, e.g., eddies, are well reconstructed. The evolution of an eddy located in the northwest Paciˉc Ocean is traced by using the assimilation model. Model results show that during three months the eddy migrated southwestward for about 6 degrees before merging into the Kuroshio. The three dimensional structure of this eddy on 12 August 2001 is further analyzed. The strength of this warm, cyclonic eddy decreases with the increase of depth. The eddy shows di?erent horizontal patterns at di?erent layers, and the SLA and temperature ˉelds agree with each other well. This study suggests that this kind of data assimilation is economic and reliable for eddy reconstruction, and can be used as a promising technique in further studies of ocean eddies as well as other ˉne circulation structures.     

A critical analysis of processes governing the nutrient profiles in the ocean

We have theoretically considered the problem of interpretation of nutrient profiles in the upper ocean (100–2000m). We compare the experimental depth profiles of nonconservative tracers, both stable and radioactive, with solutions of one-dimensional steady state transport equations of increasing complexity including situations not encountered in the real oceans. Apart from gaining insight into the nutrient transport processes, this analysis is useful in offering a way to obtain operational estimates of depth dependent/independent eddy diffusivity and dissolution fluxes in the ocean. These parameters are essential for estimating new production, total production and burial of carbon in the sediments.