Nutrients cause grassland biomass to outpace herbivory.

Human activities are transforming grassland biomass via changing climate, elemental nutrients, and herbivory. Theory predicts that food-limited herbivores will consume any additional biomass stimulated by nutrient inputs ('consumer-controlled'). Alternatively, nutrient supply is predicted to increase biomass where herbivores alter community composition or are limited by factors other than food ('resource-controlled'). Using an experiment replicated in 58 grasslands spanning six continents, we show that nutrient addition and vertebrate herbivore exclusion each caused sustained increases in aboveground live biomass over a decade, but consumer control was weak. However, at sites with high vertebrate grazing intensity or domestic livestock, herbivores consumed the additional fertilization-induced biomass, supporting the consumer-controlled prediction. Herbivores most effectively reduced the additional live biomass at sites with low precipitation or high ambient soil nitrogen. Overall, these experimental results suggest that grassland biomass will outstrip wild herbivore control as human activities increase elemental nutrient supply, with widespread consequences for grazing and fire risk.

Propagation of acoustic waves through saturated porous media.

The homogenization approach to wave propagation through saturated porous media is extended in order to include the compressibility of the interstitial fluid and the existence of several connected pore components which may or not percolate. The necessary theoretical developments are summarized and the Christoffel equation whose solutions provide the wave velocities is presented. Some analytical developments are proposed for isotropic media. Finally, a systematic application to a synthetic porous medium illustrates the methodology and its results.

Soil net nitrogen mineralisation across global grasslands.

Soil nitrogen mineralisation (Nmin), the conversion of organic into inorganic N, is important for productivity and nutrient cycling. The balance between mineralisation and immobilisation (net Nmin) varies with soil properties and climate. However, because most global-scale assessments of net Nmin are laboratory-based, its regulation under field-conditions and implications for real-world soil functioning remain uncertain. Here, we explore the drivers of realised (field) and potential (laboratory) soil net Nmin across 30 grasslands worldwide. We find that realised Nmin is largely explained by temperature of the wettest quarter, microbial biomass, clay content and bulk density. Potential Nmin only weakly correlates with realised Nmin, but contributes to explain realised net Nmin when combined with soil and climatic variables. We provide novel insights of global realised soil net Nmin and show that potential soil net Nmin data available in the literature could be parameterised with soil and climate data to better predict realised Nmin.

Sensitivity of global soil carbon stocks to combined nutrient enrichment.

Soil stores approximately twice as much carbon as the atmosphere and fluctuations in the size of the soil carbon pool directly influence climate conditions. We used the Nutrient Network global change experiment to examine how anthropogenic nutrient enrichment might influence grassland soil carbon storage at a global scale. In isolation, enrichment of nitrogen and phosphorous had minimal impacts on soil carbon storage. However, when these nutrients were added in combination with potassium and micronutrients, soil carbon stocks changed considerably, with an average increase of 0.04 KgCm-2  year-1 (standard deviation 0.18 KgCm-2  year-1 ). These effects did not correlate with changes in primary productivity, suggesting that soil carbon decomposition may have been restricted. Although nutrient enrichment caused soil carbon gains most dry, sandy regions, considerable absolute losses of soil carbon may occur in high-latitude regions that store the majority of the world's soil carbon. These mechanistic insights into the sensitivity of grassland carbon stocks to nutrient enrichment can facilitate biochemical modelling efforts to project carbon cycling under future climate scenarios.

Fast Laplace solver approach to pore-scale permeability.

We introduce a powerful and easily implemented method to calculate the permeability of porous media at the pore scale using an approximation based on the Poiseulle equation to calculate permeability to fluid flow with a Laplace solver. The method consists of calculating the Euclidean distance map of the fluid phase to assign local conductivities and lends itself naturally to the treatment of multiscale problems. We compare with analytical solutions as well as experimental measurements and lattice Boltzmann calculations of permeability for Fontainebleau sandstone. The solver is significantly more stable than the lattice Boltzmann approach, uses less memory, and is significantly faster. Permeabilities are in excellent agreement over a wide range of porosities.

Thermal convection in three-dimensional fractured porous media.

