The ecological basis for simulating phytoplankton responses to environmental change (PROTECH)

The theory and mathematical development of a model, called PROTECH, are presented. The model simulates the dynamic responses of up to eight species of phytoplankton to environmental variability in lakes and reservoirs. PROTECH models were developed originally to fulfil a commercial, decision-support role in the management of industrial water quality, where plankton growth is an issue. The progressive refinements to the model nevertheless have a robust ecological basis. This makes PROTECH a promising tool for researching plankton community ecology. The model calculates exponents describing growth and attrition, from a base of the maximum growth rates of algal species in culture. Subject to defined thresholds, growth integrates variability in the fluxes of light and nutrients. The paper develops this philosophy and its embodiment into the structure of the model. Examples of its authenticated, validated and sensitivity-tested outputs are presented.     

A Comparison of Single and Multiple Habitat Protocols for Collecting Macroinvertebrates in Wadeable Streams1

Abstract: In 2003, we compared two benthic macroinvertebrate sampling methods that are used for rapid biological assessment of wadeable streams. A single habitat method using kick sampling in riffles and runs was compared to a multiple habitat method that sampled all available habitats in proportion of occurrence. Both methods were performed side‐by‐side at 41 sites in lower gradient streams of the Piedmont and Northern Piedmont ecoregions of the United States, where riffle habitat is less abundant. Differences in sampling methods were examined using similarity indices, two multimetric indices [the family‐level Virginia Stream Condition Index (VSCI) and the species‐level Macroinvertebrate Biotic Integrity Index (MBII)], their component metrics, and bioassessment endpoints based on each index. Index scores were highly correlated between single and multiple habitat field methods, and sampling method comparability, based on comparison of similarities between and within sampling methods, was particularly high for species level data. The VSCI scores and values of most of its component metrics were not significantly higher for one particular method, but relationships between single and multiple habitat values were highly variable for percent Ephemeroptera, percent chironomids, and percent Plecoptera and Trichoptera (Hydropsychidae excluded). A similar level of variability in the relationship was observed for the MBII and most of its metrics, but Ephemeroptera richness, percent individuals in the dominant five taxa, and Hilsenhoff Biotic Index scores all exhibited differences in values between single and multiple habitat field methods. When applied to multiple habitat samples, the MBII exhibited greater precision, higher index scores, and higher assessment categories than when applied to single habitat samples at the same sites. In streams with limited or no riffle habitats, the multiple habitat method should provide an adequate sample for biological assessment, and at sites with abundant riffle habitat, little difference would be expected between the single and multiple habitat field methods. Thus, in geographic areas with a wide variety of stream types, the multiple habitat method may be more desirable. Even so, the variability in the relationship between single and multiple habitat methods indicates that the data are not interchangeable, and we suggest that any change in sampling method should be accompanied by a recalibration of any existing assessment tool (e.g., multimetric index) with data collected using the new method, regardless of taxonomic level.     

Multiphase flow and transport in fractured clay/sand sequences.

