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.     

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.     

Ecological role of stomatopods (mantis shrimps) and potential impacts of trawling in a marine ecosystem of the southeast coast of India

Ecosystems are balanced by nature and each component in the system has a role in the sustenance of other components. A change in one component would invariably have an effect on others. Stomatopods (mantis shrimps) are common and ecologically important predatory crustaceans in tropical marine waters. The ecological role of mantis shrimps and potential impacts of trawling in a marine ecosystem were estimated using Ecopath with Ecosim (EwE) Version 5.0 software, by constructing a mass balanced Ecopath model of Parangipettai (Porto Novo) ecosystem. Based on fisheries information from the region, 17 ecological groups were defined including stomatopods. Both primary and secondary data on biomass, P/B, Q/B and diet composition were used as basic inputs. The mass balanced model gave a total system throughput of 14,756 t km⁻² year⁻¹. The gross efficiency of 0.000942 indicated higher contribution of lower food chain groups in the fishery though the mean trophic level was 3.08. The immature and developing stage of the ecosystem was indicated by the ratio of total primary production and total respiration (1.832) and the net system production (2643.30 t km⁻² year⁻¹). Key indices (flow to detritus, net efficiency and omnivory index), split mortality rates and mixed trophic impact of different ecological groups were obtained from the model. A flow diagram was constructed to illustrate the trophic interactions, which explained the biomass flows in the ecosystem with reference to stomatopods. Two temporal simulations were made, with 10 year durations in the mass balanced Ecopath model by using ecosim routine incorporated in EwE software. The effect of decrease in biomass of stomatopods in the ecosystem was well defined, in the first run with increase in stomatopod fishing mortality, and the group showed a high positive impact on benthopelagic fish biomass increase (129%). The simulation with increase in trawling efforts resulted in the biomass decline of different ecological groups as elasmobranchs to 1%, stomatopods to 2%, crabs and lobsters to 36%, cephalopods to 63%, mackerel to 78%, and shrimps to 89%. Present study warns stomatopod discards and further increase in trawling efforts in the region and it explained the need for ecosystem based fisheries management practices for the sustainability of marine fisheries.     

Estimation of annual potential evapotranspiration at regional scale based on the effect of moisture on soil respiration

Potential evapotranspiration (PET) is an important component of water cycle. For traditional models derived from the principle of aerodynamics and the surface energy balance, its calculation always includes many parameters, such as net radiation, water vapor pressure, air temperature and wind speed. We found that it can be acquired in an easier way in specific regions. In this study, a new PET model (PETP model) derived from two empirical models of soil respiration was evaluated using the Penman-Monteith equation as a standard method. The results indicate that the PETP model estimation concur with the Penman-Monteith equation in sites where annual precipitation ranges from 717.71 mm to 1727.37 mm (R² = 0.68, p = 0.0002), but show large discrepancies in all sites (R² = 0.07, p = 0.1280). Then we applied our PETP model at the global scale to the regions with precipitation higher than 700 mm using 2.5° CMAP data to obtain the annual PET for 2006. As expected, the spatial pattern is satisfactory overall, with the highest PET values distributed in the lower latitudes or coastal regions, and with an average of 1292.60 ± 540.15 mm year⁻¹. This PETP model provides a convenient approach to estimate PET at regional scales.     

Modeling net primary production of a fast-growing forest using a light use efficiency model

As interest grows in the quantification of global carbon cycles, Light Use Efficiency (LUE) model predictions of the forest net primary production (NPP) are being developed at an accelerating rate. Such models can provide useful predictions at large scales, but evaluating their performance has been difficult. In this study, a remote sensing-based LUE model was established to estimate forest NPP. Using the forest inventory data (FID) from the regional forest inventory survey in China and established allometric biomass equations, we calculated the biomass, the biomass increment, and the NPP of Eucalyptus urophylla (E. urophylla) plantation plots in the forestry jurisdiction of the Leizhou Forestry Bureau, Southern China. The FID-based NPP and the NPP from LUE model predictions were then compared to each other. Results show that the NPP from model predictions at a spatial resolution of 30 m × 30 m varied from 0 to 265 gC/(m² month) and showed regional differences. In addition, the stand age had variable effects on the average individual biomass of the E. urophylla plantation plots. The average individual biomass of the young and mid-age forests increased exponentially and logarithmically with the stand age (R² = 0.9178 and R² = 0.8683), respectively. For young and mid-age E. urophylla plantation plots, the LUE model-predicted NPP was fairly consistent with the FID-based NPP, but the model predictions of the NPP were higher than the estimates from FID. Through the analysis of the causes of uncertainty and the possible reasons for the discrepancy between the model-based NPP and FID-based NPP, the FID-derived estimates provided a foundation for model evaluation.     

