Comparison of ERA5 turbulent fluxes at the air-sea interface with measurements from a wave-following platform

Turbulent fluxes at the air-sea interface were estimated with data collected in 2011-2020 with a low-profile platform named OCARINA during eight experiments in five regions:2011,2015,and 2016 in the Iroise Sea;2012 in the tropical Atlantic;2014 in the Chilie-Peru upwelling;2017 and 2018 in the Mediterranean Sea,and 2018 and 2020 in Barbados.The observations were carried out with moderate winds(2-10 m s^(-1))and average wave heights of 1.5 m.In this study,the authors used the fluxes calculated by the bulk method using OCARINA-sampled data as the input.These data can validate the fluxes estimated from ERAS reanalysis data.The OCARINA and ERA5 data were taken concomitantly.To do this,the authors established an algorithm to extract the OCARINA data as closely as possible to the reanalysis data in time and position.The measurements of the OCARINA platform can conclude on the relevance of the widely used reanalysis data.

Evaluating Parameterizations for Turbulent Fluxes over the Landfast Sea-Ice Surface in Prydz Bay, Antarctica

It is crucial to appropriately determine turbulent fluxes in numerical models.Using data collected in East Antarctica from 8 April to 26 November 2016,this study evaluates parameterization schemes for turbulent fluxes over the landfast seaice surface in five numerical models.The Community Noah Land Surface Model with Multi-Parameterizations Options(Noah_mp)best replicates the turbulent momentum flux,while the Beijing Climate System Model(BCC_CSM)produces the optimum sensible and latent heat fluxes.In particular,two critical issues of parameterization schemes,stability functions and roughness lengths,are investigated.Sensitivity tests indicate that roughness lengths play a decisive role in model performance.Based on the observed turbulent fluxes,roughness lengths over the landfast sea-ice surface are calculated.The results,which can provide a basis for setting up model parameters,reveal that the dynamic roughness length(z0m)increases with the increase of frictional velocity(u*)when u*≤0.4 m s^(−1) and fluctuates around 10^(−3 )m when u*>0.4 m s^(−1);thermal roughness length(z0t)is linearly related to the temperature gradient between air and sea-ice surface(ΔT)with a relation of lg(z0t)=−0.29ΔT−3.86;and the mean water vapor roughness length(z0q)in the specific humidity gradient(Δq)range ofΔq≤−0.6 g kg^(−1) is 10^(−6) m,3.5 times smaller than that in the range ofΔq˃−0.6 g kg^(−1).

Estimation of Turbulent Fluxes Using the Flux-Variance Method over an Alpine Meadow Surface in the Eastern Tibetan Plateau

The flux-variance similarity relation and the vertical transfer of scalars exhibit dissimilarity over different types of surfaces, resulting in different parameterization approaches of relative transport efficiency among scalars to estimate turbulent fluxes using the flux-variance method. We investigated these issues using eddycovariance measurements over an open, homogeneous and flat grassland in the eastern Tibetan Plateau in summer under intermediate hydrological conditions during rainy season. In unstable conditions, the temperature, water vapor, and CO2 followed the flux-variance similarity relation, but did not show in precisely the same way due to different roles （active or passive） of these scalars. Similarity constants of temperature, water vapor and CO2 were found to be 1.12, 1.19 and 1.17, respectively. Heat transportation was more eft% cient than water vapor and CO2. Based on the estimated sensible heat flux, five parameterization methods of relative transport efficiency of heat to water vapor and CO2 were examined to estimate latent heat and CO2 fluxes. The strategy of local determination of flux-variance similarity relation is recommended for the estimation of latent heat and CO2 fluxes. This approach is better for representing the averaged relative transport efficiency, and technically easier to apply, compared to other more complex ones.

Analysis of Turbulence Characteristics over the Northern Tibetan Plateau Area

