A methodological framework for integrating computational fluid dynamics and ecological models applied to juvenile freshwater mussel dispersal in the Upper Mississippi River

Interdisciplinary research in hydraulics and ecology for river management and restoration must integrate processes that occur over a wide range of spatial and temporal scales, which presents a challenge to ecohydraulics modelers. Computational fluid dynamics (CFD) models are being more widely used to determine flow fields for ecohydraulics applications. In the Upper Mississippi River (UMR), the mussel dynamics model was developed as a tool for management and conservation of freshwater mussels (Unionidae), which are benthic organisms, imperiled in North America, that are inextricably linked with the hydraulics of river flow. We updated the juvenile dispersal component of the mussel dynamics model by using stochastic Lagrangian particle tracking in a three dimensional flow field output from CFD models of reaches in the UMR. We developed a methodological framework to integrate hydrodynamic data with the mussel dynamics model, and we demonstrate the use of the juvenile dispersal model employed within the methodological framework in two reaches of the UMR. The method was used to test the hypothesis that impoundment affects the relationship of some hydraulic parameters with juvenile settling distribution. Simulation results were consistent with this hypothesis, and the relationships of bed shear stress and Froude number with juvenile settling were altered by impoundment most likely through effects on local hydraulics. The methodological framework is robust, integrates Eulerian and Lagrangian reference frameworks, and incorporates processes over a wide range of temporal and spatial scales, from watershed scale hydrologic processes (decades), to reach scale (km) processes that occur over hours or days, and turbulent processes on spatial scales of meter to millimeter and times scales of seconds. The methods are presently being used to assess the impacts of pre- and early post-settlement processes on mussel distributions, including the effects of bed shear stress, and the sensitivity of the location of the host fish when juveniles excyst, on juvenile settling distribution.     

Is habitat amount important for biodiversity in rocky shore systems? A study of South African mussel assemblages

Habitat-forming species on rocky shores are often subject to high levels of exploitation, but the effects of subsequent habitat loss and fragmentation on associated species and the ecosystem as a whole are poorly understood. In this study, the effects of habitat amount on the fauna associated with mussel beds were investigated, testing for the existence of threshold effects at small landscape scales. Specifically, the relationships between mussel or algal habitat amount and: associated biodiversity, associated macrofaunal abundance and density of mussel recruits were studied at three sites (Kidd’s Beach, Kayser’s Beach and Kini Bay) on the southern and south-eastern coasts of South Africa. Samples, including mussel-associated macrofauna, of 10 × 10 cm were taken from areas with 100 % mussel cover (Perna perna or a combination of P. perna and Mytilus galloprovincialis) at each site. The amount of habitat provided by mussels and algae surrounding the sampled areas was thereafter determined at the 4.0 m² scale. A number of significant positive relationships were found between the amount of surrounding mussel habitat and the abundances of several taxa (Anthozoa, Malacostraca and Nemertea). Likewise, there were positive relationships between the amount of surrounding algal habitat and total animal abundance as well as abundance of mussel recruits at one site, Kini Bay. In contrast, abundance of mussel recruits showed a significant negative relationship with the amount of mussel habitat at Kayser’s Beach. Significant negative relationships were also detected between the amount of mussel habitat and species richness and total abundance at Kidd’s Beach, and between amount of mussel habitat and the abundance of many taxa (Bivalvia, Gastropoda, Maxillopoda, Ophiuroidea, Polychaeta and Pycnogonida) at all three sites. No threshold effects were found, nor were significant relationships consistent across the investigated sites. The results indicate that the surrounding landscape is important in shaping the structure of communities associated with these mussel beds, with significant effects of the amount of surrounding habitat per se. The strength and the direction of habitat effects vary, however, between shores and probably with the scale of observation as well as with the studied dependent variables (e.g. diversity, abundance, mussel recruitment, species identity), indicating the complexity of the processes structuring macrofaunal communities on these shores.     

Influence of grain size on species-habitat models