A three-isotope approach to disentangling the diet of a generalist consumer: the yellow-legged gull in northwest Spain

The widespread omnivory of consumers and the trophic complexity of marine ecosystems make it difficult to infer the feeding ecology of species. The use of stable isotopic analysis plays a crucial role in elucidating trophic interactions. Here we analysed δ¹⁵N, δ¹³C and δ³⁴S in chick feathers, and we used a Bayesian triple-isotope mixing model to reconstruct the diet of a generalist predator, the yellow-legged gull (Larus michahellis) that breeds in the coastal upwelling area off northwest mainland Spain. The mixing model indicated that although chicks from all colonies were fed with a high percentage of fish, there are geographical differences in their diets. While chicks from northern colonies consume higher percentages of earthworms, refuse constitutes a more important source in the diet of chicks from western colonies. The three-isotope mixing model revealed a heterogeneity in foraging habitats that would not have been apparent if only two stable isotopes had been analysed. Moreover, our work highlights the potential of adding δ³⁴S for distinguishing not only between terrestrial and marine prey, but also between different marine species such as fish, crabs and mussels.     

Interpreting spatial heterogeneity of crop yield with a process model and remote sensing

A process-based crop growth model (Vegetation Interface Processes (VIP) model) is used to estimate crop yield with remote sensing over the North China Plain. Spatial pattern of the key parameter—maximum catalytic capacity of Rubisco (V_cmax) for assimilation is retrieved from Normalized Difference of Vegetation Index (NDVI) from Terra-MODIS and statistical yield records. The regional simulation shows that the agreements between the simulated winter wheat yields and census data at county-level are quite well with R² being 0.41-0.50 during 2001-2005. Spatial variability of photosynthetic capacity and yield in irrigated regions depend greatly on nitrogen input. Due to the heavy soil salinity, the photosynthetic capacity and yield in coastal region is less than 50 μmol C m⁻² s⁻¹ and 3000 kg ha⁻¹, respectively, which are much lower than that in non-salinized region, 84.5 μmol C m⁻² s⁻¹ and 5700 kg ha⁻¹. The predicted yield for irrigated wheat ranges from 4000 to 7800 kg ha⁻¹, which is significantly larger than that of rainfed, 1500-3000 kg ha⁻¹. According to the path coefficient analysis, nitrogen significantly affects yield, by which water exerts noticeably indirect influences on yield. The effect of water on yield is regulated, to a certain extent, by crop photosynthetic capacity and nitrogen application. It is believed that photosynthetic parameters retrieved from remote sensing are reliable for regional production prediction with a process-based model.     

Trophic ecology of large predatory reef fishes: energy pathways, trophic level, and implications for fisheries in a changing climate

Large predatory fishes are disproportionately targeted by reef fisheries, but little is known about their trophic ecology, which inhibits understanding of community dynamics and the potential effects of climate change. In this study, stable isotope analyses were used to infer trophic ecology of a guild of large predatory fishes that are targeted by fisheries on the Great Barrier Reef, Australia. Each of four focal predators (Plectropomus leopardus, Plectropomus maculatus, Lethrinus miniatus and Lutjanus carponotatus) was found to have a distinct isotopic signature in terms of δ¹³C and δ¹⁵N. A two-source mixing model (benthic reef-based versus pelagic) indicated that P. leopardus and L. miniatus derive the majority (72 and 62 %, respectively) of their production from planktonic sources, while P. maculatus and L. carponotatus derive the majority (89 and 74 %, respectively) of their production from benthic reef-based sources. This indicates that planktonic production is important for sustaining key species in reef fisheries and highlights the need for a whole-ecosystem approach to fisheries management. Unexpectedly, there was little isotopic niche overlap between three of four focal predators, suggesting that inter-specific competition for prey may be low or absent. δ¹⁵Nitrogen indicated that the closely related P. leopardus and P. maculatus are apex predators (trophic level > 4), while δ¹³C indicated that each species has a different diet and degree of trophic specialisation. In view of these divergent trophic ecologies, each of the four focal predators (and the associated fisheries) are anticipated to be differentially affected by climate-induced disturbances. Thus, the results presented herein provide a useful starting point for precautionary management of exploited predator populations in a changing climate.     

