An energy budget agent-based model of earthworm populations and its application to study the effects of pesticides

Earthworms are important organisms in soil communities and so are used as model organisms in environmental risk assessments of chemicals. However current risk assessments of soil invertebrates are based on short-term laboratory studies, of limited ecological relevance, supplemented if necessary by site-specific field trials, which sometimes are challenging to apply across the whole agricultural landscape. Here, we investigate whether population responses to environmental stressors and pesticide exposure can be accurately predicted by combining energy budget and agent-based models (ABMs), based on knowledge of how individuals respond to their local circumstances. A simple energy budget model was implemented within each earthworm Eisenia fetida in the ABM, based on a priori parameter estimates. From broadly accepted physiological principles, simple algorithms specify how energy acquisition and expenditure drive life cycle processes. Each individual allocates energy between maintenance, growth and/or reproduction under varying conditions of food density, soil temperature and soil moisture. When simulating published experiments, good model fits were obtained to experimental data on individual growth, reproduction and starvation. Using the energy budget model as a platform we developed methods to identify which of the physiological parameters in the energy budget model (rates of ingestion, maintenance, growth or reproduction) are primarily affected by pesticide applications, producing four hypotheses about how toxicity acts. We tested these hypotheses by comparing model outputs with published toxicity data on the effects of copper oxychloride and chlorpyrifos on E. fetida. Both growth and reproduction were directly affected in experiments in which sufficient food was provided, whilst maintenance was targeted under food limitation. Although we only incorporate toxic effects at the individual level we show how ABMs can readily extrapolate to larger scales by providing good model fits to field population data. The ability of the presented model to fit the available field and laboratory data for E. fetida demonstrates the promise of the agent-based approach in ecology, by showing how biological knowledge can be used to make ecological inferences. Further work is required to extend the approach to populations of more ecologically relevant species studied at the field scale. Such a model could help extrapolate from laboratory to field conditions and from one set of field conditions to another or from species to species.

Implementation of a new urban energy budget scheme into MetUM. Part II: Validation against observations and model intercomparison

In the first part of this article, we introduced a new urban surface scheme, the Met Office – Reading Urban Surface Exchange Scheme (MORUSES), into the Met Office Unified Model (MetUM) and compared its impact on the surface fluxes with respect to the current urban scheme. In this second part, we aim to analyze further the reasons behind the differences. This analysis is conducted by a comparison of the performance of the two schemes against observations and against a third model, the Single Column Reading Urban model (SCRUM). The key differences between the three models lie in how each model incorporates the heat stored in the urban fabric and how the surface-energy balance is coupled to the underlying substrate. The comparison of the models with observations from Mexico City reveals that the performance of MORUSES is improved if roof insulation is included by minimizing the roof thickness. A comparison of MORUSES and SCRUM reveals that, once insulation is included within MORUSES, these two models perform equally well against the observations overall, but that there are differences in the details of the simulations at the roof and canyon level. These differences are attributed to the different representations of the heat-storage term, specifically differences in the dominant frequencies captured by the urban canopy and substrate, between the models. These results strongly suggest a need for an urban model intercomparison exercise. Copyright © 2010 Royal Meteorological Society and Crown Copyright

The ocean's gravitational potential energy budget in a coupled climate model

This study examines, in a unified fashion, the budgets of ocean gravitational potential energy (GPE) and available gravitational potential energy (AGPE) in the control simulation of the coupled atmosphere–ocean general circulation model HadCM3. Only AGPE can be converted into kinetic energy by adiabatic processes. Diapycnal mixing supplies GPE, but not AGPE, whereas the reverse is true of the combined effect of surface buoyancy forcing and convection. Mixing and buoyancy forcing, thus, play complementary roles in sustaining the large scale circulation. However, the largest globally integrated source of GPE is resolved advection (+0.57 TW) and the largest sink is through parameterized eddy transports (-0.82 TW). The effect of these adiabatic processes on AGPE is identical to their effect on GPE, except for perturbations to both budgets due to numerical leakage exacerbated by non-linearities in the equation of state.

