Modeling pesticide fate in a small tidal estuary

The exposure analysis modeling system (EXAMS), a pesticide fate model developed by the U.S. Environmental Protection Agency, was modified to model the fate of the herbicides atrazine and metolachlor in a small tidally dominated estuary (Bath Creek) in North Carolina, USA where freshwater inflow accounts for only 3% of the total flow. The modifications simulated the changes that occur during the tidal cycle in the estuary, scenarios that are not possible with the original EXAMS model. Two models were created within EXAMS, a steady-state model and a time-variant tidally driven model. The steady-state model accounted for tidal flushing by simply altering freshwater input to yield an estuary residence time equal to that measured in Bath Creek. The tidal EXAMS model explicitly incorporated tidal flushing by modifying the EXAMS code to allow for temporal changes in estuary physical attributes (e.g., volume). The models were validated with empirical measurements of atrazine and metolachlor concentrations in the estuary shortly after herbicide application in nearby fields and immediately following a rain event. Both models provided excellent agreement with measured concentrations. The steady-state EXAMS model accurately predicted atrazine concentrations in the middle of the estuary over the first 3 days and under-predicted metolachlor by a factor of 2–3. The time-variant, tidally driven EXAMS model accurately predicted the rise and plateau of both herbicides over the 6-day measurement period. We have demonstrated the ability of these modified EXAMS models to be useful in predicting pesticide fate and exposure in small tidal estuaries. This is a significant improvement and expansion of the application of EXAMS, and given the wide use of EXAMS for surface water quality modeling by both researchers and regulators and the ability of EXAMS to interface with terrestrial models (e.g., pesticide root zone model) and bioaccumulation models, we now have an easily-accessible and widely accepted means of modeling chemical fate in estuaries.     

Large Eddy Simulation of turbulence generated by a weak breaking tidal bore

A tidal bore is a natural and fragile phenomenon, which is of great importance for the ecology of an estuary. The bore development is closely linked with the tidal range and the river mouth shape, and its existence is sensitive to any small change in boundary conditions. Despite their ecological and cultural value, little is known on the flow field, turbulent mixing and sediment motion beneath tidal bores. Indeed, some striking features can be highlighted in two-dimensional simulations, such as large velocity fluctuations and flow recirculation structures. Using Large Eddy Simulation method, the numerical results emphasised the complicated turbulent structures and their unsteadiness under a tidal bore.     

Mechanisms mediating plant distributions across estuarine landscapes in a low-latitude tidal estuary

Understanding of how plant communities are organized and will respond to global changes requires an understanding of how plant species respond to multiple environmental gradients. We examined the mechanisms mediating the distribution patterns of tidal marsh plants along an estuarine gradient in Georgia (USA) using a combination of field transplant experiments and monitoring. Our results could not be fully explained by the "competition-to-stress hypothesis" (the current paradigm explaining plant distributions across estuarine landscapes). This hypothesis states that the upstream limits of plant distributions are determined by competition, and the downstream limits by abiotic stress. We found that competition was generally strong in freshwater and brackish marshes, and that conditions in brackish and salt marshes were stressful to freshwater marsh plants, results consistent with the competition-to-stress hypothesis. Four other aspects of our results, however, were not explained by the competition-to-stress hypothesis. First, several halophytes found the freshwater habitat stressful and performed best (in the absence of competition) in brackish or salt marshes. Second, the upstream distribution of one species was determined by the combination of both abiotic and biotic (competition) factors. Third, marsh productivity (estimated by standing biomass) was a better predictor of relative biotic interaction intensity (RII) than was salinity or flooding, suggesting that productivity is a better indicator of plant stress than salinity or flooding gradients. Fourth, facilitation played a role in mediating the distribution patterns of some plants. Our results illustrate that even apparently simple abiotic gradients can encompass surprisingly complex processes mediating plant distributions.     

Lunar influence on distribution of a calanoid copepod in the water column of a shallow,temperate estuary

