The REFLEX project: Comparing different algorithms and implementations for the inversion of a terrestrial ecosystem model against eddy covariance data

We describe a model-data fusion (MDF) inter-comparison project (REFLEX), which compared various algorithms for estimating carbon (C) model parameters consistent with both measured carbon fluxes and states and a simple C model. Participants were provided with the model and with both synthetic net ecosystem exchange (NEE) of CO2 and leaf area index (LAI) data, generated from the model with added noise, and observed NEE and LAI data from two eddy covariance sites. Participants endeavoured to estimate model parameters and states consistent with the model for all cases over the two years for which data were provided, and generate predictions for one additional year without observations. Nine participants contributed results using Metropolis algorithms, Kalman filters and a genetic algorithm. For the synthetic data case, parameter estimates compared well with the true values. The results of the analyses indicated that parameters linked directly to gross primary production (GPP) and ecosystem respiration, such as those related to foliage allocation and turnover, or temperature sensitivity of heterotrophic respiration, were best constrained and characterised. Poorly estimated parameters were those related to the allocation to and turnover of fine root/wood pools. Estimates of confidence intervals varied among algorithms, but several algorithms successfully located the true values of annual fluxes from synthetic experiments within relatively narrow 90% confidence intervals, achieving >80% success rate and mean NEE confidence intervals <110 gC m−2 year−1 for the synthetic case. Annual C flux estimates generated by participants generally agreed with gap-filling approaches using half-hourly data. The estimation of ecosystem respiration and GPP through MDF agreed well with outputs from partitioning studies using half-hourly data. Confidence limits on annual NEE increased by an average of 88% in the prediction year compared to the previous year, when data were available. Confidence intervals on annual NEE increased by 30% when observed data were used instead of synthetic data, reflecting and quantifying the addition of model error. Finally, our analyses indicated that incorporating additional constraints, using data on C pools (wood, soil and fine roots) would help to reduce uncertainties for model parameters poorly served by eddy covariance data.

Occurrence of Poecilochirus austroasiaticus (Acari: Parasitidae) in forensic autopsies and its application on postmortem interval estimation

Despite the fact that mites were used at the dawn of forensic entomology to elucidate the postmortem interval, their use in current cases remains quite low for procedural reasons such as inadequate taxonomic knowledge. A special interest is focused on the phoretic stages of some mite species, because the phoront-host specificity allows us to deduce in many occasions the presence of the carrier (usually Diptera or Coleoptera) although it has not been seen in the sampling performed in situ or in the autopsy room. In this article, we describe two cases where Poecilochirus austroasiaticus Vitzthum (Acari: Parasitidae) was sampled in the autopsy room. In the first case, we could sample the host, Thanatophilus ruficornis (Küster) (Coleoptera: Silphidae), which was still carrying phoretic stages of the mite on the body. That attachment allowed, by observing starvation/feeding periods as a function of the digestive tract filling, the establishment of chronological cycles of phoretic behavior, showing maximum peaks of phoronts during arrival and departure from the corpse and the lowest values in the phase of host feeding. From the sarcosaprophagous fauna, we were able to determine in this case a minimum postmortem interval of 10 days. In the second case, we found no Silphidae at the place where the corpse was found or at the autopsy, but a postmortem interval of 13 days could be established by the high specificity of this interspecific relationship and the departure from the corpse of this family of Coleoptera.

Accounting for the economic risk, caused by variation in disease severity, in fungicide dose decisions: exemplified for Mycosphaerella graminicola on winter wheat

A method is presented to calculate economic optimum fungicide doses accounting for the risk-aversion of growers responding to variability in disease severity between crops. Simple dose-response and disease-yield loss functions are used to estimate net disease-related costs (fungicide cost, plus disease-induced yield loss) as a function of dose and untreated severity. With fairly general assumptions about the shapes of the probability distribution of disease severity and the other functions involved, we show that a choice of fungicide dose which minimises net costs on average across seasons results in occasional large net costs caused by inadequate control in high disease seasons. This may be unacceptable to a grower with limited capital. A risk-averse grower can choose to reduce the size and frequency of such losses by applying a higher dose as insurance. For example, a grower may decide to accept ‘high loss’ years one year in ten or one year in twenty (i.e. specifying a proportion of years in which disease severity and net costs will be above a specified level). Our analysis shows that taking into account disease severity variation and risk-aversion will usually increase the dose applied by an economically rational grower. The analysis is illustrated with data on septoria tritici leaf blotch of wheat caused by Mycosphaerella graminicola. Observations from untreated field plots at sites across England over three years were used to estimate the probability distribution of disease severities at mid-grain filling. In the absence of a fully reliable disease forecasting scheme, reducing the frequency of ‘high loss’ years requires substantially higher doses to be applied to all crops. Disease resistant cultivars reduce both the optimal dose at all levels of risk and the disease-related costs at all doses.

