Modeling environmental effects on the size-structured energy flow through marine ecosystems. Part 2: Simulations

Numerical simulations using a physiologically-based model of marine ecosystem size spectrum are conducted to study the influence of primary production and temperature on energy flux through marine ecosystems. In stable environmental conditions, the model converges toward a stationary linear log–log size-spectrum. In very productive ecosystems, the model predicts that small size classes are depleted by predation, leading to a curved size-spectrum.It is shown that the absolute level of primary production does not affect the slope of the stationary size-spectrum but has a nonlinear effect on its intercept and hence on the total biomass of consumer organisms (the carrying capacity). Three domains are distinguished: at low primary production, total biomass is independent from production changes because loss processes dominate dissipative processes (biological work); at high production, ecosystem biomass is proportional to primary production because dissipation dominates losses; an intermediate transition domain characterizes mid-production ecosystems. Our results enlighten the paradox of the very high ecosystem biomass/primary production ratios which are observed in poor oceanic regions. Thus, maximal dissipation (least action and low ecosystem biomass/primary production ratios) is reached at high primary production levels when the ecosystem is efficient in transferring energy from small sizes to large sizes. Conversely, least dissipation (most action and high ecosystem biomass/primary production ratios) characterizes the simulated ecosystem at low primary production levels when it is not efficient in dissipating energy.Increasing temperature causes enhanced predation mortality and decreases the intercept of the stationary size spectrum, i.e., the total ecosystem biomass. Total biomass varies as the inverse of the Arrhenius coefficient in the loss domain. This approximation is no longer true in the dissipation domain where nonlinear dissipation processes dominate over linear loss processes. Our results suggest that in a global warming context, at constant primary production, a 2–4 °C warming would lead to a 20–43% decrease of ecosystem biomass in oligotrophic regions and to a 15–32% decrease of biomass in eutrophic regions.Oscillations of primary production or temperature induce waves which propagate along the size-spectrum and which amplify until a “resonant range” which depends on the period of the environmental oscillations. Small organisms oscillate in phase with producers and are bottom-up controlled by primary production oscillations. In the “resonant range”, prey and predators oscillate out of phase with alternating periods of top-down and bottom-up controls. Large organisms are not influenced by bottom-up effects of high frequency phytoplankton variability or by oscillations of temperature.     

Stable carbon and nitrogen isotopes of sinking particles in the eastern bransfield strait (Antarctica)

A time-series sediment trap was deployed at 1,034 m water depth in the eastern Bransfield Strait for a complete year from December 25, 1998 to December 24, 1999. About 99% of total mass flux was trapped during an austral summer, showing distinct seasonal variation. Biogenic particles (biogenic opal, particulate organic carbon, and calcium carbonate) account for about two thirds of annual total mass flux (49.2 g m^-2), among which biogenic opal flux is the most dominant (42% of the total flux). A positive relationship (except January) between biogenic opal and total organic carbon fluxes suggests that these two variables were coupled, due to the surface-water production (mainly diatoms). The relatively low δ¹³C values of settling particles result from effects on C-fixation processes at low temperature and the high CO² availability to phytoplankton. The correspondingly low δ¹⁵N values are due to intense and steady input of nitrates into surface waters, reflecting an unlikely nitrate isotope fractionation by degree of surface-water production. The δ¹⁵N and δ¹³C values of sinking particles increased from the beginning to the end of a presumed phytoplankton bloom, except for anomalous δ¹⁵N values. Krill and the zooplankton fecal pellets, the most important carriers of sinking particles, may have contributed gradually to the increasing δ¹³C values towards the unproductive period through the biomodification of the δ¹³C values in the food web, respiring preferentially and selectively¹²C atoms. Correspondingly, the increasing δ¹⁵N values in the intermediate-water trap are likely associated with a switch in source from diatom aggregates to some remains of zooplankton, because organic matter dominated by diatom may be more liable and prone to remineralization, leading to greater isotopic alteration. In particular, the tendency for abnormally high δ¹⁵N values in February seems to be enigmatic. A specific species dominancy during the production may be suggested as a possible and speculative reason.     

Seasonal changes in food uptake by the sea cucumber Apostichopus japonicus in a farm pond: Evidence from C and N stable isotopes

This study investigated the seasonal changes in carbon (C) and nitrogen (N) stable isotope values of several typical food sources of Apostichopus japonicus in a farm pond, including particulate organic matter (POM), macroalgae, benthic microalgae and animals such as nematode and copepod. The stable isotope technique was used to quantify relative contributions of various sources to the food uptake by A. japonicus. The results showed that significant changes occurred in the C and N stable isotope values of sea cucumber food sources due to the seasonality of micro-or macroalgae prosperity and the fluctuation of environmental conditions. The sea cucumber A. japonicus exhibited corresponding alterations in feeding strategy in response to the changes in food conditions. Calculation with a stable isotope mixing model showed that macroalgae was the principal food source for A. japonicus throughout the 1-yr investigation, with the relative contribution averaging 28.1% - 63.2%. The relative contributions of other food sources such as copepod and nematode, POM, benthic microalgae to the total food uptake by sea cucumber averaged 22.6% - 39.1%, 6.3% - 22.2%, 2.8% - 6.5%, and 2.8% - 4.2%, respectively. Together these results indicated that the seasonal changes in food sources led to the obvious temporal differences in the relative contribution of various food sources utilized by A. japonicus. Such findings provide the basic scientific information for improving the aquaculture techniques of A. japonicus, particularly for optimizing the food environment of A. japonicus culture in farm ponds.     