Thermal convection is numerically computed in three-dimensional (3D) fluid saturated isotropically fractured porous media. Fractures are randomly inserted as two-dimensional (2D) convex polygons. Flow is governed by Darcy's 2D and 3D laws in the fractures and in the porous medium, respectively; exchanges take place between these two structures. Results for unfractured porous media are in agreement with known theoretical predictions. The influence of parameters such as the fracture aperture (or fracture transmissivity) and the fracture density on the heat released by the whole system is studied for Rayleigh numbers up to 150 in cubic boxes with closed-top conditions. Then, fractured media are compared to homogeneous porous media with the same macroscopic properties. Three major results could be derived from this study. The behavior of the system, in terms of heat release, is determined as a function of fracture density and fracture transmissivity. First, the increase in the output flux with fracture density is linear over the range of fracture density tested. Second, the increase in output flux as a function of fracture transmissivity shows the importance of percolation. Third, results show that the effective approach is not always valid, and that the mismatch between the full calculations and the effective medium approach depends on the fracture density in a crucial way.

Percolation in three-dimensional fracture networks for arbitrary size and shape distributions.

The percolation threshold of fracture networks is investigated by extensive direct numerical simulations. The fractures are randomly located and oriented in three-dimensional space. A very wide range of regular, irregular, and random fracture shapes is considered, in monodisperse or polydisperse networks containing fractures with different shapes and/or sizes. The results are rationalized in terms of a dimensionless density. A simple model involving a new shape factor is proposed, which accounts very efficiently for the influence of the fracture shape. It applies with very good accuracy in monodisperse or moderately polydisperse networks, and provides a good first estimation in other situations. A polydispersity index is shown to control the need for a correction, and the corrective term is modelled for the investigated size distributions.

Three-scale analysis of the permeability of a natural shale.

The macroscopic permeability of a natural shale is determined by using structural measurements on three different scales. Transmission electron microscopy yields two-dimensional (2D) images with pixels smaller than 1 nm; these images are used to reconstruct 3D nanostructures. Three-dimensional focused ion beam-scanning electron microscopy (5.95- to 8.48-nm voxel size) provides 3D mesoscale pores of limited relative volume (1.71-5.9%). Micro-computed tomography (700-nm voxel size) provides information on the mineralogy of the shale, including the pores on this scale which do not percolate; synthetic 3D media are derived on the macroscopic scale by a training image technique. Permeability of the nanoscale, of the mesoscale structures and of their superposition is determined by solving the Stokes equation and this enables us to estimate the permeabilities of the 700-nm voxels located within the clay matrix. Finally, the Darcy equation is solved on synthetic 3D macroscale media to obtain the macroscopic permeability which is found in good agreement with experimental results obtained on the centimetric scale.

Two-scale analysis of a tight gas sandstone.

Tight gas sandstones are low porosity media, with a very small permeability (i.e., below 1 mD). Their porosity is below 10%, and it is mainly composed of fine noncemented microcracks, which are present between neighboring quartz grains. While empirical models of permeability are available, their predictions, which do not compare well with macroscopic measurements, are not reliable to assess gas well productivity. The purpose of this work is to compare the permeability measured on centimetric plugs to predictions based on pore structure data. Two macroscopic measurements are performed, namely dry gas permeability and mercury intrusion porosimetry (MIP), together with a series of local measurements including focused ion beam and scanning electron microscopy (FIB-SEM), x-ray computed microtomography (CMT), and standard two-dimensional (2D) SEM. Numerical modeling is performed by combining analyses on two scales, namely the microcrack network scale (given by 2D SEM) and the individual 3D microcrack scale (given by either FIB-SEM or CMT). The network permeability is calculated by means of techniques developed for fracture networks. This permeability is proportional to the microcrack transmissivity, which is determined by solving the Stokes equation in the microcracks measured by FIB-SEM or CMT. Good correlation with experimental permeability values is only found when using transmissivity from 3D CMT data.

Phonon-Modulated Magnetic Interactions and Spin Tomonaga-Luttinger Liquid in the p-Orbital Antiferromagnet CsO2.

The magnetic response of antiferromagnetic CsO2, coming from the p-orbital S=1/2 spins of anionic O2(-) molecules, is followed by 133Cs nuclear magnetic resonance across the structural phase transition occurring at T(s1)=61 K on cooling. Above T(s1), where spins form a square magnetic lattice, we observe a huge, nonmonotonic temperature dependence of the exchange coupling originating from thermal librations of O2(-) molecules. Below T(s1), where antiferromagnetic spin chains are formed as a result of p-orbital ordering, we observe a spin Tomonaga-Luttinger-liquid behavior of spin dynamics. These two interesting phenomena, which provide rare simple manifestations of the coupling between spin, lattice, and orbital degrees of freedom, establish CsO2 as a model system for molecular solids.

A multi-locus approach resolves the phylogenetic relationships of the Simulium asakoae and Simulium ceylonicum species groups in Malaysia: evidence for distinct evolutionary lineages.