Chapter 4: Simulation of reaction injection moulding

Reaction Injection Moulding is a technology that enables the rapid production of complex plastic parts directly from a mixture of two reactive materials of low viscosity. The reactants are mixed in specific quantities and injected into a mould. This process allows large complex parts to be produced without the need for high clamping pressures. This chapter explores the simulation of the complex processes involved in reaction injection moulding. The reaction processes mean that the dynamics of the material in the mould are in constant evolution and an effective model which takes full account of these changing dynamics is introduced and incorporated in to finite element procedures, which are able to provide a complete simulation of the cycle of mould filling and subsequent curing.

Non-ergodicity of inviscid two-dimensional flow on a beta-plane and on the surface of a rotating sphere

It is shown that, for a sufficiently large value of β, two-dimensional flow on a doubly-periodic beta-plane cannot be ergodic (phase-space filling) on the phase-space surface of constant energy and enstrophy. A corresponding result holds for flow on the surface of a rotating sphere, for a sufficiently rapid rotation rate Ω. This implies that the higher-order, non-quadratic invariants are exerting a significant influence on the statistical evolution of the flow. The proof relies on the existence of a finite-amplitude Liapunov stability theorem for zonally symmetric basic states with a non-vanishing absolute-vorticity gradient. When the domain size is much larger than the size of a typical eddy, then a sufficient condition for non-ergodicity is that the wave steepness ε < 1, where ε = 2[surd radical]2Z/βU in the planar case and $\epsilon = 2^{\frac{1}{4}} a^{\frac{5}{2}}Z^{\frac{7}{4}}/\Omega U^{\frac{5}{2}}$ in the spherical case, and where Z is the enstrophy, U the r.m.s. velocity, and a the radius of the sphere. This result may help to explain why numerical simulations of unforced beta-plane turbulence (in which ε decreases in time) seem to evolve into a non-ergodic regime at large scales.

Tracing troposphere-to-stratosphere transport above a mid-latitude deep convective system

Within the project SPURT (trace gas measurements in the tropopause region) a variety of trace gases have been measured in situ in order to investigate the role of dynamical and chemical processes in the extra-tropical tropopause region. In this paper we report on a flight on 10 November 2001 leading from Hohn, Germany (52�N) to Faro, Portugal (37�N) through a strongly developed deep stratospheric intrusion. This streamer was associated with a large convective system over the western Mediterranean with potentially significant troposphere-to-stratosphere transport. Along major parts of the flight we measured unexpectedly high NOy mixing ratios. Also H2O mixing ratios were significantly higher than stratospheric background levels confirming the extraordinary chemical signature of the probed air masses in the interior of the streamer. Backward trajectories encompassing the streamer enable to analyze the origin and physical characteristics of the air masses and to trace troposphere-to-stratosphere transport. Near the western flank of the streamer features caused by long range transport, such as tropospheric filaments characterized by sudden drops in the O3 and NOy mixing ratios and enhanced CO and H2O can be reconstructed in great detail using the reverse domain filling technique. These filaments indicate a high potential for subsequent mixing with the stratospheric air. At the south-western edge of the streamer a strong gradient in the NOy and the O3 mixing ratios coincides very well with a sharp gradient in potential vorticity in the ECMWF fields. In contrast, in the interior of the streamer the observed highly elevated NOy and H2O mixing ratios up to a potential temperature level of 365K and potential vorticity values of maximum 10 PVU cannot be explained in terms of resolved troposphere-to-stratosphere transport along the backward trajectories. Also mesoscale simulations with a High Resolution Model reveal no direct evidence for convective H2O injection up to this level. Elevated H2O mixing ratios in the ECMWF and HRM are seen only up to about tropopause height at 340 hPa and 270 hPa, respectively, well below flight altitude of about 200 hPa. However, forward tracing of the convective influence as identified by satellite brightness temperature measurements and counts of lightning strokes shows that during this part of the flight the aircraft was closely following the border of an air mass which was heavily impacted by convective activity over Spain and Algeria. This is evidence that deep convection at mid-latitudes may have a large impact on the tracer distribution of the lowermost stratosphere reaching well above the thunderstorms anvils as claimed by recent studies using cloud-resolving models.