A numerical model (Queen's University Multi-Phase Flow Simulator, QUMPFS) was used to assess the rate of trichloroethylene (TCE) dense, non-aqueous phase liquid (DNAPL) migration through fractured clay, with special attention focused on the influence of interbedded sand lenses. The presence of these sand lenses was found to increase the time required for the non-wetting phase to migrate through the full 30 m vertical extent of the clay sequence from a few days to several years. Applied vertical hydraulic gradients were found to be moderately influential in systems consisting solely of fractured clays, yet one of the dominant factors controlling speed of vertical migration when sand lenses were present. Larger displacement pressure of the sands relative to that of the fractures leads to slower DNAPL migration rates, due to the delays that occur during build-up of capillary pressures. Dissolution of DNAPL and subsequent matrix diffusion of the aqueous phase has little effect on the rate of DNAPL migration through systems consisting of fractured clay only, yet slows the rate of migration in systems containing sand lenses. In all cases examined, the rate of DNAPL loading to the lower aquifer far exceeded the rate of aqueous phase mass loading. It was also found that DNAPL reaches the lower aquifer at approximately the same time as the aqueous phase plumes even for systems experiencing downward groundwater flow due to the attenuation of the aqueous phase through matrix diffusion.     

Using Multiple Chemical Criteria for Critical Loads of Acidity in Maritime Regions

Critical loads have become a well-establishedpart of the work programme of the UnitedNations Economic Commission for Europe (UNECE) Conventionon Long-Range Transboundary Air Pollution (CLRTAP). Thelinking of ecosystem response to deposition level is thecentral principal of the critical loads approach. Foreach ecosystem, a biological indicator is chosen, asuitable chemical criterion selected and a criticalchemical limit assigned. The Bc:Al ratio is the mostwidely used chemical criterion for setting criticalloads. However, critical loads based on this criterionare very sensitive to marine deposition. In regions whichreceive high depositions of marine-derived base cations,such as the Republic of Ireland (ROI) and the UnitedKingdom (UK), critical loads based on the Bc:Al criterionwill inevitably be high. Therefore, it is proposed thatcritical loads are estimated using multiple chemicalcriteria with appropriate critical limits to protect thechosen biological indicators. The range of publishedchemical criteria have been applied to the ROI and theUK. The chemical criterion corresponding to the mostsensitive critical load have been mapped and thecontribution of each to the final maps investigated. Thesimulations indicate that the most sensitive criteria forsetting critical loads are based on specifying critical Hor Al concentrations. However, the choice of critical limits andmodel parameters will ultimately effect the criticalloads. Therefore, it is important that appropriate criticallimits are chosen to protect the biological indicator andreceptor ecosystem from long-term damage.     

The Role of Mineral Weathering Rate Determinations in Generating Uncertainties in the Calculation of Critical Loads of Acidity and their Exceedance

Current legislation within Europe aimed at limitingecosystem damage resulting from inputs of atmosphericpollution is based on the critical load concept. Mineralweathering rates are central to the calculation ofcritical loads (acceptable levels) of acid deposition.The authors have undertaken a number of studies whichillustrate the complications and limitations inherent inpredicting mineral weathering rates and the implicationswhich these have for critical loads calculations andmapping. Calculated weathering rates and critical loadsfor two acid-sensitive parent materials (greywackes andgranites) are presented and are used to illustrate theimpact that uncertainty can have on critical loadexceedances. The results have obvious implications forportraying the uncertainties of critical loads to policy makers.     

The Link Between the Exceedance of Acidity Critical Loads for Freshwaters,Current Chemical Status and Biological Damage: A Re-Interpretation

Critical loads have for several years been employed bypolicymakers to aid in the development of strategies for aciddeposition abatement. They provide an effects-based approachwhereby an acid deposition flux greater than the critical load(known as critical load exceedance) implies that long-termharmful effects on a selected target organism will occur.Implicit in this approach are two assumptions: first, theexceedance of a critical load will harm the target organism,and second, the severity of biological impact is related to themagnitude of exceedance. However, static models give noindication of when the predicted damage might occur. One suchmodel, the Steady-State Water Chemistry (SSWC) model, employs aseries of empirical relationships to derive the pre-industrial,baseline leaching rate of base cations from measured waterchemistry using the so-called `F-factor'. The SSWC model setsthe critical load relative to pre-industrial base cationleaching (a permanent buffer of acid deposition) and a selectedacid neutralizing capacity (ANC) value which corresponds with aknown likelihood of damage to a biological target organism.Here we interpret the meaning of critical load exceedance as aprediction of steady-state ANC, and explore the relationshipbetween exceedance of the critical load and current chemistry. We demonstrate that a critical loadexceedance with the SSWC model does not necessarily indicatethat the critical chemical threshold (zero ANC) has alreadybeen crossed, and there may be no correlation betweenexceedance and biological status. A reformulation of the SSWCmodel is proposed which provides a direct link between currentdeposition and current chemical conditions, and is thereforemore likely to indicate current biological damage. Thereformulation illustrates the discrepancy between currentchemical status and that predicted by the SSWC model atsteady-state, which is a function of the `F-factor'.     

Influence of soil hydrological pathways on stream aluminium chemistry at Llyn Brianne, mid-wales

Storm runoff in afforested catchments at Llyn Brianne is acidic and Al-bearing. At baseflows, stream water is well-buffered with low Al levels. This paper presents the results of a study into how hydrological pathways account for these variations in stream-water chemistry. The investigation was carried out in the LI1 catchment; a 0.4-ha subcatchment covered by stagnohumic gley soils was monitored between October 1988 and September 1989. An instrumented hill-slope was established to identify the hydrological pathways that control the hydrochemistry of storm runoff draining from the subcatchment. Perched watertables developed in the surface horizons of the soil during storm episodes and produced lateral flow above the impeding subsoil. This near-surface flow path was responsible for generating acid, Al-rich storm runoff. Some water drained vertically through the soil profile into the underlying slope drift; seepage from groundwater in the drift sustained baseflows. Buffering reactions in the groundwater zone reduced the acidity and Al levels of baseflows. These hydrochemical characteristics are likely to be representative of other areas of stagnohumic gley soils, which cover 19% of the LI1 catchment: these soils may therefore provide a substantial source of acid, Al-bearing storm runoff in LI1 and similar afforested catchments.     

The scope and severity of white-nose syndrome on hibernating bats in North America

Assessing the scope and severity of threats is necessary for evaluating impacts on populations to inform conservation planning. Quantitative threat assessment often requires monitoring programs that provide reliable data over relevant spatial and temporal scales, yet such programs can be difficult to justify until there is an apparent stressor. Leveraging efforts of wildlife management agencies to record winter counts of hibernating bats, we collated data for 5 species from over 200 sites across 27 U.S. states and 2 Canadian provinces from 1995 to 2018 to determine the impact of white-nose syndrome (WNS), a deadly disease of hibernating bats. We estimated declines of winter counts of bat colonies at sites where the invasive fungus that causes WNS (Pseudogymnoascus destructans) had been detected to assess the threat impact of WNS. Three species undergoing species status assessment by the U.S. Fish and Wildlife Service (Myotis septentrionalis, Myotis lucifugus, and Perimyotis subflavus) declined by more than 90%, which warrants classifying the severity of the WNS threat as extreme based on criteria used by NatureServe. The scope of the WNS threat as defined by NatureServe criteria was large (36% of Myotis lucifugus range) to pervasive (79% of Myotis septentrionalis range) for these species. Declines for 2 other species (Myotis sodalis and Eptesicus fuscus) were less severe but still qualified as moderate to serious based on NatureServe criteria. Data-sharing across jurisdictions provided a comprehensive evaluation of scope and severity of the threat of WNS and indicated regional differences that can inform response efforts at international, national, and state or provincial jurisdictions. We assessed the threat impact of an emerging infectious disease by uniting monitoring efforts across jurisdictional boundaries and demonstrated the importance of coordinated monitoring programs, such as the North American Bat Monitoring Program (NABat), for data-driven conservation assessments and planning.