Eutrophication and salinization of urban and rural kettle lakes in Kalamazoo and Barry Counties,Michigan, USA

Salinization and eutrophication caused by runoff of road salt and nutrients was assessed in three kettle lakes, two (Woods and Asylum Lakes) located in urban Kalamazoo, MI, and one (Brewster Lake) in rural Hastings, MI. Profiles of dissolved O₂, conductivity, pH, and temperature were measured in situ, at half meter intervals. Water samples were collected at discrete depth intervals of 1 m and analyzed for Fe(II), Mn(II), ammonium, alkalinity, Cl^?, Na, Mg, K and Ca. Results of this study indicate that all three lakes are eutrophic with anoxic bottom waters. Conductivity was much greater, and Cl^? levels were more than 100 times greater, in the two urban lakes compared to the rural lake, demonstrating the significant impact of road salt deicers on urban lake water chemistry.     

Biogenic carbon flows through the planktonic food web of the Amundsen Gulf (Arctic Ocean): A synthesis of field measurements and inverse modeling analyses

Major pathways of biogenic carbon (C) flow are resolved for the planktonic food web of the flaw lead polynya system of the Amundsen Gulf (southeast Beaufort Sea, Arctic Ocean) in spring-summer 2008. This period was relevant to study the effect of climate change on Arctic marine ecosystems as it was characterized by unusually low ice cover and warm sea surface temperature. Our synthesis relied on a mass balance estimate of gross primary production (GPP) of 52.5 ± 12.5 g C m⁻² calculated using the drawdown of nitrate and dissolved inorganic C, and a seasonal f-ratio of 0.64. Based on chlorophyll a biomass, we estimated that GPP was dominated by phytoplankton (93.6%) over ice algae (6.4%) and by large cells (>5 μm, 67.6%) over small cells (<5 μm, 32.4%). Ancillary in situ data on bacterial production, zooplankton biomass and respiration, herbivory, bacterivory, vertical particle fluxes, pools of particulate and dissolved organic carbon (POC, DOC), net community production (NCP), as well as selected variables from the literature were used to evaluate the fate of size-fractionated GPP in the ecosystem. The structure and functioning of the planktonic food web was elucidated through inverse analysis using the mean GPP and the 95% confidence limits of every other field measurement as lower and upper constraints. The model computed a net primary production of 49.2 g C m⁻², which was directly channeled toward dominant calanoid copepods (i.e. Calanus hyperboreus 20%, Calanus glacialis 10%, and Metridia longa 10%), other mesozooplankton (12%), microzooplankton (14%), detrital POC (18%), and DOC (16%). Bacteria required 29.9 g C m⁻², a demand met entirely by the DOC derived from local biological activities. The ultimate C outflow comprised respiration fluxes (82% of the initial GPP), a small sedimentation (3%), and a modest residual C flow (15%) resulting from NCP, dilution and accumulation. The sinking C flux at the model limit depth (395 m) supplied 60% of the estimated benthic C demand (2.8 g C m⁻²), suggesting that the benthos relied partly on other C sources within the bottom boundary layer to fuel its activity. In summary, our results illustrate that the ongoing decline in Arctic sea ice promotes the growth of pelagic communities in the Amundsen Gulf, which benefited from a ∼80% increase in GPP in spring-summer 2008 when compared to 2004 – a year of average ice conditions and relatively low GPP. However, 53% of the secondary production was generated within the microbial food web, the net ecological efficiency of zooplankton populations was not particularly high (13.4%), and the quantity of biogenic C available for trophic export remained low (6.6 g C m⁻²). Hence it is unlikely that the increase in lower food web productivity, such as the one observed in our study, could support new harvestable fishery resources in the offshore Beaufort Sea domain.     

Seasonal variation and tidal influences on estuarine use by bottlenose dolphins (Tursiops aduncus)