A multi-locus approach was used to examine the DNA sequences of 10 nominal species of blackfly in the Simulium subgenus Gomphostilbia (Diptera: Simuliidae) in Malaysia. Molecular data were acquired from partial DNA sequences of the mitochondria-encoded cytochrome c oxidase subunit I (COI), 12S rRNA and 16S rRNA genes, and the nuclear-encoded 18S rRNA and 28S rRNA genes. No single gene, nor the concatenated gene set, resolved all species or all relationships. However, all morphologically established species were supported by at least one gene. The multi-locus sequence analysis revealed two distinct evolutionary lineages, conforming to the morphotaxonomically recognized Simulium asakoae and Simulium ceylonicum species groups.

Cryptic biodiversity in the cytogenome of bird-biting blackflies in North Africa.

Bird-biting blackflies in the Simulium (Eusimulium) aureum group (Diptera: Simuliidae) are widespread vectors of Leucocytozoon and Trypanosoma parasites. The polytene chromosomes of 619 larvae of the three nominal members of the S. aureum group in North Africa were evaluated cytogenetically for cryptic biodiversity. Seven chromosomal segregates were discovered among 29 populations in Algeria and Morocco. This diversity was based primarily on two chromosomal inversions, which have assumed unique roles in different lineages, including sex linkage, fixation, loss and autosomal polymorphism. Reproductive isolation was demonstrated for six of the seven segregates, doubling the number of species known in the area. Four species were linked with existing names: (a) Simulium mellah Giudicelli & Bouzidi, which is known only from North African high-salinity habitats; (b) Simulium petricolum (Rivosecchi), which is tentatively conspecific with continental European populations; (c) Simulium rubzovianum (Sherban) and its synonym Simulium latinum (Rubtsov), which is widely distributed from North Africa across Europe into Western Asia, and (d) Simulium velutinum (Santos Abreu) and its new synonym Simulium tenerificum Crosskey, which is restricted to North Africa and the Canary Islands. Of the remaining entities, two are new species precinctive to North Africa and one, known only from Morocco, is of undetermined taxonomic status.

Transport properties of a Bentheim sandstone under deformation.

The mechanical and transport properties of a Bentheim sandstone are studied both experimentally and numerically. Three classical classes of loads are applied to a sample whose permeability is measured. The elasticity and the Stokes equations are discretized on unstructured tetrahedral meshes which precisely follow the deformations of the sample. Numerical results are presented, discussed, and compared to the available experimental data.

k=0 magnetic structure and absence of ferroelectricity in SmFeO3.

SmFeO3 has attracted considerable attention very recently due to its reported multiferroic properties above room temperature. We have performed powder and single crystal neutron diffraction as well as complementary polarization dependent soft X-ray absorption spectroscopy measurements on floating-zone grown SmFeO3 single crystals in order to determine its magnetic structure. We found a k=0 G-type collinear antiferromagnetic structure that is not compatible with inverse Dzyaloshinskii-Moriya interaction driven ferroelectricity. While the structural data reveal a clear sign for magneto-elastic coupling at the Néel-temperature of ∼675 K, the dielectric measurements remain silent as far as ferroelectricity is concerned.

Application of level-set method for deposition in three-dimensional reconstructed porous media.

Three-dimensional (3D) porous structures which are usually discretized by voxels can also be discretized by the level-set method (LSM) and flow, reactive transport, and structure evolution can be modeled. The determination of the solid-liquid interface is detailed as well as the discretization of the governing equations. Comparisons between a one-dimensional analytical solution and LSM are conducted for validation. Deposition in 3D reconstructed media is studied under various flow and reaction conditions. The evolution of the structure is explored locally by means of the pore geometry and globally by means of the permeability and the porosity. The difference between the voxel method and LSM is discussed during the investigations.

Solid matrix partition by fracture networks.

The geometrical properties of the matrix blocks formed by a random fracture network are investigated numerically, for a wide range of fracture shapes and for fracture densities ranging from the dilute limit to well above the threshold where the material is entirely partitioned into finite blocks. The main block characteristics are the density and volume fraction, the mean volume and surface area, and their number of faces. In the dilute limit, general expressions for these characteristics are obtained, which provide a good approximation of the numerical data for any fracture shape. In the dense regime, most properties are governed by power laws, which involve two fitted exponents independent of the fracture shape. The shape factors identified in the dilute limit remain relevant for dense networks and can be used to formulate a general model for the block characteristics, valid up to the total matrix fracturation. The transition density when this occurs is determined. It can also be used to account for the fracture shape effects in a very simple and fairly accurate general model. Beyond the transition density, the block characteristics converge as expected toward those in the space tesselation by infinite planes.

