Trophic relationships in the recently impounded Bagré reservoir in Burkina Faso

The trophic dynamics of Bagré reservoir which has been recently impounded in Burkina Faso was based on the data collected during 1997–1998 period using the Ecopath model and software. Total fish biomass is 22.63 t km⁻² and mainly represents trophic levels (TLs) 2 and 3. The trophic food chain is relatively long and the overall transfer efficiency is quite low. Grazing foodweb based on primary producers is prominent in the reservoir ecosystem and detritus plays a less significant role. Seasonal and long-term variations in water quality have significant influences on the lower TLs clearly showing a bottom-up functioning of the ecosystem. Environmental degradations, such as siltation occurring in the lake, suggest possible risks in limiting ecosystem productivity.     

Mass balanced and dynamic simulations of trophic models of kelp ecosystems near the Mejillones Peninsula of northern Chile (SE Pacific): Comparative network structure and assessment of harvest strategies

Mass balanced trophic models for kelp ecosystems which include subsystems dominated by Macrocystis integrifolia, Lessonia trabeculata and areas of barren ground (BG) were constructed for subtidal areas near the Mejillones Peninsula (SE Pacific), Chile. Information on biomass, P/B ratios, catches, food spectrum, consumption and dynamics of commercial and non-commercial populations was obtained and examined using Ecopath with Ecosim software analyses. The biomass of blades of L. trabeculata and M. integrifolia represented the compartments most relevant to the subsystems studied. Within the herbivores, the sea urchin Tetrapigus niger was dominant, followed by the snails Turritella sp. and Tegula sp. The fishes Pinguipes chilensis and Cheilodactylus variegatus were the dominant predators, followed by the asteroids Heliaster helianthus and Meyenaster gelatinosus. The highest system throughput (72,512 g wet weight m⁻² year⁻¹) was calculated for the subsystem dominated by M. integrifolia. The mean trophic level of the catch ranged from 1.1 (subsystem dominated by L. trabeculata) to 1.3 (subsystem dominated by M. integrifolia) to 3.2 (barren ground subsystem), showing that harvesting in each system was concentrated either on primary producers (blades of kelp species) or top predator fishes. Although the values for the Relative Ascendency (A/C) fluctuated from 36.5 to 45%, suggesting that all the systems were immature, the subsystem dominated by M. integrifolia emerged as the least resistant to external disturbances (e.g. fisheries). This result agreed with the high value of the system recovery time (SRT) for the M. integrifolia subsystem as a response to combined fisheries scenarios. The results obtained using mixed trophic impact (MTI) and Ecosim [increasing the fishing mortality F_i by 4×] showed that in most of the cases the predictions had the same qualitative tendencies. One of the most important results obtained in this study was that exploitation of kelp blades as an alternative strategy to harvesting the whole plants appeared to be ecologically sustainable, since harvesting the blades propagated only small effects on the entire subsystem. The fish P. chilensis may be considered as a top predator species with a strong top-down control since an increase in its fishing mortality in the subsystem dominated by M. integrifolia produced a high SRT value, and the F_MSY was less than the originally entered F_i in Ecopath. Based on the results obtained, it was concluded that the trophic mass balanced models and simulated management scenarios offered good possibilities for the planning of interventions and manipulations or the planning of more sustainable management strategies in highly disturbed natural systems.     

Assessing the carrying capacity of tilapia in an intertidal mangrove-based polyculture system of Pearl River Delta, China

A trophic model of an intertidal mangrove-based polyculture system in Pearl River Delta, China, was constructed using the Ecopath with Ecosim software. This polyculture system was chosen since it is the first integrated multi-trophic aquaculture (IMTA) system that was constructed on the basis of mangrove planting in China. The energy flows, ecosystem property, and carrying capacity of tilapia in the polyculture system were analyzed and evaluated. The results show the trophic level of 1.00 for primary producers and detritus to 2.85 for grass carp. The geometric mean of the trophic transfer efficiencies was 7.0%, with 7.2% from detritus and 6.8% from primary producers within the system. The ecosystem property indices show that this polyculture system has a high value of total primary production/total respiration (TPP/TR) and total primary production/total biomass (TPP/TB), together with low Finn's cycled index (FCI), Finn's mean path length (FML), and connectance index (CI), indicating that this system is at a development stage according to Odum's theory. The principal fish cultured in the system is tilapia, and mixed trophic impacts (MTI) show that tilapia has a marked impact on most compartments in this system, and the carrying capacity was found to be a tilapia culture biomass of 5.8 t ha⁻¹ in the system.     

