A contour-advective semi-lagrangian numerical algorithm for simulating fine-scale conservative dynamical fields

This paper describes a novel numerical algorithm for simulating the evolution of fine-scale conservative fields in layer-wise two-dimensional flows, the most important examples of which are the earth's atmosphere and oceans. the algorithm combines two radically different algorithms, one Lagrangian and the other Eulerian, to achieve an unexpected gain in computational efficiency. The algorithm is demonstrated for multi-layer quasi-geostrophic flow, and results are presented for a simulation of a tilted stratospheric polar vortex and of nearly-inviscid quasi-geostrophic turbulence. the turbulence results contradict previous arguments and simulation results that have suggested an ultimate two-dimensional, vertically-coherent character of the flow. Ongoing extensions of the algorithm to the generally ageostrophic flows characteristic of planetary fluid dynamics are outlined.

Fine-scale temporal characterization of trends in soil water dissolved organic carbon and potential drivers

Long-term monitoring of surface water quality has shown increasing concentrations of Dissolved Organic Carbon (DOC) across a large part of the Northern Hemisphere. Several drivers have been implicated including climate change, land management change, nitrogen and sulphur deposition and CO2 enrichment. Analysis of stream water data, supported by evidence from laboratory studies, indicates that an effect of declining sulphur deposition on catchment soil chemistry is likely to be the primary mechanism, but there are relatively few long term soil water chemistry records in the UK with which to investigate this, and other, hypotheses directly. In this paper, we assess temporal relationships between soil solution chemistry and parameters that have been argued to regulate DOC production and, using a unique set of co-located measurements of weather and bulk deposition and soil solution chemistry provided by the UK Environmental Change Network and the Intensive Forest Monitoring Level II Network . We used statistical non-linear trend analysis to investigate these relationships at 5 forested and 4 non-forested sites from 1993 to 2011. Most trends in soil solution DOC concentration were found to be non-linear. Significant increases in DOC occurred mostly prior to 2005. The magnitude and sign of the trends was associated qualitatively with changes in acid deposition, the presence/absence of a forest canopy, soil depth and soil properties. The strongest increases in DOC were seen in acidic forest soils and were most clearly linked to declining anthropogenic acid deposition, while DOC trends at some sites with westerly locations appeared to have been influenced by shorter-term hydrological variation. The results indicate that widespread DOC increases in surface waters observed elsewhere, are most likely dominated by enhanced mobilization of DOC in surficial organic horizons, rather than changes in the soil water chemistry of deeper horizons. While trends in DOC concentrations in surface horizons have flattened out in recent years, further increases may be expected as soil chemistry continues to adjust to declining inputs of acidity.

Predicting mesh density for adaptive modelling of the global atmosphere

The shallow water equations are solved using a mesh of polygons on the sphere, which adapts infrequently to the predicted future solution. Infrequent mesh adaptation reduces the cost of adaptation and load-balancing and will thus allow for more accurate mapping on adaptation. We simulate the growth of a barotropically unstable jet adapting the mesh every 12 h. Using an adaptation criterion based largely on the gradient of the vorticity leads to a mesh with around 20 per cent of the cells of a uniform mesh that gives equivalent results. This is a similar proportion to previous studies of the same test case with mesh adaptation every 1–20 min. The prediction of the mesh density involves solving the shallow water equations on a coarse mesh in advance of the locally refined mesh in order to estimate where features requiring higher resolution will grow, decay or move to. The adaptation criterion consists of two parts: that resolved on the coarse mesh, and that which is not resolved and so is passively advected on the coarse mesh. This combination leads to a balance between resolving features controlled by the large-scale dynamics and maintaining fine-scale features.

Earthworms produce granules of intricately zoned calcite

Earthworms of the family Lumbricidae, which includes many common species, produce and secrete up to millimeter-sized calcite granules, and the intricate fine-scale zoning of their constituent crystals is unique for a biomineral. Granule calcite is produced by crystallization of amorphous calcium carbonate (ACC) that initially precipitates within the earthworm calciferous glands, then forms protogranules by accretion on quartz grain cores. Crystallization of ACC is mediated by migrating fluid films and is largely complete within 24 11 of ACC production and before granules leave the earthworm. Variations in the density of defects formed as a byproduct of trace element incorporation during calcite crystall growth have generated zoning that can be resolved by cathodoluminescence imaging at ultraviolet to blue wavelengths and using the novel technique of scanning electron microscope charge contrast imaging. Mapping of calcite crystal orientations by electron backscatter diffraction reveals an approximate radial fabric to the granules that reflects crystal growth from internal nucleation sites toward their margins. The survival within granules of ACC inclusions for months after they enter soils indicates that they crystallize only within the earthworm and in the presence of fluids containing biochemical catalysts. The earthworm probably promotes crystallization of ACC in order to prevent remobilization of the calcium carbonate by dissolution. Calcite granules vividly illustrate the role of transient precursors in biomineralization, but the underlying question of why earth-worms produce granules in volumes sufficient to have a measurable impact on soil carbon cycling remains to be answered.

