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.     

Vertical Distribution of Radiation and Energy Balance Partitioning Within and Above a Lodgepole Pine Stand Recovering from a Recent Insect Attack

The current outbreak of mountain pine beetle (MPB) that started in the late 1990s in British Columbia, Canada, is the largest ever recorded in the north American native habitat of the beetle. The killing of trees is expected to change the vertical distribution of net radiation ( $Q^*$ Q ? ) and the partitioning of latent ( $Q_\mathrm{E}$ Q E ) and sensible ( $Q_\mathrm{H}$ Q H ) heat fluxes in the different layers of an attacked forest canopy. During an intensive observation period in the summer of 2010, eddy-covariance flux and radiation measurements were made at seven heights from ground level up to 1.34 times the canopy height in an MPB-attacked open-canopy forest stand $(\hbox {leaf area index} = 0.55~\mathrm{{m}}^{2}\ \mathrm{{m}}^{-2})$ ( leaf area index = 0.55 m 2 m - 2 ) in the interior of British Columbia, Canada. The lodgepole pine dominated stand with a rich secondary structure (trees and understorey not killed by the beetle) was first attacked by the MPB in 2003 and received no management. In this study, the vertical distribution of the energy balance components and their sources and sinks were analyzed and energy balance closure (EBC) was determined for various levels within the canopy. The low stand density resulted in approximately 60 % of the shortwave irradiance and 50 % of the daily total $Q^*$ Q ? reaching the ground. Flux divergence calculations indicated relatively strong sources of latent heat at the ground and where the secondary structure was located. Only very weak sources of latent heat were found in the upper part of the canopy, which was mainly occupied by dead lodgepole pine trees. $Q_\mathrm{H}$ Q H was the dominant term throughout the canopy, and the Bowen ratio ( $Q_\mathrm{H}/Q_\mathrm{E}$ Q H / Q E ) increased with height in the canopy. Soil heat flux ( $Q_\mathrm{G}$ Q G ) accounted for approximately 4 % of $Q^*$ Q ? . Sensible heat storage in the air ( $\Delta Q_\mathrm{S,H}$ Δ Q S , H ) was the largest of the energy balance storage components in the upper canopy during daytime, while in the lower canopy sensible heat storage in the boles ( $\Delta Q_\mathrm{S,B}$ Δ Q S , B ) and biochemical energy storage ( $\Delta Q_\mathrm{S,C}$ Δ Q S , C ) were the largest terms. $\Delta Q_\mathrm{S,H}$ Δ Q S , H was almost constant from the bottom to above the canopy. $\Delta Q_\mathrm{S,C}$ Δ Q S , C , $\Delta Q_\mathrm{S,B}$ Δ Q S , B and latent heat storage in the air ( $\Delta Q_\mathrm{S,E}$ Δ Q S , E ) varied more than $\Delta Q_\mathrm{S,H}$ Δ Q S , H throughout the canopy. During daytime, energy balance closure was high in and above the upper canopy, and in the lowest canopy level. However, where the secondary structure was most abundant, ${\textit{EBC}} \le 66\,\%$ EBC ≤ 66 % . During nighttime, the storage terms together with $Q_\mathrm{G}$ Q G made up the largest part of the energy balance, while $Q_\mathrm{H}$ Q H and $Q_\mathrm{E}$ Q E were relatively small. These radiation and energy balance measurements in an insect-attacked forest highlight the role of secondary structure in the recovery of attacked stands.     

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.     

Cross-Validation of Open-Path and Closed-Path Eddy-Covariance Techniques for Observing Methane Fluxes

