Air-sea carbon-dioxide flux estimated by eddy covariance method from a buoy observation

Precise measurements of the CO2 gas transfer across the air-sea interface provide a better under- standing of the global carbon cycle. The air-sea CO2 fluxes are obtained by the eddy covariance method and the bulk method from a buoy observation in the northern Huanghai sea. The effects of buoy motion on flux calculated by the eddy covariance method are demonstrated. The research shows that a motion correction can improve the correlation coefficient between the C02 fluxes esti- mated from two different levels. Without the CO2-H20 cross-correlation correction which is termed as PKT correction, the air-sea CO2 fluxes estimated by eddy covariance method using the motion corrected data are nearly an order of magnitude larger than those estimated by the bulk method. After the CO2-H20 cross-correlation correction, some eddy covariance CO2 fluxes indeed become closer to the bulk CO2 flux, whereas some are overcorrected which are in response to small water vapor flux.

Multi-level CO_(2) fluxes over Beijing megacity with the eddy covariance method

本文基于北京325米气象塔在47,140,和280米三层高度的5年涡动相关观测资料,研究了城市下垫面与大气间的CO_(2)交换过程.由于北京市2011年开始实行工作日汽车尾号限行,140米高度CO_(2)通量的年增长率由2008-2010年的7.8%降低到2010-2012年的2.3%.140米高度通量源区内植被比例最小且人口密度最大,因此140米高度的5年平均CO_(2)通量年总量(6.41 kg C m^(−2) yr^(−1)(大于47米)5.78 kg C m^(−2) yr^(−1)(和280米)3.99 kg C m^(−2) yr^(−1)).在年尺度上,北京汽车总保有量和总人口是最重要的CO_(2)通量控制因子.CO_(2)通量随风向的变化主要与风向对应的通量源区内下垫面土地利用方式有关.三层高度的夏季CO_(2)通量均与道路的比例呈正相关关系.47,140,和280米的决定系数分别为0.69,0.57,和0.54(P<0.05).植被比例的下降,会导致CO_(2)年总量上升,两者存在近似于指数的关系.城市人口密度的上升会引起CO_(2)年总量上升.

Eddy Covariance Tilt Corrections over a Coastal Mountain Area in South-east China:Significance for Near-Surface Turbulence Characteristics

Turbulence characteristics of an atmospheric surface layer over a coastal mountain area were investigated under different coordinate frames. Performances of three methods of coordinate rotation： double rotation （DR）, triple rotation （TR）, and classic planar-fit rotation （PF） were examined in terms of correction of eddy covariance flux. Using the commonly used DR and TR methods, unreasonable rotation angles are encountered at low wind speeds and cause significant run-to-run errors of some turbulence characteristics. The PF method rotates the coordinate system to an ensemble-averaged plane, and shows large tilt error due to an inaccurate fit plane over variable terrain slopes. In this paper, we propose another coordinate rotation scheme. The observational data were separated into two groups according to wind direction. The PF method was adapted to find an ensemble-averaged streamline plane for each group of hourly runs with wind speed exceeding 1.0 m s-1. Then, the coordinate systems were rotated to their respective best- fit planes for all available hourly observations. We call this the PF10 method. The implications of tilt corrections for the turbulence characteristics are discussed with a focus on integral turbulence characteristics, the spectra of wind-velocity components, and sensible heat and momentum fluxes under various atmospheric stabilities. Our results show that the adapted application of PF provides greatly improved estimates of integral turbulence characteristics in complex terrain and maintains data quality. The comparisons of the sensible heat fluxes for four coordinate rotation methods to fluxes before correction indicate that the PF10 scheme is the best to preserve consistency between fluxes.

Eddy covariance measurements of water vapor and energy flux over a lake in the Badain Jaran Desert,China

