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, 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.     

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.     

Numerical simulation for impacts of hydrodynamic conditions on algae growth in Chongqing Section of Jialing River, China

Hydrodynamic conditions are important factors for planktonic algae growth, through introducing two parameters which express the optimal velocity and the velocity range for planktonic algae growth, a new velocity factor was put forward for the formula of growth rate. Therefore, the two-dimensional unsteady ecological dynamic model for algae growth was established to analyze the effects of hydrodynamic conditions on algae growth in Chongqing Reach of Jialing River in China. The temporal and spatial distribution of Chlorophyll-a (Chl-a) concentration was simulated numerically for various water levels, under climate conditions in period of high frequency for algae blooms of Three Gorges Reservoir and nutrition status at present in the research reach. The corresponding locations and areas of likely algae blooms were analyzed and forecasted. The results showed that about 0.04 m s⁻¹ was the optimal velocity for algae growth, and the occurrence of algae blooms in large scale is almost impossible because of relatively high water flow velocity for Jialing River.     

Estimating Plecoglossus altivelis altivelis migration using a mass balance model expressed by hydrological distribution parameters in a major limpid river basin in Japan

We examined the influence of several hydrological and meteorological parameters on the migratory movements of ayu Plecoglossus altivelis altivelis in central Japan. When comprehensively evaluating rivers and ayu behaviour on a catchment scale, the subjects of analysis typically include human activities and hydrological and meteorological phenomena. However, limiting analyses to such factors may be too restrictive when human activities are being conducted. Accordingly, we incorporated a biological viewpoint into the evaluation method, analysing hydrological data (river discharge, river water temperature, sea water temperature) to determine watershed characteristics and examining the relationship between these characteristics and the habitat conditions of ayu. Then we constructed a numerical model for ayu migratory runs that incorporated ayu ecology and watershed characteristics. Analyses of ayu movements from a lower estuarine dam demonstrated that downstream displacements were associated with high water flows of more than 200 m³ s⁻¹ at the beginning of summer. We conclude that it is important to consider the effects of environmental parameters on the movements of different fish species to understand the causes of spatial variation in fish distribution in lowland rivers.     

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.     

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.     

Calculation of the unit emergy value of water in an Italian watershed

Emergy is a thermodynamics-based entity that enables the implementation of a holistic environmental accounting system. It contributes to identify and measure all the inputs (energy and matter) supporting a given system, expressed in a common unit, namely solar emergy joule (sej). The emergy per unit product (called unit emergy value, UEV), is a measure of the environmental cost of a given resource. It is specific of the system/process and gives information on the dynamics, components and functioning of it. This paper presents the emergy evaluation of water resources within the watershed of the river Sieve, located in the Province of Florence (Italy). Along the river, an artificial basin has been created by means of a dam to preserve water quantity and quality, and to protect the Florentine area from dangerous floods and inundations. Different UEVs of water can be identified along the course of the river, especially upstream and downstream of the dam. These values quantify both the environmental and human efforts made to ensure and regulate the presence of water at different points of the river. The UEV of water flowing in the river increases from 1.35 × 10⁵ sej/g upstream, to 5.80 × 10⁵ sej/g downstream of the dam, depending mainly on man-made infrastructure. Along the watershed, three different systems of extraction, purification and distribution of water have been chosen on the basis of their dimension, type and location. UEVs of water distributed and the emergy investment necessary to implement different water management strategies are presented. The value of water purified and distributed decreases from 2.00 × 10⁶ sej/g for the smallest plant in the mountainous area, to 1.72 × 10⁶ sej/g for the largest plant, in the city of Florence, depending on production efficiency.     

Estimating soil carbon turnover using radiocarbon data: A case-study for European Russia

