Impact of Megacity Shanghai on the Urban Heat-Island Effects over the Downstream City Kunshan

The impact of upstream urbanization on the enhanced urban heat-island (UHI) effects between Shanghai and Kunshan is investigated by analyzing seven years of surface observations and results from mesoscale model simulations. The observational analysis indicates that, under easterly and westerly winds, the temperature difference between Shanghai and Kunshan increases with wind speed when the wind speed \(<\) 5 m s \(^{-1}\) . The Weather Research and Forecasting (WRF) numerical model, coupled with a one-layer urban canopy model (UCM), is used to examine the UHI structure and upstream effects by replacing the urban surface of Shanghai and/or Kunshan with cropland. The WRF/UCM modelling system is capable of reproducing the surface temperature and wind field reasonably well. The simulated urban canopy wind speed is a better representation of the near-surface wind speed than is the 10-m wind speed at the centre of Shanghai. Without the urban landscape of Shanghai, the surface air temperature over downstream Kunshan would decrease by 0.2–0.4  \(^{\circ }\) C in the afternoon and 0.4–0.6  \(^{\circ }\) C in the evening. In the simulation with the urban landscape of Shanghai, a shallow cold layer is found above the UHI, with a minimum temperature of about \(-0.2\) to \(-\) 0.5  \(^{\circ }\) C during the afternoon hours. Strong horizontal divergence is found in this cold layer. The easterly breeze over Shanghai is strengthened at the surface by strong UHI effects, but weakened at upper levels. With the appearance of the urban landscape specific humidity decreases by 0.5–1 g kg \(^{-1}\) within the urban area because of the waterproof property of an urban surface. On the other hand, the upper-level specific humidity is increased because of water vapour transferred by the strong upward vertical motions.     

A Semi-Analytical Approach for Parametrization of the Obukhov Stability Parameter in the Unstable Atmospheric Surface Layer

A semi-analytical scheme is proposed to parametrize the Obukhov stability parameter \(\zeta \) (= \(z/L\) ; \(z\) is the height above the ground and \(L\) is the Obukhov length) in terms of the bulk Richardson number ( \(R_{iB}\) ) in unstable conditions within the framework of Monin–Obukhov similarity (MOS) theory. The scheme involves, (i) a solution of a cubic equation in \(\zeta \) whose coefficients depend on the gradient Richardson number ( \(R_{i}\) ), and (ii) a relationship between \(R_{i}\) and \(R_{iB}\) . The proposed scheme is applicable for a wide range (i) \(-5\le R_{iB}\le 0\) , (ii) \(0\le \hbox {ln}(z_{0}/z_{h})\le 29.0\) , and (iii) \(10\le z/z_{0}\le 10^{5}\) and performs relatively better than all other schemes in terms of accuracy in computation of surface-layer transfer coefficients. The absolute errors in computing the transfer coefficients do not exceed 7 %. The analysis presented here is found to be valid for different \(\gamma _{m}\) and \(\gamma _{h}\) appearing in the expressions of the similarity functions \(\varphi _{m}\) and \(\varphi _{h}\) (representing non-dimensional wind and temperature profiles), so long as the ratio of \(\gamma _{m}\) to \(\gamma _{h} \ge 1\) . The improved scheme can be easily employed in atmospheric modelling for a comprehensive range of \(R_{iB}\) and a variety of surfaces.     

Storm Surges in New York During Hurricane Sandy in 2012: A Verification of the Wind-Stress Tide Relation

In October 2012 Hurricane Sandy devastated New York City and its vicinity caused mainly by the storm surge, which is the water height above normal astronomical tide level. The meteorological conditions were as follows: minimum central pressure, 962 hPa, highest sustained wind speed 27.1 m s $^{-1}$ ? 1 and maximum gust 37.8 m s $^{-1}$ ? 1 . The peak storm surge was at 3.9 m and the peak storm tide at 4.4 m (which is referenced above mean lower low water). The wind-stress tide relation shows that $S=K\,V^{2}$ S = K V 2 , where $S$ S is the storm surge, $V$ V is the wind speed and $K$ K is the coefficient. It is found that with $S$ S in units of m, and $V$ V in  m s $^{-1}$ ? 1 , $K = 0.0051$ K = 0.0051 with $R^{2}= 0.91$ R 2 = 0.91 ( $R$ R is the correlation coefficient) indicating that 91 % of the total variation of the storm surge can be explained by variations in the wind stress, which is proportional to $V^{2}$ V 2 . Similar results were obtained during Hurricane Irene in 2011, which also affected the New York area. Therefore, this simple wind stress-tide relation should be useful in coastal engineering, urban planning, and emergency management.     