Methane ( ${\mathrm {CH}}_{4}$ ) fluxes observed with the eddy-covariance technique using an open-path ${\mathrm {CH}}_{4}$ analyzer and a closed-path ${\mathrm {CH}}_{4}$ analyzer in a rice paddy field were evaluated with an emphasis on the flux correction methodology. A comparison of the fluxes obtained by the analyzers revealed that both the open-path and closed-path techniques were reliable, provided that appropriate corrections were applied. For the open-path approach, the influence of fluctuations in air density and the line shape variation in laser absorption spectroscopy (hereafter, spectroscopic effect) was significant, and the relative importance of these corrections would increase when observing small ${\mathrm {CH}}_{4}$ fluxes. A new procedure proposed by Li-Cor Inc. enabled us to accurately adjust for these effects. The high-frequency loss of the open-path ${\mathrm {CH}}_{4}$ analyzer was relatively large (11 % of the uncorrected covariance) at an observation height of 2.5 m above the canopy owing to its longer physical path length, and this correction should be carefully applied before correcting for the influence of fluctuations in air density and the spectroscopic effect. Uncorrected ${\mathrm {CH}}_{4}$ fluxes observed with the closed-path analyzer were substantially underestimated (37 %) due to high-frequency loss because an undersized pump was used in the observation. Both the bandpass and transfer function approaches successfully corrected this flux loss. Careful determination of the bandpass frequency range or the transfer function and the cospectral model is required for the accurate calculation of ${\mathrm {CH}}_{4}$ fluxes with the closed-path technique.     

An Approach to Minimizing Artifacts Caused by Cross-Sensitivity in the Determination of Air–Sea \text{ CO }_{2} Flux Using the Eddy-Covariance Technique

The air–sea $\text{ CO }_{2}$ flux was measured from a research vessel in the North Yellow Sea in October 2007 using an open-path eddy-covariance technique. In 11 out of 64 samples, the normalized spectra of scalars ( $\text{ CO }_{2}$ , water vapour, and temperature) showed similarities. However, in the remaining samples, the normalized $\text{ CO }_{2}$ spectra were observed to be greater than those of water vapour and temperature at low frequencies. In this paper, the noise due to cross-sensitivity was identified through a combination of intercomparisons among the normalized spectra of three scalars and additional analyses. Upon examination, the cross-sensitivity noise appeared to be mainly present at frequencies ${<}0.8\,\text{ Hz }$ . Our analysis also suggested that the high-frequency fluctuations of $\text{ CO }_{2}$ concentration (frequency ${>}0.8\,\text{ Hz }$ ) was probably less affected by the cross-sensitivity. To circumvent the cross-sensitivity issue, the cospectrum in the high-frequency range 0.8–1.5 Hz, instead of the whole range, was used to estimate the $\text{ CO }_{2}$ flux by taking the contribution of the high frequency to the $\text{ CO }_{2}$ flux to be the same as the contribution to the water vapour flux. The estimated air–sea $\text{ CO }_{2}$ flux in the North Yellow Sea was $-0.039\,\pm \,0.048\,\text{ mg } \text{ m }^{-2}\,\text{ s }^{-1},$ a value comparable to the estimates using the inertial dissipation method and Edson’s method (Edson et al., J Geophys Res 116:C00F10, 2011).     

Investigation of the Flow Structure in Step-Up Street Canyons—Mean Flow and Turbulence Statistics

A step-up street canyon is a characteristic urban element composed of two buildings in which the height of the upwind building ( $H_\mathrm{u}$ ) is less than the height of the downwind building ( $H_\mathrm{d}$ ). Here, the effect of canyon geometry on the flow structure in isolated step-up street canyons is investigated through isothermal wind-tunnel measurements. The measurements were acquired along the vertical symmetry plane of model buildings using two-dimensional particle image velocimetry (PIV) for normal approach flow. The building-height ratios considered were: $H_\mathrm{d}/ H_\mathrm{u} \approx 3$ , and $H_\mathrm{d}/ H_\mathrm{u} \approx 1.67$ . For each building-height ratio, the along-wind lengths (L) of the upwind and downwind buildings, and the street-canyon width (S) were kept constant, with $L \approx S$ . The cross-wind widths (W) of the upwind and downwind buildings were varied uniformly from $W/S \approx 1$ through $W/S \approx 4$ , in increments of $W/S \approx 1$ . The objective of the work was to characterize the changes in the flow structure in step-up canyons as a function of W/S, for fixed L, S, and $H_\mathrm{d}/H_\mathrm{u}$ values. The results indicate that the in-canyon flow structure does not vary significantly for $H_\mathrm{d}/H_\mathrm{u} \approx 3$ for the W/S values considered. Qualitatively, for $H_\mathrm{d}/H_\mathrm{u} \approx 3$ , the upwind building behaves as an obstacle in the upwind cavity of the downwind building. In contrast, the flow patterns observed for the $H_\mathrm{d}/H_\mathrm{u} \approx 1.67$ configurations are unique and counter-intuitive, and depend strongly on building width (W/S). For $W/S \approx 1$ and $W/S \approx 2$ , the effect of lateral flow into the canyon is so prominent that even the mean flow patterns are highly ambiguous. For $W/S \approx 3$ and 4, the flow along the vertical symmetry plane is more shielded from the lateral flow, and hence a stable counter-rotating vortex pair is observed in the canyon. In addition to these qualitative features, a quantitative analysis of the mean flow field and turbulence stress field is presented.     

