Proposing an information criterion for individual-based models developed in a pattern-oriented modelling framework

Individual-based models (IBMs) have been improved in quality and reliability in recent years with an approach called pattern-oriented modelling (POM). POM proposes guidelines to develop models reproducing multiple patterns observed on the field and to test systematically how well the IBMs reproduce them. POM studies used generally traditional methods of goodness of fit such as the sum of squares evaluation or ad hoc comparisons of fitting errors and variations. Model selection, however, can be a rigorous statistical approach based on information theory and information criteria such as the Akaike's information criterion (AIC) or the deviance information criterion (DIC). So far, it has not been tried to link POM to these rigorous techniques. The main problems to achieve that are: (a) the difficulty to have likelihood functions for IBMs’ parameters and (b) the possibility to obtain posterior distributions of IBMs’ parameters given the patterns to reproduce. In a first part, this paper answers problem (a) by proposing and explaining how to calculate a deviance measure (POMDEV) for models developed in a context of POM. And while answering the second problem, a second part of the paper proposes an information criterion for model selection in a POM context (the pattern-oriented modelling information criterion: POMIC). This criterion does not yet have the same theoretical foundation as, e.g., AIC, but uses formal analogies to the DIC. In a third part POMIC is tested with a modelling exercise. This exercise shows the potential of POMIC to use multiple patterns for selecting among multiple potential submodels and eventually select the most parsimonious and well fitting model version. We conclude that POMIC, although being a heuristically derived approach, can greatly improve the POM framework.     

Spatio-temporal dissolved oxygen dynamics in the Orbetello lagoon by fuzzy pattern recognition

This paper models the dissolved oxygen (DO) dynamics in the Orbetello lagoon as a function of the physico-chemical and ecological system variables, including the submerged vegetation, nutrients, and hydrodynamics. It should be viewed as the concluding sequel to a previous paper describing the dynamics of the lagoon ecosystem [Giusti, E., Marsili-Libelli, S., 2006. An integrated model for the Orbetello lagoon ecosystem, Ecol. Model. 196, 379–394] by introducing the missing DO dynamics. The model considers the oxygen demand originating from the decay of carbonaceous and nitrogenous compounds, as well as photosynthesis and natural reaeration by winds and currents as the oxygen producing processes. With a fixed-parameter set the model could accurately reproduce each single circadian DO cycle, but in the long run it failed to extend this fit and could not accommodate the large DO fluctuations induced by the seasonal variability. In order to enhance the model flexibility, a fuzzy pattern recognition algorithm was designed to classify the circadian DO patterns into four typical behaviours, related to the season, and estimate the corresponding parameters, with the overall model output being a fuzzy combination of these sets. The paper discusses several methods to patch the parameter sets and compares their performance in tracking long-term DO variations. A final assessment of the model validity is obtained by incorporating the whole DO dynamics (model, fuzzy pattern recognition and parameter combination) into the general lagoon model and producing a consistently correct series of DO daily distributions over a yearly cycle. Thus the paper contains both a practical and a methodological aspect. The practical one is the linking of all the lagoon dynamics to the dissolved oxygen kinetics in order to clarify to what extent macroalgae and macrophytes influence the oxygen balance. The methodological aspect consists of extending the validity of short-term models to long time-horizons through a patching technique supported by fuzzy pattern recognition.     

The regional economic impact of reducing greenhouse gas emissions: Western Australia

This study analyses the general-equilibrium impacts of an international climate change response policy on the economy of Western Australia (WA), one of the most mining-based and energy-intensive states of Australia. It finds that emissions would fall by up to 11% from the base level in WA. However, such environmental benefits emanate at some costs to the state economy; in terms of foregone gross state product, the costs are up to 3% of the base level. Indeed, the actual costs and benefits depend on the precise design of the climate change response policy as well as on the other policies within which it operates. For example, when emission quota permits are sold to industries and no tradeable carbon credits (i.e. credits for the carbon sequestrated in Kyoto forests) are granted, emissions decline by about 8% and GSP falls by about 3% of the base levels. If carbon credits are tradeable, however, the environmental benefits could be increased and the GSP cost could be reduced substantially. Also, the reduced economic activity caused by emission abatement results in a modest fall in net government revenue, despite the additional revenue from permit sales in some cases. Accordingly, government’s fiscal package surrounding the emission permits would influence the emission abatement impacts on the economy. With regard to the effects on the structure of the state economy, the oil and gas industry suffers only a slight contraction but the energy-supplying sector as a whole contracts substantially. It is therefore not surprising that the impacts on the WA economy of curbing emissions by energy and transport industries alone are quite significant when compared to those resulted from all industries’ compliance with the abatement scheme. It needs to be noted that the model projections analysed in the paper are based on simplifying assumptions and tentative scenarios, and hence should be viewed with caution and not be understood as unconditional forecasts.     