Exploring the surface energy exchange between atmosphere and water bodies is essential to gain a quantitative understanding of regional climate change, especially for the lakes in the desert. In this study, measurements of energy flux and water vapor were performed over a lake in the Badain Jaran Desert, China from March 2012 to March 2013. The studied lake had about a 2-month frozen period （December and January） and a 10-month open-water period （February-November）. Latent heat flux （LE） and sensible heat flux （Hs） acquired using the eddy covariance technique were argued by measurements of long＇wave and shortwave radiation. Both fluxes of longwave and shortwave radiation showed seasonal dynamics and daily fluctuations during the study period. The reflected solar radiation was much higher in winter than in other seasons. LE exhibited diurnal and seasonal variations. On a daily scale, LE was low in the morning and peaked in the afternoon. From spring （April） to winter （January）, the diurnal amplitude of LE decreased slowly. LE was the dominant heat flux throughout the year and consumed most of the energy from the lake. Generally speaking, LE was mostly affected by changes in the ambient wind speed, while Hs was primarily affected by the product of water-air temperature difference and wind speed. The diurnal LE and Hs were negatively correlated in the open-water period. The variations in Hs and LE over the lake were differed from those on the nearby land surface. The mean evaporation rate on the lake was about 4.0 mm/d over the entire year, and the cumulative annual evaporation rate was 1445 mm/a. The cumulative annual evaporation was 10 times larger than the cumulative annual precipitation. Furthermore, the average evaporation rates over the frozen period and open-water period were approximately 0.6 and 5.0 mm/d, respectively. These results can be used to analyze the water balance and quantify the source of lake water in the Badain Jaran Desert.

A large eddy simulation of flows around an underwater vehicle model using an immersed boundary method

A large eddy simulation （LES） of the flows around an underwater vehicle model at intermediate Reynolds numbers is performed. The underwater vehicle model is taken as the DARPA SUBOFF with full appendages, where the Reynolds number based on the hull length is 1.0x 105, An immersed boundary method based on the moving-least-squares reconstruction is used to handle the complex geometric boundaries. The adaptive mesh refinement is utilized to resolve the flows near the hull, The parallel scalabilities of the flow solver are tested on meshes with the number of cells varying from 50 million to 3.2 billion, The parallel solver reaches nearly linear scalability for the flows around the underwater vehicle model, The present simulation captures the essential features of the vortex structures near the hull and in the wake, Both of the time-averaged pressure coefficients and srreamwise velocity profiles obtained from the LES are consistent with the characteristics of the flows pass an appended axisymmetric body. The code efficiency and its correct predictions on flow features allow us to perform the full-scale simulations on tens of thousands of cores with billions of grid points for higher-Reynolds-number flows around the underwater vehicles.

A Re-examination of Density Effects in Eddy Covariance Measurements of CO_2 Fluxes

Corrections of density effects resulting from air-parcel expansion/compression are important in interpreting eddy covariance fluxes of water vapor and CO2 when open-path systems are used. To account for these effects, mean vertical velocity and perturbation of the density of dry air are two critical parameters in treating those physical processes responsible for density variations. Based on various underlying assumptions, different studies have obtained different formulas for the mean vertical velocity and perturbation of the density of dry air, leading to a number of approaches to correct density effects. In this study, we re-examine physical processes related to different assumptions that are made to formulate the density effects. Specifically, we re-examine the assumptions of a zero dry air flux and a zero moist air flux in the surface layer, used for treating density variations, and their implications for correcting density effects. It is found that physical processes in relation to the assumption of a zero dry air flux account for the influence of dry air expansion/compression on density variations. Meanwhile, physical processes in relation to the assumption of a zero moist air flux account for the influence of moist air expansion/compression on density variations. In this study, we also re-examine mixing ratio issues. Our results indicate that the assumption of a zero dry air flux favors the use of the mixing ratio relative to dry air, while the assumption of a zero moist air flux favors the use of the mixing ratio relative to the total moist air. Additionally, we compare different formula for the mean vertical velocity, generated by air-parcel expansion/compression, and for density effect corrections using eddy covariance data measured over three boreal ecosystems.