Assessing biomass gains from marsh restoration in Delaware Bay using Ecopath with Ecosim

The Delaware Bay ecosystem has been the focus of extensive habitat restoration efforts to offset finfish losses due to mortality associated with power plant water intake. As a result, a 45 km² or a 3% increase in total marsh area was achieved by 1996-1997 through the restoration efforts of the Public Service Enterprise Group (PSEG). To quantify the impact of restoration efforts on system productivity, an Ecopath with Ecosim model was constructed that represented all major components of the ecosystem. The model consisted of 47 functional groups including: 27 fish species, 5 invertebrate groups, 4 multi-species benthic groups, 6 multi-species fish groups, 3 plankton groups, 1 shorebird group and 1 marine mammal group. Biomass, abundance, catch, and demographic data were obtained from the literature or from individual stock assessments conducted for principal ecosystem components. A base Ecosim model was fitted to time series of key species in the Bay representing the period 1966-2003. To access the gains from marsh restoration, model simulations reflecting no restoration were conducted to estimate the productivity that would have been lost if restoration efforts had not occurred. The results indicated that restoration increased total ecosystem biomass by 47.7 t km⁻² year⁻¹. Simulations indicated increased biomasses across a wide range of species including important forage and commercially important species. The marsh restoration also significantly impacted ecosystem structure increasing the ratio of production-to-respiration, increasing system path length and decreasing the ratio of production-to-biomass.     

Response of balanced network models to large-scale perturbation: Implications for evaluating the role of small pelagics in the Gulf of Maine

Exploring the response of an ecosystem, and subsequent tradeoffs among its biological community, to human perturbations remains a key challenge for the implementation of an ecosystem approaches to fisheries (EAF). To address this and related issues, we developed two network (or energy budget) models, Ecopath and Econetwrk, for the Gulf of Maine ecosystem. These models included 31 network “nodes” or biomass state variables across a broad range of trophic levels, with the present emphasis to particularly elucidate the role of small pelagics. After initial network balancing, various perturbation scenarios were evaluated to explore how potential changes to different fish, fisheries and lower trophic levels can affect model outputs. Categorically across all scenarios and interpretations thereof, there was minimal change at the second trophic levels and most of the “rebalancing” after a perturbation occurred via alteration of the diet matrix. Yet the model results from perturbations to a balanced energy budget fall into one of three categories. First, some model results were intuitive and in obvious agreement with established ecological and fishing theory. Second, some model results were counter-intuitive upon initial observation, seemingly contradictory to known ecological and fishing theory; but upon further examination the results were explainable given the constraints of an equilibrium energy budget. Finally, some results were counter-intuitive and difficult to reconcile with theory or further examination of equilibrium constraints. A detailed accounting of biomass flows for example scenarios explores some of the non-intuitive results more rigorously. Collectively these results imply a need to carefully track biomass flows and results of any given perturbation and to critically evaluate the conditions under which a new equilibrium is obtained for these types of models, which has implications for dynamic simulations based off of them. Given these caveats, the role of small pelagics as a prominent component of this ecosystem remains a robust conclusion. We discuss how one might use this approach in the context of further developing an EAF, recognizing that a more holistic, integrated perspective will be required as we continue to evaluate tradeoffs among marine biological communities.     

Calculating ecological carrying capacity of shellfish aquaculture using mass-balance modeling: Narragansett Bay, Rhode Island