Moonlight had a significant influence on the distribution of the copepodPseudodiaptomus hessei in the water column. Copepod samples were taken between March 1986 and March 1987 in the Sundays River estuary. South Africa. Late copepodid stages and adults showed vertical differences in abundance on moonlit nights, with more animals in bottom waters, but not on dark nights. Numbers of adults present in the water column were markedly lower on full moon nights compared with new moon nights. The strongest decrease in abundance was displayed by ovigerous females. These observed distribution patterns support published predation hypotheses.     

Animal-sediment relations and community analysis of a Florida estuary

Five species assemblages of the intertidal infauna of Old Tampa Bay, Florida, USA are identified. Two assemblages are judged to constitute distinct communities, while a third is shown to be an interdigitation of the two communities. Dominance by on species is the prevalent pattern within the assemblages. Numbers of deposit feeders are found to be inversely correlated to that of filter feeders, and both trophic types are found to be correlated to the sediment parameters of median grain size, sorting and skewness. Three transects with three stations each were established along the south side of Courtney Campbell Causeway in Tampa Bay, Florida. A faunal sample (0.4 m²), a sediment sample, and a water sample were taken at each station in September, December, and March of 1968/1969. Sediment samples were wet-seived. Animal samples were reduced to numbers of organisms and biomass per species. Trellis diagrams and correlation tests were generated. Support is shown for the trophic group-amensalims hypothesis, however, the silt-clay fraction is apparently of lesser importance to deposit feeders in Florida sediments than in Buzzards Bay sediments. An attempt is made to relate an analysis of the optimal grain size for filter feeders to the trophic group-amensalism hypothesis. A view of communities as abstractions from continua is more realistic than communities as discrete units.     

Long-term variation of organochlorine residues and assemblages of epibenthic organisms in a shallow North Florida (USA) estuary