High resolution remote sensing data facilitate the use of small-scale habitat features such as trees or hedges in the analysis of species-habitat relationships. Such data potentially enable more accurate species-habitat mapping than lower resolution data. Here, for the first time, we systematically investigated this hypothesis by altering the spatial resolution from 1 m up to 1000 m grain size in species-habitat models of 13 bird species. The study area covered the Nidda river catchment in central Germany, a large heterogeneous landscape of 1620 km². A high resolution habitat map of the area was converted to coarser spatial and thematic resolutions in seven steps. We investigated how model performance responded to grain size, and we compared the differential effects of spatial resolution and thematic resolution on model performance. Explained deviance (D²) of the bird models generally decreased with coarser spatial resolution of the data, although it did not decrease monotonically in all species. On average across all species, model D² decreased from 41.5 at 1 m grain size to 15.9 at 1000 m grain size. Ten species were best modelled at 1 m, two species at 3 m and one species at 32 m grain size. Model performance degraded continuously with increasing grain size, both in habitat generalist and habitat specialist bird species, and was systematically lower in habitat generalists. The higher model performance observed at finer grain sizes was most likely caused by the combination of three factors: (1) high spatial accuracy of bird records and (2) a more precise location and delineation of habitat features and, (3) to a lesser degree, by more habitat types differentiated in maps of finer resolution. We conclude that higher spatial and thematic resolution data can be essential for deriving accurate predictions on bird distribution patterns from species-habitat models. Especially for bird species that are sensitive to specific land-use types or to small-scaled habitat features, a grain size of 1-3 m seems most promising.     

Assessing spawning ground hydraulic suitability for Chinese sturgeon (Acipenser sinensis) from horizontal mean vorticity in Yangtze River

The Chinese sturgeon, Acipenser sinensis, is an anadromous protected species. Its migration pattern has been blocked since the construction of Gezhouba dam and the length of the natural spawning site reduced to less than 7 km. However, the fish eventually established an alternative spawning ground in the narrow downstream reach of Gezhouba dam. In this article, we applied Delft3D-Flow model to simulate the hydraulic suitability of the spawning ground downstream Gezhouba dam. Horizontal mean vorticity was used to assess the hydraulic environment of spawning ground. The flow field state was determined through model simulation and field-measured data used to validate the model. The computational method was improved by calculating absolute horizontal mean vorticity from estimates the literature. The final vorticity was determined from the simulation output and its distribution pattern retrieved. The horizontal mean vorticity range was 0.71–4.61 10⁻³ s⁻¹ for the entire spawning grounds, with egg mass field upper limit of 1.0 × 10⁻³ s⁻¹. Vorticity strength selection of Chinese sturgeon spawning can enhance our understanding of egg fertilization rate, hence the protection of fertilized eggs. Furthermore, the results of the study would add to existing scientific database for spawning ground hydraulic environment protection, especially in the ecological regulation drive of the Three Gorges Reservoir.     

Two-dimensional habitat modeling of Chinese sturgeon spawning sites

Since the construction of the Gezhouba Dam in the 1980s, the number of Chinese sturgeon in the Yangtze River has been rapidly declining. The Gezhouba Dam has cutoff the migration path of these sturgeon, resulting in an overall reduction of suitable sturgeon habitat. This paper describes a habitat suitability index model that is used to evaluate the impacts of the Gezhouba Dam and Three Gorges Project on Chinese sturgeon spawning sites. Based on research concerning the reproduction characteristics of Chinese sturgeon, ten major ecological factors that influence reproduction were analyzed, including: water temperature, velocity, water depth, substrate, suspended sediment concentration, and the amount of egg predatory fish. The suitability index (SI) curves based on these ten ecological factors were obtained, and a habitat suitability function was developed. A two-dimensional mathematical model was also created to simulate and predict physical habitat situation (such as hydraulic, sediment, and substrate) of the Chinese sturgeon. By coupling the habitat suitability function and a two-dimensional mathematical model, a habitat suitability index model for Chinese sturgeon was established. The habitat suitability index model was validated by comparing measured data with predictions from the model. These comparisons showed that the computed results agreed well with the measured results, and the high calculated habitat suitability index (HSI) corresponded to high measured quantity of eggs per unit (1000 m³) discharge (CPUE_d). The calculated habitat suitability index for Chinese sturgeon also showed that the habitat suitability index was better in 1999, before the impoundment of the Three Gorges Project, compared with the habitat suitability in 2003. Simulation results of different discharges from Gezhouba Dam predicted that flow discharges between 10,000 and 30,000 m³/s were most suitable for sturgeon spawning.     

Modeling amphibian energetics, habitat suitability, and movements of western toads, Anaxyrus (=Bufo) boreas, across present and future landscapes