A kinetic energy budget and internal instabilities in the Fine Resolution Antarctic Model

An energy analysis of the Fine Resolution Antarctic Model (FRAM) reveals the instability processes in the model. The main source of time-mean kinetic energy is the wind stress and the main sink is transfer to mean potential energy. The wind forcing thus helps maintain the density structure. Transient motions result from internal instabilities of the Bow rather than seasonal variations of the forcing. Baroclinic instability is found to be an important mechanism in FRAM. The highest values of available potential energy are found in the western boundary regions as well as in the Antarctic Circumpolar Current (ACC) region. All subregions with predominantly zonal flow are found to be baroclinically unstable. The observed deficit of eddy kinetic energy in FRAM occurs as a result of the high lateral friction, which decreases the growth rates of the most unstable waves. This high friction is required for the numerical stability of the model and can only be made smaller by using a finer horizontal resolution. A grid spacing of at least 10-15 km would be required to resolve the most unstable waves in the southern part of the domain. Barotropic instability is also found to be important for the total domain balance. The inverse transfer (that is, transfer from eddy to mean kinetic energy) does not occur anywhere, except in very localized tight jets in the ACC. The open boundary condition at the northern edge of the model domain does not represent a significant source or sink of eddy variability. However, a large exchange between internal and external mode energies is found to occur. It is still unclear how these boundary conditions affect the dynamics of adjacent regions.

Factors controlling the surface energy budget over snow and ice

Polar Regions are an energy sink of the Earth system, as the Sun rays do not reach the Poles for half of the year, and hit them only at very low angles for the other half of the year. In summer, solar radiation is the dominant energy source for the Polar areas, therefore even small changes in the surface albedo strongly affect the surface energy balance and, thus, the speed and amount of snow and ice melting. In winter, the main heat sources for the atmosphere are the cyclones approaching from lower latitudes, and the atmosphere-surface heat transfer takes place through turbulent mixing and longwave radiation, the latter dominated by clouds. The aim of this thesis is to improve the knowledge about the surface and atmospheric processes that control the surface energy budget over snow and ice, with particular focus on albedo during the spring and summer seasons, on horizontal advection of heat, cloud longwave forcing, and turbulent mixing during the winter season. The critical importance of a correct albedo representation in models is illustrated through the analysis of the causes for the errors in the surface and near-surface air temperature produced in a short-range numerical weather forecast by the HIRLAM model. Then, the daily and seasonal variability of snow and ice albedo have been examined by analysing field measurements of albedo, carried out in different environments. On the basis of the data analysis, simple albedo parameterizations have been derived, which can be implemented into thermodynamic sea ice models, as well as numerical weather prediction and climate models. Field measurements of radiation and turbulent fluxes over the Bay of Bothnia (Baltic Sea) also allowed examining the impact of a large albedo change during the melting season on surface energy and ice mass budgets. When high contrasts in surface albedo are present, as in the case of snow covered areas next to open water, the effect of the surface albedo heterogeneity on the downwelling solar irradiance under overcast condition is very significant, although it is usually not accounted for in single column radiative transfer calculations. To account for this effect, an effective albedo parameterization based on three-dimensional Monte Carlo radiative transfer calculations has been developed. To test a potentially relevant application of the effective albedo parameterization, its performance in the ground-based retrieval of cloud optical depth was illustrated. Finally, the factors causing the large variations of the surface and near-surface temperatures over the Central Arctic during winter were examined. The relative importance of cloud radiative forcing, turbulent mixing, and lateral heat advection on the Arctic surface temperature were quantified through the analysis of direct observations from Russian drifting ice stations, with the lateral heat advection calculated from reanalysis products.

Wave breaking on turbulent energy budget in the ocean surface mixed layer