Shift in foraging grounds and diet broadening during ontogeny in southern elephant seals from Kerguelen Islands

Highly dimorphic species like southern elephant seals (Mirounga leonina, SES hereafter) frequently exhibit resource partitioning according to sex and/or age classes. We measured carbon and nitrogen stable isotopes (δ¹³C and δ¹⁵N) of 404 blood samples (136 males and 268 females from Kerguelen Islands, 49°21′S, 70°18′E) from 2004 to 2011. Assuming that the distribution of carbon isotopes (δ¹³C value) reflects the two main foraging grounds (Polar Frontal and Antarctic Zones), we quantified the proportion of SES foraging within each zone in relation with size, a proxy for their age. We found a clear shift from Polar Frontal to Antarctic waters as male SES aged, but no relation as far as females is concerned. We also observed a widening range of nitrogen isotopic (δ¹⁵N) values, suggesting that both males and females expanded their diet spectra with age. Whereas males increased their trophic level, females remained constant on average, with some adult females feeding both at lower and at higher trophic levels than juveniles.     

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.     

Trophic diversity in amphipods within a temperate eelgrass ecosystem as determined by gut contents and C and N isotope analysis

A conjoint analysis of gut contents and stable C and N isotopes was applied to determine the main food sources and feeding habits of dominant amphipods in an eelgrass bed (Zostera marina) in Gwangyang Bay, Korea. Gut content observations demonstrated that, while Gammaropsis japonicus and Jassa slatteryi are herbivorous, feeding on epiphytes and detritus, Pontogeneia rostrata and Monocorophium acherusicum are omnivorous, feeding on mesozooplankton fragments and detritus. Stable isotope data confirmed that epiphytes, detritus, and mesozooplankton fragments were major food sources for amphipods in the eelgrass bed. Isotopic mixing model calculations clearly showed an interspecific difference in diet composition. A high isotopic dissimilarity between amphipod taxa demonstrated interspecific trophic diversity, reflecting their herbivorous (G. japonicus and J. slatteryi) and omnivorous (P. rostrata and M. acherusicum) feeding habits and confirmed the detrivorous feeding habits of caprellids. Such trophic diversity at interspecific level of the amphipod species indicates that they use different food resources within their microhabitats and play species-specific functional roles as mediators in trophic pathways from producers to higher-level consumers of the eelgrass ecosystem. Finally, our findings suggest that information on the species-specific trophic ecology of amphipods is needed to better understand their potential role in the trophic dynamics and carbon flow of seagrass bed ecosystems.     

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.     

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.     

Development of TRIPLEX-Management model for simulating the response of forest growth to pre-commercial thinning

In order to simulate forest growth response to pre-commercial thinning (PCT), TRIPLEX1.0 - a process-based model designed to predict forest growth as well as carbon (C) and nitrogen (N) dynamics - was modified and improved to also simulate managed forest ecosystem thinning practices. A three-parameter Weibull distribution model was integrated to simulate thinning treatments within the newly developed TRIPLEX-Management model. The thinning intensity component within the model allows users to simulate thinning treatments by applying basal area, stand density and volume to quantify thinning intensity. Natural mortality decreased following thinning due to an increase in growing space for residual stems. Predicted litterfall pools also increased after thinning events took place. The TRIPLEX-Management model was tested against published observational data for Jack Pine (Pinus banksiana Lamb.) stands subjected to PCT in Northwestern Ontario, Canada. The coefficients of determination (R²) between the predicted and observed variables including stand density, mean DBH (diameter at breast height), the quadratic mean DBH, total volume and merchantable volume as well as belowground, aboveground, and total biomass ranged from 0.50 to 0.88 (n = 20, P < 0.001) with the exception of mean tree height (R² = 0.25, n = 20, P < 0.05). Overall, the Willmott index of agreement between predicted and observed variables ranged from 0.97 to 1.00. Results show that the TRIPLEX-Management model is generally capable of simulating growth response to PCT for Jack Pine stands.     

Formulation and testing of a novel river nitrification model