Increasing growth in the aquaculture industry demands ecosystem-based techniques for management if that growth is to be ecologically sustainable and promote equity among users of the ecosystems in which it occurs. Models of carrying capacity can be used to responsibly limit the growth of aquaculture in increasingly crowded coastal areas. Narragansett Bay, Rhode Island, USA is one such crowded coastal region experiencing a rapid increase in bivalve aquaculture. An ecosystem mass-balance model was used to calculate the ecological carrying capacity of bivalve aquaculture. Cultured oyster biomass is currently at 0.47 t km⁻² and could be increased 625 times without exceeding the ecological carrying capacity of 297 t km⁻². This translates to approximately 38,950 t of harvested cultured oysters annually which is 4 times the total estimated annual harvest of finfish. This potential for growth is due to the high primary productivity and large energy throughput to detritus of this ecosystem. Shellfish aquaculture has potential for continued growth and is unlikely to become food limited due, in part, to the large detritus pool.     

Estimating California ecosystem carbon change using process model and land cover disturbance data: 1951-2000

Land use change, natural disturbance, and climate change directly alter ecosystem productivity and carbon stock level. The estimation of ecosystem carbon dynamics depends on the quality of land cover change data and the effectiveness of the ecosystem models that represent the vegetation growth processes and disturbance effects. We used the Integrated Biosphere Simulator (IBIS) and a set of 30- to 60-m resolution fire and land cover change data to examine the carbon changes of California's forests, shrublands, and grasslands. Simulation results indicate that during 1951-2000, the net primary productivity (NPP) increased by 7%, from 72.2 to 77.1 Tg C yr⁻¹ (1 teragram = 10¹² g), mainly due to CO₂ fertilization, since the climate hardly changed during this period. Similarly, heterotrophic respiration increased by 5%, from 69.4 to 73.1 Tg C yr⁻¹, mainly due to increased forest soil carbon and temperature. Net ecosystem production (NEP) was highly variable in the 50-year period but on average equalled 3.0 Tg C yr⁻¹ (total of 149 Tg C). As with NEP, the net biome production (NBP) was also highly variable but averaged −0.55 Tg C yr⁻¹ (total of -27.3 Tg C) because NBP in the 1980s was very low (-5.34 Tg C yr⁻¹). During the study period, a total of 126 Tg carbon were removed by logging and land use change, and 50 Tg carbon were directly removed by wildland fires. For carbon pools, the estimated total living upper canopy (tree) biomass decreased from 928 to 834 Tg C, and the understory (including shrub and grass) biomass increased from 59 to 63 Tg C. Soil carbon and dead biomass carbon increased from 1136 to 1197 Tg C.Our analyses suggest that both natural and human processes have significant influence on the carbon change in California. During 1951-2000, climate interannual variability was the key driving force for the large interannual changes of ecosystem carbon source and sink at the state level, while logging and fire were the dominant driving forces for carbon balances in several specific ecoregions. From a long-term perspective, CO₂ fertilization plays a key role in maintaining higher NPP. However, our study shows that the increase in C sequestration by CO₂ fertilization is largely offset by logging/land use change and wildland fires.     

Spatial and mass balanced trophic models of La Rinconada Marine Reserve (SE Pacific coast), a protected benthic ecosystem: Management strategy assessment

Steady-state, dynamic, and spatial models were constructed for the benthic system of La Rinconada Marine Reserve off northern Chile (SE Pacific coast). We examined data on biomass, P/B ratios, catches, food spectrum, consumption, and the dynamics of commercial and non-commercial populations using three theoretical frameworks: Ecopath, Ecosim, and Ecospace. The biomass of the scallop Argopecten purpuratus and the clam Tagelus dombeii were the most relevant compartments of the studied ecosystem. Among the carnivores, the functional crab group Cancer spp. was the most relevant. The Rhodophyta was the dominant macroalga compartment of the system. The results obtained using mixed trophic impacts (MTI) showed that the predatory snail Thais chocolata propagated higher magnitudes of direct and indirect effects on the other species or functional groups. The sea star Luidia magallanica and Rhodophyta had the least effects on the remaining compartments. According to the Ecosim estimates (increasing mortality by fishing), the scallop A. purpuratus had the highest impact on the other compartments. The Ecospace model showed similar qualitative and quantitative effects for changes in biomass under three different exploitation scenarios (by subsystems and globally). Nevertheless, the greatest changes were provoked by using the top-down control and the vulnerabilities estimated by Ecosim. System recovery times were highest with increased mortality of the asteroid L. magallanica and the carnivorous snail T. chocolata, suggesting that the sea star could be considered to be a top predator with a top-down control. The F_MSY estimated for the scallop A. purpuratus was close to the F_i originally entered in Ecopath, limiting the design and execution of an exploitation plan within ecologically sustainable boundaries. The situation was different (F_MSY ? F_i) for the other commercial species, making possible multi-species exploitation programs. The Ecospace trophic-spatially explicit model shows a similar pattern of direct and indirect effects generated when exerting exploitation separately by subsystems. Therefore, habitat rotation of fisheries is not justified.     