Intercomparison of Methods for the Simultaneous Estimation of Zero-Plane Displacement and Aerodynamic Roughness Length from Single-Level Eddy-Covariance Data

We applied three approaches to estimate the zero-plane displacement $d$ through the aerodynamic measurement height $z$ (with $z = z_{m}- d$ and $z_{m}$ being the measurement height above the surface), and the aerodynamic roughness length $z_{0}$ , from single-level eddy covariance data. Two approaches (one iterative and one regression-based) were based on the universal function in the logarithmic wind profile and yielded an inherently simultaneous estimation of both $d$ and $z_{0}$ . The third approach was based on flux–variance similarity, where estimation of $d$ and consecutive estimation of $z_{0}$ are independent steps. Each approach was further divided into two methods differing either with respect to the solution technique (profile approaches) or with respect to the variable (variance of vertical wind and temperature, respectively). All methods were applied to measurements above a large, growing wheat field where a uniform canopy height and its frequent monitoring provided plausibility limits for the resulting estimates of time-variant $d$ and $z_{0}$ . After applying, for each approach, a specific data filtering that accounted for the range of conditions (e.g. stability) for which it is valid, five of the six methods were able to describe the temporal changes of roughness parameters associated with crop growth and harvest, and four of them agreed on $d$ to within 0.3 m most of the time. Application of the same methods to measurements with a more heterogeneous footprint consisting of fully-grown sugarbeet and a varying contribution of adjacent harvested fields exhibited a plausible dependence of the roughness parameters on the sugarbeet fraction. It also revealed that the methods producing the largest outliers can differ between site conditions and stability. We therefore conclude that when determining $d$ for canopies with unknown properties from single-level measurements, as is increasingly done, it is important to compare the results of a number of methods rather than rely on a single one. An ensemble average or median of the results, possibly after elimination of methods that produce outliers, can help to yield more robust estimates. The estimates of $z_{0}$ were almost exclusively physically plausible, although $d$ was considered unknown and estimated simultaneously with the methods and results described above.     

Air–Sea \mathrm{CO}_{2} Gas Transfer Velocity in a Shallow Estuary

The air–sea transfer velocity of $\mathrm{CO}_{2}\, (k_{\mathrm{CO}_{2}})$ was investigated in a shallow estuary in March to July 2012, using eddy-covariance measurements of $\mathrm{CO}_{2}$ fluxes and measured air–sea $\mathrm{CO}_{2}$ partial-pressure differences. A data evaluation method that eliminates data by nine rejection criteria in order to heighten parametrization certainty is proposed. We tested the data evaluation method by comparing two datasets: one derived using quality criteria related solely to the eddy-covariance method, and the other derived using quality criteria based on both eddy-covariance and cospectral peak methods. The best parametrization of transfer velocity normalized to a Schmidt number of 600 $(k_{600})$ was determined to be: $k_{600} = 0.3\,{U_{10}}^{2.5}$ where $U_{10}$ is the wind speed in m $\mathrm{s}^{-1}$ at 10 m; $k_{600}$ is based on $\mathrm{CO}_{2}$ fluxes calculated by the eddy-covariance method and including the cospectral peak method criteria. At low wind speeds, the transfer velocity in the shallow water estuary was lower than in other coastal waters, possibly a symptom of low tidal amplitude leading to low intensity water turbulence. High transfer velocities were recorded above wind speeds of 5 m $\mathrm{s}^{-1}$ , believed to be caused by early-breaking waves and the large fetch (6.5 km) of the estuary. These findings indicate that turbulence in both air and water influences the transfer velocity.     