Biogenic nitric oxide emissions from soils: impact on NOx and ozone over West Africa during AMMA (African Monsoon Multidisciplinary Analysis) - observational study

Chemical and meteorological parameters measured on board the Facility for Airborne Atmospheric Measurements (FAAM) BAe 146 Atmospheric Research Aircraft during the African Monsoon Multidisciplinary Analysis (AMMA) campaign are presented to show the impact of NOx emissions from recently wetted soils in West Africa. NO emissions from soils have been previously observed in many geographical areas with different types of soil/vegetation cover during small scale studies and have been inferred at large scales from satellite measurements of NOx. This study is the first dedicated to showing the emissions of NOx at an intermediate scale between local surface sites and continental satellite measurements. The measurements reveal pronounced mesoscale variations in NOx concentrations closely linked to spatial patterns of antecedent rainfall. Fluxes required to maintain the NOx concentrations observed by the BAe-146 in a number of cases studies and for a range of assumed OH concentrations (1×106 to 1×107 molecules cm−3) are calculated to be in the range 8.4 to 36.1 ng N m−2 s−1. These values are comparable to the range of fluxes from 0.5 to 28 ng N m−2 s−1 reported from small scale field studies in a variety of non-nutrient rich tropical and sub-tropical locations reported in the review of Davidson and Kingerlee (1997). The fluxes calculated in the present study have been scaled up to cover the area of the Sahel bounded by 10 to 20 N and 10 E to 20 W giving an estimated emission of 0.03 to 0.30 Tg N from this area for July and August 2006. The observed chemical data also suggest that the NOx emitted from soils is taking part in ozone formation as ozone concentrations exhibit similar fine scale structure to the NOx, with enhancements over the wet soils. Such variability can not be explained on the basis of transport from other areas. Delon et al. (2008) is a companion paper to this one which models the impact of soil NOx emissions on the NOx and ozone concentration over West Africa during AMMA. It employs an artificial neural network to define the emissions of NOx from soils, integrated into a coupled chemistry-dynamics model. The results are compared to the observed data presented in this paper. Here we compare fluxes deduced from the observed data with the model-derived values from Delon et al. (2008).

Regional variation in cortical bone properties from broiler fowl - a first look

1. We used microhardness testing as a probe for fine-scale regional variation in the mechanical performance of bone and present data showing the extent of regional variation in the femora and humeri of 7-week-old broiler birds. 2. Ash content of dry bone was broadly correlated with microhardness, although there is some evidence that the relationship linking the two differs between the femur and the humerus. 3. Regional variations in the properties of bone from poultry are widely overlooked in the literature. Awareness of them is vital and existing measures of bone 'strength' may be misleading if local variation in properties is not taken into account when exploring the effects of nutrition and husbandry practices on bone mechanical performance.

Constraint release in entangled binary blends of linear polymers: a molecular dynamics study

We present extensive molecular dynamics simulations of the dynamics of diluted long probe chains entangled with a matrix of shorter chains. The chain lengths of both components are above the entanglement strand length, and the ratio of their lengths is varied over a wide range to cover the crossover from the chain reptation regime to tube Rouse motion regime of the long probe chains. Reducing the matrix chain length results in a faster decay of the dynamic structure factor of the probe chains, in good agreement with recent neutron spin echo experiments. The diffusion of the long chains, measured by the mean square displacements of the monomers and the centers of mass of the chains, demonstrates a systematic speed-up relative to the pure reptation behavior expected for monodisperse melts of sufficiently long polymers. On the other hand, the diffusion of the matrix chains is only weakly perturbed by the diluted long probe chains. The simulation results are qualitatively consistent with the theoretical predictions based on constraint release Rouse model, but a detailed comparison reveals the existence of a broad distribution of the disentanglement rates, which is partly confirmed by an analysis of the packing and diffusion of the matrix chains in the tube region of the probe chains. A coarse-grained simulation model based on the tube Rouse motion model with incorporation of the probability distribution of the tube segment jump rates is developed and shows results qualitatively consistent with the fine scale molecular dynamics simulations. However, we observe a breakdown in the tube Rouse model when the short chain length is decreased to around N-S = 80, which is roughly 3.5 times the entanglement spacing N-e(P) = 23. The location of this transition may be sensitive to the chain bending potential used in our simulations.