Effective conservation of amphibian populations requires the prediction of how amphibians use and move through a landscape. Amphibians are closely coupled to their physical environment. Thus an approach that uses the physiological attributes of amphibians, together with knowledge of their natural history, should be helpful. We used Niche Mapper™ to model the known movements and habitat use patterns of a population of Western toads (Anaxyrus (=Bufo) boreas) occupying forested habitats in southeastern Idaho. Niche Mapper uses first principles of environmental biophysics to combine features of topography, climate, land cover, and animal features to model microclimates and animal physiology and behavior across landscapes. Niche Mapper reproduced core body temperatures (T_c) and evaporation rates of live toads with average errors of 1.6 ± 0.4 °C and 0.8 ± 0.2 g/h, respectively. For four different habitat types, it reproduced similar mid-summer daily temperature patterns as those measured in the field and calculated evaporation rates (g/h) with an average error rate of 7.2 ± 5.5%. Sensitivity analyses indicate these errors do not significantly affect estimates of food consumption or activity. Using Niche Mapper we predicted the daily habitats used by free-ranging toads; our accuracy for female toads was greater than for male toads (74.2 ± 6.8% and 53.6 ± 15.8%, respectively), reflecting the stronger patterns of habitat selection among females. Using these changing to construct a cost surface, we also reconstructed movement paths that were consistent with field observations. The effect of climate warming on toads depends on the interaction of temperature and atmospheric moisture. If climate change occurs as predicted, results from Niche Mapper suggests that climate warming will increase the physiological cost of landscapes thereby limiting the activity for toads in different habitats.     

Application of least-cost path model to identify a giant panda dispersal corridor network after the Wenchuan earthquake—Case study of Wolong Nature Reserve in China

With growing levels of human-activity and frequent natural disturbances throughout the world, it is increasingly important that both research and management efforts take into account the widespread landscape fragmentation and its consequences for biodiversity conservation. The magnitude 5.12 Wenchuan earthquake in China caused dramatic impacts on giant panda (Ailuropoda melanoleuca) habitat in the nature reserves within Minshan and Qionglai mountains. With the combined stresses of the natural disaster and the extensive human activities during postquake reconstruction, giant panda habitat in this region may become more fragmented in the future. In order to preserve the giant panda population after the earthquake and protect the species against habitat fragmentation, this article explores a method of identifying giant panda migration corridors involving habitat suitability assessments and a least-cost path model. Focusing on postquake Wolong Nature Reserve, our results demonstrate that it contains 430.3 km² of suitable habitat (21.1% of total area), 463.8 km² of marginally suitable habitat (22.8%) and 1141.9 km² of unsuitable habitat (50.1%). We further show that several giant panda dispersal corridors exist in the reserve, including four corridor groups that cross the provincial highway and five corridors that do not intersect areas of human activity. This study will contribute to management and conservation efforts in Wolong Nature Reserve and beyond after the Wenchuan earthquake.     

Modeling ecosystem responses to prescribed fires in a phosphorus-enriched Everglades wetland: I. Phosphorus dynamics and cattail recovery

We have developed and applied a process-based model, the Wetland Ecosystem Model (WEM), to evaluate the effects of a prescribed fire on the phosphorus (P) dynamics and cattail (Typha domingensis) growth in a P-enriched area in the Florida Everglades. The WEM couples major ecosystem processes including carbon (C), nitrogen (N) and P biogeochemical cycles, plant growth, hydrology, and fire disturbance. The model is used to assess the effects of a prescribed fire on P dynamics and cattail growth through dynamic interaction among four modules: fire, water chemistry, soil, and vegetation. The simulation results are in agreement with observed data including cattail above- and belowground biomass and dead mass, P concentration in surface-water, pore-water, and soil, and soil and water temperature. Cattail aboveground biomass reached the unburned level one year after burn; belowground biomass recovered to unburned level one and half years after the fire, however, dead mass did not completely reach unburned level two years after fires. The fire increased water and soil temperatures in the short term, while indirectly increasing the sensitivity of water and soil temperature post-fire response to air temperature by altering the energy exchange between air and water through a canopy gap created by fire. The fire also altered the P dynamics in surface-water and pore-water. A post-fire P pulse that lasted for less than one month was observed in surface-water. A similar P pulse, but in a small magnitude and a longer duration, was also observed in the pore-water total phosphorus (TP), and then came back to normal level after approximately three months. No significant changes in soil TP was observed during the study period. Meanwhile, no significant changes in water nutrients were observed downstream of the study plot. This finding indicated that the P-enriched wetlands in Everglades act as a buffer in regulating the P concentration in surface-water. Our study showed that the distance of fire effects on a 300 m × 300 m plot was less than 300 m downstream. Sensitivity analysis identified that the air temperature and hydrological conditions are two important driving factors which may alter the cattail community dynamics in response to prescribed fires. Similar to the filed studies, this study provided evidences that fire played an important role in managing plant growth and P dynamics in the Florida Everglades.     