Very-Large-Scale Motions in the Atmospheric Boundary Layer Educed by Snapshot Proper Orthogonal Decomposition

Large-eddy simulations of the atmospheric boundary layer (ABL) under a wide range of stabilities are conducted to educe very-large-scale motions and then to study their dynamics and how they are influenced by buoyancy. Preliminary flow visualizations suggest that smaller-scale motions that resemble hairpins are embedded in much larger scale streamwise meandering rolls. Using simulations that represent more than 150 h of physical time, many snapshots in the \(xy\) -, \(yz\) - and \(xz\) -planes are then collected to perform snapshot proper orthogonal decomposition and further investigate the large structures. These analyses confirm that large streamwise rolls that share several features with the very-large-scale motions observed in laboratory studies arise as the dominant modes under most stabilities, but the effect of the surface kinematic buoyancy flux on the energy content of these dominant modes is very significant. The first two modes in the \(yz\) -plane in the neutral case contain up to 3 % of the total turbulent kinetic energy; they also have a vertical tilt angle in the \(yz\) -plane of about 0 to 30 \(^\circ \) due to the turning effect associated with the Coriolis force. Unstable cases also feature streamwise rolls, but in the convective ABL they are strengthened by rising plumes in between them, with two to four rolls spanning the whole domain in the first few modes; the Coriolis effect is much weaker in the unstable ABL. These rolls are no longer the dominant modes under stable conditions where the first mode is observed to contain sheet-like motions with high turbulent kinetic energy. Using these proper orthogonal decomposition modes, we are also able to extract the vertical velocity fields corresponding to individual modes and then to correlate them with the horizontal velocity or temperature fields to obtain the momentum and heat flux carried by individual modes. Structurally, the fluxes are explained by the topology of their corresponding modes. However, the fraction of the fluxes produced by the modes is invariably smaller than the fraction of energy they contain, particularly under stable conditions where the first modes are found to perform weak counter-gradient fluxes.     

The Effect of Surface Heterogeneity on the Structure Parameters of Temperature and Specific Humidity: A Large-Eddy Simulation Case Study for the LITFASS-2003 Experiment

We conduct a high-resolution large-eddy simulation (LES) case study in order to investigate the effects of surface heterogeneity on the (local) structure parameters of potential temperature \(C_T^2\) and specific humidity \(C_q^2\) in the convective boundary layer (CBL). The kilometre-scale heterogeneous land-use distribution as observed during the LITFASS-2003 experiment was prescribed at the surface of the LES model in order to simulate a realistic CBL development from the early morning until early afternoon. The surface patches are irregularly distributed and represent different land-use types that exhibit different roughness conditions as well as near-surface fluxes of sensible and latent heat. In the analysis, particular attention is given to the Monin–Obukhov similarity theory (MOST) relationships and local free convection (LFC) scaling for structure parameters in the surface layer, relating \(C_T^2\) and \(C_q^2\) to the surface fluxes of sensible and latent heat, respectively. Moreover we study possible effects of surface heterogeneity on scintillometer measurements that are usually performed in the surface layer. The LES data show that the local structure parameters reflect the surface heterogeneity pattern up to heights of 100–200 m. The assumption of a blending height, i.e. the height above the surface where the surface heterogeneity pattern is no longer visible in the structure parameters, is studied by means of a two-dimensional correlation analysis. We show that no such blending height is found at typical heights of scintillometer measurements for the studied case. Moreover, \(C_q^2\) does not follow MOST, which is ascribed to the entrainment of dry air at the top of the boundary layer. The application of MOST and LFC scaling to elevated \(C_T^2\) data still gives reliable estimates of the surface sensible heat flux. We show, however, that this flux, derived from scintillometer data, is only representative of the footprint area of the scintillometer, whose size depends strongly on the synoptic conditions.     

Scalar Flux–Gradient Relationships Under Unstable Conditions over Water in Coastal Regions