Based on CAMP/Tibet [Coordinated Enhanced Observing Period （CEOP） Asia-Australia Monsoon Project （CAMP） on the Tibetan Plateau] turbulent data collected at the Bujiao （B J） site of the Nagqu area, the turbulent structure and transportation characteristics in the near surface layer during summer are analyzed. The main results show that the relationship between the normalized standard deviation of 3D wind speed and stability satisfies the similarity law under both unstable and stable stratifications. The relations of normalized standard deviation of temperature and specific humidity to stability only obey the ＂-1/3 power law＂ under unstable conditions. In the case of stable stratifications, their relations to stability are dispersing. The sensible heat dominates in the dry period, while in the wet period, the latent heat is larger than the sensible heat.

Characteristics of land-atmosphere energy and turbulentfluxes over the plateau steppe in central Tibetan Plateau

The land-atmosphere energy and turbulence exchange is key to understanding land surface processes on the Tibetan Plateau(TP). Using observed data for Aug. 4 to Dec. 3, 2012 from the Bujiao observation point(BJ) of the Nagqu Plateau Climate and Environment Station(NPCE-BJ), different characteristics of the energy flux during the Asian summer monsoon(ASM) season and post-monsoon period were analyzed. This study outlines the impact of the ASM on energy fluxes in the central TP. It also demonstrates that the surface energy closure rate during the ASM season is higher than that of the post-monsoon period. Footprint modeling shows the distribution of data quality assessments(QA) and quality controls(QC) surrounding the observation point. The measured turbulent flux data at the NPCE-BJ site were highly representative of the target land-use type. The target surface contributed more to the fluxes under unstable conditions than under stable conditions. The main wind directions(180° and 210°) with the highest data density showed flux contributions reaching 100%, even under stable conditions. The lowest flux contributions were found in sectors with low data density, e.g., 90.4% in the 360° sector under stable conditions during the ASM season. Lastly, a surface energy water balance(SEWAB) model was used to gap-fill any absent or corrected turbulence data. The potential simulation error was also explored in this study. The Nash-Sutcliffe model efficiency coefficients(NSEs) of the observed fluxes with the SEWAB model runs were 0.78 for sensible heat flux and 0.63 for latent heat flux during the ASM season, but unrealistic values of-0.9 for latent heat flux during the post-monsoon period.

Evaluation of monthly turbulent heat fluxes from WHOI analysis and NCEP reanalysis in the tropical Atlantic

The biases and their sources in monthly turbulent heat fluxes from the Woods Hole Oceanographic Institution （WHOI） analysis, and the National Centers for Environmental Prediction-National Center for Atmospheric Research reanalyses 1 and 2 （NCEPI and NCEP2） are checked in the climatically representative regions in the tropical Atlantic using the fluxes from the Southampton Oceanographic Centre （SOC） and the pilot research moored array in the tropical Atlantic （PIRATA） as references. For the WHOI analysis, the biases in turbulent heat fluxes mainly exist in equatorial regions which are due to the overestimated sea surface temperature and the underestimated 2 m air humidity. For the NCEP2 reanalysis, the maximum biases, about （40±5） W/m^2, exist in southeast and northeast trade wind regions, which are mainly caused by the flux algorithm used because the biases in wind speed and air-sea humidity difference are relatively small. In the equatorial regions, the flux biases in the NCEP2 derived from both flux-related basic variables and algorithm are equally large. Although the estimations of time series trends and air-sea humidity difference of the NCEPI are improved greatly in the NCEP2, the biases of latent heat flux in the NCEP2 are about 20 W/m^2 greater than those from the NCEP1 in the trade wind regions. The result shows that the climatologies and monthly variabilities of the turbulent heat fluxes from the WHOI are more accurate than those from the NCEP1 and NCEP2 in the tropical Atlantic, especially on outside of the equatorial regions.

A Note on Reviving the Goddard Satellite-Based Surface Turbulent Fluxes (GSSTF) Dataset