Inter-annual variability of ecosystem production in boreal jack pine forests (1975-2004) estimated from tree-ring data using CBM-CFS3

We describe and apply a method of using tree-ring data and an ecosystem model to reconstruct past annual rates of ecosystem production. Annual data on merchantable wood volume increment and mortality obtained by dendrochronological stand reconstruction were used as input to the Carbon Budget Model of the Canadian Forest Sector (CBM-CFS3) to estimate net ecosystem production (NEP), net primary production (NPP), and heterotrophic respiration (Rh) annually from 1975 to 2004 at 10 boreal jack pine (Pinus banksiana Lamb.) stands in Saskatchewan and Manitoba, Canada. From 1975 (when sites aged 41-60 years) to 2004 (when they aged 70-89 years), all sites were moderate C sinks except during some warmer than average years where estimated Rh increased. Across all sites and years, estimated annual NEP averaged 57 g Cm⁻² yr⁻¹ (range −31 to 176 g Cm⁻² yr⁻¹), NPP 244 g Cm⁻² yr⁻¹ (147-376 g Cm⁻² yr⁻¹), and Rh 187 g Cm⁻² yr⁻¹ (124-270 g Cm⁻² yr⁻¹). Across all sites, NPP was related to stand age and density, which are proxies for successional changes in leaf area. Regionally, warm spring temperature increased NPP and defoliation by jack pine budworm 1 year previously reduced NPP. Our estimates of NPP, Rh, and NEP were plausible when compared to regional eddy covariance and carbon stock measurements. Inter-annual variability in ecosystem productivity contributes uncertainty to inventory-based assessments of regional forest C budgets that use yield curves predicting averaged growth over time. Our method could expand the spatial and temporal coverage of annual forest productivity estimates, providing additional data for the development of empirical models accounting for factors not presently considered by these models.     

EcoTroph: Modelling marine ecosystem functioning and impact of fishing

EcoTroph (ET) is a model articulated around the idea that the functioning of aquatic ecosystems may be viewed as a biomass flow moving from lower to higher trophic levels, due to predation and ontogenetic processes. Thus, we show that the ecosystem biomass present at a given trophic level may be estimated from two simple equations, one describing biomass flow, the other their kinetics (which quantifies the velocity of biomass transfers towards top predators). The flow kinetic of prey partly depends on the abundance of their predators, and a top-down equation expressing this is included in the model. Based on these relationships, we simulated the impact on a virtual ecosystem of various exploitation patterns. Specifically, we show that the EcoTroph approach is able to mimic the effects of increased fishing effort on ecosystem biomass expected from theory. Particularly, the model exhibits complex patterns observed in field data, notably cascading effects and ‘fishing down the food web’. EcoTroph also provides diagnostic tools for examining the relationships between catch and fishing effort at the ecosystem scale and the effects of strong top-down controls and fast-flow kinetics on ecosystems resilience. Finally, a dynamic version of the model is derived from the steady-state version, thus allowing simulations of time series of ecosystem biomass and catches. Using this dynamic model, we explore the propagation of environmental variability in the food web, and illustrated how exploitation can induce a decrease of ecosystem stability. The potential for applying EcoTroph to specific ecosystems, based on field data, and similarities between EcoTroph and Ecopath with Ecosim (EwE) are finally discussed.     