The scalar flux–gradient relationships of temperature ( $\phi _{T}$ ? T ) and specific humidity ( $\phi _{q}$ ? q ) under unstable conditions are investigated using eddy-covariance measurements of air–sea turbulent fluxes and vertical profiles of temperature and specific humidity collected from a marine meteorological platform. The gradients of temperature and specific humidity are obtained from measurements at five heights above the sea surface using the log-square fitting method and the simpler first-order approximation method. The two methods yield similar results. The proposed flux–gradient relationships $\phi _{T}$ ? T and $\phi _{q}$ ? q covers a wide range of instability: the stability parameter $\zeta $ ζ ranges from $-$ ? 0.1 to $-$ ? 50. The functional form of the proposed flux–gradient relationships is an interpolation between the Businger–Dyer relation and the free convection relation, which includes the “ $-$ ? 1/2” and “ $-$ ? 1/3” scaling laws at two different stability regimes. The widely used COARE 3.0 algorithm, which is an interpolation between the integrals of the Businger–Dyer and the free convection relations, is also evaluated and compared. The analysis and comparisons show that both schemes generate reasonable values of $\phi _{q}$ ? q in the whole unstable regime. The COARE 3.0 algorithm, however, overestimates $\phi _{T}$ ? T values under very unstable conditions. The errors in the flux–gradient relationships induced by the random errors in the turbulence measurements are assessed. When the random errors are taken into account, the observations agree with predictions of various schemes fairly well, implying that the dominant transport mechanism is adequately captured by the Monin–Obukhov similarity theory. The study also shows that $\phi _{q}$ ? q is significantly ${>}\phi _{T}$ > ? T under unstable conditions and that the ratio $\phi _{q}/\phi _{T}$ ? q / ? T increases with $-\zeta $ ? ζ . The ratio of $\phi _{q}$ ? q to $\phi _{T}$ ? T and the ratio of turbulent transport efficiencies of heat and water vapour ( $R_{wT}/R_{wq}$ R wT / R wq ) suggest that heat is transported more efficiently than water vapour under unstable conditions.     

Advances in Air–Sea \hbox {CO}_2 Flux Measurement by Eddy Correlation

Eddy-correlation measurements of the oceanic \(\hbox {CO}_2\) flux are useful for the development and validation of air–sea gas exchange models and for analysis of the marine carbon cycle. Results from more than a decade of published work and from two recent field programs illustrate the principal interferences from water vapour and motion, demonstrating experimental approaches for improving measurement precision and accuracy. Water vapour cross-sensitivity is the greatest source of error for \(\hbox {CO}_2\) flux measurements using infrared gas analyzers, often leading to a ten-fold bias in the measured \(\hbox {CO}_2\) flux. Much of this error is not related to optical contamination, as previously supposed. While various correction schemes have been demonstrated, the use of an air dryer and closed-path analyzer is the most effective way to eliminate this interference. This approach also obviates density corrections described by Webb et al. (Q J R Meteorol 106:85–100, 1980). Signal lag and frequency response are a concern with closed-path systems, but periodic gas pulses at the inlet tip provide for precise determination of lag time and frequency attenuation. Flux attenuation corrections are shown to be \(<\) 5 % for a cavity ring-down analyzer (CRDS) and dryer with a 60-m inlet line. The estimated flux detection limit for the CRDS analyzer and dryer is a factor of ten better than for IRGAs sampling moist air. While ship-motion interference is apparent with all analyzers tested in this study, decorrelation or regression methods are effective in removing most of this bias from IRGA measurements and may also be applicable to the CRDS.     

Evolution of Stratocumulus Over Land: Comparison of Ground and Aircraft Observations with Numerical Weather Prediction Simulations