A spatial model of white sturgeon rearing habitat in the lower Columbia River,USA

Concerns over the potential effects of in-water placement of dredged materials prompted us to develop a GIS-based model that characterizes in a spatially explicit manner white sturgeon Acipenser transmontanus rearing habitat in the lower Columbia River, USA. The spatial model was developed using water depth, riverbed slope and roughness, fish positions collected in 2002, and Mahalanobis distance (D²). We created a habitat suitability map by identifying a Mahalanobis distance under which >50% of white sturgeon locations occurred in 2002 (i.e., high-probability habitat). White sturgeon preferred relatively moderate to high water depths, and low to moderate riverbed slope and roughness values. The eigenvectors indicated that riverbed slope and roughness were slightly more important than water depth, but all three variables were important. We estimated the impacts that fill might have on sturgeon habitat by simulating the addition of fill to the thalweg, in 3-m increments, and recomputing Mahalanobis distances. Channel filling simulations revealed that up to 9 m of fill would have little impact on high-probability habitat, but 12 and 15 m of fill resulted in habitat declines of ∼12% and ∼45%, respectively. This is the first spatially explicit predictive model of white sturgeon rearing habitat in the lower Columbia River, and the first to quantitatively predict the impacts of dredging operations on sturgeon habitat. Future research should consider whether water velocity improves the accuracy and specificity of the model, and to assess its applicability to other areas in the Columbia River.     

Modelling habitat overlap among sympatric mesocarnivores in southern Illinois,USA

Few researchers have developed large-scale habitat models for sympatric carnivore species. We created habitat models for red foxes (Vulpes vulpes), coyotes (Canis latrans) and bobcats (Lynx rufus) in southern Illinois, USA, using the Penrose distance statistic, remotely sensed landscape data, and sighting location data within a GIS. Our objectives were to quantify and spatially model potential habitat differences among species. Habitat variables were quantified for 1-km² buffered areas around mesocarnivore sighting locations. Following variable reduction procedures, five habitat variables (percentage of grassland patches, interspersion–juxtaposition of forest patches, mean fractal dimension of wetland patches and the landscape, and road density) were used for analysis. Only one variable differed (P < 0.05) between red fox and coyote sighting areas (road density) and bobcat and coyote sighting areas (mean fractal dimension of the landscape). However, all five variables differed between red fox and bobcat sighting areas, indicating considerable differences in habitat affiliation between this pair-group. Compared to bobcats, red fox sightings were affiliated with more grassland cover and larger grassland patches, higher road densities, lower interspersion and juxtaposition of forest patches, and lower mean fractal dimension of wetland patches. These differences can be explained by different life history requirements relative to specific cover types. We then used the Penrose distance statistic to create habitat models for red foxes and bobcats, respectively, based on the five-variable dataset. An independent set of sighting locations were used to validate these models; model fit was good with 65% of mesocarnivore locations within the top 50% of Penrose distance values. In general, red foxes were affiliated with mixtures of agricultural and grassland cover, whereas bobcats were associated with a combination of grassland, wetland, and forest cover. The greatest habitat overlap between red foxes and bobcats was found at the interface between forested areas and more open cover types. Our study provides insight into habitat overlap among sympatric mesocarnivores, and the distance-based modelling approach we used has numerous applications for modelling wildlife–habitat relationships over large scales.     

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.     

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.     

Carbon sequestration in soils of SW-Germany as affected by agricultural management—Calibration of the EPIC model for regional simulations