Hyperresolution global land surface modeling: meeting a grand challenge for monitoring Earth's terrestrial water

Monitoring Earth's terrestrial water conditions is critically important to many hydrological applications such as global food production; assessing water resources sustainability; and flood, drought, and climate change prediction. These needs have motivated the development of pilot monitoring and prediction systems for terrestrial hydrologic and vegetative states, but to date only at the rather coarse spatial resolutions (∼10–100 km) over continental to global domains. Adequately addressing critical water cycle science questions and applications requires systems that are implemented globally at much higher resolutions, on the order of 1 km, resolutions referred to as hyperresolution in the context of global land surface models. This opinion paper sets forth the needs and benefits for a system that would monitor and predict the Earth's terrestrial water, energy, and biogeochemical cycles. We discuss six major challenges in developing a system: improved representation of surface‐subsurface interactions due to fine‐scale topography and vegetation; improved representation of land‐atmospheric interactions and resulting spatial information on soil moisture and evapotranspiration; inclusion of water quality as part of the biogeochemical cycle; representation of human impacts from water management; utilizing massively parallel computer systems and recent computational advances in solving hyperresolution models that will have up to 109 unknowns; and developing the required in situ and remote sensing global data sets. We deem the development of a global hyperresolution model for monitoring the terrestrial water, energy, and biogeochemical cycles a “grand challenge” to the community, and we call upon the international hydrologic community and the hydrological science support infrastructure to endorse the effort.

Reactions to peripheral image motion using a head/eye platform

The authors demonstrate four real-time reactive responses to movement in everyday scenes using an active head/eye platform. They first describe the design and realization of a high-bandwidth four-degree-of-freedom head/eye platform and visual feedback loop for the exploration of motion processing within active vision. The vision system divides processing into two scales and two broad functions. At a coarse, quasi-peripheral scale, detection and segmentation of new motion occurs across the whole image, and at fine scale, tracking of already detected motion takes place within a foveal region. Several simple coarse scale motion sensors which run concurrently at 25 Hz with latencies around 100 ms are detailed. The use of these sensors are discussed to drive the following real-time responses: (1) head/eye saccades to moving regions of interest; (2) a panic response to looming motion; (3) an opto-kinetic response to continuous motion across the image and (4) smooth pursuit of a moving target using motion alone.

Absolute high spectral resolution measurements of surface solar radiation for detection of water vapour continuum absorption

Solar-pointing Fourier transform infrared (FTIR) spectroscopy offers the capability to measure both the fine scale and broadband spectral structure of atmospheric transmission simultaneously across wide spectral regions. It is therefore suited to the study of both water vapour monomer and continuum absorption behaviours. However, in order to properly address this issue, it is necessary to radiatively calibrate the FTIR instrument response. A solar-pointing high-resolution FTIR spectrometer was deployed as part of the ‘Continuum Absorption by Visible and Infrared radiation and its Atmospheric Relevance’ (CAVIAR) consortium project. This paper describes the radiative calibration process using an ultra-high-temperature blackbody and the consideration of the related influence factors. The result is a radiatively calibrated measurement of the solar irradiation at the ground across the IR region from 2000 to 10 000 cm−1 with an uncertainty of between 3.3 and 5.9 per cent. This measurement is shown to be in good general agreement with a radiative-transfer model. The results from the CAVIAR field measurements are being used in ongoing studies of atmospheric absorbers, in particular the water vapour continuum.

The integrated WRF/urban modelling system: development, evaluation, and applications to urban environmental problems

To bridge the gaps between traditional mesoscale modelling and microscale modelling, the National Center for Atmospheric Research, in collaboration with other agencies and research groups, has developed an integrated urban modelling system coupled to the weather research and forecasting (WRF) model as a community tool to address urban environmental issues. The core of this WRF/urban modelling system consists of the following: (1) three methods with different degrees of freedom to parameterize urban surface processes, ranging from a simple bulk parameterization to a sophisticated multi-layer urban canopy model with an indoor–outdoor exchange sub-model that directly interacts with the atmospheric boundary layer, (2) coupling to fine-scale computational fluid dynamic Reynolds-averaged Navier–Stokes and Large-Eddy simulation models for transport and dispersion (T&D) applications, (3) procedures to incorporate high-resolution urban land use, building morphology, and anthropogenic heating data using the National Urban Database and Access Portal Tool (NUDAPT), and (4) an urbanized high-resolution land data assimilation system. This paper provides an overview of this modelling system; addresses the daunting challenges of initializing the coupled WRF/urban model and of specifying the potentially vast number of parameters required to execute the WRF/urban model; explores the model sensitivity to these urban parameters; and evaluates the ability of WRF/urban to capture urban heat islands, complex boundary-layer structures aloft, and urban plume T&D for several major metropolitan regions. Recent applications of this modelling system illustrate its promising utility, as a regional climate-modelling tool, to investigate impacts of future urbanization on regional meteorological conditions and on air quality under future climate change scenarios. Copyright © 2010 Royal Meteorological Society