Numerical Simulation of Sloshing Using the MPS-FSI Method with Large Eddy Simulation

A numerical model has been developed to study sloshing of turbulent flow in a tank with elastic baffles. The Moving-Particle Semi-implicit method(MPS) is a kind of meshless Lagrangian calculation method. The large eddy simulation(LES) approach is employed to model the turbulence by using the Smagorinsky Sub-Particle Scale(SPS)closure model. This paper uses MPS-FSI method with LES to simulate the interaction between free surface flow and a thin elastic baffle in sloshing. Then, the numerical model is validated, and the numerical solution has good agreement with experimental data for sloshing in a tank with elastic baffles. Furthermore, under external excitations,the MPS is applied to viscous laminar flow and turbulent flow, with both the deformation of elastic baffles and the wave height of the free surface are compared with each other. Besides, the impact pressure with/without baffles and wave height of free surface are investigated and discussed in detail. Finally, preliminary simulations are carried out in the damage problem of elastic baffles, taking the advantage of the MPS-FSI method in computations of the fluid–structure interaction with large deformation.

Variations of mesoscale eddy SST fronts based on an automatic detection method in the northern South China Sea

SST fronts at the mesoscale eddy edge(ME fronts)were investigated from 2007–2017 in the northern South China Sea(NSCS)based on an automatic method using satellite sea level anomaly(SLA)and SST data.The relative probabilities between the number of anticyclonic/cyclonic ME fronts(AEF/CEF)and the number of anticyclones/cyclones reached 20%.The northeastern and southwestern parts of these anticyclones had more fronts than the northwestern and southeastern parts,although CEFs were nearly equally distributed in all directions.The number of ME fronts had remarkable seasonal variations,while the eddy kinetic energy(EKE)showed no seasonal variations.The total EKE at the ME fronts was three times of that within the MEs,and it was much stronger in AEFs than in CEFs.The interannual variability in the number of ME fronts and EKE had no significant correlation with the El Ni?o-Southern Oscillation(ENSO)index.Possible mechanisms of ME fronts were discussed,but the contributions of mesoscale eddies to SST fronts need to be quantified in future studies.

Wall-modeling for large-eddy simulation of flows around an axisymmetric body using the diffuse-interface immersed boundary method

A novel method is proposed to combine the wall-modeled large-eddy simulation(LES) with the diffuse-interface direct-forcing immersed boundary(IB) method.The new developments in this method include:(i) the momentum equation is integrated along the wall-normal direction to link the tangential component of the effective body force for the IB method to the wall shear stress predicted by the wall model;(ii) a set of Lagrangian points near the wall are introduced to compute the normal component of the effective body force for the IB method by reconstructing the normal component of the velocity. This novel method will be a classical direct-forcing IB method if the grid is fine enough to resolve the flow near the wall. The method is used to simulate the flows around the DARPA SUBOFF model. The results obtained are well comparable to the measured experimental data and wall-resolved LES results.

Applicability of an eddy covariance system based on a close-path quantum cascade laser spectrometer for measuring nitrous oxide fluxes from subtropical vegetable fields

亚热带蔬菜地是氧化亚氮(N_2O)的重要排放源,本研究首次采用基于闭路量子级联激光吸收光谱仪的涡动相关(QCLAS-EC)法观测亚热带蔬菜地秋冬季非施肥阶段的N_2O排放通量,以评估QCLAS-EC法测量亚热带蔬菜地N_2O通量的适用性。结果表明,QCLAS-EC观测系统在野外条件下能长期稳定运行,可观测到N_2O排放的季节变化趋势,其检测限为18.5μg N m^(-2) h^(-1)(95%置信水平),试验期间获得的97.5%的N_2O通量大于此检测限,表明QCLAS-EC观测系统可有效并准确测量亚热带蔬菜地的N_2O排放通量,仪器本身的灵敏度不是通量检测的限制因素。

Sensitivity of Near Real-time MODIS Gross Primary Productivity in Terrestrial Forests Based on Eddy Covariance Measurements