Eco-hydrological modelling in a highly regulated lowland catchment to find measures for improving water quality

Water quality modelling in the meso-scale Rhin catchment in the German federal state Brandenburg was done (1) to answer some specific questions concerning identification of point and diffuse sources of nutrient pollution in the catchment, (2) to assess the influences of possible climate and land use changes on water quantity and quality and (3) to evaluate potential measures to be done in order to achieve a “good ecological status” of the river and its lakes as required by the Water Framework Directive (WFD).The Rhin catchment is a typical highly regulated lowland river basin in Northern Germany. The regulations complicate water quantity and quality modelling in the catchment. The research was done by using the eco-hydrological model SWIM (Soil and Water Integrated Model), which simulates water and nutrient fluxes in soil and vegetation, as well as transport of water and nutrients to and within the river network. The modelling period was from 1981 until 2005. After calibrating the hydrological processes at different gauges within the basin with satisfactory results, water quality (nitrogen and phosphorus) modelling was done taking into account the emissions of different point sources (sewage treatment plants, etc.) and identifying the amount of diffuse pollution caused mainly by agriculture.For suggesting some feasible measures to improve water quality and to reduce diffuse pollution considering possible climate and land use changes, different reasonable scenarios were applied in consultation with the Environmental Agency of Brandenburg (LUA). The study revealed that the amount of water discharge has significant influence on the concentration of nutrients in the river network, and that nitrogen pollution, caused mainly by diffuse sources, could be notably reduced by application of agricultural measures, whereas the pollution by phosphorus could be diminished most effectively by the reduction of point source emissions.     

Modelling of stem biomass accumulation in Pinus densiflora seedlings exposed to aqueous-phase OH radicals generating mist

The present study aimed to evaluate the short-term biomass accumulation of forest trees exposed to wet acidic depositions. A hierarchical Bayesian model of tree growth was developed based on the data of a short-term experiment in which 2-year-old Japanese red pine (Pinus densiflora Sieb. et Zucc.) seedlings were exposed to aqueous phase OH radicals generated by an iron-oxalate-H₂O₂ mist (a pseudo-polluted dew) over two growing periods. We conducted a statistical comparison of tree growth between the control and pollution treatment groups by using the growth model incorporated the random effects due to the unknown characteristics of each seedling. The variability among seedlings is expressed in this model by the posterior probabilistic distributions of unobserved dry weight of a stem cohort before exposure treatment. The analysis of the effects of pollution treatment on the stem growth revealed that this treatment decreases the biomass allocation in the current year stems. However, the effects on the relative growth rate of pre-existing stems were unclear. Based on these results, we can speculate that in a polluted environment, the short-term growth of the young stems in the seedlings inhibited by pollution treatment, thereby resulting in the slowdown of long-term biomass accumulation. This can explain the patterns observed in the declining Japanese red pine forests that are subjected to OH-generating dews in the polluted area of western Japan.     

Consequences of changed wetness on riverine nitrogen – human impact on retention vs. natural climatic variability

The HBV-N model was used for a scenario analysis of changes in nitrogen retention and transport caused by alterations of wetness due to land drainage, lowering of lakes, building of dams and climatic variability in a river basin in south-central Sweden (1885–1994). In general, dams were situated in locations more favourable for retention, compared to the lowered lakes. Rather modest conversions of water bodies only changed nitrogen transport by about 3%. The 180-times-larger increase of (mainly) tile-drained agricultural land had, according to simulations, increased the nitrogen transport by 17%, due to reduced retention. However, compared to human-induced alteration of the landscape N retention, the choice of 10-year periods of climatological data had the overriding effect on the calculated nitrogen transport. Weather-induced variations resulted in a 13% difference in nitrogen retention between various 10-year periods. When the model was driven by climatological data from the driest 10-year period (1905–1914), the estimated average annual load was only half of that obtained with climatological data from the wettest 10-year period (1975–1984). Electronic Publication     