Spatially distributed modeling of the long-term carbon balance of a boreal landscape

Spatially and temporally distributed information on the sizes of biomass carbon (C) pools (BCPs) and soil C pools (SCPs) is vital for improving our understanding of biosphere-atmosphere C fluxes. Because the sizes of C pools result from the integrated effects of primary production, age-effects, changes in climate, atmospheric CO₂ concentration, N deposition, and disturbances, a modeling scheme that interactively considers these processes is important. We used the InTEC model, driven by various spatio-temporal datasets to simulate the long-term C-balance in a boreal landscape in eastern Canada. Our results suggested that in this boreal landscape, mature coniferous stands had stabilized their productivity and fluctuated as a weak C-sink or C-source depending on the interannual variations in hydrometeorological factors. Disturbed deciduous stands were larger C-sinks (NEP₂₀₀₄ = 150 gC m⁻² yr⁻¹) than undisturbed coniferous stands (e.g. NEP₂₀₀₄ = 8 gC m⁻² yr⁻¹). Wetlands had lower NPP but showed temporally consistent C accumulation patterns. The simulated spatio-temporal patterns of BCPs and SCPs were unique and reflected the integrated effects of climate, plant growth and atmospheric chemistry besides the inherent properties of the C pool themselves. The simulated BCPs and SCPs generally compared well with the biometric estimates (BCPs: r = 0.86, SCPs: r = 0.84). The largest BCP biases were found in recently disturbed stands and the largest SCP biases were seen in locations where moss necro-masses were abundant. Reconstructing C pools and C fluxes in the ecosystem in such a spatio-temporal manner could help reduce the uncertainties in our understanding of terrestrial C-cycle.     

Structural changes in lake functioning induced from nutrient loading and climate variability

Climate variability is increasingly recognized as an important regulatory factor, capable of influencing the structural properties of aquatic ecosystems. Lakes appear to be particularly sensitive to the ecological impacts of climate variability, and several long time series have shown a close coupling between climate, lake thermal properties and individual organism physiology, population abundance, community structure, and food web dynamics. Thus, understanding the complex interplay among meteorological forcing, hydrological variability, and ecosystem functioning is essential for improving the credibility of model-based water resources/fisheries management. Our objective herein is to examine the relative importance of the ecological mechanisms underlying plankton seasonal variability in Lake Washington, Washington State (USA), over a 35-year period (1964–1998). Our analysis is founded upon an intermediate complexity plankton model that is used to reproduce the limiting nutrient (phosphate)–phytoplankton–zooplankton–detritus (particulate phosphorus) dynamics in the lake. Model parameterization is based on a Bayesian calibration scheme that offers insights into the degree of information the data contain about model inputs and allows obtaining predictions along with uncertainty bounds for modeled output variables. The model accurately reproduces the key seasonal planktonic patterns in Lake Washington and provides realistic estimates of predictive uncertainty for water quality variables of environmental management interest. A principal component analysis of the annual estimates of the underlying ecological processes highlighted the significant role of the phosphorus recycling stemming from the zooplankton excretion on the planktonic food web variability. We also identified a moderately significant signature of the local climatic conditions (air temperature) on phytoplankton growth (r = 0.41), herbivorous grazing (r = 0.38), and detritus mineralization (r = 0.39). Our study seeks linkages with the conceptual food web model proposed by Hampton et al. [Hampton, S.E., Scheuerell, M.D., Schindler, D.E., 2006b. Coalescence in the Lake Washington story: interaction strengths in a planktonic food web. Limnol. Oceanogr. 51, 2042–2051.] to emphasize the “bottom-up” control of the Lake Washington plankton phenology. The posterior predictive distributions of the plankton model are also used to assess the exceedance frequency and confidence of compliance with total phosphorus (15 μg L⁻¹) and chlorophyll a (4 μg L⁻¹) threshold levels during the summer-stratified period in Lake Washington. Finally, we conclude by underscoring the importance of explicitly acknowledging the uncertainty in ecological forecasts to the management of freshwater ecosystems under a changing global environment.     