As an important product of Moderate Resolution Imaging Spectroradiometer(MODIS), MOD17A2 provides dramatic improvements in our ability to accurately and continuously monitor global terrestrial primary production, which is also significant in effort to advance scientific research and eco-environmental management. Over the past decades, forests have moderated climate change by sequestrating about one-quarter of the carbon emitted by human activities through fossil fuels burning and land use/land cover change. Thus, the carbon uptake by forests reduces the rate at which carbon accumulates in the atmosphere. However, the sensitivity of near real-time MODIS gross primary productivity(GPP) product is directly constrained by uncertainties in the modeling process, especially in complicated forest ecosystems. Although there have been plenty of studies to verify MODIS GPP with ground-based measurements using the eddy covariance(EC) technique, few have comprehensively validated the performance of MODIS estimates(Collection 5) across diverse forest types. Therefore, the present study examined the degree of correspondence between MODIS-derived GPP and EC-measured GPP at seasonal and interannual time scales for the main forest ecosystems, including evergreen broadleaf forest(EBF), evergreen needleleaf forest(ENF), deciduous broadleaf forest(DBF), and mixed forest(MF) relying on 16 flux towers with a total of 68 site-year datasets. Overall, site-specific evaluation of multi-year mean annual GPP estimates indicates that the current MOD17A2 product works highly effectively for MF and DBF, moderately effectively for ENF, and ineffectively for EBF. Except for tropical forest, MODIS estimates could capture the broad trends of GPP at 8-day time scale for all other sites surveyed. On the annual time scale, the best performance was observed in MF, followed by ENF, DBF, and EBF. Trend analyses also revealed the poor performance of MODIS GPP product in EBF and DBF. Thus, improvements in the sensitivity of MOD17A2 to forest productivity require continued efforts.

New Method for 3D Transient Eddy Current Field Calculation and Its Application in Magneto-Acoustic Tomography

A new method of 3D transient eddy current field calculation is proposed. The Maxwell equations with time component elimination （METCE） are derived under the assumption of magnetic quasi static approximation, especially for the sample of low conductivity. Based on METCE, we deduce a more efficient reconstruction algorithm of a 3D transient eddy current field. The computational burden is greatly reduced through the new algorithm, and the computational efficiency is improved. This new algorithm decompounds the space-time variables into two individual variables. The idea is to solve the spatial vector component firstly, and then multiply it by the corresponded time component. The iterative methods based on METCE are introduced to recover the distribution of conductivity in magneto-acoustic tomography. The reconstructed images of conductivity are consistent with the original distribution, which validate the new method.

Eddy covariance measurements of turbulent fluxes in the surf zone

Turbulent eddies play a critical role in oceanic flows. Direct measurements of turbulent eddy fluxes beneath the sea surface were taken to study the direction of flux-carrying eddies as a means of supplementing our understanding of vertical fluxes exchange processes and their relationship to tides. The observations were made at 32 Hz at a water depth of ~1.5 m near the coast of Sanya, China, using an eddy covariance system, which mainly consists of an acoustic doppler velocimeter(ADV) and a fast temperature sensor. The cospectra-fit method-an established semi-empirical model of boundary layer turbulence to the measured turbulent cospectra at frequencies below those of surface gravity waves-was used in the presence of surface gravity waves to quantify the turbulent eddy fluxes(including turbulent heat flux and Reynolds stress). As much as 87% of the total turbulent stress and 88% of the total turbulent heat flux were determined as being at band frequencies below those of surface gravity waves. Both the turbulent heat flux and Reynolds stress showed a daily successive variation;the former peaked during the low tide period and the later peaked during the ebb tide period.Estimation of roll-off wavenumbers, k0, and roll-off wavelengths, λ0(where λ0=2π/k0), which were estimated as the horizontal length scales of the dominant flux-carrying turbulent eddies, indicated that the λ0 of the turbulent heat flux was approximately double that of the Reynolds stress. Wavelet analysis showed that both the turbulent heat flux and the Reynolds stress have a close relationship to the semi-diurnal and diurnal tides, and therefore indicate the energy that is transported from tides to turbulence.