As an important physical process at the air-sea interface, wave movement and breaking have a significant effect on the ocean surface mixed layer (OSML). When breaking waves occur at the ocean surface, turbulent kinetic energy (TKE) is input downwards, and a sublayer is formed near the surface and turbulence vertical mixing is intensively enhanced. A one-dimensional ocean model including the Mellor-Yamada level 2.5 turbulence closure equations was employed in our research on variations in turbulent energy budget within OSML. The influence of wave breaking could be introduced into the model by modifying an existing surface boundary condition of the TKE equation and specifying its input. The vertical diffusion and dissipation of TKE were effectively enhanced in the sublayer when wave breaking was considered. Turbulent energy dissipated in the sublayer was about 92.0% of the total depth-integrated dissipated TKE, which is twice higher than that of non-wave breaking. The shear production of TKE decreased by 3.5% because the mean flow fields tended to be uniform due to wave-enhanced turbulent mixing. As a result, a new local equilibrium between diffusion and dissipation of TKE was reached in the wave-enhanced layer. Below the sublayer, the local equilibrium between shear production and dissipation of TKE agreed with the conclusion drawn from the classical law-of-the-wall (Craig and Banner, 1994).

Parameterisation of bivalve functional traits for mechanistic eco-physiological dynamic energy budget (DEB) models

Mechanistic models such as those based on dynamic energy budget (DEB) theory are emergent ecomechanics tools to investigate the extent of fitness in organisms through changes in life history traits as explained by bioenergetic principles. The rapid growth in interest around this approach originates from the mechanistic characteristics of DEB, which are based on a number of rules dictating the use of mass and energy flow through organisms. One apparent bottleneck in DEB applications comes from the estimations of DEB parameters which are based on mathematical and statistical methods (covariation method). The parameterisation process begins with the knowledge of some functional traits of a target organism (e. g. embryo, sexual maturity and ultimate body size, feeding and assimilation rates, maintenance costs), identified from the literature or laboratory experiments. However, considering the prominent role of the mechanistic approach in ecology, the reduction of possible uncertainties is an important objective. We propose a revaluation of the laboratory procedures commonly used in ecological studies to estimate DEB parameters in marine bivalves. Our experimental organism was Brachidontes pharaonis. We supported our proposal with a validation exercise which compared life history traits as obtained by DEBs (implemented with parameters obtained using classical laboratory methods) with the actual set of species traits obtained in the field. Correspondence between the 2 approaches was very high (>95%) with respect to estimating both size and fitness. Our results demonstrate a good agreement between field data and model output for the effect of temperature and food density on age-size curve, maximum body size and total gamete production per life span. The mechanistic approach is a promising method of providing accurate predictions in a world that is under in creasing anthropogenic pressure.

Implementation of a new urban energy budget scheme in the MetUM. Part I: Description and idealized simulations

This paper describes the formulation of a new urban scheme, MORUSES (Met Office–Reading Urban Surface Exchange Scheme) for use in the Met Office Unified Model. The implementation of the new scheme ensures that (1) the new scheme offers more flexibility in the parametrization of the building properties, and hence provides a more realistic representation of the fluxes; (2) the bulk outputs are in satisfactory agreement with previous observational studies; and (3) the impact of the new scheme on the energy balance fluxes is similar to the impact of the current urban scheme when set up to mimic it. As well as having a better physical basis, MORUSES also gains in flexibility in applications and adaptations to different urban materials as well as urban planning. The new scheme represents the urban area as a composition of two tiles, a canyon and a roof, using a simple 2D geometry. Sensitivity analysis to canyon geometry and thickness of the roof canopy emphasizes the gain in flexibility captured by the new scheme. Copyright © 2010 Royal Meteorological Society and Crown Copyright

Modeling study of the aspect ratio influence on urban canopy energy fluxes with a modified wall-canyon energy budget scheme