Turnover rates of soil carbon for 20 soil types typical for a 3.7 million km² area of European Russia were estimated based on ¹⁴C data. The rates are corrected for bomb radiocarbon which strongly affects the topsoil ¹⁴C balance. The approach is applied for carbon stored in the organic and mineral layers of the upper 1 m of the soil profile. The turnover rates of carbon in the upper 20 cm are relatively high for forest soils (0.16–0.78% year⁻¹), intermediate for tundra soils (0.25% year⁻¹), and low for grassland soils (0.02–0.08% year⁻¹) with the exception of southern Chernozems (0.32% year⁻¹). In the soil layer of 20–100 cm depth, the turnover rates were much lower for all soil types (0.01–0.06% year⁻¹) except for peat bog soils of the southern taiga (0.14% year⁻¹). Combined with a map of soil type distribution and a dataset of several hundred soil carbon profiles, the method provides annual fluxes for the slowest components of soil carbon assuming that the latter is in equilibrium with climate and vegetation cover. The estimated carbon flux from the soil is highest for forest soils (12–147 gC/(m² year)), intermediate for tundra soils (33 gC/(m² year)), and lowest for grassland soils (1–26 gC/(m² year)). The approach does not distinguish active and recalcitrant carbon fractions and this explains the low turnover rates in the top layer. Since changes in soil types will follow changes in climate and land cover, we suggest that pedogenesis is an important factor influencing the future dynamics of soil carbon fluxes. Up to now, both the effect of soil type changes and the clear evidence from ¹⁴C measurements that most soil organic carbon has a millennial time scale, are basically neglected in the global carbon cycle models used for projections of atmospheric CO₂ in 21st century and beyond.     

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.     

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.     

Assessing the effect of long-term crop cultivation on distribution of Cd in the root zone

Trace elements such as cadmium (Cd) may be inadvertently added to cropland soils through application of fertilizers, irrigation water, and other amendments. These toxic trace elements pose a potentially threat to soil quality and, through the food chain transfer, to human health. A generalized soil trace element mass balance model that accounts for the interactive processes governing the reactions of trace elements in soils, and consequently removed with crop harvest and leaching out of the soil profile with irrigation water was developed in this research. The model conceptually approximates the mechanisms and kinetics of a real field cropland system. The model was used to evaluate the long-term cultivation on distribution of Cd in California croplands. Under typical California cropping practices, Cd added into the soils accumulated primarily in the plow layer while the Cd content below the plow layer was barely affected. After 100 years of continuous cultivation, the soil Cd content of the plow layer increases from the background level 0.22 mg kg⁻¹ to 0.40 mg kg⁻¹. The accumulation of Cd in the plow layer is in proportion to the external inputs and is affected by the soil and plant characteristics, and management practices. The model can be used to evaluate the environmental fates of other toxic element in soils with case specific parameters.     

Resistance and re-organization of an ecosystem in response to biological invasion: Some hypotheses

Using a dynamic model of Lake Chozas developed by Marchi et al. (2011), we tested three hypotheses about recovery of the indigenous community and water quality after radical changes caused by introduction of an invasive allochthonous crayfish, Procambarus clarkii:

1.

Can the lake resist the pressure of an invasive species, like P. clarkii, by adaptation?

2.

Can the ecosystem recover when all the crayfish are removed and low phosphorus concentrations persist in inflow water?

3.

Does the simulated recovery of submerged vegetation occur at a total phosphorus concentration below 100 mg TP m⁻³, as estimated by Scheffer et al. (1993), Scheffer (1997), Jeppesen et al. (1998) and Zhang et al. (2003)?

We obtained the following answers:

1.

Lake Chozas can at least partly resist by adaptation. A combination of possible parameter changes could lead to a significant increase in eco-exergy.

2.

Removal of the phosphorus represented by crayfish (by harvesting) implies complete recovery of the lake and its eco-exergy, albeit not necessarily with the same organisms having the same properties.

3.

The expected hysteresis created by introduction and harvesting of crayfish is observed under the following conditions: phytoplankton dominance at total phosphorus ≥ about 200-250 mg TP m⁻³ and submerged vegetation returns at total phosphorus < 100 mg TP m⁻³.

     

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.     

Increasing accuracy of dispersal kernels in grid-based population models

Dispersal kernels in grid-based population models specify the proportion, distance and direction of movements within the model landscape. Spatial errors in dispersal kernels can have large compounding effects on model accuracy. Circular Gaussian and Laplacian dispersal kernels at a range of spatial resolutions were investigated, and methods for minimizing errors caused by the discretizing process were explored. Kernels of progressively smaller sizes relative to the landscape grid size were calculated using cell-integration and cell-center methods. These kernels were convolved repeatedly, and the final distribution was compared with a reference analytical solution. For large Gaussian kernels (σ > 10 cells), the total kernel error was <10⁻¹¹ compared to analytical results. Using an invasion model that tracked the time a population took to reach a defined goal, the discrete model results were comparable to the analytical reference. With Gaussian kernels that had σ ≤ 0.12 using the cell integration method, or σ ≤ 0.22 using the cell center method, the kernel error was greater than 10%, which resulted in invasion times that were orders of magnitude different than theoretical results. A goal-seeking routine was developed to adjust the kernels to minimize overall error. With this, corrections for small kernels were found that decreased overall kernel error to <10⁻¹¹ and invasion time error to <5%.     

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.     

