Turbulent Pressure Statistics in the Atmospheric Boundary Layer from Large-Eddy Simulation

We use large-eddy simulation (LES) to study the turbulent pressure field in atmospheric boundary layers with free convection, forced convection, and stable stratification. We use the Poisson equation for pressure to represent the pressure field as the sum of mean-shear, turbulence–turbulence, subfilter-scale, Coriolis, and buoyancy contributions. We isolate these contributions and study them separately. We find that in the energy-containing range in the free-convection case the turbulence–turbulence pressure dominates over the entire boundary layer. That part dominates also up to midlayer in the forced-convection case; above that the mean-shear pressure dominates. In the stable case the mean-shear pressure dominates over the entire boundary layer.We find evidence of an inertial subrange in the pressure spectrum in the free and forced-convection cases; it is dominated by the turbulence–turbulence pressure and has a three-dimensional spectral constant of about 4.0. This agrees well with quasi-Gaussian predictions but is a factor of 2 less than recent results from direct numerical simulations at moderate Reynolds numbers. Measurements of the inertial subrange pressure spectral constant at high Reynolds numbers, which might now be possible, would be most useful.     

Influence of heterogeneous air entry pressure on large scale unsaturated flow in porous media

The paper presents numerical simulations of water infiltration in unsaturated porous media containing coarse-textured inclusions embedded in fine-textured background material. The calculations are performed using the two-phase model for water and air flow and a simplified model known as the Richards equation. It is shown that the Richards equation cannot correctly describe flow in the presence of heterogeneities. However, its performance can be improved by introducing appropriately defined effective capillary and permeability functions, representing largescale behaviour of the heterogeneous medium.     

Simulating the impact of the large-scale circulation on the 2-m temperature and precipitation climatology

The impact of the simulated large-scale atmospheric circulation on the regional climate is examined using the Weather Research and Forecasting (WRF) model as a regional climate model. The purpose is to understand the potential need for interior grid nudging for dynamical downscaling of global climate model (GCM) output for air quality applications under a changing climate. In this study we downscale the NCEP-Department of Energy Atmospheric Model Intercomparison Project (AMIP-II) Reanalysis using three continuous 20-year WRF simulations: one simulation without interior grid nudging and two using different interior grid nudging methods. The biases in 2-m temperature and precipitation for the simulation without interior grid nudging are unreasonably large with respect to the North American Regional Reanalysis (NARR) over the eastern half of the contiguous United States (CONUS) during the summer when air quality concerns are most relevant. This study examines how these differences arise from errors in predicting the large-scale atmospheric circulation. It is demonstrated that the Bermuda high, which strongly influences the regional climate for much of the eastern half of the CONUS during the summer, is poorly simulated without interior grid nudging. In particular, two summers when the Bermuda high was west (1993) and east (2003) of its climatological position are chosen to illustrate problems in the large-scale atmospheric circulation anomalies. For both summers, WRF without interior grid nudging fails to simulate the placement of the upper-level anticyclonic (1993) and cyclonic (2003) circulation anomalies. The displacement of the large-scale circulation impacts the lower atmosphere moisture transport and precipitable water, affecting the convective environment and precipitation. Using interior grid nudging improves the large-scale circulation aloft and moisture transport/precipitable water anomalies, thereby improving the simulated 2-m temperature and precipitation. The results demonstrate that constraining the RCM to the large-scale features in the driving fields improves the overall accuracy of the simulated regional climate, and suggest that in the absence of such a constraint, the RCM will likely misrepresent important large-scale shifts in the atmospheric circulation under a future climate.     

Institutional Challenges for National Groundwater Governance: Policies and Issues

Understanding the issues surrounding groundwater governance is a precondition for developing policy recommendations for both national and transboundary groundwater governance. This review discusses groundwater attributes relevant to the design of governance systems and provides a systematic review of current national groundwater governance differentiated by various policy instruments. The synthesis of both resource system characteristics and experience with policy instruments allows us to conceptualize institutional aspects of groundwater governance. This leads to six institutional aspects: (1) voluntary compliance; (2) tradition and mental models; (3) administrative responsibility and bureaucratic inertia; (4) conflict resolution mechanisms; (5) political economy; and (6) information deficits. Each of these issues embodies institutional challenges for national and international policy implementation.     

Simulation of Meteorological Fields Within and Above Urban and Rural Canopies with a Mesoscale Model