Molecular systematics of the Simulium jenningsi species group (Diptera: Simuliidae), with three new fast-evolving nuclear genes for phylogenetic inference.

A molecular phylogeny was inferred for the 22 nominal species of black flies in the Simulium jenningsi species group, which includes major pests of humans and livestock in North America. Females are structurally monomorphic, presenting a problem for identification of the pests. For each species, we sequenced approximately two kilobases from the mitochondrial genome (ND2, Cox I, proximal one-half of Cox II) and about six kilobases from the nuclear genome (ca. 2 kilobases each from 3 rapidly evolving nuclear genes: big zinc finger [BZF], "5-intron gene" [5intG], and elongation complex protein 1 [ECP1]) and analyzed them phylogenetically using maximum likelihood and Bayesian methods. The three nuclear loci have not previously been used in phylogenetic studies. The mitochondrial region recovered 6 group members as monophyletic. BZF, 5intG, and ECP1 sequences each permitted identification of 13 species and recovered the S. fibrinflatum and S. taxodium subgroups. Simulium aranti Stone and Snoddy and S. luggeri Nicholson and Mickel were consistently recovered at the base of the group. Simulium ozarkense Moulton and Adler, S. dixiense Stone and Snoddy, S. krebsorum Moulton and Adler, and S. haysi Stone and Snoddy branched off before two well-supported sister groups of the remaining species. This remainder consisted of species occupying slow, sandy lowland streams-S. definitum Moulton and Adler, S. jonesi Stone and Snoddy, and the S. taxodium subgroup (S. taxodium Snoddy and Beshear, S. chlorum Moulton and Adler, S. confusum Moulton and Adler, and S. lakei Snoddy)-as sister to two clades of species inhabiting swift, rocky upland streams-the S. fibrinflatum subgroup (S. fibrinflatum Twinn, S. notiale Stone and Snoddy, and S. snowi Stone and Snoddy) and a clade comprised of S. anchistinum Moulton and Adler, S. jenningsi Malloch, and S. nyssa Stone and Snoddy, plus species having cocoons without anterolateral apertures (S. infenestrum Moulton and Adler, S. podostemi Snoddy, S. penobscotense Snoddy and Bauer, and S. remissum Moulton and Adler). Simulium snowi Stone and Snoddy is here considered a synonym of S. notiale Stone and Snoddy. Trees inferred from BZF and 5intG were largely concordant with those from ECP1, but slightly less resolved. Combining mitochondrial and nuclear data sets did not greatly improve the performance of the ECP1 data set alone. We, therefore, propose ECP1 as the gold standard for identification of members of the S. jenningsi group. Maximum likelihood analysis of combined sequences from all three nuclear genes, with three morphological constraints imposed, yielded a tree proposed as the best hypothesis of relationships among group members, based on all available data.

Interaction between a fracture network and a cubic cavity.

The intersection between a network of polygonal fractures and a cubic cavity is numerically studied. Several probabilities are defined and particular attention is paid to the probabilities of intersection or not of the percolating cluster with the cavity; they depend on the size of the domain, on the fracture density, and on the relative size of the fractures and of the cavity. These probabilities are extrapolated to infinite domains. Analytical approximations are proposed which are in good agreement with the numerical data for sufficiently large densities. Some extensions of practical and theoretical interests are given in the concluding remarks.

Microscale simulation and numerical upscaling of a reactive flow in a plane channel.

A bimolecular homogeneous irreversible reaction of the kind A+B→C is simulated in a plane channel as a base example of reactive transport processes taking place at the microscale within porous and/or fractured media. The numerical study explores the way microscale processes embedded in dimensionless quantities such as Péclet (Pe) and Damköhler (Da) numbers propagate to upscaled coefficients describing effective system dynamics. The microscale evolution of the reactant concentrations is obtained through a particle-based numerical method which has been specifically tailored to the considered problem. Key results include a complete documentation of the process evolution for a wide range of Pe and Da, in terms of the global reaction rate, space-time distribution of reactants, and local mixing features leading to characterization of effective reaction and dispersion coefficients governing a section-averaged upscaled model of the system. The robustness of previously presented theoretical analyses concerning closures of volume-averaged (upscaled) formulations is assessed. The work elucidates the dependence of the effective dispersion and reactive parameters on the microscale mixing and reactive species evolution. Our results identify the role played by Da and Pe on the occurrence of incomplete mixing of reactants, which affects the features of the reactive transport scenario.