Forecasting of low cloud continues to challenge numerical weather prediction. With this in mind, surface and airborne observations were made over East Anglia, UK, during March 2011 to investigate stratus and stratocumulus advecting from the sea over land. Four surface sites were deployed at various distances inland aligned approximately along the flow. In situ data include cloud-droplet measurements from an aircraft operating off the coast and a tethered balloon 100 km inland. Comparisons of thermodynamic and cloud properties are made with Met Office operational model simulations at horizontal resolutions of 4 and 1.5 km. The clouds contained droplet concentrations up to 600 cm \(^{-3}\) within polluted outflow off Europe. These measurements were compared to three different model schemes for predicting droplet concentration: two of them perform well at low to moderate concentrations but asymptote to 375 cm \(^{-3}\) . Microwave radiometers at the ground sites retrieved liquid water paths that reduced with distance inland and were generally below 200 g m \(^{-2}\) . The modelled water path performs well upstream but more erratically far inland. Comparisons of thermodynamic profiles are made within both Lagrangian and Eulerian frameworks and show the model predicted changes in equivalent potential temperature generally within 1 K, with occasional errors of 2 K or more. The modelled cloud-top temperatures were in good agreement with the observations down to \(-\) 7  \(^{\circ }\) C, but the magnitude of the temperature inversion, although good at times, was too small by on average 1.6 K. The different simulations produced different cloud-top water contents due to a combination of resolution and scientific upgrades to the model, but they generally underestimate the amount of cloud water. Major changes, such as the mesoscale temporary cloud breaks on 2 March 2011 and the complete clearance on 4 March, were seemingly predicted by the model for the correct reasons.     

Impacts of Mixing Processes in Nocturnal Atmospheric Boundary Layer on Urban Ozone Concentrations

A number of open questions remain regarding the role of low-level jets (LLJs) and nocturnal mixing processes in the buildup of tropospheric ozone. The prevalence of southerly winds and LLJs in the U.S. Southern Great Plains during summer makes this region an ideal site for investigating the structure of the nocturnal boundary layer and its impacts on urban air quality. Ozone $(\mathrm{O}_{3})$ and nitrogen oxide concentrations measured at regulatory monitoring sites in the Oklahoma City (OKC) area and simulations with the Weather Research and Forecasting with Chemistry (WRF/Chem) model were analyzed to show how the nocturnal LLJ moderates boundary-layer mixing processes and air quality. Datasets collected during the Joint Urban 2003 campaign, which took place in July 2003 in OKC, provided detailed information about nocturnal boundary-layer structure and dynamics. In general, ${\mathrm{O}_{3}}$ time series show the expected behavior that urban ${\mathrm{O}_{3}}$ concentrations decrease at night due to nitrogen oxide titration reactions, but elevated ${\mathrm{O}_{3}}$ concentrations and secondary ${\mathrm{O}_{3}}$ peaks are also seen quite frequently after sunset. LLJs developed on most nights during the study period and were associated with strong vertical wind shear, which affected the boundary-layer stability and structure. Near-surface ${\mathrm{O}_{3}}$ concentrations are higher during less stable nights when active mixing persists throughout the night. The WRF/Chem model results agree well with the observations and further demonstrate the role of LLJs in moderating nocturnal mixing processes and air quality. The highest nocturnal ${\mathrm{O}_{3}}$ concentrations are linked to a strong LLJ that promotes both nocturnal long-range transport and persistent downward mixing of ${\mathrm{O}_{3}}$ from the residual layer to the surface.     

High Resolution Simulation of the Variability of Surface Energy Balance Fluxes Across Central London with Urban Zones for Energy Partitioning

The parameterization of surface heat-flux variability in urban areas relies on adequate representation of surface characteristics. Given the horizontal resolutions (e.g. $\approx $ 0.1–1 km) currently used in numerical weather prediction (NWP) models, properties of the urban surface (e.g. vegetated/built surfaces, street-canyon geometries) often have large spatial variability. Here, a new approach based on Urban Zones to characterize Energy partitioning (UZE) is tested within a NWP model (Weather Research and Forecasting model; WRF v3.2.1) for Greater London. The urban land-surface scheme is the Noah/Single-Layer Urban Canopy Model (SLUCM). Detailed surface information (horizontal resolution 1 km) in central London shows that the UZE offers better characterization of surface properties and their variability compared to default WRF-SLUCM input parameters. In situ observations of the surface energy fluxes and near-surface meteorological variables are used to select the radiation and turbulence parameterization schemes and to evaluate the land-surface scheme and choice of surface parameters. For radiative fluxes, improved performance (e.g. $>$ 25 W m $^{-2}$ root-mean-square error reduction for the net radiation) is attained with UZE parameters compared to the WRF v3.2.1 default for all three methods from the simplest to the most detailed. The UZE-based spatial fluxes reproduce a priori expectations of greater energy storage and less evaporation in the dense city centre compared to the residential surroundings. Problems in Noah/SLUCM partitioning of energy between the daytime turbulent fluxes are identified with the overestimation of the turbulent sensible heat and underestimation of the turbulent latent heat fluxes.     