Integrative trophic network assessments of a lentic ecosystem by key ecological approaches of water chemistry, trophic guilds, and ecosystem health assessments along with an ECOPATH model

The objective of this study was to describe the trophic structure and energy flow in a lentic ecosystem in South Korea. Physicochemical water conditions were evaluated along with the reservoir ecosystem health using a multimetric IBI model. Nutrient analyses of the reservoir showed a nutrient rich and hypereutrophic system. Guild analysis revealed that tolerant and omnivorous species dominated the ecosystem. Tolerant fish, as a proportion of the number of individuals, were associated (R² > 0.90, p < 0.01) with TN and TP, the key indicators of trophic state in lentic ecosystems. The mean Reservoir Ecosystem Health Assessment (REHA) score was 19.3 during the study, which was judged as in ‘fair to poor’ condition. A trophic analysis of the reservoir estimated by the ECOPATH model shows that most activity in terms of energy flow occurred in the lower part of the trophic web, where there was intensive use of primary producers as a food source. Consequently, of the 10 consumer groups, nine fell within trophic levels <2.8. Trophic levels (TL) estimated from the weighted average of prey trophic levels varied from 1.0 for phytoplankton, macrophytes, and detritus to 3.25 for the top predator, Pseudobagrus fulvidraco. Our integrated approach to trophic network analysis may provide a key tool for determining the effects of nutrient influx on energy flow pathways in lentic ecosystems.     

Database-driven models of the world's Large Marine Ecosystems

We present a new methodology for database-driven ecosystem model generation and apply the methodology to the world's 66 currently defined Large Marine Ecosystems. The method relies on a large number of spatial and temporal databases, including FishBase, SeaLifeBase, as well as several other databases developed notably as part of the Sea Around Us project. The models are formulated using the freely available Ecopath with Ecosim (EwE) modeling approach and software. We tune the models by fitting to available time series data, but recognize that the models represent only a first-generation of database-driven ecosystem models. We use the models to obtain a first estimate of fish biomass in the world's LMEs. The biggest hurdles at present to further model development and validation are insufficient time series trend information, and data on spatial fishing effort.     

Modeling an exploited rocky coastal ecosystem: Bahia Tortugas, Mexico

A trophic structure model of the rocky coastal ecosystem in Bahia Tortugas, Mexico was constructed using Ecopath software to represent the main biomass flows in the system. Data for the model came from field observations (biomass estimates, stomach contents, and ecological observations for sea snails, abalones, lobster, some demersal finfishes, and macroalgae) carried out through ten field trips from 2006 to 2008. The results provide a snapshot of how the ecosystem operates. The model considers 23 functional groups. The total system throughput was 553 t/km²/year, 57% corresponds to internal consumption, 28% to respiration, 14% becomes detritus, and only 1% is removed through commercial fishing. The model suggests that even for exploited populations, predation and competition are heavier stresses than current fishing effort; however, because spiny lobster showed the second highest keystoneness’ index value, increasing fishing pressure on this group could strongly impact the entire ecosystem. We believe that this model has the potential to support management by allowing the exploration of the potential impacts of different fishing decisions at ecosystem level.     

Food web structure and trophic interactions of the tropical highland lake Hayq, Ethiopia

We generated a mass-balance model to figure out the food web structure and trophic interactions of the major functional groups of the Ethiopian highland Lake Hayq. Moreover, the study lay down a baseline data for future ecosystem-based investigations and management activities. Extensive data collection has been taken place between October 2007 and May 2009. Ecotrophic efficiency (EE) of several functional groups including phytoplankton (0.8) and detritus (0.85) was high indicating the utilization of the groups within the system. However, the EE of Mesocyclops (0.03) and Thermocyclops (0.30) was very low implying these resources were rather a ‘sink’ in the trophic hierarchy. Flows based on aggregated trophic level sensu Lindeman revealed the importance of both phytoplankton and detritus to higher trophic levels. The computed average transfer efficiency of 11.5% for the first four trophic levels was within the range for highly efficient African lakes. The primary production to respiration (P/R) ratio (1.05) of Lake Hayq indicates the maturity of the ecosystem. We also modeled the food-web by excluding Tilapia and reduced phytoplankton biomass to get insight into the mass balance before Tilapia was introduced. The analysis resulted in a lower system omnivory index (SOI = 0.016) and a reduced P/R ratio (0.13) that described the lake as immature ecosystem, suggesting the introduction of Tilapia might have contributed to the maturity of the lake. Tilapia in Lake Hayq filled an ecological empty niche of pelagic planktivores, and contributed for the better transfer efficiency observed from primary production to fish yield.     