Inspection of aluminum alloys by a multi-frequency eddy current method

The paper proposes an experimental method of material inspection,which is based on digital processing of multi-frequency eddy current measurement data.The influences of various factors(conductivity,the gap between the sample surface and the sensor,the thickness of the sample) on the obtained hodographs are examined by taking the aluminum alloys for example,and the possibility of separation of various factors is analyzed.The results obtained are indicative of how much promise the proposed method offers for the inspection and testing of products made of aluminum alloys.

A, ■-Ω METHOD FOR 3-D EDDY CURRENT ANALYSIS

After the field equations and the snonumuoo conditions between the interfaces for 3D eddy current problems Under various gauges were discussed, it was pointed cut in this paper that using the magnetic vector potential A. the electric scalar potential and Coulomb gauge △ .A = 0 in eddy current regions and using the magntetic scalar potential Ω in the non-conducting regions are more suitable. All field equations, the boundary conditions, the interface continuity conditions and the corresponding variational principle of this method are also given

Eddy Current Analyses by Domain Decomposition Method Using Double-Double Precision

A matrix equation solved in an eddy current analysis,??-??method based on a domain decomposition method becomes a complex symmetric system.In general,iterative method is used as the solver.Convergence of iterative method in an interface problem is improved by increasing an accuracy of a solution of an iterative method of a subdomain problem.However,it is difficult to improve the convergence by using a small convergence criterion in the subdomain problem.Therefore,authors propose a method to introduce double-double precision into the interface problem and the subdomain problem.This proposed method improves the convergence of the interface problem.In this paper,first,we describe proposed method.Second,we confirm validity of the method by using Team Workshop Problem 7,standard model for eddy current analysis.Finally,we show effectiveness of the method from two numerical results.

基于REddyProc的干旱区枣林通量数据插补及能量平衡分析

【目的】获取完整且有效的枣林通量数据,分析生态系统能量平衡,为评估干旱区枣林生态系统与大气间的能量和物质交换提供理论依据。【方法】选择R语言REddyProc包插补涡度相关法测量的通量数据,采用交叉验证法和能量平衡闭合度评价插补数据质量,并修正热储存项,分析2018、2019年新疆阿克苏地区枣林能量平衡闭合度与能量变化趋势。【结果】1)REddyProc包插补通量数据效果较好,交叉验证的误差统计参数为回归系数b=0.86~0.99、决定系数R^(2)=0.86~0.95、一致性指数d=0.96~0.98、模拟效率EF=0.84~0.92、均方根误差与观测值标准差比率RSR=0.28~0.40;2)经REddyProc包插补后,2018、2019年干旱区枣林能量平衡闭合度分别为73.45%、73.11%,有效能量和可利用能量的决定系数均为0.97;加入热储存项后,能量闭合度分别提高3.72%、2.75%,达77.17%、75.86%,增幅较小;3)各能量分项(净辐射、潜热通量、显热通量、土壤热通量)的日均变化规律在生育期和休眠期相似,日均变化均呈以净辐射变化规律为基础的单峰变化。【结论】干旱区枣林全年能量平衡闭合度符合ChinaFlux范围,热储存项对改善枣林能量平衡闭合度有一定影响,可为研究枣林能量和物质交换提供理论依据。

Finite Element Method Based on Equivalent Magnetic Energy Method for Computation of 2D Nonlinear Eddy Current Field

In this paper, the finite element method using vector potential in applications to 2D nonlinear eddy current field is discussed. The authors use the equivalent magnetic energy method to deal with magnetization curve of ferromagnetic material,and present the formulation of 2D nonlinear eddy current field.With this method the authors analyze the eddy current field in an induction ladle furnace and the force distribution in the charge (molten metal),and plot the corresponding curves.