Towards end-to-end models for investigating the effects of climate and fishing in marine ecosystems

End-to-end models that represent ecosystem components from primary producers to top predators, linked through trophic interactions and affected by the abiotic environment, are expected to provide valuable tools for assessing the effects of climate change and fishing on ecosystem dynamics. Here, we review the main process-based approaches used for marine ecosystem modelling, focusing on the extent of the food web modelled, the forcing factors considered, the trophic processes represented, as well as the potential use and further development of the models. We consider models of a subset of the food web, models which represent the first attempts to couple low and high trophic levels, integrated models of the whole ecosystem, and size spectrum models. Comparisons within and among these groups of models highlight the preferential use of functional groups at low trophic levels and species at higher trophic levels and the different ways in which the models account for abiotic processes. The model comparisons also highlight the importance of choosing an appropriate spatial dimension for representing organism dynamics. Many of the reviewed models could be extended by adding components and by ensuring that the full life cycles of species components are represented, but end-to-end models should provide full coverage of ecosystem components, the integration of physical and biological processes at different scales and two-way interactions between ecosystem components. We suggest that this is best achieved by coupling models, but there are very few existing cases where the coupling supports true two-way interaction. The advantages of coupling models are that the extent of discretization and representation can be targeted to the part of the food web being considered, making their development time- and cost-effective. Processes such as predation can be coupled to allow the propagation of forcing factors effects up and down the food web. However, there needs to be a stronger focus on enabling two-way interaction, carefully selecting the key functional groups and species, reconciling different time and space scales and the methods of converting between energy, nutrients and mass.     

Status of butyltin pollution along the coasts of western Japan in 2001, 11 years after partial restrictions on the usage of tributyltin

Restrictions on the use of tributyltin (TBT) in aquaculture and on boats in coastal regions, except for ocean-going vessels, have been in place in Japan since 1990 due to the strong toxic effects of TBT on marine organisms. However, TBT pollution along the Japanese coasts has been reported after this legislation was enacted. In order to elucidate the present status of contamination by butyltin (BT) compounds, we measured the levels of BTs [monobutyltin (MBT), dibutyltin (DBT) and TBT] in seawater and Caprella spp. samples obtained from the western part of Seto Inland Sea, Uwa Sea and Uranouchi Bay in western Japan during March to September, 2001. Butyltins were detected in more than 90% of the seawater samples (n = 59), with average concentrations of 8.2 ± 9.2 (SD) ng MBT L⁻¹, 3.3 ± 3.0 ng DBT L⁻¹ and 9.0 ± 7.0 ng TBT L⁻¹. Among 41 stations situated on coastal lines, a sufficient number of Caprella organisms for chemical analysis could be collected from only 16 stations. The butyltin concentrations in seawater and Caprella samples from Uwa Sea and Uranouchi Bay, in which the dominant industry in both waters is aquaculture, showed significantly higher than or no significant differences from those samples from the western part of Seto Inland Sea, a major heavy-industry area in Japan. As the TBT concentration in seawater increased, the number of Caprella organisms collected decreased. The mean TBT concentration among the seawater samples was above the estimated lowest observable effect concentration (LOEC) that reduces the survival rate of Caprella danilevskii. Thus, the present study indicates that TBT is still a potential ecological hazard to the survival of marine invertebrates inhabiting coasts along western Japan, even 11 years after the partial ban on TBT usage was implemented.     

Laboratory model tests for cyclic load-induced settlement of a strip foundation on a clayey soil

Laboratory model test results are presented for the cyclic load-induced settlement of a strip foundation supported by a saturated clayey soil. In performing the tests, the foundation was initially subjected to an allowable static load, after which a cyclic load with a frequency of one cycle per second was superimposed on it. The magnitudes of the static load and the amplitude of the cyclic load were varied. Based on the model test results, relationships for the foundation settlement and intensities of the static and cyclic loads are presented.     

Estimation of asymptotic variances of quantiles for the generalized logistic distribution

In this study, the parameter estimations for the 3-parameter generalized logistic (GL) distribution are presented based on the methods of moments (MOM), maximum likelihood (ML), and probability weighted moments (PWM). The asymptotic variances of the MOM, ML, and PWM quantile estimators for the GL distribution are expressed as functions of the sample size, return period, and parameters. A Monte Carlo simulation was performed to verify the derived expressions for variances and covariances between parameters and to evaluate the applicability of the derived asymptotic variances of quantiles for the MOM, ML and PWM methods. The simulation results generally show good agreement with the analytical results estimated from the asymptotic variances of parameters and quantiles when the shape parameter (β) of the GL distribution is between −0.10 and 0.10 for the MOM method and between −0.25 and 0.45 for the ML and PWM methods, respectively. In addition, the actual sample variances and the root mean square error (RMSE) of asymptotic variances of quantiles for various sample sizes, return periods, and shape parameters were presented. In order to evaluate the applicability of the estimation methods to real data and to compare the values of estimated parameter, quantiles, and confidence intervals based on each parameter estimation method, the GL distribution was fitted to the 24-h annual maximum rainfall data at Pohang, Korea.