Accurate simulation of air quality at neighbourhood scales (on order of 1-km horizontal grid spacing) requires detailed meteorological fields inside the roughness sub-layer (RSL). Since the assumptions of the roughness approach, used by most of the mesoscale models, are unsatisfactory at this scale, a detailed urban and rural canopy parameterisation, called DA-SM2-U, is developed inside the Penn State/NCAR Mesoscale Model (MM5) to simulate the meteorological fields within and above the urban and rural canopies. DA-SM2-U uses the drag-force approach to represent the dynamic and turbulent effects of the buildings and vegetation, and a modified version of the soil model SM2-U, called SM2-U(3D), to represent the thermodynamic effects of the canopy elements. The turbulence length scale is also modified inside the canopies. SM2-U(3D) assesses the sensible and latent heat fluxes from rural and urban surfaces in each of the computational layers inside the canopies by considering the shadowing effect, the radiative trapping by the street canyons, and the storage heat flux by the artificial surfaces. DA-SM2-U is tested during one simulated day above the city of Philadelphia, U.S.A. It is shown that DA-SM2-U is capable of simulating the important features observed in the urban and rural RSL, as seen in the vertical profiles of the shear stress, turbulent kinetic energy budget components, eddy diffusivity, potential air temperature, and specific humidity. Within the canopies, DA-SM2-U simulates the decrease of the wind speed inside the dense canopies, the skirting of the flow around the canopy blocks, warmer air inside the vegetation canopy than above open areas during the night and conversely during the day, and constantly warmer air inside the urban canopy. The comparison with measurements shows that the surface air temperature above rural and urban areas is better simulated by DA-SM2-U than by the `standard version' of MM5.     

IACT observations of gamma-ray bursts: prospects for the Cherenkov Telescope Array

Gamma rays at rest frame energies as high as 90 GeV have been reported from gamma-ray bursts (GRBs) by the Fermi Large Area Telescope (LAT). There is considerable hope that a confirmed GRB detection will be possible with the upcoming Cherenkov Telescope Array (CTA), which will have a larger effective area and better low-energy sensitivity than current-generation imaging atmospheric Cherenkov telescopes (IACTs). To estimate the likelihood of such a detection, we have developed a phenomenological model for GRB emission between 1 GeV and 1 TeV that is motivated by the high-energy GRB detections of Fermi-LAT, and allows us to extrapolate the statistics of GRBs seen by lower energy instruments such as the Swift-BAT and BATSE on the Compton Gamma-ray Observatory. We show a number of statistics for detected GRBs, and describe how the detectability of GRBs with CTA could vary based on a number of parameters, such as the typical observation delay between the burst onset and the start of ground observations. We also consider the possibility of using GBM on Fermi as a finder of GRBs for rapid ground follow-up. While the uncertainty of GBM localization is problematic, the small field-of-view for IACTs can potentially be overcome by scanning over the GBM error region. Overall, our results indicate that CTA should be able to detect one GRB every 20–30 months with our baseline instrument model, assuming consistently rapid pursuit of GRB alerts, and provided that spectral breaks below ~100 GeV are not a common feature of the bright GRB population. With a more optimistic instrument model, the detection rate can be as high as 1 to 2 GRBs per year.     

Numerical modelling of floodplain hydraulics and suspended sediment transport and deposition at the event scale in the middle river Elbe,Germany

The objective of this case study was to calibrate and verify detailed transport model of sediment in a 4‐kilometre stretch of the middle Elbe floodplains in Germany. The hydraulic RMA‐2 model and the SED2d‐WES sediment transport model were used. These models were calibrated and validated by detailed measurement of the surface water elevations, the velocities at six profiles, and the suspended sediment concentration and deposition (by means of 10 sediment traps). The flow was modelled for three steady‐state discharges. The surface water elevations were calculated to an accuracy of less than 5 cm compared to measurements. The differences between the calculated and measured velocities were with one exception smaller than 0.2 m/s (measured range 0.1…?1.0 m/s). An average sediment input of 35 g/(m² d) was calculated for the flood event studied. The highest calculated sedimentation rates of 700 g/(m² d) (dry density 90 kg/m³) occurred in quiescent zones and abandoned channels. Twenty‐five percent of the deposited sediment settled in the quiescent zones (which only account for 13% of the area). The most sensitive parameters of the sediment transport model were the settling velocity and critical shear stress. The modelling techniques used allowed sediment deposition on the floodplains of the Elbe to be realistically depicted.     

An evaluation of neutral and convective planetary boundary-layer parameterizations relative to large eddy simulations

This paper compares a number of one-dimensional closure models for the planetary boundary layer (PBL) that are currently in use in large-scale atmospheric models. Using the results of a large-eddy simulation (LES) model as the standard of comparison, the PBL models are evaluated over a range of stratifications from free convective to neutral and a range of surface shear stresses. Capping inversion strengths for the convective cases range from weakly to strongly capped. Six prototypical PBL models are evaluated in this study, which focuses on the accuracy of the boundary-layer fluxes of momentum, heat, and two passive scalars. One scalar mimics humidity and the other is a top-down scalar entrained into the boundary layer from above. A set of measures based on the layer-averaged differences of these fluxes from the LES solutions is developed. In addition to the methodological framework and suite of LES solutions, the main result of the evaluation is the recognition that all of the examined PBL parameterizations have difficulty reproducing the entrainment at the top of the PBL, as given by the LES, in most parameter regimes. Some of the PBL models are relatively accurate in their entrainment flux in a subset of parameter regimes. The sensitivity of the PBL models to vertical resolution is explored, and substantive differences are observed in the performance of the PBL models, relative to LES, at low resolution typical of large scale atmospheric models.     

Paleobiological Implications of the Isotopic Signatures (C,N) of Fossil Mammal Collagen in Scladina Cave (Sclayn, Belgium)

An isotopic investigation of upper Pleistocene mammal bones and teeth from Scladina cave (Sclayn, Belgium) demonstrated the very good quality of collagen preservation. A preliminary screening of the samples used the amount of nitrogen in whole bone and dentine in order to estimate the preserved amount of collagen before starting the extraction process. The isotopic abundances of fossil specimens from still-extant species are consistent with their trophic position. Moreover, the¹⁵N isotopic abundance is higher in dentine than in bone in bears and hyenas, a phenomenon already observed in modern specimens. These results demonstrate that the isotopic compositions of samples from Scladina cave can be interpreted in ecological terms. Mammoths exhibit a high¹⁵N isotopic abundance relative to other herbivores, as was the case in Siberian and Alaskan samples. These results suggest distinctive dietary adaptations in herbivores living in the mammoth steppe. Cave bears are clearly isotopically different from coeval brown bears, suggesting an ecological separation between species, with a pure vegetarian diet for cave bear and an omnivorous diet for brown bear.