Bulk Mixing and Decoupling of the Nocturnal Stable Boundary Layer Characterized Using a Ubiquitous Natural Tracer

Vertical mixing of the nocturnal stable boundary layer (SBL) over a complex land surface is investigated for a range of stabilities, using a decoupling index ( $0 < D_{rb} < 1$ ) based on the 2–50 m bulk gradient of the ubiquitous natural trace gas radon-222. The relationship between $D_{rb}$ and the bulk Richardson number ( $R_{ib}$ ) exhibits three broad regions: (1) a well-mixed region ( $D_{rb} \approx 0.05$ ) in weakly stable conditions ( $R_{ib} < 0.03$ ); (2) a steeply increasing region ( $0.05 < D_{rb} < 0.9$ ) for “transitional” stabilities ( $0.03 < R_{ib} < 1$ ); and (3) a decoupled region ( $D_{rb} \approx 0.9$ –1.0) in very stable conditions ( $R_{ib} > 1$ ). $D_{rb}$ exhibits a large variability within individual $R_{ib}$ bins, however, due to a range of competing processes influencing bulk mixing under different conditions. To explore these processes in $R_{ib}$ – $D_{rb}$ space, we perform a bivariate analysis of the bulk thermodynamic gradients, various indicators of external influences, and key turbulence quantities at 10 and 50 m. Strong and consistent patterns are found, and five distinct regions in $R_{ib}$ – $D_{rb}$ space are identified and associated with archetypal stable boundary-layer regimes. Results demonstrate that the introduction of a scalar decoupling index yields valuable information about turbulent mixing in the SBL that cannot be gained directly from a single bulk thermodynamic stability parameter. A significant part of the high variability observed in turbulence statistics during very stable conditions is attributable to changes in the degree of decoupling of the SBL from the residual layer above. When examined in $R_{ib}$ – $D_{rb}$ space, it is seen that very different turbulence regimes can occur for the same value of $R_{ib}$ , depending on the particular combination of values for the bulk temperature gradient and wind shear, together with external factors. Extremely low turbulent variances and fluxes are found at 50 m height when $R_{ib} > 1$ and $D_{rb} \approx 1$ (fully decoupled). These “quiescent” cases tend to occur when geostrophic forcing is very weak and subsidence is present, but are not associated with the largest bulk temperature gradients. Humidity and net radiation data indicate the presence of low cloud, patchy fog or dew, any of which may aid decoupling in these cases by preventing temperature gradients from increasing sufficiently to favour gravity wave activity. The largest temperature gradients in our dataset are actually associated with smaller values of the decoupling index ( $D_{rb} < 0.7$ ), indicating the presence of mixing. Strong evidence is seen from enhanced turbulence levels, fluxes and submeso activity at 50 m, as well as high temperature variances and heat flux intermittencies at 10 m, suggesting this region of the $R_{ib}$ – $D_{rb}$ distribution can be identified as a top-down mixing regime. This may indicate an important role for gravity waves and other wave-like phenomena in providing the energy required for sporadic mixing at this complex terrain site.     

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.     

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.     

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.     

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.     

Improvements in Sensible Heat-Flux Parametrization in the High-Resolution Regional Model (HRM) Through the Modified Treatment of the Roughness Length for Heat