The nitrification process in many river water quality models has been approximated by a simple first order dependency on the water column ammonia concentration, while the benthic contribution has routinely been neglected. In this study a mathematical framework was developed for sediment bed nitrification based on mass transfer theory and Monod bacterial growth kinetics. The model describes ammonia transport across the boundary layer and consumption within the biofilm to quantify the overall nitrification flux. Model results suggest that nitrification is usually controlled by the boundary layer thickness, and we estimated a nitrification velocity range between 0.14 and 0.97 m d⁻¹, assuming typical boundary thicknesses of 0.1–1.0 mm. These ranges compared favorably with reported literature values, including our own measurements. The model was applied to several river systems of different depths where nitrification rates and river depths were available. Assuming that nitrification is exclusively a benthic process, the average velocity of all the rivers evaluated was 0.85 m d⁻¹ (r² = 0.72).     

Determination of δ13C and δ15N and trophic fractionation in jellyfish: implications for food web ecology

Application of stable isotope analysis (SIA) in jellyfish allows definition of trophic patterns not detectable using gut content analysis alone, but analytical protocols require standardization to avoid bias in interpreting isotopic data. We determined δ¹³C and δ¹⁵N in Aurelia sp. from the northern Gulf of Mexico (30°00′N, 89°00′W–30°24′N, 88°00′W) to define differences in stable isotope composition between body parts and whole body, the effect of lipid extraction on δ¹³C in tissues, and fractionation values from medusa to prey. The isotopic composition of bell and whole Aurelia sp. was not different. The increase in δ¹³C values after lipid removal suggested a correction is needed. To aid future analyses, we derived a correction equation from empirical data for jellyfish samples. Laboratory feeding experiments indicated medusae increased +4 ‰ in δ¹³C and +0.1 ‰ in δ¹⁵N compared to their diet. These results suggest protocols commonly applied for other species may be inaccurate to define Aurelia sp. trophic ecology. Because Aurelia spp. are commonly found in marine ecosystems, accurately defining their trophic role by use of SIA has implications for understanding marine food webs worldwide.     

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.     

Local calibration of coliforms parameters of water quality problem at Igapó I Lake, Londrina, Paraná, Brazil

This article presents results concerning the local calibration of the transport parameters (longitudinal and transversal diffusions and decay coefficient) for a two-dimensional problem of water quality at Igapó I Lake, located in Londrina, Paraná, Brazil, using fecal coliforms as an indicator of water quality. The simulation of fecal coliforms concentrations all over the water body is conducted by means of a structured discretization of the geometry of Igapó I Lake, together with the finite difference and finite element methods. By using the velocity field, modeled by the Navier-Stokes and Poisson equations, the flow of fecal coliforms is described by means of a transport model, which considers advective and diffusive processes, as well as a process of fecal coliforms decay. In the checkpoint, the longitudinal and transversal diffusion coefficients and the coliforms decay coefficient that best fitted the value of the fecal coliforms concentration were D_x = D_y = 0.001 m²/h and k = 0.5 d⁻¹ = 0.02083 h⁻¹. A qualitative and quantitative analysis of the numerical simulations conducted in function of the diffusion coefficients and of the coliforms decay parameter provided a better understanding of the local water quality at Igapó I Lake.     

Simulating soil organic matter with CQESTR (v. 2.0): Model description and validation against long-term experiments across North America

Soil carbon (C) models are important tools for examining complex interactions between climate, crop and soil management practices, and to evaluate the long-term effects of management practices on C-storage potential in soils. CQESTR is a process-based carbon balance model that relates crop residue additions and crop and soil management to soil organic matter (SOM) accretion or loss. This model was developed for national use in U.S and calibrated initially in the Pacific Northwest. Our objectives were: (i) to revise the model, making it more applicable for wider geographic areas including potential international application, by modifying the thermal effect and incorporating soil texture and drainage effects, and (ii) to recalibrate and validate it for an extended range of soil properties and climate conditions. The current version of CQESTR (v. 2.0) is presented with the algorithms necessary to simulate SOM at field scale. Input data for SOM calculation include crop rotation, aboveground and belowground biomass additions, tillage, weather, and the nitrogen content of crop residues and any organic amendments. The model was validated with long-term data from across North America. Regression analysis of 306 pairs of predicted and measured SOM data under diverse climate, soil texture and drainage classes, and agronomic practices at 13 agricultural sites having a range of SOM (7.3–57.9 g SOM kg⁻¹), resulted in a linear relationship with an r² of 0.95 (P < 0.0001) and a 95% confidence interval of 4.3 g SOM kg⁻¹. Using the same data the version 1.0 of CQESTR had an r² of 0.71 with a 95% confidence interval of 5.5 g SOM kg⁻¹. The model can be used as a tool to predict and evaluate SOM changes from various management practices and offers the potential to estimate C accretion required for C credits.