The influence of the aspect ratio (building height/street canyon width) and the mean building height of cities on local energy fluxes and temperatures is studied by means of an Urban Canopy Model (UCM) coupled with a one-dimensional second-order turbulence closure model. The UCM presented is similar to the Town Energy Balance (TEB) model in most of its features but differs in a few important aspects. In particular, the street canyon walls are treated separately which leads to a different budget of radiation within the street canyon walls. The UCM has been calibrated using observations of incoming global and diffuse solar radiation, incoming long-wave radiation and air temperature at a site in So Paulo, Brazil. Sensitivity studies with various aspect ratios have been performed to assess their impact on urban temperatures and energy fluxes at the top of the canopy layer. In these simulations, it is assumed that the anthropogenic heat flux and latent heat fluxes are negligible. Results show that the simulated net radiation and sensible heat fluxes at the top of the canopy decrease and the stored heat increases as the aspect ratio increases. The simulated air temperature follows the behavior of the sensible heat flux. (C) 2010 Elsevier Ltd. All rights reserved.

Habitat characteristics and energy budget of the ocean quahog (Arctica islandica)

We compared lifetime and population energy budgets of the extraordinary long-lived ocean quahog Arctica islandica from 6 different sites - the Norwegian coast, Kattegat, Kiel Bay, White Sea, German Bight, and off northeast Iceland - covering a temperature and salinity gradient of 4-10°C (annual mean) and 25-34, respectively. Based on von Bertalanffy growth models and size-mass relationships, we computed organic matter production of body (PSB) and of shell (PSS), whereas gonad production (PG) was estimated from the seasonal cycle in mass. Respiration (R) was computed by a model driven by body mass, temperature, and site. A. islandica populations differed distinctly in maximum life span (40 y in Kiel Bay to 197 y in Iceland), but less in growth performance (phi' ranged from 2.41 in the White Sea to 2.65 in Kattegat). Individual lifetime energy throughput, as approximated by assimilation, was highest in Iceland (43,730 kJ) and lowest in the White Sea (313 kJ). Net growth efficiency ranged between 0.251 and 0.348, whereas lifetime energy investment distinctly shifted from somatic to gonad production with increasing life span; PS/PG decreased from 0.362 (Kiel Bay, 40 y) to 0.031 (Iceland, 197 y). Population annual energy budgets were derived from individual budgets and estimates of population mortality rate (0.035/y in Iceland to 0.173/y in Kiel Bay). Relationships between budget ratios were similar on the population level, albeit with more emphasis on somatic production; PS/ PG ranged from 0.196 (Iceland) to 2.728 (White Sea), and P/B ranged from 0.203-0.285/y. Life span is the principal determinant of the relationship between budget parameters, whereas temperature affects net growth efficiency only. In the White Sea population, both growth performance and net growth efficiency of A. islandica were lowest. We presume that low temperature combined with low salinity represent a particularly stressful environment for this species.

Global macrozoobenthos production and energy budget data base V150731, with link to database

I developed a new model for estimating annual production-to-biomass ratio P/B and production P of macrobenthic populations in marine and freshwater habitats. Self-learning artificial neural networks (ANN) were used to model the relationships between P/B and twenty easy-to-measure abiotic and biotic parameters in 1252 data sets of population production. Based on log-transformed data, the final predictive model estimates log(P/B) with reasonable accuracy and precision (r2 = 0.801; residual mean square RMS = 0.083). Body mass and water temperature contributed most to the explanatory power of the model. However, as with all least squares models using nonlinearly transformed data, back-transformation to natural scale introduces a bias in the model predictions, i.e., an underestimation of P/B (and P). When estimating production of assemblages of populations by adding up population estimates, accuracy decreases but precision increases with the number of populations in the assemblage.

On the measurement of the surface energy budget over a land surface during the summer monsoon

The measurement of surface energy balance over a land surface in an open area in Bangalore is reported. Measurements of all variables needed to calculate the surface energy balance on time scales longer than a week are made. Components of radiative fluxes are measured while sensible and latent heat fluxes are based on the bulk method using measurements made at two levels on a micrometeorological tower of 10 m height. The bulk flux formulation is verified by comparing its fluxes with direct fluxes using sonic anemometer data sampled at 10 Hz. Soil temperature is measured at 4 depths. Data have been continuously collected for over 6 months covering pre-monsoon and monsoon periods during the year 2006. The study first addresses the issue of getting the fluxes accurately. It is shown that water vapour measurements are the most crucial. A bias of 0.25% in relative humidity, which is well above the normal accuracy assumed the manufacturers but achievable in the field using a combination of laboratory calibration and field intercomparisons, results in about 20 W m(-2) change in the latent heat flux on the seasonal time scale. When seen on the seasonal time scale, the net longwave radiation is the largest energy loss term at the experimental site. The seasonal variation in the energy sink term is small compared to that in the energy source term.

Solubilization of n-decane into n-decyl-polyoxyethylene surfactants: A cavity forming free energy based model

The present paper reports the results of a theoretical study of the forces and factors driving the solubilization of n-alkane solubilizates into the micellar core of some non-ionic surfactants, based on a micellar model which includes the cavity forming free energy as a component of micellization. The solubilizate is n-decane and the non-ionic surfactants considered are n-decyl-polyoxyethylene surfactants. The extent of solubilization, i.e. the mole fraction of the solubilizate within the core has been calculated. The results indicate that the incorporated solubilizate has more translational and rotational degrees of freedom as compared to those of the tail parts of the surfactants present in the core. This drives the total free energy of aggregation after solubilization into a more favourable direction. The results are in fair agreement with the experimental results.

Growth and energy budget of F-2 'all-fish' growth hormone gene transgenic common carp

The growth and energy budget for F-2 'all-fish' growth hormone gene transgenic common carp Cyprinus carpio of two body sizes were investigated at 29.2 degrees C for 21 days. Specific growth rate, feed intake, feed efficiency, digestibility coefficients of dry matter and protein, gross energy intake (I-E), and the proportion of I-E utilized for heat production (H-E) were significantly higher in the transgenics than in the controls. The proportion of I-E directed to waste products [faecal energy (F-E) and excretory energy loss (Z(E) + U-E) where Z(E) is through the gills and U-E through the kidney], and the proportion of metabolizable energy (M-E) for recovered energy (R-E) were significantly lower in the transgenics than in the controls. The average energy budget equation of transgenic fish was as follows: 100 I-E = 19.3 F-E + 6.0 (Z(E) + U-E) + 45.2 H-E + 29.5 R-E or 100 M-E = 60.5 H-E + 39.5 R-E. The average energy budget equation of the controls was: 100 I-E = 25.2 F-E + 7.4 (Z(E) + U-E) + 35.5 H-E + 31.9 R-E or 100 M-E = 52.7 H-E + 47.3 R-E. These findings indicate that the high growth rate of 'all-fish' transgenic common carp relative to their non-transgenic counterparts was due to their increased feed intake, reduced lose of waste productions and improved feed efficiency. The benefit of the increased energy intake by transgenic fish, however, was diminished by their increased metabolism.

Effect of ration on the growth and energy budget of Chinese longsnout catfish, Leiocassis longirostris Gunther

The effect of ration on growth and energy budget of Chinese longsnout catfish was investigated in a growth trial. Fish of initial body weight of 6.5 g were fed at six ration levels (RLs): starvation, 0.8%, 1.6%, 2.4%, 3.2% of body weight per day, and apparent satiation for 8 weeks. Fish were weighed biweekly to adjust feed amount. The results showed that specific growth rate in wet weight, protein and energy increased logarithmically with increased RLs. The relationship of specific growth rate in wet weight (SGRw, % day(-1)) and RL (%) was a decelerating curve: SGRw=-0.62+3.10 Ln(RL+1). The energy budget equation at satiation was: 100 IE=12.94 FE+5.50(ZE+UE)+40.07 HE+41.49 RE, where IE, FE, (ZE+UE), HE, RE are food energy, faecal energy, excretory energy, heat production and recovered energy respectively. Body composition was slightly but significantly affected by ration size except for protein content. The most efficient ration based on the relationship between RL and feed efficiency ratio in energy (FERe) was 1.8% of body weight per day.