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.     

Diagnostic Analysis of the Quasi-Balanced Flow of a Mesoscale Vortex During the 12 June 2008 Guangxi Rainstorm

By using the high-resolution observation data and MM5 model simulation data, the analysis on the 12 June 2008 Guangxi flash-flood rainstorm shows that the associated major mesoscale weather system of this event is a quasi-stationary mesoscale vortex, which resulted from the interaction between the midlatitude synoptic-scale waves in the westerly belt and the low-latitude warm-moist flow under the terrain effect. The genesis, development, and movement of the mesoscale vortex have significant impacts on the intensity and persistence of the severe precipitation from the Guangxi flash-flood rainstorm. This vortex is characterized by the coexistence of strong vorticity and divergence with the same order of magnitude. Well organized, deep, and moist convection was observed for a long period of time, and was produced by the interaction between the mesoscale vortex and the gravity waves. The latter was generated by the terrain effect and the ageostrophic effect of high winds in the low-level jet. According to the quasi-balanced dynamical theory, quasi-balanced flow must have existed in the mesoscale motions with both divergent and rotational winds. Thus, based on the diagnosis of the quasi-balanced flow, the PV-ωinversion method is employed to analyze the organized moist convection. The results show that 50%-70% of the vertical circulation in the rainstorm areas was quasi-balanced, so the quasi-balanced flow could well reflect features of the strong vertical motions associated with the coexistence of vorticity and divergence during this event.     

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.     

Early Warning Signals for Regime Transition in the Stable Boundary Layer: A Model Study

The evening transition is investigated in an idealized model for the nocturnal boundary layer. From earlier studies it is known that the nocturnal boundary layer may manifest itself in two distinct regimes, depending on the ambient synoptic conditions: strong-wind or overcast conditions typically lead to weakly stable, turbulent nights; clear-sky and weak-wind conditions, on the other hand, lead to very stable, weakly turbulent conditions. Previously, the dynamical behaviour near the transition between these regimes was investigated in an idealized setting, relying on Monin–Obukhov (MO) similarity to describe turbulent transport. Here, we investigate a similar set-up, using direct numerical simulation; in contrast to MO-based models, this type of simulation does not need to rely on turbulence closure assumptions. We show that previous predictions are verified, but now independent of turbulence parametrizations. Also, it appears that a regime shift to the very stable state is signaled in advance by specific changes in the dynamics of the turbulent boundary layer. Here, we show how these changes may be used to infer a quantitative estimate of the transition point from the weakly stable boundary layer to the very stable boundary layer. In addition, it is shown that the idealized, nocturnal boundary-layer system shares important similarities with generic non-linear dynamical systems that exhibit critical transitions. Therefore, the presence of other, generic early warning signals is tested as well. Indeed, indications are found that such signals are present in stably stratified turbulent flows.     

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.     

Fractional Flow Speed-Up from Porous Windbreaks for Enhanced Wind-Turbine Power

The potential for porous windbreaks to enhance wind-turbine power production is studied using linearized theory and wind-tunnel experiments. Results suggest that windbreaks have the potential to substantially increase power production, while lowering mean shear, and leading to negligible changes in turbulence intensity. The fractional increase in turbine power output is found to vary roughly linearly with windbreak height, where a windbreak 10% the height of the turbine hub increases power by around 10%. Wind-tunnel experiments with a windbreak imposed beneath a turbulent boundary layer show the linearized predictions to be in good agreement with particle-image-velocimetry data. Power measurements from a model turbine further corroborate predictions in power increase. Moreover, the wake of the windbreak showed a significant interaction with the turbine wake, which may inform windbreak use in large wind farms. Power measurements from a second turbine downwind of the first with its own windbreak show that the net effect for multiple turbines is dependent on windbreak height.     

Influence of the Surf Zone on the Marine Aerosol Concentration in a Coastal Area

Sea-salt aerosol concentrations in the coastal zone are assessed with the numerical aerosol-transport model MACMod that applies separate aerosol source functions for open ocean and the surf zone near the sea–land transition. Numerical simulations of the aerosol concentration as a function of offshore distance from the surf zone compare favourably with experimental data obtained during a surf-zone aerosol experiment in Duck, North Carolina in autumn 2007. Based on numerical simulations, the effect of variations in aerosol production (source strength) and transport conditions (wind speed, air–sea temperature difference), we show that the surf-zone aerosols are replaced by aerosols generated over the open ocean as the airmass advects out to sea. The contribution from the surf-generated aerosol is significant during high wind speeds and high wave events, and is significant up to 30 km away from the production zone. At low wind speeds, the oceanic component dominates, except within 1–5 km of the surf zone. Similar results are obtained for onshore flow, where no further sea-salt aerosol production occurs as the airmass advects out over land. The oceanic aerosols that are well-mixed throughout the boundary layer are then more efficiently transported inland than are the surf-generated aerosols, which are confined to the first few tens of metres above the surface, and are therefore also more susceptible to the type of surface (trees or grass) that determines the deposition velocity.