Assessing the economic feasibility of regional deep saline aquifer CO2 injection and storage: A geomechanics-based workflow applied to the Rose Run sandstone in Eastern Ohio, USA

Typical top-down regional assessments of CO₂ storage feasibility are sufficient for determining the maximum volumetric capacity of deep saline aquifers. However, they do not reflect the regional economic feasibility of storage. This is controlled, in part, by the number and type of injection wells that are necessary to achieve regional CO₂ storage goals. In contrast, the geomechanics-based assessment workflow that we present in this paper follows a bottom-up approach for evaluating regional deep saline aquifer CO₂ storage feasibility. The CO₂ storage capacity of an aquifer is a function of its porous volume as well as its CO₂ injectivity. For a saline aquifer to be considered feasible in this assessment it must be able to store a specified amount of CO₂ at a reasonable cost per ton of CO₂. The proposed assessment workflow has seven steps that include (1) defining the storage project and goals, (2) characterizing the geology and developing a geomechanical model of the aquifer, (3) constructing 3D aquifer models, (4) simulating CO₂ injection, (5,6) evaluating CO₂ injection and storage feasibility (with and without injection well stimulation), and (7) determining whether it is economically feasible to proceed with the storage project. The workflow was applied to a case study of the Rose Run sandstone aquifer in the Eastern Ohio River Valley, USA. We found that it is feasible in this region to inject 113 Mt CO₂/year for 30 years at an associated well cost of less than US $1.31/t CO₂, but only if injectivity enhancement techniques such as hydraulic fracturing and injection induced micro-seismicity are implemented.     

Long-term geochemical evolution of the near field repository: insights from reactive transport modelling and experimental evidences

The KBS-3 underground nuclear waste repository concept designed by the Swedish Nuclear Fuel and Waste Management Co. (SKB) includes a bentonite buffer barrier surrounding the copper canisters and the iron insert where spent nuclear fuel will be placed. Bentonite is also part of the backfill material used to seal the access and deposition tunnels of the repository. The bentonite barrier has three main safety functions: to ensure the physical stability of the canister, to retard the intrusion of groundwater to the canisters, and in case of canister failure, to retard the migration of radionuclides to the geosphere. Laboratory experiments (< 10 years long) have provided evidence of the control exerted by accessory minerals and clay surfaces on the pore water chemistry. The evolution of the pore water chemistry will be a primordial factor on the long-term stability of the bentonite barrier, which is a key issue in the safety assessments of the KBS-3 concept.In this work we aim to study the long-term geochemical evolution of bentonite and its pore water in the evolving geochemical environment due to climate change. In order to do this, reactive transport simulations are used to predict the interaction between groundwater and bentonite which is simulated following two different pathways: (1) groundwater flow through the backfill in the deposition tunnels, eventually reaching the top of the deposition hole, and (2) direct connection between groundwater and bentonite rings through fractures in the granite crosscutting the deposition hole. The influence of changes in climate has been tested using three different waters interacting with the bentonite: present-day groundwater, water derived from ice melting, and deep-seated brine. Two commercial bentonites have been considered as buffer material, MX-80 and Deponit CA-N, and one natural clay (Friedland type) for the backfill. They show differences in the composition of the exchangeable cations and in the accessory mineral content. Results from the simulations indicate that pore water chemistry is controlled by the equilibrium with the accessory minerals, especially carbonates. pH is buffered by precipitation/dissolution of calcite and dolomite, when present. The equilibrium of these minerals is deeply influenced by gypsum dissolution and cation exchange reactions in the smectite interlayer. If carbonate minerals are initially absent in bentonite, pH is then controlled by surface acidity reactions in the hydroxyl groups at the edge sites of the clay fraction, although its buffering capacity is not as strong as the equilibrium with carbonate minerals. The redox capacity of the bentonite pore water system is mainly controlled by Fe(II)-bearing minerals (pyrite and siderite). Changes in the groundwater composition lead to variations in the cation exchange occupancy, and dissolution–precipitation of carbonate minerals and gypsum. The most significant changes in the evolution of the system are predicted when ice-melting water, which is highly diluted and alkaline, enters into the system. In this case, the dissolution of carbonate minerals is enhanced, increasing pH in the bentonite pore water. Moreover, a rapid change in the population of exchange sites in the smectite is expected due to the replacement of Na for Ca.