Frequency of Boundary-Layer-Top Fluctuations in Convective and Stable Conditions Using Laser Remote Sensing

The planetary boundary-layer (PBL) height is determined with high temporal and altitude resolution from lidar backscatter profiles. Then, the frequencies of daytime thermal updrafts and downdrafts and of nighttime gravity waves are obtained applying a fast Fourier transform on the temporal fluctuation of the PBL height. The principal frequency components of each spectrum are related to the dominant processes occurring at the daytime and nighttime PBL top. Two groups of cases are selected for the study: one group combines daytime cases, measured in weak horizontal wind conditions and dominated by convection. The cases show higher updraft and downdraft frequencies for the shallow, convective boundary layer and lower frequencies for a deep PBL. For cases characterized by strong horizontal winds, the frequencies directly depend on the wind speed. The temporal variation of the PBL height is determined also in the likely presence of lee waves. For nighttime cases, the main frequency components in the spectra do not show a real correlation with the nocturnal PBL height. Altitude fluctuations of the top of the nocturnal boundary layer are observed even though the boundary layer is statically stable. These oscillations are associated with the wind shear effect and with buoyancy waves at the PBL top.     

Estimating Variances of Horizontal Wind Fluctuations in Stable Conditions

Information concerning the average wind speed and the variances of lateral and longitudinal wind velocity fluctuations is required by dispersion models to characterise turbulence in the atmospheric boundary layer. When the winds are weak, the scalar average wind speed and the vector average wind speed need to be clearly distinguished and both lateral and longitudinal wind velocity fluctuations assume equal importance in dispersion calculations. We examine commonly-used methods of estimating these variances from wind-speed and wind-direction statistics measured separately, for example, by a cup anemometer and a wind vane, and evaluate the implied relationship between the scalar and vector wind speeds, using measurements taken under low-wind stable conditions. We highlight several inconsistencies inherent in the existing formulations and show that the widely-used assumption that the lateral velocity variance is equal to the longitudinal velocity variance is not necessarily true. We derive improved relations for the two variances, and although data under stable stratification are considered for comparison, our analysis is applicable more generally.     

Variability of ozone in the marine boundary layer of the equatorial Pacific Ocean

This study examines the processes controlling the diurnal variability of ozone (O₃) in the marine boundary layer of the Kwajalein Atoll, Republic of the Marshall Islands (latitude 8° 43′ N, longitude 167° 44′ E), during July to September 1999. At the study site, situated in the equatorial Pacific Ocean, O₃ mixing ratios remained low, with an overall average of 9–10 parts per billion on a volume basis (ppbv) and a standard deviation of 2.5 ppbv. In the absence of convective storms, daily O₃ mixing ratios decreased after sunrise and reached minimum during the afternoon in response to photochemical reactions. The peak-to-peak amplitude of O₃ diurnal variation was approximately 1–3 ppbv. During the daytime, O₃ photolysis, hydroperoxyl radicals, hydroxyl radicals, and bromine atoms contributed to the destruction of O₃, which explained the observed minimum O₃ levels observed in the afternoon. The entrainment of O₃-richer air from the free troposphere to the local marine boundary layer provided a recovery mechanism of surface O₃ mixing ratio with a transport rate of 0.04 to 0.2 ppbv per hour during nighttime. In the presence of convection, downward transport of O₃-richer tropospheric air increased surface O₃ mixing ratios by 3–12 ppbv. The magnitude of O₃ increase due to moist convection was lower than that observed over the continent (as high as 20–30 ppbv). Differences were ascribed to the higher O₃ levels in the continental troposphere and weaker convection over the ocean. Present results suggest that moist convection plays a role in surface-level O₃ dynamics in the tropical marine boundary layer.     

Political orientation moderates Americans’ beliefs and concern about climate change

I offer some theoretical insights to help us better understand the moderator effect of political orientation that Larry Hamilton and others have found in recent years. Reflexive modernization theory highlights an emerging tension between those who direct attention to the negative consequences of industrial capitalism such as climate change (e.g., the scientific community and environmental organizations) and those who defend the economic system against such critiques (e.g., the conservative movement). Political divisions in the American public increasingly map onto these societal divisions between critics and defenders of the industrial capitalist order—especially for the issue of climate change. This alignment is facilitated by increased polarization among political elites and balkanization of the news media. Strong evidence of the moderator effect is consistent with the expectations of information processing theory and elite cues hypothesis from political science. Recent empirical findings in political psychology and neuroscience also seem pertinent for explaining this moderator effect. I end by outlining a few implications for climate change research and communication.     