A 4-year study (1972–1976) determined long-term trends of organochlorine residues (DDT, DDE, DDD, PCB's, mirex) and trawl-susceptible organisms in a shallow, river-dominated estuary in North Florida (Apalachicola Bay, USA). Moderate levels of such compounds were found in various species prior to the restricted use of DDT in 1972. A subsequent precipitous decline in organochlorine besidues was attributed to decreased upland usage, major flushing of the river basin in early 1973, and various factors associated with estuarine function. No mirex was found in sediments or aquatic organisms. Apparently, the half-life of organochlorines is relatively short in this bay system. Various statistical methods were used to test the relationships of different physico-chemical and biological parameters. During the 4-year study period, seasonal river flow fluctuations dominated water color, turbidity, salinity, nutrients (NO₃), chlorophyll a, and the temporal succession of fishes in the bay. Certain long-term trends of fish associations were noted; relative dominance of key fish species declined and stabilized while bay-wide species richness and diversity increased with time. Qualitative changes in species representation determined the long-term pattern of community variability. This was consonant with a distinctive fish fauna during the first year of sampling. The bay anchovy Anchoa mitchilli was dominant during 9 of the first 12 monts of the project; this influenced the time-related changes in community indices. Temporally clustered fish associations reflected the importance of river flow in the estuarine environment. Direct correlation of fish distribution with the rapid disappearance of organochlorine compounds was complicated by aperiodic natural phenomena such as storms and river fluctuations. Population and community trends appeared consistent with other studies showing similar patterns of dominance of stress-resistant fish populations and related changes in community parameters. In this case, the relatively predictable annual succession of fish associations allowed an appraisal of key forcing functions. Due to the high level of seasonal and annual biological variability in this estuary, there were some problems in the application of linear statistical models to the data base. Although the long-term trend of relative species representation is useful as an index of stress, new techniques are needed to analyze extensive field data so that functions such as trophic interactions are included in the estimation of causal relationships. There are indications that such effects could be related to the impact of organochlorine compounds on estuarine systems.     

Assesment of the response of a shallow macrotidal estuary to changes in hydrological and wastewater inputs through numerical modelling

The aim of this study was to investigate the response to short-term changes in river freshwater discharges and in nutrients loadings (mainly from the treatment of urban wastewater), of the shallow macrotidal Urdaibai estuary (north of Spain), by using numerical tools. A two-dimensional hydrodynamic model and a water quality model were applied to the estuary, in order to better use it as a prediction tool in the study of the effects of variations in hydrodynamic conditions and in waste water inputs. The model was calibrated and verified using data measured under different hydrological conditions (spring and summer). A model calibration was carried out with field data measured during the summer, while the model validation was conducted for spring conditions. The calibration process allowed the model parameter definition, while the model validation permitted the verification of the calibrated parameters under different environmental conditions. The model results were in reasonable agreement with field measurements, in both the calibration and the validation phases. The model showed a significant decrease in phytoplankton concentration with river input increase. A study on the effects of nutrient input reduction from the Gernika Waste Water Treatment Plant (WWTP) was conducted. It showed a decrease in phytoplankton concentration with decreasing levels of nutrient discharge. This reduction was more pronounced in conjunction with the highest river discharge. In that case, a 50% decrease for the elimination of the WWTP discharge was observed.     

Factors influencing community structure and distribution of different life-cycle categories of fishes in shallow waters of a large Australian estuary

Fish were collected at regular intervals over 5 yr (February 1977 to December 1981) from ten shallow-water sites located throughout the lower, middle and upper regions of the large Swan Estuary in temperate southwestern Australia. Analysis of the catch data showed that the total number of species and total density of fishes were both influenced to a greater extent by site and season within the estuary than by year. The number of species and density of fishes within the whole system were greatest during the summer and autumn, when salinities and temperatures were at a maximum, and declined with distance from the estuary mouth. This reflects the trends shown by marine species, which comprise many species that occur only occasionally in the estuary (marine stragglers) and others which enter estuaries regularly and in considerable numbers (marine estuarine-opportunists). The density of marine estuarine-opportunists was also correlated with temperature, reflecting the tendency for the species of this category to congregate in the shallows during the summer and autumn. The similar seasonal aggregations of the single anadromous species and representatives of species that complete their whole life cycle in the estuary were frequently related to spawning. The density of the estuarine category was correlated neither with salinity nor distance from the estuary mouth. The number of estuarine species was also not correlated with distance from the estuary mouth. The density of freshwater species was inversely correlated with salinity and positively correlated with distance from the estuary mouth. The composition of the fish fauna changed progressively through the estuary, with that of the lower estuary being the most discrete. The composition also changed seasonally, particularly in the upper estuary where, during the winter and spring, the volume of freshwater discharge increased greatly and as a consequence the salinity declined markedly. The species diagnostic of the lower estuary were generally marine estuarine-opportunists, whereas those of the upper estuary typically belonged to either the estuarine or anadromous categories. The marine estuarine-opportunistMugil cephalus was, however, also one of the diagnostic species in the upper estuary during the winter and spring.Please address all correspondence and requests for reprints to I. C. Potter at Murdoch University     

Bioavailability and accumulation of trace elements in soils and plants of a highly contaminated estuary (Domingo Rubio tidal channel,SW Spain)

The Domingo Rubio tidal channel (Palos de la Frontera, Huelva, Spain) is an estuary located in the mouth of the Tinto River. The estuary is affected by different sources of pollution (waters of the Tinto River, contaminated with trace elements from the Iberian Pyrite belt, and effluent from the Huelva chemical industrial area). Soil and the most frequent plant species were collected in 2004 and 2006 at six different locations on the estuary. In general, N-Kjedahl, Total Organic Carbon values, salinity and contamination (total trace elements up to 1,000 mg kg⁻¹ As, 6 mg kg⁻¹ Cd, 2,500 mg kg⁻¹ Cu, 1,900 mg kg⁻¹ Pb and 1,300 mg kg⁻¹ Zn) tended to increase downstream of the tidal channel. Soil biochemical properties were not negatively affected either by the high salinity or by trace element contamination. Despite the high values of the trace elements, analysed plant samples showed that Cu was the only metal that could be a serious risk for the food chain.     

Testing the Stochastic Dynamic Methodology (StDM) as a management tool in a shallow temperate estuary of south Europe (Mondego, Portugal)

A long-term monitoring program has been carried out since the early 1990s in the Mondego estuary, on Portugal's west coast, which is presently under heavy human pressure. In this shallow warm-temperate estuary, a significant macroalgal proliferation has been observed, which is a clear sign of nutrient enrichment. As a result of competition with algae, the extension of the seagrass meadows (mainly Zostera noltii) has been reduced. The present paper examined the applicability of a holistic Stochastic Dynamic Methodology (StDM) in predicting the tendencies of trophic key-components (macrophytes, macroalgae, benthic macroinvertebrate and wading birds) as a response to the changes in estuarine environmental conditions. The StDM is a sequential modelling process developed in order to predict the ecological status of changed ecosystems, from which management strategies can be designed. The data used in the dynamic model construction included true gradients of environmental changes and was sampled from January 1993 to September 1995 and from December 1998 to December 2005. The dynamic model developed was preceded by a conventional multivariate statistical procedure performed to discriminate the significant relationships between the selected ecological components. The model validation was based on independent data collected from January 1996 to January 1997 and from February 1999 to April 2000 for all the state variables considered. Overall, the simulation results are encouraging since they seem to demonstrate the StDM reliability in capturing the trophic dynamics of the studied estuary, by predicting the behavioural pattern for the most part of the components selected, with a focus on the Zostera noltii meadows recovery after the implementation of important management measures.     

Eddy covariance CO2 flux measurements in nocturnal conditions: an analysis of the problem

A detailed analysis of the various processes at work in stable boundary layers was made. It pointed out that two main mechanisms may affect eddy covariance measurements in stable conditions and that their impacts were different. On one hand, intermittent turbulence produces strongly nonstationary events during which the validity of turbulent transport and storage measurements is uncertain. On the other hand, during breeze and drainage flow events, significant advection takes place and competes with turbulent flux and storage. Intermittent turbulence questions both the ability of eddy covariance systems to adequately capture turbulent flux and storage and the representativeness of the measurements. Ability of the systems to capture the fluxes could be improved by adapting the averaging time period or the high pass filtering characteristics. However, beyond this, the question of representativeness of the flux remains open as the flux measured during an intermittent turbulence event represents not only the source term, but also the removal of CO2 that built up in the control volume and that cannot be simply related to the source term. In these conditions, the u* discrimination is likely to be insufficient and should be completed with a stationarity criterion. Further research should allow determining better selection criteria. Advection occurs mainly in presence of flows associated with topographical slopes (drainage flows) or with land use changes (breezes). Direct advection measurements were performed at several sites, but the results were shown to be strongly site dependent. A classification based on the general flow pattern and on the source intensity evolution along streamlines was proposed here. Five different patterns were identified that helped to classify the different observations. The classification was found to be a fairly good fit for the observations. This could serve as a tool to better understand and quantify the fluxes at sites subjected to repeatable patterns.     

Multiple determinants of community structure in shallow marsh habitats,Cape Fear River estuary,North Carolina,USA

The nekton of tidal creeks was studied at 17 sampling localities from September 1977 through August 1978, in the Cape Fear River estuary, North Carolina, USA. Prior to these dates, collections were made at 9 stations beginning in January 1977; these data were used to supplement conclusions drawn from the larger effort. Species recruited from the ocean utilized marsh habitats only temporarily and dominated the catches with over 70% of the total abundance. Their distribution was influenced by salinity gradients and to a lesser extent by substrate characteristics. In addition, temporal habitat partitioning with associated size differences of related species played an important role in structuring marsh nekton communities. A clearly defined ecotone was associated with the mesohaline-polyhaline transition zone, in slainities between 14 and 21 S. Numerous marine stenohaline forms were restricted to salinities above 16 S, thus increasing species richness in high salinity marshes. Despite differences in freshwater flows in 1977 and 1978, major features of the various marsh communities (species associations and relative abundances) exhibited little change throughout the Cape Fear estuary, indicating that these communities were relatively persistent in time. Standing crops for ocean-spawned species at the end of the growing season indicated that considerable annual export in the form of living biomass of fish and shellfish takes place from the marshes. Since most individuals of these species return to the ocean in the fall, an important energy link between the marshes and nearshore marine environment is demonstrated.     

高分辨气相色谱/高分辨质谱法测定南极样品中有机氯农药

建立了同位素稀释-高分辨气相色谱/高分辨质谱法(HRGC/HRMS)测定南极土壤、苔藓和地衣样品中23种有机氯农药的分析方法.样品经冷冻干燥、研磨处理后用正己烷∶二氯甲烷(1∶1,V∶V)混合溶剂进行加速溶剂萃取(ASE),萃取液经硅胶-氧化铝层析柱和C18小柱净化后,进HRGC/HRMS检测分析.样品中目标物定量采用平均相对响应因子法,6点标准曲线响应因子的相对标准偏差(RSD)≤20%,方法的线性范围为0.4—800μg·L-1,回收率在62%—101%之间.实际样品分析结果表明,23种OCPs的加标回收率为40%—100%,在土壤、苔藓和地衣样品中的检出限(LODs)分别为0.024—5.01、0.2—12.2、0.020—13.7 pg·g-1,可以满足南极环境样品中有机氯农药的检测分析.     

Blooms of the dinoflagellate Gyrodinium aureolum in a Long Island estuary:

A dense dinoflagellate bloom of Gyrodinium aureolum Hulburt in a shallow temperate zone estuary was monitored during the summers of 1982 and 1983. The bloom was typically extremely localized, its densest part exceeding 1000g chlorophyll a liter^-1 (2x10⁴ cells ml^-1). The bloom began at temperatures between 24.5° and 27°C, existed at as high as 30°C and terminated when water temperature dropped to between 19° and 22°C. The highest specific growth rate measured was 0.90d^-1 (1.3 divisions d^-1) and near the termination of the bloom decreased to 0.28d^-1 (0.4 divisions d^-1). A diel vertical migration of the bloom was observed. A box model analysis, based on division rates, vertical migration and water circulation patterns, indicated that the bloom must move downward at the estuary mouth to maintain itself in the estuary, either by means of a convergence system or by downward swimming. High growth rate, low grazing pressure, and a stratified water column are proposed to stimulate bloom formation. Decreasing growth rate appeared to reduce the intensity of the bloom and finally allowed its disappearance by estuarine flushing and mixing.Communicated by J. M. Shick, Orono     

Catastrophic-like shifts in shallow Turkish lakes: a modeling approach

A generalized logistic model (GLM) was developed to determine occurrence of submerged macrophytes in shallow Lake Eymir, and the model was tested independently on the upstream shallow Lake Mogan using the data collected fortnightly from both lakes during 1997–2002. The independent variables included concentrations of chlorophyll-a (chl-a), suspended solids (SS) and total phosphorus (TP), Secchi disc transparency and z-scores of water levels. The dependent variable was the binary index of submerged plant occurrence. We used bootstrapping to determine the maximum number of epochs to train the model and to execute training when the corrected average cross entropy (c-index) leveled off. The model predicted that SS concentration, z-scores of water levels and TP concentration were the most important variables for determining occurrence of submerged plants. Sensitivity analyses showed that the probability of submerged plant occurrence followed a strong hysterisis response to varying water levels and the concentrations of SS and TP, with the break points being ±50 cm, 12–17 mg l⁻¹ and 200–300 μg l⁻¹, respectively. This observed sensitivity was in accordance with the alternative stable states hypothesis of shallow lakes. For occurrence of submerged plants, chlorophyll-a concentration and Secchi disc transparency had low significance. This was in concert with both relevances of input variables and the field results. The model gave a good definition of the system since the c-index and corrected c-index on the training data were high (0.970 and 0.963, respectively). Testing the model on Lake Mogan produced a c-index of 0.815 with around 80% of the cases being correctly classified. This showed that the model had a high ability to generalize over a spatially independent test set; therefore, it had a great reliability as well. In addition, the predictive power of the model was indeed very high. Consequently, the model captured the relationships between the input and output variables successfully and consistently with alternative stable states hypothesis.     

Local extinction of a foundation species in a hypoxic estuary: integrating individuals to ecosystem

We integrated across individual, population, community, and ecosystem levels to understand the impact of environmental stress by tracking the foundation species Mytilus edulis in the hypoxic estuary Narragansett Bay, Rhode Island, USA. Our initial surveys revealed that the mussels occurred in nine extensive (2-28 ha) dense (814-9943 individuals/m2) subtidal reefs that attracted a diverse suite of predators (sea stars, crabs, gastropods). Hypoxia occurred in the summer of 2001, and a mussel transplant experiment revealed overall reduced growth rates of individuals, and higher mortality rates among larger mussels. At the population level, large decreases in densities and cover of mussels were correlated with dissolved oxygen concentrations, leading to extinction at one site and reductions of over an order of magnitude at others. Within one year, seven of the eight remaining populations were edged to extinction, and the previously extinct population was recolonized. At the community level, a predator exclusion experiment indicated that predation was an unimportant source of mussel mortality during the hypoxic period, in part due to the emigration of sea stars, as predicted by the Consumer Stress Model. However, mussels were too intolerant to hypoxia to have a net benefit from the predation refuge. The seasonal (summer) occurrence of hypoxia allowed sea stars to return following a lag, as predicted by a stress return time model, and the resumption of predation contributed to the subsequent extinction of mussel populations. At the ecosystem level, the initial filtration rate of the mussel reefs was estimated at 134.6 x 10(6) m3/d, equivalent to filtering the volume of the bay 1.3 times during the 26-d average residence time. That function was reduced by >75% following hypoxia. The effect of hypoxia on each level of organization had consequences at others. For example, size-specific mortality and decreased growth of individuals, and reduced filtration capacity of reefs, indicated a loss of the ability of mussels to entrain planktonic productivity and potential to control future eutrophication and hypoxia. Our study quantified patterns of loss and identified pathways within an integrative framework of feedbacks, summarized in a conceptual model that is applicable to similar foundation species subjected to environmental stress.     

Automatic sensitivity analysis of a finite volume model for two-dimensional shallow water flows

Given a numerical model for solving two-dimensional shallow water equations, we are interested in the robustness of the simulation by identifying the rate of change of the water depths and discharges with respect to a change in the bottom friction coefficients. Such a sensitivity analysis can be carried out by computing the corresponding derivatives. Automatic differentiation (AD) is an efficient numerical method, free of approximation errors, to evaluate derivatives of the objective function specified by the computer program, Rubar20 for example. In this paper AD software tool Tapenade is used to compute forward derivatives. Numerical tests were done to show the robustness of the model and to demonstrate the efficiency of these AD-derivatives.     

Verification of the plume rise/dispersion model USPR: plume rise for single stack emissions

Despite proliferation of the use of air pollution models for regulatory application, major discrepancies still occur between models and also between models and observations, especially when oversimplistic models are used. The problem of predicting plume rise (and subsequently ground level concentrations) from a single source is evaluated here in terms of an integral plume rise and dispersion model (USPR) which encompasses both bouyant rise and turbulent spreading; thus avoiding the problems of the concatenation of separate plume rise and dispersion models. The wide range of validity of the USPR model is demonstrated is terms of plume rise by comparison with the highly buoyant GCOS and Kincaid plumes as well as with dense effluents. It is also shown to be in agreement with Briggs' two-thirds law when the restrictions applicable to the latter model are imposed.