Agriculture and the life histories of Mauritius kestrels

Habitat modification for agriculture is one of the greatest current threats to global biodiversity. Studies show large-scale population declines and short-term demographic impacts, but knowledge of the long-term effects of agriculture on individuals remains poor. This thesis examines the short- and long-term impact of agriculture on a reintroduced population of the Mauritius kestrel Falco punctatus, a tropical forest-dwelling raptor endemic to the island of Mauritius, that also utilises agricultural habitats. This population is a particularly appropriate model system, because complete life history data exists for individuals over a 22-year period, alongside detailed habitat and climate data. Agriculture has a short-term detrimental effect on Mauritius kestrel breeding success by exacerbating the seasonal decline in fledgling production. This is partly driven by the habitat-specific composition of the prey community that kestrels exploit to feed their chicks. The fledglings from agricultural territories tend to recruit in agricultural territories. This is largely due to poor natal dispersal and fine-scale spatial autocorrelation in the habitat matrix. Breeders do not respond to agriculture in the breeding territory by dispersing, unless the pair bond is broken. Therefore, individuals originating in agricultural territories tend to recruit, and remain in, agricultural territories throughout their lives. In addition to this, females from agricultural natal territories have shorter lifespans, schedule their peak reproductive output earlier in life, and exhibit more rapid senescence than non-agricultural females. The combination of this long-term effect and the adult experience of agriculture imposed by life history and environmental constraints, leads to a lower mean lifetime reproductive rate compared to females originating in non-agricultural habitats. These results demonstrate that agriculture experienced in early life has a lifelong effect on individuals. The effects can persist in time and space, with potentially delayed effects on population dynamics. These findings are important for understanding species’ responses to agricultural expansion.

Detection of coherent airstreams using cluster analysis: application to an extratropical cyclone

Flow in geophysical fluids is commonly summarized by coherent streams, for example conveyor belt flows in extratropical cyclones or jet streaks in the upper troposphere. Typically, parcel trajectories are calculated from the flow field and subjective thresholds are used to distinguish coherent streams of interest. This methodology contribution develops a more objective approach to distinguish coherent airstreams within extratropical cyclones. Agglomerative clustering is applied to trajectories along with a method to identify the optimal number of cluster classes. The methodology is applied to trajectories associated with the low-level jets of a well-studied extratropical cyclone. For computational efficiency, a constraint that trajectories must pass through these jet regions is applied prior to clustering; the partitioning into different airstreams is then performed by the agglomerative clustering. It is demonstrated that the methodology can identify the salient flow structures of cyclones: the warm and cold conveyor belts. A test focusing on the airstreams terminating at the tip of the bent-back front further demonstrates the success of the method in that it can distinguish fine-scale flow structure such as descending sting jet airstreams.

Spatial structuring of arbuscular mycorrhizal communities in benchmark and modified temperate eucalypt woodlands

Arbuscular mycorrhizal fungi (AMF) are crucial to the functioning of the plant–soil system, but little is known about the spatial structuring of AMF communities across landscapes modified by agriculture. AMF community composition was characterized across four sites in the highly cleared south-western Australian wheatbelt that were originally dominated by forb-rich eucalypt woodlands. Environmentally induced spatial structuring in AMF composition was examined at four scales: the regional scale associated with location, the site scale associated with past management (benchmark woodlands with no agricultural management history, livestock grazing, recent revegetation), the patch scale associated with trees and canopy gaps, and the fine scale associated with the herbaceous plant species beneath which soils were sourced. Field-collected soils were cultured in trap pots; then, AMF composition was determined by identifying spores and through ITS1 sequencing. Structuring was strongest at site scales, where composition was strongly related to prior management and associated changes in soil phosphorus. The two fields were dominated by the genera Funneliformis and Paraglomus, with little convergence back to woodland composition after revegetation. The two benchmark woodlands were characterized by Ambispora gerdemannii and taxa from Gigasporaceae. Their AMF communities were strongly structured at patch scales associated with trees and gaps, in turn most strongly related to soil N. By contrast, there were few patterns at fine scales related to different herbaceous plant species, or at regional scales associated with the 175 km distance between benchmark woodlands. Important areas for future investigation are to identify the circumstances in which recolonization by woodland AMF may be limited by fungal propagule availability, reduced plant diversity and/or altered chemistry in agricultural soils.