Characterization of Oscillatory Motions in the Stable Atmosphere of a Deep Valley

In a valley sheltered from strong synoptic effects, the dynamics of the valley atmosphere at night is dominated by katabatic winds. In a stably stratified atmosphere, these winds undergo temporal oscillations, whose frequency is given by $N \sin {\alpha }$ N sin α for an infinitely long slope of constant slope angle $\alpha $ α , $N$ N being the buoyancy frequency. Such an unsteady flow in a stably stratified atmosphere may also generate internal gravity waves (IGWs). The numerical study by Chemel et al. (Meteorol Atmos Phys 203:187–194, 2009) showed that, in the stable atmosphere of a deep valley, the oscillatory motions associated with the IGWs generated by katabatic winds are distinct from those of the katabatic winds. The IGW frequency was found to be independent of $\alpha $ α and about $0.8N$ 0.8 N . Their study did not consider the effects of the background stratification and valley geometry on these results. The present work extends this study by investigating those effects for a wide range of stratifications and slope angles, through numerical simulations for a deep valley. The two oscillatory systems are reproduced in the simulations. The frequency of the oscillations of the katabatic winds is found to be equal to $N$ N times the sine of the maximum slope angle. Remarkably, the IGW frequency is found to also vary as $C_\mathrm{w}N$ C w N , with $C_\mathrm{w}$ C w in the range $0.7$ 0.7 – $0.95$ 0.95 . These values for $C_\mathrm{w}$ C w are similar to those reported for IGWs radiated by any turbulent field with no dominant frequency component. Results suggest that the IGW wavelength is controlled by the valley depth.     

Multi-Scale Sensible Heat Fluxes in the Suburban Environment from Large-Aperture Scintillometry and Eddy Covariance

Sensible heat fluxes ( \(Q_{H}\) ) are determined using scintillometry and eddy covariance over a suburban area. Two large-aperture scintillometers provide spatially integrated fluxes across path lengths of 2.8 and 5.5 km over Swindon, UK. The shorter scintillometer path spans newly built residential areas and has an approximate source area of 2–4 \(\text {km}^{2}\) , whilst the long path extends from the rural outskirts to the town centre and has a source area of around 5–10 \(\text {km}^{2}\) . These large-scale heat fluxes are compared with local-scale eddy-covariance measurements. Clear seasonal trends are revealed by the long duration of this dataset and variability in monthly \(Q_{H}\) is related to the meteorological conditions. At shorter time scales the response of \(Q_{H}\) to solar radiation often gives rise to close agreement between the measurements, but during times of rapidly changing cloud cover spatial differences in the net radiation ( \(Q^{*}\) ) coincide with greater differences between heat fluxes. For clear days \(Q_{H}\) lags \(Q^{*}\) , thus the ratio of \(Q_{H}\) to \(Q^{*}\) increases throughout the day. In summer the observed energy partitioning is related to the vegetation fraction through use of a footprint model. The results demonstrate the value of scintillometry for integrating surface heterogeneity and offer improved understanding of the influence of anthropogenic materials on surface-atmosphere interactions.     

Equilibrium Atmospheric Boundary-Layer Flow: Computational Fluid Dynamics Simulation with Balanced Forces