Accurate sea surface flux measurements are crucial for understanding the global water and energy cycles. The oceanic evaporation, which is a major component of the global oceanic fresh water flux, is useful for predicting oceanic circulation and transport. The global Goddard Satellite-based Surface Turbulent Fluxes Version-2 （GSSTF2; July 1987–December 2000） dateset that was o？cially released in 2001 has been widely used by scientific community for global energy and water cycle research, and regional and short period data analyses. We have recently been funded by NASA to resume processing the GSSTF dataset with an objective of continually producing a uniform dataset of sea surface turbulent fluxes, derived from remote sensing data. The dataset is to be reprocessed and brought up-to-date （GSSTF2b） using improved input datasets such as a recently upgraded NCEP/DOE sea surface temperature reanalysis, and an upgraded surface wind and microwave brightness temperature V6 dataset （Version 6） from the Special Sensor Microwave Imager （SSM/I） produced by Remote Sensing Systems （RSS）. A second new product （GSSTF3） is further proposed with a finer temporal （12-h） and spatial （0.25° × 0.25°） resolution. GSSTF2b （July 1987–December 2008） and GSSTF3 （July 1999–December 2009） will be released for the research community to use by late 2009 and early 2011, respectively.

Seasonal variability of turbulent heat fluxes in the tropical Atlantic Ocean based on WHOI flux product

The mean seasonal variability of turbulent heat fluxes in the tropical Atlantic Ocean is examined using the Woods Hole Oceanographic Institution （WHOI） flux product. The most turbulent heat fluxes occur during winter seasons in the two hemispheres, whose centers are located at 10° -20°N and 5° 15°S respectively. In climatological ITCZ, the turbulent heat fluxes are the greatest from June to August, and in equatorial cold tongue the turbulent heat fluxes are the greatest from March to May. Seasonal variability of sensible heat flux is smaller than that of latent heat flux and mainly is dominated by the variations of air-sea temperature difference. In the region with larger climatological mean wind speed （air-sea humidity difference）, the variations of air-sea humidity difference （wind speed） dominate the variability of latent heat flux. The characteristics of turbulent heat flux yielded from theory analysis and WHOI dataset is consistent in physics which turns out that WHOI＇ s flux data are pretty reliable in the tropical Atlantic Ocean.

Effect of Counter-Gradient in the Computation of Turbulent Fluxes of Heat and Moisture in a Regional Model

The counter-gradient terms in the computations of turbulent fluxes of heat and moisture have been included in the PBL parameterization of a regional model for monsoon prediction. Results show that inclusion of counter-gradient terms has a marginal impact in the prediction of large scale monsoon circulation and rainfall rates.

Observation Research of the Turbulent Fluxes of Momentum, Sensible Heat and Latent Heat over the West Pacific Tropical Ocean Area

This paper describes results of the fluxes of momentum , sensible heat and latent heat for the West Pacific Tropical Ocean Area ( 127 ° E - 150 ° E , 5 ° N -3 ° S ). The data were collected by the small tethered balloon sounding system over this ocean area including 6 continuous stations (140 ° E. 0 ° ; 145 ° E, 0 ° ; 150 ° E, 0 ° ; 140° E, 5 ° N; 145 ° E, 5° N and 150 ° E, 5 ° N) from 11 October to 15 December, 1986 . These fluxes were calculated by the semiempirical flux-profile relationships of Monin-Obukhov similarity theory using these observed data. The results show that for this tropical ocean area the drag coefficient CD is equal to (1.53 ± 0.25) × 10 3 and the daily mean latent flux Hl is greater than its daily mean sensible flux HV by a factor of about 9.

Estimations of Land Surface Characteristic Parameters and Turbulent Heat Fluxes over the Tibetan Plateau Based on FY-4A/AGRI Data

Accurate estimates of land surface characteristic parameters and turbulent heat fluxes play an important role in the understanding of land-atmosphere interaction. In this study, Fengyun-4A (FY-4A) Advanced Geostationary Radiation Imager (AGRI) satellite data and the China Land Data Assimilation System (CLDAS) meteorological forcing dataset CLDAS-V2.0 were applied for the retrieval of broadband albedo, land surface temperature (LST), radiation flux components, and turbulent heat fluxes over the Tibetan Plateau (TP). The FY-4A/AGRI and CLDAS-V2.0 data from 12 March 2018 to 30 April 2018 were first used to estimate the hourly turbulent heat fluxes over the TP. The time series data of in-situ measurements from the Tibetan Observation and Research Platform were divided into two halves-one for developing retrieval algorithms for broadband albedo and LST based on FY-4A, and the other for the cross validation. Results show the root-mean-square errors (RMSEs) of the FY-4A retrieved broadband albedo and LST were 0.0309 and 3.85 K, respectively, which verifies the applicability of the retrieval method. The RMSEs of the downwelling/upwelling shortwave radiation flux and downwelling/upwelling longwave radiation flux were 138.87/32.78 W m^(−2) and 51.55/17.92 W m^(−2), respectively, and the RMSEs of net radiation flux, sensible heat flux, and latent heat flux were 58.88 W m^(−2), 82.56 W m^(−2) and 72.46 W m^(−2), respectively. The spatial distributions and diurnal variations of LST and turbulent heat fluxes were further analyzed in detail.