In order to show that dolphins use estuary habitats differently depending on the season and tidal state, possibly in response to prey distribution, temperature, risk of stranding and accessibility, Indo-Pacific bottlenose dolphins (Tursiops aduncus) were observed year-round during a 3-year study in the Clarence River estuary (CR) and Richmond River estuary (RR) in northern New South Wales, Australia. Peak dolphin sightings occurred during the spring season and one or 2 h prior to high tide. The spatial distribution of the dolphins in each estuary was analysed using the distance in kilometres that the dolphins travelled upstream with seasons and tidal phase as determinants. A General Linear Model showed that in the CR the dolphin spatial distribution in the estuary was not determined by season (F = 0.434, df = 3, P = 0.729) but was by tidal phase (F = 9.943, df = 3, P < 0.001) and the interaction between season and tidal phase (F = 3.398, df = 9, P < 0.002). However, in the RR the spatial distribution of the dolphin use of the estuary was not determined by either season (F = 1.647, df = 3, P = 0.194) or tidal phase (F = 0.302, df = 3, P = 0.824). In the CR, the spatial distribution of the dolphins was largest on high and flood tides. This pattern of spatial distribution may occur because the CR is a relatively shallow estuary and this increased spatial distribution may reflect a lower stranding risk and an increase in accessibility of shallow areas during periods of higher tide. These areas could also provide access to their preferred prey items of sea mullet (Mugil cephalus) and sand whiting (Sillago ciliata).     

Understanding the effect of single-fracture heterogeneity from single-well injection-withdrawal (SWIW) tests

The single-well injection-withdrawal (SWIW) tracer test is a method used to estimate the tracer retardation properties of a fracture or fracture zone. The effects of single-fracture aperture heterogeneity on SWIW-test tracer breakthrough curves are examined by numerical modelling. The effects of the matrix diffusion and sorption are accounted for by using a particle tracking method through the addition of a time delay added to the advective transport time. For a given diffusion and sorption property (P _m) value and for a heterogeneous fracture, the peak concentration is larger compared to a homogeneous fracture. The cumulative breakthrough curve for a heterogeneous fracture is similar to that for a homogeneous fracture and a less sorptive/diffusive tracer. It is demonstrated that the fracture area that meets the flowing water, the specific flow-wetted surface (sFWS) of the fracture, can be determined by matching the observed breakthrough curve for a heterogeneous fracture to that for a homogeneous fracture with an equivalent property parameter. SWIW tests are also simulated with a regional pressure gradient present. The results point to the possibility of distinguishing the effect of the regional pressure gradient from that of diffusion through the use of multiple tracers with different P _m values.     

Environmental Determinants of Emergent Macrophyte Vegetation in Pacific Northwest Estuarine Tidal Wetlands

We investigated whether within wetland environmental conditions or surrounding land cover measured at multiple scales were more influential in structuring regional vegetation patterns in estuarine tidal wetlands in the Pacific Northwest, USA. Surrounding land cover was characterized at the 100, 250, and 1,000 m, and watershed buffer scales. Vegetation communities were characterized by high species richness, lack of monotypic zonation, and paucity of invasive species. The number of species per site ranged between 4 and 20 (mean?±?standard deviation?=?10.2?±?3.1). Sites supported a high richness (mean richness of native species 8.7?±?2.8) and abundance of native macrophytes (mean relative abundance 85 %?±?19 %). Vegetation assemblages were dominated by a mix of grasses, sedges, and herbs with Sarcocornia pacifica and Distichlis spicata being common at sites in the oceanic zone of the estuary and Carex lyngbyei and Agrostis stolonifera being common at the fresher sites throughout the study area. The vegetation community was most strongly correlated with salinity and land cover within close proximity to the study site and less so with land cover variables at the watershed scale. Total species richness and richness of native species were negatively correlated with the amount of wetland in the buffer at all scales, while abundance of invasive species was significantly correlated to within wetland factors, including salinity and dissolved phosphorus concentrations. Landscape factors related to anthropogenic disturbances were only important at the 100-m buffer scale, with anthropogenic disturbances further from the wetland not being influential in shaping the vegetation assemblage. Our research suggests that the traditional paradigms of tidal wetland vegetation structure and environmental determinants developed in east coast US tidal wetlands might not hold true for Pacific Northwest wetlands due to their unique chemical and physical factors, necessitating further detailed study of these systems.     

In situ Raman observation of the crystallization in NaNO3–Na2SO4–H2O solution droplets