Climate change links fate of glaciers and an endemic alpine invertebrate

Climate warming in the mid- to high-latitudes and high-elevation mountainous regions is occurring more rapidly than anywhere else on Earth, causing extensive loss of glaciers and snowpack. However, little is known about the effects of climate change on alpine stream biota, especially invertebrates. Here, we show a strong linkage between regional climate change and the fundamental niche of a rare aquatic invertebrate—the meltwater stonefly Lednia tumana—endemic to Waterton-Glacier International Peace Park, Canada and USA. L. tumana has been petitioned for listing under the U.S. Endangered Species Act due to climate-change-induced glacier loss, yet little is known on specifically how climate impacts may threaten this rare species and many other enigmatic alpine aquatic species worldwide. During 14 years of research, we documented that L. tumana inhabits a narrow distribution, restricted to short sections (~500 m) of cold, alpine streams directly below glaciers, permanent snowfields, and springs. Our simulation models suggest that climate change threatens the potential future distribution of these sensitive habitats and the persistence of L. tumana through the loss of glaciers and snowfields. Mountaintop aquatic invertebrates are ideal early warning indicators of climate warming in mountain ecosystems. Research on alpine invertebrates is urgently needed to avoid extinctions and ecosystem change.     

How well do integrated assessment models simulate climate change?

Integrated assessment models (IAMs) are regularly used to evaluate different policies of future emissions reductions. Since the global costs associated with these policies are immense, it is vital that the uncertainties in IAMs are quantified and understood. We first demonstrate the significant spread in the climate system and carbon cycle components of several contemporary IAMs. We then examine these components in more detail to understand the causes of differences, comparing the results with more complex climate models and earth system models (ESMs), where available. Our results show that in most cases the outcomes of IAMs are within the range of the outcomes of complex models, but differences are large enough to matter for policy advice. There are areas where IAMs would benefit from improvements (e.g. climate sensitivity, inertia in climate response, carbon cycle feedbacks). In some cases, additional climate model experiments are needed to be able to tune some of these improvements. This will require better communication between the IAM and ESM development communities.     

Relation of the Second Type Thermal Helicity to Precipitation of Landfalling Typhoons: A Case Study of Typhoon Talim

This study utilized the MM5 mesoscale model to simulate the landfalling process of Typhoon Talim. The simulated typhoon track, weather patterns, and rainfall process are consistent with the observation. Using the simulation results, the relation of the second type thermal helicity (H ₂) to rainfall caused by the landfalling typhoon Talim was analyzed. The results show that H ₂ could well indicate the heavy inland rainfall but it did not perform as well as the helicity in predicting rainfall during the beginning stage of the typhoon landfall. In particular, H ₂ was highly correlated with rainfall of Talim at 1-h lead time. For 1–5-h lead time, it also had a higher correlation with rainfall than the helicity did, and thus showing a better potential in forecasting rainfall intensification. Further analyses have shown that when Talim was in the beginning stage of landfall, 1) the 850–200-hPa vertical wind shear around the Talim center was quite small (about 5 m s⁻¹); 2) the highest rainfall was to the right of the Talim track and in the area with a 300-km radius around the Talim center, exhibiting no obvious relation to low-level temperature advection, low-level air convergence, and upper-level divergence; 3) the low-level relative vorticity reflected the rainfall change quite well, which was the main reason why helicity had a better performance than H ₂ in this period. However, after Talim moved inland further, 1) it weakened gradually and was increasingly affected by the northern trough; 2) the vertical wind shear was enhanced as well; 3) the left side of the down vertical wind shear lay in the Lushan and Dabieshan mountain area, which could have contributed to triggering a secondary vertical circulation, helping to produce the heavy rainfall over there; hence, H ₂ showed a better capacity to reflect the rainfall change during this stage.     

Numerical Considerations for Lagrangian Stochastic Dispersion Models: Eliminating Rogue Trajectories,and the Importance of Numerical Accuracy