Surface Turbulent Flux Measurements over the Loess Plateau for a Semi-Arid Climate Change Study

In order to provide high quality data for climate change studies, the data quality of turbulent flux measurements at the station of SACOL （Semi-Arid Climate ＆ Environment Observatory of Lanzhou University）, which is located on a semi-arid grassland over the Loess Plateau in China, has been analyzed in detail. The effects of different procedures of the flux corrections on CO2, momentum, and latent and sensible heat fluxes were assessed. The result showed that coordinate rotation has a great influence on the momentum flux but little on scalar fluxes. For coordinate rotation using the planar fit method, different regression planes should be determined for different wind direction sectors due to the heterogeneous nature of the ground surface. Sonic temperature correction decreased the sensible heat flux by about 9%, while WPL correction （correction for density fluctuations） increased the latent heat flux by about 10%. WPL correction is also particularly important for CO2 fluxes. Other procedures of flux corrections, such as the time delay correction and frequency response correction, do not significantly influence the turbulent fluxes. Furthermore, quality tests on stationarity and turbulence development conditions were discussed. Parameterizations of integral turbulent characteristics （ITC） were tested and a specific parameterization scheme was provided for SACOL. The ITC test on turbulence development conditions was suggested to be applied only for the vertical velocity. The combined results of the quality tests showed that about 62%-65% of the total data were of high quality for the latent heat flux and CO2 flux, and as much as about 76% for the sensible heat flux. For the momentum flux, however, only about 35% of the data were of high quality.

荒漠戈壁下垫面表面动量和感热湍流通量参数化研究(英文)

The parameterization of surface turbulent fluxes over the Gobi Desert in arid regions is studied by using rationally screened observational data. First, the characteristics of Monin-Obukhov similarity functions are analyzed and their empirical formulae are fitted. The results show that fitted curves of changes of similarity functions of wind speed and temperature with stability parameter differ little from the typical empirical curves and are within the ranges of scatter of the empirical curves, but their values in the neutral condition arc different from the typical values to some extent. Furthermore, average values of momentum and scalar (sensible heat) roughness lengths as well as changes of scalar roughness length with friction velocity are determined by utilizing the data. It is found that the average values of scalar roughness length are about one order smaller than that of the momentum roughness length and decrease with increasing friction velocity, but they are evidently larger than their theoretically forecasted values.

Turbulent Variance Characteristics of Temperature and Humidity over a Non-uniform Land Surface for an Agricultural Ecosystem in China

This paper describes the application of the variance method for flux estimation over a mixed agricultural region in China. Eddy covariance and flux variance measurements were conducted in a near-surface layer over a non-uniform land surface in the central plain of China from 7 June to 20 July 2002. During this period, the mean canopy height was about 0.50 m. The study site consisted of grass （10% of area）, beans （15%）, corn （15%） and rice （60%）. Under unstable conditions, the standard deviations of temperature and water vapor density （normalized by appropriate scaling parameters）, observed by a single instrument, followed the Monin-Obukhov similarity theory. The similarity constants for heat （CT） and water vapor （Cq） were 1.09 and 1.49, respectively. In comparison with direct measurements using eddy covariance techniques, the flux variance method, on average, underestimated sensible heat flux by 21% and latent heat flux by 24%, which may be attributed to the fact that the observed slight deviations （20% or 30% at most） of the similarity ＂constants＂ may be within the expected range of variation of a single instrument from the generally-valid relations.

Simulation of Snow Processes Beneath a Boreal Scots Pine Canopy