Toward vicarious calibration of microwave remote-sensing satellites in arid environments

The Soil Moisture and Ocean Salinity (SMOS) satellite marks the commencement of dedicated global surface soil moisture missions, and the first mission to make passive microwave observations at L-band. On-orbit calibration is an essential part of the instrument calibration strategy, but on-board beam-filling targets are not practical for such large apertures. Therefore, areas to serve as vicarious calibration targets need to be identified. Such sites can only be identified through field experiments including both in situ and airborne measurements. For this purpose, two field experiments were performed in central Australia. Three areas are studied as follows: 1) Lake Eyre, a typically dry salt lake; 2) Wirrangula Hill, with sparse vegetation and a dense cover of surface rock; and 3) Simpson Desert, characterized by dry sand dunes. Of those sites, only Wirrangula Hill and the Simpson Desert are found to be potentially suitable targets, as they have a spatial variation in brightness temperatures of <4 K under normal conditions. However, some limitations are observed for the Simpson Desert, where a bias of 15 K in vertical and 20 K in horizontal polarization exists between model predictions and observations, suggesting a lack of understanding of the underlying physics in this environment. Subsequent comparison with model predictions indicates a SMOS bias of 5 K in vertical and 11 K in horizontal polarization, and an unbiased root mean square difference of 10 K in both polarizations for Wirrangula Hill. Most importantly, the SMOS observations show that the brightness temperature evolution is dominated by regular seasonal patterns and that precipitation events have only little impact.

Pesticide-handling practices of smallholder coffee farmers in Eastern Jamaica

Pesticide use among smallholder coffee producers in Jamaica has been associated with significant occupational health effects. Research on pesticide handling practices, however, has been scarce, especially in eastern Jamaica. This explorative study aims at filling this gap and provides a first basis to develop effective interventions to promote a safer pesticide use. A random sample of 81 coffee farmers was surveyed. The majority of farmers reported to suffer from at least one health symptom associated with pesticide handling, but safety practices were scarcely adopted. There was also the risk that other household members and the wider local community are exposed to pesticides. The lack of training on pesticide management, the role of health services and the cost for protective equipment seemed to be the most significant factors that influence current pesticide handling practices in eastern Jamaica. Further research is recommended to develop a systemic understanding of farmer’s behaviour to provide a more solid basis for the development of future intervention programmes.

Hydrology and geomorphology of the Upper White Nile lakes and their relevance for water resources management in the Nile basin

Remote sensing data and digital elevation models were utilized to extract the catchment hydrological parameters and to delineate storage areas for the Ugandan Equatorial Lakes region. Available rainfall/discharge data are integrated with these morphometric data to construct a hydrological model that simulates the water balance of the different interconnected basins and enables the impact of potential management options to be examined. The total annual discharges of the basins are generally very low (less than 7% of the total annual rainfall). The basin of the shallow (5 m deep) Lake Kioga makes only a minor hydrological contribution compared with other Equatorial Lakes, because most of the overflow from Lake Victoria basin into Lake Kioga is lost by evaporation and evapotranspiration. The discharge from Lake Kioga could be significantly increased by draining the swamps through dredging and deepening certain channel reaches. Development of hydropower dams on the Equatorial Lakes will have an adverse impact on the annual water discharge downstream, including the occasional reduction of flow required for filling up to designed storage capacities and permanently increasing the surface areas of water that is exposed to evaporation. On the basis of modelling studies, alternative sites are proposed for hydropower development and water storage schemes

Concurrent recordings of bladder afferents from multiple nerves using a microfabricated PDMS microchannel electrode array