Several double salts have been detected in building materials and most of these salts are incongruently soluble compounds. In contrast to single salts, however, no systematic investigations of the crystallization behavior and deleterious effects of incongruently soluble double salts exist. To assess the damage potential of these salts, a systematic investigation of their highly complex behavior is desirable. This paper deals with the crystallization behavior of various solids in the ternary mixed NaNO₃–Na₂SO₄ system including the formation of the double salt darapskite, Na₃NO₃SO₄·H₂O. The crystallization sequence during droplet evaporation experiments at room conditions was determined using Raman and polarization microscopy. The basic idea of this research is to use deviations of the crystallization sequence of a salt or a mixed salt solution from the equilibrium pathway as an indicator to detect the degree of supersaturation. The observed crystallization pathway includes the formation of the metastable phases Na₂SO₄(III), Na₂SO₄(V) and darapskite. The experimental observations are discussed on the basis of the NaNO₃–Na₂SO₄–H₂O phase diagram and the results provide evidence for crystal growth from highly supersaturated solutions in both systems. If the crystals growing under these conditions are confined, these supersaturations result in substantial crystallization pressures.     

On the statistical distribution of seismic velocities in Earth's deep mantle

The existence of uncoupled shear (S) and compression (P) wave velocity variations in Earth's mantle is a characteristic that might only be explained by the presence of significant chemical and/or phase heterogeneity, with important implications for the dynamics and evolution of Earth's interior. While making a one-to-one comparison between tomographic models for P and S velocity (V_P and V_S) variations for a particular geographic region is ill-posed, their global statistical distributions reveal several robust characteristics indicative of the nature of uncoupled V_P and V_S in the deep mantle. We find that all of the V_P and V_S model distributions at a given depth are Gaussian-like throughout the lowermost mantle. However, a distinct low velocity feature is present in V_S distributions below ≈ 2200 km depth that is not present or is relatively weak in V_P models. The presence of anomalously low V_S material cannot be explained as an artifact, nor can the absence of a similarly strong feature in P models be ascribed to under-resolution. We propose that this feature can be partly explained by laterally variable occurrences of post-perovskite (pPv) lenses in the D″ layer, however, the persistence of significantly slow V_S regions at heights up to ≈ 700 km or more above the core–mantle boundary is likely to be incompatible with a pPv origin and might only be explained by the presence of a laterally discontinuous layer of chemically distinct material and/or some other kind of phase heterogeneity. There also exist significant discrepancies between tomographic models with respect to the width of the distributions as well as differences between the modeled peak values. We propose a scheme for comparison between different seismic models in which the widths of the dominant features in their statistical distributions is exploited.     

Uniwavelength lidar sensitivity to spherical aerosol microphysical properties for the interpretation of Lagrangian stratospheric observations

The determination of stratospheric particle microphysical properties from multiwavelength lidar, including Rayleigh and/or Raman detection, has been widely investigated. However, most lidar systems are uniwavelength operating at 532 nm. Although the information content of such lidar data is too limited to allow the retrieval of the full size distribution, the coupling of two or more uniwavelength lidar measurements probing the same moving air parcel may provide some meaningful size information. Within the ORACLE-O3 IPY project, the coordination of several ground-based lidars and the CALIPSO (Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation) space-borne lidar is planned during measurement campaigns called MATCH-PSC (Polar Stratospheric Clouds). While probing the same moving air masses, the evolution of the measured backscatter coefficient (BC) should reflect the variation of particles microphysical properties. A sensitivity study of 532 nm lidar particle backscatter to variations of particles size distribution parameters is carried out. For simplicity, the particles are assumed to be spherical (liquid) particles and the size distribution is represented with a unimodal log-normal distribution. Each of the four microphysical parameters (i.e. log-normal size distribution parameters, refractive index) are analysed separately, while the three others are remained set to constant reference values. Overall, the BC behaviour is not affected by the initial values taken as references. The total concentration (N₀) is the parameter to which BC is least sensitive, whereas it is most sensitive to the refractive index (m). A 2% variation of m induces a 15% variation of the lidar BC, while the uncertainty on the BC retrieval can also reach 15%. This result underlines the importance of having both an accurate lidar inversion method and a good knowledge of the temperature for size distribution retrieval techniques. The standard deviation (σ) is the second parameter to which BC is most sensitive to. Yet, the impact of m and σ on BC variations is limited by the realistic range of their variations. The mean radius (r_m) of the size distribution is thus the key parameter for BC, as it can vary several-fold. BC is most sensitive to the presence of large particles. The sensitivity of BC to r_m and σ variations increases when the initial size distributions are characterized by low r_m and large σ. This makes lidar more suitable to detect particles growing on background aerosols than on volcanic aerosols.