A physically-based multi-layer snow model Snow-Atmosphere-Soil-Transfer scheme （SAST） and a land surface model Biosphere-Atmosphere Transfer Scheme （BATS） were employed to investigate how boreal forests influence snow accumulation and ablation under the canopy. Mass balance and energetics of snow beneath a Scots pine canopy in Finland at different stages of the 2003-2004 and 2004-2005 snow seasons are analyzed. For the fairly dense Scots pine forest, drop-off of the canopy-intercepted snow contributes, in some cases, twice as much to the underlying snowpack as the direct throughfall of snow. During early winter snow melting, downward turbulent sensible and condensation heat fluxes play a dominant role together with downward net longwave radiation. In the final stage of snow ablation in middle spring, downward net all- wave radiation dominates the snow melting. Although the downward sensible heat flux is comparable to the net solar radiation during this period, evaporative cooling of the melting snow surface makes the turbulent heat flux weaker than net radiation. Sensitivities of snow processes to leaf area index （LAI） indicate that a denser canopy speeds up early winter snowmelt, but also suppresses melting later in the snow season. Higher LAI increases the interception of snowfall, therefore reduces snow accumulation under the canopy during the snow season; this effect and the enhancement of downward longwave radiation by denser foliage outweighs the increased attenuation of solar radiation, resulting in earlier snow ablation under a denser canopy. The difference in sensitivities to LAI in two snow seasons implies that the impact of canopy density on the underlying snowpack is modulated by interannual variations of climate regimes.

Discrepancy and Applicability of Various Similarity Functions in Flux Calculations Under Stable Conditions

Based on data obtained during the Hualhe River Basin Experiment （HUBEX） in 1999, this study intends to detect the quantitative discrepancies in the momentum （τ0）, sensible heat （H0） and latent heat （E0） fluxes among six sets of similarity functions with the aerodynamic method. It also aims to clarify the applicability of the functions under stable conditions. The relative discrepancy was studied with the normalized transfer coefficients for τ0, H0 and E0, namely CD, CH and CQ, respectively. Except for one set of functions that adopted a rather small von Kármán＇s constant （0.365）, the relative discrepancy in τ0 among the other functions was less than 10%, while that in H0（E0） sometimes reached 25% when the bulk Richardson number （R/B） was less than 0.07. The absolute discrepancy in the fluxes was studied with statistical computations. Among the six sets of functions, the discrepancy in τ0, H0 and E0 sometimes reached 0.03 kg m^-1 s^-2, 4 W m^-2 and 10 W m^-2, respectively, and the discrepancy in the energy balance ratio sometimes exceeded 0.1. Furthermore, when RiB exceeded the critical value （Ric） for a specific set of functions, no fluxes could be derived with the functions. It is therefore suggested that RiB be compared with Ric before computing the fluxes if RiB is less than Ric. Finally, two sets of nonlinear similarity functions are recommended, due to their unlimited applicability in terms of RiB.

非均匀地表类型有效粗糙度的计算和特征分析

A new method is proposed for the effective roughness length (ERL) in heterogeneous terrain based on the principle of equalisation of momentum or heat fluxes calculated by the drag coefficient parameterization scheme used in the ECMWF numerical model. Compared with the area-weighted logarithmically averaged ERL (drag coefficient), the newly calculated ERL (drag coefficient) is about 40% (16%) larger with a roughness step of 2.3. These differences reach their maximum values when the ratio of smooth to rough surface is 60% to 40%. Since the determination by this method is not sensitive to the atmospheric stratification, it is suitable for use in climate models.

Seasonal variations of the near surfacelayer parameters over the Antarctic ice sheet in Princess Elizabeth Land,East Antarctica