When Lagrangian stochastic models for turbulent dispersion are applied to complex atmospheric flows, some type of ad hoc intervention is almost always necessary to eliminate unphysical behaviour in the numerical solution. Here we discuss numerical strategies for solving the non-linear Langevin-based particle velocity evolution equation that eliminate such unphysical behaviour in both Reynolds-averaged and large-eddy simulation applications. Extremely large or ‘rogue’ particle velocities are caused when the numerical integration scheme becomes unstable. Such instabilities can be eliminated by using a sufficiently small integration timestep, or in cases where the required timestep is unrealistically small, an unconditionally stable implicit integration scheme can be used. When the generalized anisotropic turbulence model is used, it is critical that the input velocity covariance tensor be realizable, otherwise unphysical behaviour can become problematic regardless of the integration scheme or size of the timestep. A method is presented to ensure realizability, and thus eliminate such behaviour. It was also found that the numerical accuracy of the integration scheme determined the degree to which the second law of thermodynamics or ‘well-mixed condition’ was satisfied. Perhaps more importantly, it also determined the degree to which modelled Eulerian particle velocity statistics matched the specified Eulerian distributions (which is the ultimate goal of the numerical solution). It is recommended that future models be verified by not only checking the well-mixed condition, but perhaps more importantly by checking that computed Eulerian statistics match the Eulerian statistics specified as inputs.     

Conditionally Averaged Large-Scale Motions in the Neutral Atmospheric Boundary Layer: Insights for Aeolian Processes

Aeolian erosion of flat, arid landscapes is induced (and sustained) by the aerodynamic surface stress imposed by flow in the atmospheric surface layer. Conceptual models typically indicate that sediment mass flux, Q (via saltation or drift), scales with imposed aerodynamic stress raised to some exponent, n, where \(n > 1\). This scaling demonstrates the importance of turbulent fluctuations in driving aeolian processes. In order to illustrate the importance of surface-stress intermittency in aeolian processes, and to elucidate the role of turbulence, conditional averaging predicated on aerodynamic surface stress has been used within large-eddy simulation of atmospheric boundary-layer flow over an arid, flat landscape. The conditional-sampling thresholds are defined based on probability distribution functions of surface stress. The simulations have been performed for a computational domain with \(\approx 25 H\) streamwise extent, where H is the prescribed depth of the neutrally-stratified boundary layer. Thus, the full hierarchy of spatial scales are captured, from surface-layer turbulence to large- and very-large-scale outer-layer coherent motions. Spectrograms are used to support this argument, and also to illustrate how turbulent energy is distributed across wavelengths with elevation. Conditional averaging provides an ensemble-mean visualization of flow structures responsible for erosion ‘events’. Results indicate that surface-stress peaks are associated with the passage of inclined, high-momentum regions flanked by adjacent low-momentum regions. Fluid in the interfacial shear layers between these adjacent quasi-uniform momentum regions exhibits high streamwise and vertical vorticity.     

The Impact of Land-Surface Parameter Properties and Resolution on the Simulated Cloud-Topped Atmospheric Boundary Layer

A method to simulate characteristics of wind speed in the boundary layer of tropical cyclones in an idealized manner is developed and evaluated. The method can be used in a single-column modelling set-up with a planetary boundary-layer parametrization, or within large-eddy simulations (LES). The key step is to include terms in the horizontal velocity equations representing advection and centrifugal acceleration in tropical cyclones that occurs on scales larger than the domain size. Compared to other recently developed methods, which require two input parameters (a reference wind speed, and radius from the centre of a tropical cyclone) this new method also requires a third input parameter: the radial gradient of reference wind speed. With the new method, simulated wind profiles are similar to composite profiles from dropsonde observations; in contrast, a classic Ekman-type method tends to overpredict inflow-layer depth and magnitude, and two recently developed methods for tropical cyclone environments tend to overpredict near-surface wind speed. When used in LES, the new technique produces vertical profiles of total turbulent stress and estimated eddy viscosity that are similar to values determined from low-level aircraft flights in tropical cyclones. Temporal spectra from LES produce an inertial subrange for frequencies \(\gtrsim \)0.1 Hz, but only when the horizontal grid spacing \(\lesssim \)20 m.