Dealing with microtopography of an ombrotrophic bog for simulating ecosystem-level CO2 exchanges

Peatlands contain approximately 25% of the global soil carbon (C), despite covering only 3% of the earth's land surface. In order to evaluate the role of peatlands in global C cycling, models of ecosystem biogeochemistry are required, but peatland ecosystems present a number of unique challenges, particularly how to deal with the large variability that occurs at scales of one to several metres. In models, spatial variability is considered either explicitly for each individual unit and the outputs averaged, referred to as flux upscaling, or implicitly by weighting model parameters by the fractional occurrence of the individual units, referred to as parameter upscaling. The advantage of parameter upscaling is that it is much more computationally efficient: a requirement for hemispheric scale simulations. In this study we determined the differences between modelling a raised bog peatland with hummock-hollow microtopography using flux and parameter upscaling. We used the McGill Wetland Model (MWM), a process-based ecosystem C model for peatlands, configured for hummocks and hollows separately and then a weighted mixture of both. The simulated output based on flux and parameter upscaling was compared with eddy-covariance tower measurements. We found that net ecosystem production (NEP) for hollows was much larger than that for hummocks because total ecosystem respiration (TER) for hummocks was greater while gross primary production (GPP) did not differ significantly between the two topographic features. However, despite differences in components of NEP between hummocks and hollows, there was no statistically significant difference between the NEP based on flux and parameter upscaling using the MWM. Both flux and parameter upscaling show equivalent capability to capture the magnitude, direction, seasonality and inter-annual variability. The root-mean-square-errors (RMSE) are 0.66, 0.45, and 0.49 g C m⁻² day⁻¹, respectively for GPP, TER and NEP based on the flux upscaling, while 0.67, 0.44, and 0.48 g C m⁻² day⁻¹, respectively based on the parameter upscaling. The degree of agreement (d*) is 0.96, 0.97, and 0.88, respectively for GPP, TER and NEP based on the flux upscaling, while 0.96, 0.97, and 0.89, respectively based on the parameter upscaling. This result suggests that differences in processes caused by peatland microtopography scale linearly, which means an ecosystem-level model set-up (i.e. parameter upscaling scheme), is sufficient to simulate the C cycling.     

Modeling carbon sequestration under zero-tillage at the regional scale. II. The influence of crop rotation and soil type

For policy decisions with respect to CO₂-mitigation measures in the agricultural sector, national and regional estimations of the efficiency of such measures are required. The conversion of ploughed cropland to zero-tillage is discussed as an option to reduce CO₂ emissions and promises at the same time effective soil and water conservation. Based on the upscaling of simulation results with the soil and land resources information system SLISYS-BW, estimations of CO₂-mitigation rates in relation to crop rotations and soil type have been made for the state of Baden-Württemberg (Germany). The results indicate considerable differences in the CO₂-mitigation rates between crop rotations ranging from 0.48 to 0.03 Mg C ha⁻¹ a⁻¹ for winter cereals–spring cereals–rape rotations and winter cereals–spring cereals–corn silage rotations, respectively. The efficiency of the crop rotations is strongly related to the total carbon input and in particular the amount of crop residues. Among the considered soil types, highest CO₂-mitigation rates are associated with Cumulic Anthrosols (0.62 Mg C ha⁻¹ a⁻¹) and the lowest with Gleysols (−0.01 Mg C ha⁻¹ a⁻¹). An agricultural extensification scenario with conventional plowing but conversion of the presently applied intensive crop rotations to a clover–clover–winter cereals rotation indicated a CO₂-mitigation potential of 466 Gg C a⁻¹. However, the present high market prices for cereals and increasing demand for energy production from biomass encourages an intensification of the agricultural production and an excessive removal of biomass which in future will seriously reduce the potential for carbon sequestration on cropland.     