Analysis of sensible heat flux （ Qh ）, latent heat flux （ Qe ）, Richardson number （Ri） ,bulk transport coefficient （Cd） and katabatic windsare presented by using the meteorological data in the near surface layer from an automatic weather station （AWS） in Princess Elizabeth Land, East Antarctica ice sheet and the data of corresponding period at Zhongshan station in 2002. It shows that annual mean air temperature at LGB69 is -25.6℃, which is 16.4℃ lower than that at Zhongshan, where the elevation is lower and located on the coast. The temperature lapse rate is about 1.0℃/110 m for the initial from coast to inland. The turbulence heat flux at LGB69 displays obvious seasonal variations with the average sensible heat flux -17.9 W/m^2 and latent heat flux -0.9 W/m^2. The intensity （Qh ＋ Qe ） of coolling source is - 18.8 W/m^2 meaning the snow surface layer obtains heat from atmosphere. The near surface atmosphere is near-neutral stratified with bulk transport coefficients （Cd） around 2.8 ×10^-3 ,and it is near constant when the wind speed higher than 8 m/s. The speed and the frequency of easterly Katabatic winds at LGB69 were higher than that at Zhongshan Station.

Numerical simulation of the climate effect of high-altitude lakes on the Tibetan Plateau

Lakes regulate the water and heat exchange between the ground and the atmosphere on different temporal and spatial scales. However, studies of the lake effect in the high-altitude Tibetan Plateau(TP) rarely have been performed until recently, and little attention has been paid to modelling of frozen lakes. In this study, the Weather Research and Forecasting Model(WRF v. 3.6.1) is employed to conduct three numerical experiments in the Ngoring Lake Basin(the original experiment, an experiment with a tuned model, and a no-lake experiment) to investigate the influences of parameter optimization on the lake simulation and of the high-altitude lake on the regional climate. After the lake depth, the roughness lengths, and initial surface temperature are corrected in the model, the simulation of the air temperature is distinctly improved. In the experiment using a tuned model, the simulated sensible-heat flux(H) is clearly improved, especially during periods of ice melting(from late spring to early summer) and freezing(late fall). The improvement of latent-heat flux(LE) is mainly manifested by the sharp increase in the correlation coefficient between simulation and observation, whereas the improvement in the average value is small. The optimization of initial surface temperature shows the most prominent effect in the first year and distinctly weakens after a freezing period. After the lakes become grassland in the model, the daytime temperature clearly increases during the freezing and melting periods; but the nocturnal cooling appears in other stages, especially from September to October. The annual mean H increases by 6.4 times in the regions of the Ngoring Lake and the Gyaring Lake, and the LE declines by 56.2%. The sum of H and LE increases from 71.2 W/m2(with lake) to 84.6 W/m2(no lake). For the entire simulation region, the sum of H and LE also increases slightly. After the lakes are removed, the air temperature increases significantly from June to September over the area corresponding to the two lakes, and an abnormal convergence field appears; at the same time, the precipitation clearly increases over the two lakes and surrounding areas.

A Sensitivity Study of Arctic Ice-Ocean Heat Exchange to the Three-Equation Boundary Condition Parametrization in CICE6

In this study,we perform a stand-alone sensitivity study using the Los Alamos Sea ice model version 6(CICE6)to investigate the model sensitivity to two Ice-Ocean(IO)boundary condition approaches.One is the two-equation approach that treats the freezing temperature as a function of the ocean mixed layer(ML)salinity,using two equations to parametrize the IO heat exchanges.Another approach uses the salinity of the IO interface to define the actual freezing temperature,so an equation describing the salt flux at the IO interface is added to the two-equation approach,forming the so-called three-equation approach.We focus on the impact of the three-equation boundary condition on the IO heat exchange and associated basal melt/growth of the sea ice in the Arctic Ocean.Compared with the two-equation simulation,our three-equation simulation shows a reduced oceanic turbulent heat flux,weakened basal melt,increased ice thickness,and reduced sea surface temperature(SST)in the Arctic.These impacts occur mainly at the ice edge regions and manifest themselves in summer.Furthermore,in August,we observed a downward turbulent heat flux from the ice to the ocean ML in two of our three-equation sensitivity runs with a constant heat transfer coefficient(0.006),which caused heat divergence and congelation at the ice bottom.Additionally,the influence of different combinations of heat/salt transfer coefficients and thermal conductivity in the three-equation approach on the model simulated results is assessed.The results presented in this study can provide insight into sea ice model sensitivity to the three-equation IO boundary condition for coupling the CICE6 to climate models.