Global emissions trading allows for agricultural measures to be accounted for the carbon sequestration in soils. The Environmental Policy Integrated Climate (EPIC) model was tested for central European site conditions by means of agricultural extensification scenarios. Results of soil and management analyses of different management systems (cultivation with mouldboard plough, reduced tillage, and grassland/fallow establishment) on 13 representative sites in the German State Baden-Württemberg were used to calibrate the EPIC model. Calibration results were compared to those of the Intergovernmental Panel on Climate Change (IPCC) prognosis tool. The first calibration step included adjustments in (a) N depositions, (b) N₂-fixation by bacteria during fallow, and (c) nutrient content of organic fertilisers according to regional values. The mixing efficiency of implements used for reduced tillage and four crop parameters were adapted to site conditions as a second step of the iterative calibration process, which should optimise the agreement between measured and simulated humus changes. Thus, general rules were obtained for the calibration of EPIC for different criteria and regions. EPIC simulated an average increase of +0.341 Mg humus-C ha⁻¹ a⁻¹ for on average 11.3 years of reduced tillage compared to land cultivated with mouldboard plough during the same time scale. Field measurements revealed an average increase of +0.343 Mg C ha⁻¹ a⁻¹ and the IPCC prognosis tool +0.345 Mg C ha⁻¹ a⁻¹. EPIC simulated an average increase of +1.253 Mg C ha⁻¹ a⁻¹ for on average 10.6 years of grassland/fallow establishment compared to an average increase of +1.342 Mg humus-C ha⁻¹ a⁻¹ measured by field measurements and +1.254 Mg C ha⁻¹ a⁻¹ according to the IPCC prognosis tool. The comparison of simulated and measured humus C stocks was r² ≥ 0.825 for all treatments. However, on some sites deviations between simulated and measured results were considerable. The result for the simulation of yields was similar. In 49% of the cases the simulated yields differed from the surveyed ones by more than 20%. Some explanations could be found by qualitative cause analyses. Yet, for quantitative analyses the available information from farmers was not sufficient. Altogether EPIC is able to represent the expected changes by reduced tillage or grassland/fallow establishment acceptably under central European site conditions of south-western Germany.     

Modelling potential habitat for cougars in midwestern North America

Cougars (Puma concolor) are of considerable interest to wildlife biologists and the general public in midwestern North America, yet no researchers have modelled potential habitat in the region. We created a model of potential cougar habitat in 9 midwestern states using geospatial data, expert-opinion surveys, the analytical hierarchy process, and a GIS. About 8% of the study region contained highly favorable habitat (with favorability scores ≥75%) for cougars; the states of Arkansas (19%) and Missouri (16%) contained the highest proportions of potentially favorable habitat. We identified 6 large (≥2500 km² in size), contiguous areas of highly favorable habitat for cougars. Model testing procedures indicated a valid model when compared to an independent set of cougar locations, a null dataset, and similar studies. Our model is useful as a planning tool to proactively address future human-cougar conflicts should cougars re-colonize the Midwest via subadult dispersal.     

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.     

Effects of mussel aquaculture on the nitrogen cycle and benthic communities in Kenepuru Sound,Marlborough Sounds,New Zealand

Nitrogen pools and transformations and benthic communities at a Perna canaliculus farm and a nearby reference site without direct influence of marine farming in Kenepuru Sound, New Zealand, were compared on four dates between September 1982 and May 1983. The organic nitrogen pool in the top 12 cm sediment was 7.4 to 10.8 mol m^-2 at the mussel farm and 6.1 to 8.9 mol m^-2 at the reference site. The nitrate and nitrite pools were similar in both sediments, but the ammonium pool in the mussel farm sediment was about twice as high as in the reference sediment. In January, the sediment ammonium concentrations ranged from 418 nmol cm^-3 (surface) to 149 nmol cm^-3 (12 cm depth) at the mussel farm and from 86 to 112 nmol cm^-3 at the reference site. The molar C:N ratio of the sediment organic matter was 6.2 to 7.2 at the mussel farm and 7.9 to 10.0 at the reference site. The molar N:P ratio of the sediment organic matter was 4.3 to 7.2 and 3.3 to 6.1 at mussel farm and reference site, respectively. The total nitrogen mineralisation rate in the top 12 cm sediment ranged from 21.7 to 37.1 mmol m^-2 d^-1 at the mussel farm and from 8.5 to 25.0 mmol m^-2 d^-1 at the reference site. Ammonium excretion by mussels was about 4.7% (January) and 7.4% (May) of the combined nitrogen mineralisation by mussels and sediment. The sediment-denitrification rate was 0.7 to 6.1 mmol m^-2 d^-1 at the mussel farm and 0.1 to 0.9 mmol m^-2 d^-1 at the reference site. In January, 76 and 93% of the nitrate reduced in the sediments were denitrified at the mussel farm and reference site, respectively. The denitrification rate on the mussel lines (determined on detritus-covered mussels) was twice the mussel farm sediment-denitrification rate and 10 times the reference sediment-denitrification rate. Total denitrification at the mussel farm was 21% higher than at the reference site. The loss of nitrogen through mussel harvest and denitrification was 68% higher at the mussel farm. The surface layers of both sediments contained about 75 mg m^-2 chlorophyll a. Sediment phaeophytin levels were 52 mg m^-2 at the reference site and 137 mg m^-2 at the mussel farm. While the benthic infauna of the mussel-farm sediment consisted only of polychaete worms, the reference sediment contained also bivalve molluscs, brittle stars and crustaceans.     