We discuss the impact of the differential treatment of the roughness lengths for momentum and heat ( $z_{0\mathrm{m}}$ and $z_{0\mathrm{h}}$ ) in the flux parametrization scheme of the high-resolution regional model (HRM) for a heterogeneous terrain centred around Thiruvananthapuram, India (8.5°N, 76.9°E). The magnitudes of sensible heat flux (H) obtained from HRM simulations using the original parametrization scheme differed drastically from the concurrent in situ observations. With a view to improving the performance of this parametrization scheme, two distinct modifications are incorporated: (1) In the first method, a constant value of 100 is assigned to the $z_{0\mathrm{m}}/z_{0\mathrm{h}}$ ratio; (2) and in the second approach, this ratio is treated as a function of time. Both these modifications in the HRM model showed significant improvements in the H simulations for Thiruvananthapuram and its adjoining regions. Results obtained from the present study provide a first-ever comparison of H simulations using the modified parametrization scheme in the HRM model with in situ observations for the Indian coastal region, and suggest a differential treatment of $z_{0\mathrm{m}}$ and $z_{0\mathrm{h}}$ in the flux parametrization scheme.     

Characteristics and Mechanisms of Low-Level Jets in the Yangtze River Delta of China

A dataset obtained using a wind-profile radar located at the Yangtze River Delta in China ( $31.14^{\circ }$ N, $121.81^{\circ }$ E) in 2009 was used to investigate the characteristics and evolution of low-level jets (LLJs) along the east China coast. The study investigated the daily and seasonal structures of LLJs as well as several possible causes. A total of 1,407 1-h LLJ periods were detected based on an adaptive definition that enabled determination of four LLJ categories. The majority (77 %) of LLJs were found to have speeds $<$ 14.0 m s $^{-1}$ (maximum of 34.6 m s $^{-1})$ and occur at an average altitude below 600 m (76 % of the observed LLJs). The dominant direction of the LLJs was from the south-south-west, which accounted for nearly 32 %, with the second most common wind direction ranging from $040^{\circ }$ to $100^{\circ }$ , albeit with a number of stronger LLJs from the west-south-west. A comparison of LLJs and South-west Jets revealed that the frequencies of occurrence in summer are totally different. Results also revealed that in spring and summer, most LLJs originate from the south-south-west, whereas in autumn and winter, north-east is the dominant direction of origin. The peak heights of LLJs tended to be higher in winter than in other seasons. The horizontal wind speed and peak height of the LLJs displayed pronounced diurnal cycles. The Hilbert–Huang transform technique was applied to demonstrate that the intrinsic mode functions with a cycle of nearly 23 h at levels below 800 m, and the instantaneous amplitudes of inertial events (0.0417–0.0476 h $^{-1}$ frequencies) have large values at 300–600 m. The variations in the occurrences of LLJs suggested connections between the formation mechanisms of LLJs and the South-west Jet stream, steady occupation of synoptic-scale pressure system, and land–sea temperature contrasts.     

One Year of Surface-Based Temperature Inversions at Dome C, Antarctica

In 2005 the Study of Stable Boundary Layer Environment at Dome C (STABLEDC) experimental campaign was conducted at the plateau station of Concordia at Dome C, Antarctica. Temperature profiles measured with a microwave radiometer were used to study the characteristics of surface-based temperature inversions over the course of a year. Statistics of temperature profiles for every month are discussed; the difference between daytime and nocturnal cases observed during the summer months disappears during winter. Surface-based temperature inversions occurred in 70 % of the time during summer, and almost all of the time during winter. During winter the occurrence of warming events leads to a decrease in the temperature difference between the top and the base of the inversion (i.e. the inversion strength). The inversion strength maxima ranged between $3\,^{\circ }\mathrm{C}$ (December) and $35\,^{\circ }\mathrm{C}$ (August) corresponding to gradients of 0.1 and $0.3\,^{\circ }\mathrm{C}\, \mathrm{m}^{-1}$ , respectively. The average surface-based inversion height presents a daily cycle during the summer months with values up to 200 m in the morning hours, while it affects a layer always deeper than 100 m during the winter months. The relationships between inversion strength and the downward longwave radiative flux, absolute temperature, and wind speed are examined. The inversion strength decreases as the longwave radiation increases. A clear anti-correlation between inversion strength and near-surface temperature is evident throughout the year. During the winter, the largest inversion strength values were observed under low wind-speed conditions; in contrast, a clear dependence was not found during the summer.     