Carbon, nitrogen, oxygen and sulfide budgets in the Black Sea: A biogeochemical model of the whole water column coupling the oxic and anoxic parts

Carbon, nitrogen, oxygen and sulfide budgets are derived for the Black Sea water column from a coupled physical-biogeochemical model. The model is applied in the deep part of the sea and simulates processes over the whole water column including the anoxic layer that extends from ?115 m to the bottom (?2000 m). The biogeochemical model involves a refined representation of the Black Sea foodweb from bacteria to gelatinous carnivores. It includes notably a series of biogeochemical processes typical for oxygen deficient conditions with, for instance, bacterial respiration using different types of oxidants (i.e denitrification, sulfate reduction), the lower efficiency of detritus degradation, the ANAMMOX (ANaerobic AMMonium OXidation) process and the occurrence of particular redox reactions. The model has been calibrated and validated against all available data gathered in the Black Sea TU Ocean Base and this exercise is described in Gregoire et al. (2008). In the present paper, we focus on the biogeochemical flows produced by the model and we compare model estimations with the measurements performed during the R.V. KNORR expedition conducted in the Black Sea from April to July 1988 (Murray and the Black Sea Knorr Expedition, 1991). Model estimations of hydrogen sulfide oxidation, metal sulfide precipitation, hydrogen sulfide formation in the sediments and water column, export flux to the anoxic layer and to the sediments, denitrification, primary and bacterial production are in the range of field observations.With a simulated Gross Primary Production (GPP) of 7.9 mol C m⁻² year⁻¹ and a Community Respiration (CR) of 6.3 mol C m⁻² year⁻¹, the system is net autotrophic with a Net Community Production (NCP) of 1.6 mol C m⁻² year⁻¹. This NCP corresponds to 20% of the GPP and is exported to the anoxic layer. In order to model Particulate Organic Matter (POM) fluxes to the bottom and hydrogen sulfide profiles in agreement with in situ observations, we have to consider that the degradation of POM in anoxic conditions is less efficient that in oxygenated waters as it has often been observed (see discussion in Hedges et al., 1999). The vertical POM profile produced by the model can be fitted to the classic power function describing the oceanic carbon rate (C_R=Z^−α) using an attenuation coefficient α of 0.36 which is the value proposed for another anoxic environment (i.e. the Mexico Margin) by Devol and Hartnett (2001). Due to the lower efficiency of detritus degradation in anoxic conditions and to the aggregation of particles that enhanced the sinking, an important part of the export to the anoxic layer (i.e. 33%, 0.52 mol C m⁻² year⁻¹) escapes remineralization in the water column and reaches the sediments. Therefore, sediments are active sites of sulfide production contributing to 26% of the total sulfide production.In the upper layer, the oxygen dynamics is mainly governed by photosynthesis and respiration processes as well as by air-sea exchanges. ?71% of the oxygen produced by phytoplankton (photosynthesis+nitrate reduction) is lost through respiration, ?21% by outgasing to the atmosphere, ?5% through nitrification and only ?2% in the oxidation of reduced components (e.g. Mn²⁺, Fe²⁺, H₂S).The model estimates the amount of nitrogen lost through denitrification at 307 mmol N m⁻² year⁻¹ that can be partitioned into a loss of ?55% through the use of nitrate for the oxidation of detritus in low oxygen conditions, ?40% in the ANAMMOX process and the remaining ?5% in the oxidation of reduced substances by nitrate.In agreement with data analysis performed on long time series collected since the 1960s (Konovalov and Murray, 2001), the sulfide and nitrogen budgets established for the anoxic layer are not balanced in response to the enhanced particle fluxes induced by eutrophication: the NH₄ and H₂S concentrations increase.