Towards carrying capacity assessment for aquaculture in the Bolinao Bay, Philippines: A numerical study of tidal circulation

The increase in intensive aquaculture production in the Bolinao Bay, Philippines reached the point of harmful influence to production stock. Up to the present, there has been no estimation of aquaculture carrying capacity which is based on quantification of processes responsible for (a) water quality inside the units, (b) impact to the seabed and (c) water quality in the whole Bolinao Bay. The numerical estimation of tidal circulation, which is the most dominant part of the hydrodynamical regime in the Bolinao Bay, represents an unavoidable step in the carrying capacity determination considering points (a), (b) and (c). The hydrodynamical model we apply is a free surface, 3D finite element tidal model, forced with sea elevation dynamics at three open boundaries. It incorporates an increase in bottom stress drag coefficients in the coral areas and advanced transport corrected advection scheme. The model simulations of water flow show good agreement with measured currents in the central part of the Bolinao Bay, southern and northeastern channel, while agreement in the areas in vicinity of northern channel is not as tight. In order to provide an insight into the water exchange in aquaculture units, to support local water quality models and seabed deposition models (points (a) and (b)), the area is mapped with neap tide mean, spring tide mean and 14-day mean current velocity contours. The highest 14-day mean velocities are attained in the southern channel (>17 cm s⁻¹), while the mean velocities characteristic for northern part (<6 cm s⁻¹) and shallow areas of coral reefs (<2 cm s⁻¹) are much lower. Area-mean difference between spring tide and 14-day mean velocities is estimated to be 18.32% (std=7.31%), while difference between neap tide and 14-day mean velocities is 17.62 % (std=11.19%). To support global basin-wide water quality models (point (c)), retention of water in the bay is estimated by Lagrangian and Eulerian calculation procedure of mean residence time field. Both calculations estimate the highest (no-wind) residence times (Lan: 25.4 days, Eul: 21.03 days) in the central northern part, but most of the area (Lan: 73.11%, Eul: 79.31%) is characterized by residence time values significantly lower than 15 days. The results are readily applicable for upcoming implementation in the models of local (aquaculture units, seabed) and basin-wide (Bolinao Bay) nutrient dynamics, primary and secondary production, organic matter decomposition and oxygen dynamics.     

Rates of nitrification,distribution of nitrifying bacteria,and nitrate fluxes in different types of sediment from Danish waters

Nitrification rates, measured in different sediment types from Danish waters, are in the range of 0.3 to 1.4 mmol NO ₃ ^- m^-2 d^-1. There is no signification between sandy and muddy sediments, nor between shallow and deeper stations. The extent of nitrification is probably limited to the zone of oxygen penetration, 1.5 to 5.5 mm. There are, however, nitrifying bacteria located in the anoxic sediment layers. There relative numbers were found by measuring the nitrification potential of the sediment. These potential rates (22°C) can also be used to calculate actual rates of nitrification, by adjusting to in situ temperature and oxygen penetration. These calculated rates agree with the actual measured rates of nitrification for a wide range of sediment types and may be used for the estimation of actual nitrification rates. Nitrate flux out from the sediment/water interface is in the range of 0 to 0.7 mmol NO ₃ ^- m^-2 d^-1. There is no correlation between concentration gradients of nitrate across the sediment/water interface and the measured flux of nitrate. Approximately 50% of nitrate production is released to the water column. The remainder (0 to 0.35 mmol N m^-2 d^-1) may have been denitrified.     

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.     

Nitrification model with an inhibitory effect of sea water

A dynamic model of nitrification including bacterial growth equations was developed for the understanding of nitrification in estuaries. Sensitivity analysis based on the model revealed that not only maximum specific growth rates of nitrifiers but also yield constants and initial bacterial concentrations have great effect on the nitrification process.Batch experiments were carried out in order to investigate nitrification rates and an inhibitory effect on nitrification by sea water. The maximum specific growth rate of ammonia oxidizing bacteria in fresh water was estimated as 0.04 (h^?1) using a nonlinear regression method. The sea water content of 10%, 25%, 50% and 100% reduced the maximum specific growth rate of ammonia oxidizing bacteria to 0.027, 0.016, 0.0085 and 0.0028 (h^?1) respectively. The effect of sea water on nitrification was interpreted as noncompetitive inhibition. The computed results of inhibition of nitrification based on the noncompetitive inhibition model are presented.