Forcing relationships in steady, neutrally stratified atmospheric boundary-layer (ABL) flow are thoroughly analyzed. The ABL flow can be viewed as balanced between a forcing and a drag term. The drag term results from turbulent stress divergence, and above the ABL, both the drag and the forcing terms vanish. In computational wind engineering applications, the ABL flow is simulated not by directly specifying a forcing term in the ABL but by specifying boundary conditions for the simulation domain. Usually, these include the inflow boundary and the top boundary conditions. This ‘boundary-driven’ ABL flow is dynamically different from its real counterpart, and this is the major reason that the simulated boundary-driven ABL flow does not maintain horizontal homogeneity. Here, first a dynamical approach is proposed to develop a neutrally stratified equilibrium ABL flow. Computational fluid dynamics (CFD) software (Fluent 6.3) with the standard \(k\) – \(\varepsilon \) turbulence model is employed, and by applying a driving force profile, steady equilibrium ABL flows are simulated by the model. Profiles of wind speed and turbulent kinetic energy (TKE) derived using this approach are reasonable in comparison with the conventional logarithmic law and with observational data respectively. Secondly, the equilibrium ABL profiles apply as inflow conditions to simulate the boundary-driven ABL flow. Simulated properties between the inlet and the outlet sections across a fetch of 10 km are compared. Although profiles of wind speed, TKE, and its dissipation rate are consistently satisfactory under higher wind conditions, a deviation of TKE and its dissipation rate between the inlet and outlet are apparent (7–8 %) under lower wind-speed conditions (2 m s \(^{-1}\) at 10 m). Furthermore, the simulated surface stress systematically decreases in the downwind direction. A redistribution of the pressure field is also found in the simulation domain, which provides a different driving pattern from the realistic case in the ABL.     

Transfer Coefficients of Momentum, Heat and Water Vapour in the Atmospheric Surface Layer of a Large Freshwater Lake

In studies of lake–atmosphere interactions, the fluxes of momentum, water vapour and sensible heat are often parametrized as being proportional to the differences in wind, humidity and air temperature between the water surface and a reference height above the surface. Here, the proportionality via transfer coefficients in these relationships was investigated with the eddy-covariance method at three sites within an eddy-covariance mesonet across Lake Taihu, China. The results indicate that the transfer coefficients decreased with increasing wind speed for weak winds and approached constant values for strong winds. The presence of submerged macrophytes reduced the momentum transfer (drag) coefficient significantly. At the two sites free of submerged macrophytes, the 10-m drag coefficients under neutral stability were 1.8 $(\pm \,0.4) \times \,10^{-3}$ ( ± 0.4 ) × 10 ? 3 and $1.7\,(\pm \,0.3) \times \,10^{-3 }$ 1.7 ( ± 0.3 ) × 10 ? 3 at the wind speed of $9\,\text{ m } \text{ s }^{-1}$ 9 m s ? 1 , which are 38 and 34 % greater than the prediction by the Garratt model for the marine environment.     

Land-Surface Heterogeneity Effects in the Planetary Boundary Layer

We investigate the cumulative added value of assimilating temperature, moisture, and wind observations in the three-dimensional non-hydrostatic Fifth-Generation Pennsylvania State University/National Center for Atmospheric Research Mesoscale Model MM5 and use these forecasts to analyze the relationship between surface forcing and planetary boundary-layer (PBL) depth. A data assimilation methodology focused on the surface and the PBL, previously tested in a one-dimensional version of MM5, is applied to 29 May, 6 June, and 7 June 2002 during the International $\hbox {H}_{2}\hbox {O}$ Project over the Southern Great Plains. Model-predicted PBL depth is evaluated against PBL depth diagnosed from data across 4,800 km of airborne lidar data (flight tracks 100–300 km long). The forecast with data assimilation verifies better against observations and is thus used to investigate the environmental conditions that govern PBL depth. The spatial structure in PBL depth is found to be most affected by spatial variations in surface buoyancy flux and capping inversion strength. The spatial scales of surface flux forcing reflected in the PBL depth are found through Fourier analysis and multiresolution decomposition. Correlations are ${<}0.50$ at scales of 64 km or less and increase at larger scales for 29 May and 6 June, but on 7 June low correlations are found at all scales, possibly due to greater within-PBL wind speeds, a stronger capping inversion on this day, and clouds. The results suggest a minimum scale, a function of wind speed, below which heterogeneity in surface buoyancy fluxes is not reflected directly in PBL depth.     