Residual levels and health risk of polycyclic aromatic hydrocarbons in freshwater fishes from Lake Small Bai-Yang-Dian, Northern China

The residual levels of polycyclic aromatic hydrocarbons (PAHs) in the liver, brain, gill and muscle tissues of four common edible freshwater fish species including crucian carp, snakehead fish, grass carp and silver carp collected from Lake Small Bai-Yang-Dian in northern China were measured by GC-MS. The distribution and composition pattern of PAHs in the fish tissues, and the effects of lipid contents in fish tissues and the octanol-water partition coefficient (Kow) of PAHs congeners on them were analyzed. The human health risk of PAHs though fish consumptions was estimated. The following results were obtained: (1) The average residual levels of total PAHs (PAH₁₆) on wet weight base in the different tissues of each fish species ranged from 4.764 to 144.254 ng/g ww. The differences in the average residual levels on wet weight base for PAH₁₆ within four fish species were not statistically significant (P > 0.05); however, these within four fish tissues were statistically significant (P < 0.01). (2) There were very similar distribution patterns of PAH congeners among both the fish tissues and the fish species, as indicated by statistically significant positive interrelationships (R = 0.58-0.97, P < 0.01 or P < 0.05). Low molecular weight (LMW) PAHs predominated the distribution in the fish tissues, accounting for 89.97% of total PAHs. Phe was the most dominant component, according for 37.79% of total PAHs, followed by Ant (18.59%), Flo (12.59%), Nap (10.79%), Fla (9.82%) and Pyr (6.43%). (3) The PAHs residues and distribution in the fish tissues are dependent on both the Kow of PAH congeners and the lipid contents in the fish tissues. There was a significant positive relationship (R = 0.7116, P < 0.0001) between lipid contents and PAHs residual levels. The statistically significant negative relationships (P < 0.05) were found between LogKow and log-transformed PAHs contents on wet weight base for all fish tissues except for the muscle tissue of snakehead fish, the brain and liver tissues of crucian carp. (4) The risk levels of total PAHs were lower than 10⁻⁵ for the muscle tissues of four studied fish species and for the brain tissues of grass carp and snakehead fish; while these were higher than 10⁻⁵ for the brain tissues of crucian carp and silver carp. The risk levels of total PAHs in the liver tissues of four studied fish species except for snakehead fish exceeded 10⁻⁵ for 2-4.5 times. However, the potency equivalent concentration (PEC) of total PAHs in four studied fish tissues were still lower than the maximum permissible BaP limits for crops and baked meat and for plants in the national criterions. The distributions of PAH congeners in fish were well simulated by a level III fugacity model, especially for low molecule weight PAHs.     

Waterfowls habitat modeling: Simulation of nest site selection for the migratory Little Tern (Sterna albifrons) in the Nakdong estuary

This paper aims to find patterns in nest site selection by Little Terns Sterna albifrons, in the Nakdong estuary in South Korea. This estuary is important waterfowl stopover and breeding habitat, located in the middle of the East Asia-Australasian Flyway. The Little Tern is a common species easily observed near the seashore but their number is gradually declining around the world. We investigated their nests and eggs on a barrier islet in the Nakdong estuary during the breeding season (May to June, 2007), and a pattern for the nest site selection was identified using genetic programming (GP). The GP generated a predictive rule-set model for the number of Little Tern nests (training: R² = 0.48 and test: 0.46). The physical features of average elevation, variation of elevation, plant coverage, and average plant height were estimated to determine the influence on nest numbers for Little Tern. A series of sensitivity analyses stressed that mean elevation and vegetation played an important role in nest distribution for Little Tern. The influence of these two variables could be maximized when elevation changed moderately within the sampled quadrats. The study results are regarded as a good example of applying GP to vertebrate distribution patterning and prediction with several important advantages compared to conventional modeling techniques, and can help establish a management or restoration strategy for the species.     

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.