In this paper we present a compliant neural interface designed to record bladder afferent activity. We developed the implant's microfabrication process using multiple layers of silicone rubber and thin metal so that a gold microelectrode array is embedded within four parallel polydimethylsiloxane (PDMS) microchannels (5 mm long, 100 μm wide, 100 μm deep). Electrode impedance at 1 kHz was optimized using a reactive ion etching (RIE) step, which increased the porosity of the electrode surface. The electrodes did not deteriorate after a 3 month immersion in phosphate buffered saline (PBS) at 37 °C. Due to the unique microscopic topography of the metal film on PDMS, the electrodes are extremely compliant and can withstand handling during implantation (twisting and bending) without electrical failure. The device was transplanted acutely to anaesthetized rats, and strands of the dorsal branch of roots L6 and S1 were surgically teased and inserted in three microchannels under saline immersion to allow for simultaneous in vivo recordings in an acute setting. We utilized a tripole electrode configuration to maintain background noise low and improve the signal to noise ratio. The device could distinguish two types of afferent nerve activity related to increasing bladder filling and contraction. To our knowledge, this is the first report of multichannel recordings of bladder afferent activity.

Coupled subdaily and multiweek cycles during the lava dome erupt of Soufriere Hills Volcano, Montserrat

Observations of volcanoes extruding andesitic lava to produce lava domes often reveal cyclic behaviour. At Soufriere Hills Volcano, Montserrat, cycles with sub-daily and multi-week periods have been recognised on many occasions. These two types of cycle have been modelled separately as stick-slip magma flow at the junction between a dyke and an overlying cylindrical conduit (Costa et al. 2012), and as the filling and discharge of magma through the elastic-walled dyke (Costa et al., 2007a) respectively. Here, we couple these two models to simulate the behaviour over a period of well-observed multi-week cycles, with accompanying sub-daily cycles, from 13 May to 21 September 1997. The coupled model captures well the asymmetrical first-order behaviour: the first 40% of the multi-week cycle consists of high rates of lava extrusion during short period/high amplitude sub-daily cycles as the dyke reservoir discharges itself. The remainder of the cycle involves increasing pressurization as more magma is stored, and extrusion rate falls, followed by a gradual increase in the period of the sub-daily cycles.

A gradual depth-dependent change of connectivity features of supragranular pyramidal cells in rat barrel cortex

Recent experimental evidence suggests a finer genetic, structural and functional subdivision of the layers which form a cortical column. The classical layer II/III (LII/III) of rodent neocortex integrates ascending sensory information with contextual cortical information for behavioral read-out. We systematically investigated to which extent regular-spiking supragranular pyramidal neurons, located at different depths within the cortex, show different input-output connectivity patterns. Combining glutamate-uncaging with whole-cell recordings and biocytin filling, we revealed a novel cellular organization of LII/III: (i) “Lower LII/III” pyramidal cells receive a very strong excitatory input from lemniscal LIV and much fewer inputs from paralemniscal LVa. They project to all layers of the home column, including a feedback projection to LIV whereas transcolumnar projections are relatively sparse. (ii) “Upper LII/III” pyramidal cells also receive their strongest input from LIV, but in addition, a very strong and dense excitatory input from LVa. They project extensively to LII/III as well as LVa and Vb of their home and neighboring columns, (iii) “Middle LII/III” pyramidal cell show an intermediate connectivity phenotype that stands in many ways in-between the features described for lower versus upper LII/III. “Lower LII/III” intracolumnarly segregates and transcolumnarly integrates lemniscal information whereas “upper LII/III” seems to integrate lemniscal with paralemniscal information. This suggests a finegrained functional subdivision of the supragranular compartment containing multiple circuits without any obvious cytoarchitectonic, other structural or functional correlate of a laminar border in rodent barrel cortex.

Patterns of maximum body size evolution in Cenozoic land mammals: eco-evolutionary processes and abiotic forcing