Some Statistics of the Temperature Structure Parameter in the Convective Boundary Layer Observed by Sodar

The characteristics of the temporal and height variations of the temperature structure parameter $C_\mathrm{T}^{2}$ in strongly convective situations derived from the sodar echo-signal intensity measurements were analyzed for the first 100 m. It was corroborated that the probability density function (pdf) of the logarithm of $C_\mathrm{T}^{2}$ in the lower convective boundary layer is markedly non-Gaussian, whereas turbulence theory predicts it to be normal. It was also corroborated that the sum of two weighted Gaussians, which characterize the statistics of $C_\mathrm{T}^{2}$ within convective plumes and in their environment and the probability of plume occurrence, well approximates the observed pdfs. It was shown that the height behaviour of the arithmetic mean of $ C_\mathrm{T}^{2}$ (both total and within plumes) follows well a power law $C_\mathrm{T}^{2} (z) \sim z^{-q}$ with the exponent $q$ close to the theoretically predicted value of 4/3. But for the geometrical means of $C_\mathrm{T}^{2}$ (both total and within the plumes), $q$ is close to 1. The difference between arithmetically and geometrically averaged $C_\mathrm{T}^{2}$ profiles was analyzed. The vertical profiles of the standard deviation, skewness and kurtosis of $\hbox {ln}C_\mathrm{T}^{2}$ pdfs were analyzed to show their steady behaviour with height. The standard deviations of the logarithm of $C_\mathrm{T}^{2}$ within the plumes and between them are similar and are 1.5 times less than the total standard deviation. The estimate of the variability index $F_\mathrm{T}$ and its height behaviour were obtained, which can be useful to validate some theoretical and modelling predictions. The vertical profiles of the skewness and kurtosis show the negative asymmetry of pdfs and their flatness, respectively. The spectra of variations in $\hbox {ln}C_\mathrm{T}^{2}$ are shown to be satisfactorily fitted by the power law $f^{-\gamma } $ in the frequency range 0.02 and 0.2 Hz, with the average exponent $\approx $ 1.27  $\pm $  0.22.     

Variability of the Structure Parameters of Temperature and Humidity Observed in the Atmospheric Surface Layer Under Unstable Conditions

The structure parameters of temperature and humidity are important in scintillometry as they determine the structure parameter of the refractive index of air, the primary atmospheric variable obtained with scintillometers. In this study, we investigate the variability of the logarithm of the Monin-Obukhov-scaled structure parameters (denoted as $\log ({\widetilde{C_{s}^2}_{\mathrm {}}})$ ) of temperature and humidity. We use observations from eddy-covariance systems operated at three heights (2.5, 50, and 90 m) within the atmospheric surface layer under unstable conditions. The variability of $\log ({\widetilde{C_{s}^2}_{\mathrm {}}})$ depends on instability and on the size of the averaging window over which $\log ({\widetilde{C_{s}^2}_{\mathrm {}}})$ is calculated. If instability increases, differences in $\log ({\widetilde{C_{s}^2}_{\mathrm {}}})$ between upward motions (large $C_{s}^2$ ) and downward motions (small $C_{s}^2$ ) increase. The differences are, however, not sufficiently large to result in a bimodal probability density function. If the averaging window size increases, the variances of $\log ({\widetilde{C_{s}^2}_{\mathrm {}}})$ decrease. A linear regression of the variances of $\log ({\widetilde{C_{s}^2}_{\mathrm {}}})$ versus the averaging window size for various stability classes shows an increase of both the offset and slope (in absolute sense) with increasing instability. For temperature, data from the three heights show comparable results. For humidity, in contrast, the offset and slope are larger at 50 and 90 m than at 2.5 m. In the end we discuss how these findings could be used to assess whether observed differences in $C_{s}^2$ along a scintillometer path or aircraft flight leg are just within the range of local variability in $C_{s}^2$ or could be attributed to surface heterogeneity. This is important for the interpretation of data measured above a heterogeneous surface.