Statistics of Turbulence in the Stable Boundary Layer Affected by Gravity Waves

In order to investigate effects of interactions between turbulence and gravity waves in the stable boundary layer on similarity theory relationships, we re-examined a dataset, collected during three April nights in 1978 and in 1980 on the 300-m tower of the Boulder Atmospheric Observatory (BAO). The BAO site, located in Erie, Colorado, USA, 30 km east of the foothills of the Rocky Mountains, has been known for the frequent detection of wave activities. The considered profiles of turbulent fluxes and variances were normalized by two local, gradient-based scaling systems, and subsequently compared with similarity functions of the Richardson number, obtained based on data with no influence of gravity currents and topographical factors. The first scaling system was based on local values of the vertical velocity variance $\sigma _\mathrm{w}$ and the Brunt–Väisäla frequency $ N,$ while the second one was based on the temperature variance $\sigma _{\theta }$ and $N.$ Analysis showed some departures from the similarity functions (obtained for data with virtually no influence of mesoscale motions); nonetheless the overall dependency of dimensionless moments on the Richardson number was maintained.     

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.     

Turbulent Transport of Momentum and Scalars Above an Urban Canopy

Turbulent transport of momentum and scalars over an urban canopy is investigated using the quadrant analysis technique. High-frequency measurements are available at three levels above the urban canopy (47, 140 and 280 m). The characteristics of coherent ejection–sweep motions (flux contributions and time fractions) at the three levels are analyzed, particularly focusing on the difference between ejections and sweeps, the dissimilarity between momentum and scalars, and the dissimilarity between the different scalars (i.e., temperature, water vapour and $\hbox {CO}_{2})$ . It is found that ejections dominate momentum and scalar transfer at all three levels under unstable conditions, while sweeps are the dominant eddy motions for transporting momentum and scalars in the urban roughness sublayer under neutral and stable conditions. The flux contributions and time fractions of ejections and sweeps can be adequately captured by assuming a Gaussian joint probability density function for flow variables. However, the inequality of flux contributions from ejections and sweeps is more accurately reproduced by the third-order cumulant expansion method (CEM). The incomplete cumulant expansion method (ICEM) also works well except for $\hbox {CO}_{2}$ at 47 m where the skewness of $\hbox {CO}_{2}$ fluctuations is significantly larger than that for vertical velocity. The dissimilarity between momentum and scalar transfers is linked to the dissimilarity in the characteristics of ejection–sweep motions and is further quantified by measures of transport efficiencies. Atmospheric stability is the controlling factor for the transport efficiencies of momentum and heat, and fitted functions from the literature describe their behaviour fairly accurately. However, transport efficiencies of water vapour and $\hbox {CO}_{2}$ are less affected by the atmospheric stability. The dissimilarity among the three scalars examined in this study is linked to the active role of temperature and to the surface heterogeneity effect.     

Determination of the onset dates of the rainy season in central Amazon with equatorially antisymmetric outgoing longwave radiation

The present work provides a new methodology to determine onset dates of the rainy season (ONR) in central Amazon (CAM) using the antisymmetric in relation to the equator outgoing longwave radiation (OLR) (AOLR) data, for the 1979–2006 period. Spatial averages of the AOLR ( $\overline {AOLR} $ ) over the CAM for the ONR periods are obtained. These periods correspond to 25 pentads centered on the mean pentad for the ONR. The sign changes from positive to negative of the $\overline {AOLR} $ for the ONR periods indicate the transition from dry to wet season. Composites of several variables are done for pentads before and after the ONR dates. These composites show physically consistent features. The potential of the $\overline {AOLR} $ time series as an index for monitoring tasks is analyzed. The results here show that the $\overline {AOLR} $ for the ONR period captures the transition from dry to wet conditions in the CAM area during 2006. The advantages of this method are discussed. The new simple method proposed here seems to be efficient in determining the ONR in the CAM.