There is accumulating evidence that macroevolutionary patterns of mammal evolution during the Cenozoic follow similar trajectories on different continents. This would suggest that such patterns are strongly determined by global abiotic factors, such as climate, or by basic eco-evolutionary processes such as filling of niches by specialization. The similarity of pattern would be expected to extend to the history of individual clades. Here, we investigate the temporal distribution of maximum size observed within individual orders globally and on separate continents. While the maximum size of individual orders of large land mammals show differences and comprise several families, the times at which orders reach their maximum size over time show strong congruence, peaking in the Middle Eocene, the Oligocene and the Plio-Pleistocene. The Eocene peak occurs when global temperature and land mammal diversity are high and is best explained as a result of niche expansion rather than abiotic forcing. Since the Eocene, there is a significant correlation between maximum size frequency and global temperature proxy. The Oligocene peak is not statistically significant and may in part be due to sampling issues. The peak in the Plio-Pleistocene occurs when global temperature and land mammal diversity are low, it is statistically the most robust one and it is best explained by global cooling. We conclude that the macroevolutionary patterns observed are a result of the interplay between eco-evolutionary processes and abiotic forcing

Investigation of discrepancies in satellite rainfall estimates over Ethiopia

Tropical Applications of Meteorology Using Satellite and Ground-Based Observations (TAMSAT) rainfall estimates are used extensively across Africa for operational rainfall monitoring and food security applications; thus, regional evaluations of TAMSAT are essential to ensure its reliability. This study assesses the performance of TAMSAT rainfall estimates, along with the African Rainfall Climatology (ARC), version 2; the Tropical Rainfall Measuring Mission (TRMM) 3B42 product; and the Climate Prediction Center morphing technique (CMORPH), against a dense rain gauge network over a mountainous region of Ethiopia. Overall, TAMSAT exhibits good skill in detecting rainy events but underestimates rainfall amount, while ARC underestimates both rainfall amount and rainy event frequency. Meanwhile, TRMM consistently performs best in detecting rainy events and capturing the mean rainfall and seasonal variability, while CMORPH tends to overdetect rainy events. Moreover, the mean difference in daily rainfall between the products and rain gauges shows increasing underestimation with increasing elevation. However, the distribution in satellite–gauge differences demon- strates that although 75% of retrievals underestimate rainfall, up to 25% overestimate rainfall over all eleva- tions. Case studies using high-resolution simulations suggest underestimation in the satellite algorithms is likely due to shallow convection with warm cloud-top temperatures in addition to beam-filling effects in microwave- based retrievals from localized convective cells. The overestimation by IR-based algorithms is attributed to nonraining cirrus with cold cloud-top temperatures. These results stress the importance of understanding re- gional precipitation systems causing uncertainties in satellite rainfall estimates with a view toward using this knowledge to improve rainfall algorithms.

First-principles study of the inversion thermodynamics and electronic structure of FeM2X4 (thio)spinels (M = Cr, Mn, Co, Ni; X = O, S)

FeM2X4 spinels, where M is a transition metal and X is oxygen or sulfur, are candidate materials for spin filters, one of the key devices in spintronics. We present here a computational study of the inversion thermodynamics and the electronic structure of these (thio)spinels for M = Cr, Mn, Co, Ni, using calculations based on the density functional theory with on-site Hubbard corrections (DFT+U). The analysis of the configurational free energies shows that different behaviour is expected for the equilibrium cation distributions in these structures: FeCr2X4 and FeMn2S4 are fully normal, FeNi2X4 and FeCo2S4 are intermediate, and FeCo2O4 and FeMn2O4 are fully inverted. We have analyzed the role played by the size of the ions and by the crystal field stabilization effects in determining the equilibrium inversion degree. We also discuss how the electronic and magnetic structure of these spinels is modified by the degree of inversion, assuming that this could be varied from the equilibrium value. We have obtained electronic densities of states for the completely normal and completely inverse cation distribution of each compound. FeCr2X4, FeMn2X4, FeCo2O4 and FeNi2O4 are half-metals in the ferrimagnetic state when Fe is in tetrahedral positions. When M is filling the tetrahedral positions, the Cr-containing compounds and FeMn2O4 are half-metallic systems, while the Co and Ni spinels are insulators. The Co and Ni sulfide counterparts are metallic for any inversion degree together with the inverse FeMn2S4. Our calculations suggest that the spin filtering properties of the FeM2X4 (thio)spinels could be modified via the control of the cation distribution through variations in the synthesis conditions.