Life Histories, Salinity Zones, and Sublethal Contributions of Contaminants to Pelagic Fish Declines Illustrated with a Case Study of San Francisco Estuary, California, USA

Human effects on estuaries are often associated with major decreases in abundance of aquatic species. However, remediation priorities are difficult to identify when declines result from multiple stressors with interacting sublethal effects. The San Francisco Estuary offers a useful case study of the potential role of contaminants in declines of organisms because the waters of its delta chronically violate legal water quality standards; however, direct effects of contaminants on fish species are rarely observed. Lack of direct lethality in the field has prevented consensus that contaminants may be one of the major drivers of coincident but unexplained declines of fishes with differing life histories and habitats (anadromous, brackish, and freshwater). Our review of available evidence indicates that examining the effects of contaminants and other stressors on specific life stages in different seasons and salinity zones of the estuary is critical to identifying how several interacting stressors could contribute to a general syndrome of declines. Moreover, warming water temperatures of the magnitude projected by climate models increase metabolic rates of ectotherms, and can hasten elimination of some contaminants. However, for other pollutants, concurrent increases in respiratory rate or food intake result in higher doses per unit time without changes in the contaminant concentrations in the water. Food limitation and energetic costs of osmoregulating under altered salinities further limit the amount of energy available to fish; this energy must be redirected from growth and reproduction toward pollutant avoidance, enzymatic detoxification, or elimination. Because all of these processes require energy, bioenergetics methods are promising for evaluating effects of sublethal contaminants in the presence of other stressors, and for informing remediation. Predictive models that evaluate the direct and indirect effects of contaminants will be possible when data become available on energetic costs of exposure to contaminants given simultaneous exposure to non-contaminant stressors.     

Tectonic relationships between the Proterozoic Gawler and Willyama orogenic domains,Australia

Stratotectonic and morphotectonic data from the two principal exposed domains (pre‐Adelaidean rocks) of the Gawler sub‐province are used to characterize the Proterozoic Olarian orogeny and to distinguish its effects from those of the later Phanerozoic Delamerian orogeny.

The principal metasedimentary sequences in the Gawler domain and in the Willama domain are inferred to have been deposited in a single broad zone of early Proterozoic shallow‐water sedimentation on older (presumed Archaean) continental crust. The sequence becomes more pelitic upwards and may be interpreted as a transgressive sequence with more distal facies to the east.

Three main phases of deformation are recognized, and each phase has similar characteristics and age in both domains. D ₁ 2nd D₂ can be dated between 1850 and 1650 Ma, while D₃ appears to be about 1650–1540 Ma.

In high grade rocks, D₁ gave rise to a layer‐parallel schistosity, while D ₂ is characterized by tight folds with a high‐grade axial‐plane schistosity. The whole sub‐province was characterized by high geothermal gradients so that medium‐ to high‐grade metamorphism affected the lower parts of the succession before and during the D₁ and D₂ deformation episodes. No distinct tectonic zones can be recognized but large‐scale stratigraphic inversions (i.e. nappe tectonics) during D ₁ have been recognized only in the east of the Willyama domain. The higher parts of the stratigraphic succession are generally less deformed and exhibit only low‐grade metamorphism.

D ₃ produced relatively open, upright macroscopic folds and was characteristically associated with retrogression, but was demonstrably of pre‐Adelaidean age. The Gawler domain exhibits D ₃ structures although it lies in the platform west of the Adelaide Geosyncline and was not affected by deformation during Adelaidean sedimentation or by the subsequent Delamerian orogeny. A network of retrograde shear zones is the principal expression of post‐Olarian deformation in the Willyama domain which forms part of the basement to the Adelaide Geosyncline.

The trends of D ₂ and D ₃ folding in the two domains are similar and it is shown therefore that no large‐scale rotations of one domain relative to the other has been produced by the Delamerian orogeny. Large‐scale translations on discrete faults or on broad zones of simple shear in the basement are not easily ruled out, but if they exist, are probably largely of pre‐Adelaidean age. However, a significant relationship between Olarian structures and variable Adelaidean fold trends has been deduced.

The Olarian orogeny may have occurred in close proximity to a continental margin to the east and may thus be related to subduction processes. It differs from linear gneissic belts in Phanerozoic orogenies since it occurs in a more stable stratotectonic environment and over a wider area.     

Sediment mobility and bed armoring in the St Clair River: insights from hydrodynamic modeling

The lake levels in Lake Michigan‐Huron have recently fallen to near historical lows, as has the elevation difference between Lake Michigan‐Huron compared to Lake Erie. This decline in lake levels has the potential to cause detrimental impacts on the lake ecosystems, together with social and economic impacts on communities in the entire Great Lakes region. Results from past work suggest that morphological changes in the St Clair River, which is the only natural outlet for Lake Michigan‐Huron, could be an appreciable factor in the recent trends of lake level decline. A key research question is whether bed erosion within the river has caused an increase in water conveyance, therefore, contributed to the falling lake level. In this paper, a numerical modeling approach with field data is used to investigate the possibility of sediment movement in the St Clair River and assess the likelihood of morphological change under the current flow regime. A two‐dimensional numerical model was used to study flow structure, bed shear stress, and sediment mobility/armoring over a range of flow discharges. Boundary conditions for the numerical model were provided by detailed field measurements that included high‐resolution bathymetry and three‐dimensional flow velocities. The results indicate that, without considering other effects, under the current range of flow conditions, the shear stresses produced by the river flow are too low to transport most of the coarse bed sediment within the reach and are too low to cause substantial bed erosion or bed scour. However, the detailed maps of the bed show mobile bedforms in the upper St Clair River that are indicative of sediment transport. Relatively high shear stresses near a constriction at the upstream end of the river and at channel bends could cause local scour and deposition. Ship‐induced propeller wake erosion also is a likely cause of sediment movement in the entire reach. Other factors that may promote sediment movement, such as ice cover and dredging in the lower river, require further investigation. Copyright © 2012 John Wiley & Sons, Ltd.     

Quantifying the Morphological Print of Bedload Transport

Bedload and river morphology interact in a strong feedback manner. Bedload conditions the development of river morphology along different space and time scales; however, by concentrating the flow in preferential paths, a given morphology controls bedload for a given discharge. As bedload is a non‐linear response of shear stress, local morphology is likely to have a strong impact on bedload prediction when the shear stress is averaged over the section, as is usually done. This was investigated by comparing bedload measured in different bed morphologies (step‐pool, plane bed, riffle‐pool, braiding, and sand beds), with bedload measured in narrow flumes in the absence of any bed form, used here as a reference. The initial methodology consisted of fitting a bedload equation to the flume data. Secondly, the morphological signature of each river was studied as the distance to this referent equation. It was concluded that each morphology affects bedload in a different way. For a given average grain shear stress, the larger the river, the larger the deviation from the flume transport. Narrow streams are those morphologies that behave more like flumes; this is particularly true with flat beds, whereas results deviate from flumes to a greater extent in step‐pools. The riffle‐pool's morphology impacts bedload at different levels depending on the degree of bar development, considered here through the ratio D₈₄/D₅₀ which is used as a proxy for the local bed patchiness and morphology. In braiding rivers morphological effects are important but difficult to assess because width is dependent on transport rate. Bed morphology was found to have negligible effects in sand bed rivers where the Shields stress is usually sufficiently high to minimize the non‐linearity effects when hydraulics is averaged over the section. Copyright © 2015 John Wiley & Sons, Ltd.     

Membrane Interface Probe Protocol for Contaminants in Low‐Permeability Zones

Accurate characterization of contaminant mass in zones of low hydraulic conductivity (low k) is essential for site management because this difficult‐to‐treat mass can be a long‐term secondary source. This study developed a protocol for the membrane interface probe (MIP) as a low‐cost, rapid data‐acquisition tool for qualitatively evaluating the location and relative distribution of mass in low‐k zones. MIP operating parameters were varied systematically at high and low concentration locations at a contaminated site to evaluate the impact of the parameters on data quality relative to a detailed adjacent profile of soil concentrations. Evaluation of the relative location of maximum concentrations and the shape of the MIP vs. soil profiles led to a standard operating procedure (SOP) for the MIP to delineate contamination in low‐k zones. This includes recommendations for: (1) preferred detector (ECD for low concentration zones, PID or ECD for higher concentration zones); (2) combining downlogged and uplogged data to reduce carryover; and (3) higher carrier gas flow rate in high concentration zones. Linear regression indicated scatter in all MIP‐to‐soil comparisons, including R² values using the SOP of 0.32 in the low concentration boring and 0.49 in the high concentration boring. In contrast, a control dataset with soil‐to‐soil correlations from borings 1‐m apart exhibited an R² of ≥0.88, highlighting the uncertainty in predicting soil concentrations using MIP data. This study demonstrates that the MIP provides lower‐precision contaminant distribution and heterogeneity data compared to more intensive high‐resolution characterization methods. This is consistent with its use as a complementary screening tool.     

A Biogeochemical View of Estuarine Eutrophication: Seasonal and Spatial Trends and Correlations in the Delaware Estuary

The Delaware River and Bay Estuary is one of the major urbanized estuaries of the world. The 100-km long tidal river portion of the estuary suffered from major summer hypoxia in the past due to municipal and industrial inputs in the urban region; the estuary has seen remarkable water quality improvements from recent municipal sewage treatment upgrades. However, the estuary still has extremely high nutrient loading, which appears to not have much adverse impact. Since the biogeochemistry of the estuary has been relatively similar for the past two decades, our multiple year research database is used in this review paper to address broad spatial and seasonal patterns of conditions in the tidal river and 120 km long saline bay. Dissolved oxygen concentrations show impact from allochthonous urban inputs and meteorological forcing as well as biological influences. Nutrient concentrations, although high, do not stimulate excessive algal biomass due to light and multiple nutrient element limitations. Since the bay does not have strong persistent summer stratification, there is little potential for bottom water hypoxia. Elevated chlorophyll concentrations do not exert much influence on light attenuation since resuspended bottom inorganic sediments dominate the turbidity. Dissolved inorganic carbon and dissolved and particulate organic carbon distributions show significant variability from watershed inputs and lesser impact from urban inputs and biological processes. Ratios of dissolved and particulate carbon, nitrogen, and phosphorus help to understand watershed and urban inputs as well as autochthonous biological influences. Owing to the relatively simple geometry of the system and localized anthropogenic inputs as well as a broad spatial and seasonal database, it is possible to develop these biogeochemical trends and correlations for the Delaware Estuary. We suggest that this biogeochemical perspective allows a revised evaluation of estuarine eutrophication that should have generic value for understanding other estuarine and coastal waters.     

The size distributions of asteroid families in the SDSS Moving Object Catalog 4

Asteroid families, traditionally defined as clusters of objects in orbital parameter space, often have distinctive optical colors. We show that the separation of family members from background interlopers can be improved with the aid of SDSS colors as a qualifier for family membership. Based on an ∼88,000 object subset of the Sloan Digital Sky Survey Moving Object Catalog 4 with available proper orbital elements, we define 37 statistically robust asteroid families with at least 100 members (12 families have over 1000 members) using a simple Gaussian distribution model in both orbital and color space. The interloper rejection rate based on colors is typically ∼10% for a given orbital family definition, with four families that can be reliably isolated only with the aid of colors. About 50% of all objects in this data set belong to families, and this fraction varies from about 35% for objects brighter than an H magnitude of 13 and rises to 60% for objects fainter than this. The fraction of C-type objects in families decreases with increasing H magnitude for H>13, while the fraction of S-type objects above this limit remains effectively constant. This suggests that S-type objects require a shorter timescale for equilibrating the background and family size distributions via collisional processing. The size distribution varies significantly among families, and is typically different from size distributions for background populations. The size distributions for 15 families display a well-defined change of slope and can be modeled as a “broken” double power-law. Such “broken” size distributions are twice as likely for S-type familes than for C-type families (73% vs. 36%), and are dominated by dynamically old families. The remaining families with size distributions that can be modeled as a single power law are dominated by young families (<1 Gyr). When size distribution requires a double power-law model, the two slopes are correlated and are steeper for S-type families. No such slope-color correlation is discernible for families whose size distribution follows a single power law. For several very populous families, we find that the size distribution varies with the distance from the core in orbital-color space, such that small objects are more prevalent in the family outskirts. This “size sorting” is consistent with predictions based on the Yarkovsky effect.     

The Use of Indices to Identify Changes in Climatic Extremes

Changes in the frequencies of extremes are investigated by a variety of methods using daily temperature data from the British Isles, and monthly 5° latitude × 5° longitude grid-box temperatures over the land and marine regions of the world. The 225 year long daily Central England Temperature record shows no significant increase in very warm days in recent years but there is a marked decrease in the frequency of very cold days. Thus the rise in temperature in the last two decades is principally associated with a reduction in very cold days. Temperatures on days with particular wind circulation or pressure pattern types over the British Isles show multidecadal variations. Analyses using monthly gridded temperature data around the world since 1951 indicate that the recent rise in global surface temperatures is accompanied both by reductions in the areas affected by extremely cool temperatures and by increases in the areas with extremely warm temperatures.     

A case study of the nocturnal boundary layer over a complex terrain

A case study of the structure of the nocturnal boundary layer (NBL) over complex terrain is presented. Observations were made during the third night of Project STABLE (Weber and Kurzeja, 1991), whose main goal was to study turbulence and diffusion over the complex terrain of the Savannah River Site (SRS) near Augusta, Georgia.The passage of a mesoscale phenomenon, defined as a turbulent meso-flow (TMF) with an explanation of the nomenclature used, and a composite structure of the lowest few hundred meters over complex terrain are presented. The spatial extent of the TMF was at least 30–50 km, but the forcing is not well understood. The TMF occurred without the presence of a synoptic-scale cold front, under clear conditions, and with no discernible discontinuity in a microbarograph pressure trace. The structure of the NBL over the complex terrain at SRS differed from the expected homogeneous terrain NBL. The vertical structure exhibited dual low level wind maxima, dual inversions, and a persistent elevated turbulent layer.The persistent elevated turbulent layer, with a spatial extent of at least 30 km, was observed for the entire night. The persistent adiabatic layer may have resulted from turbulence induced by shear instability.     

Trace element partitioning between ilmenite, armalcolite and anhydrous silicate melt: Implications for the formation of lunar high-Ti mare basalts

We performed a series of experiments at high pressures and temperatures to determine the partitioning of a wide range of trace elements between ilmenite (Ilm), armalcolite (Arm) and anhydrous lunar silicate melt, to constrain geochemical models of the formation of titanium-rich melts in the Moon. Experiments were performed in graphite-lined platinum capsules at pressures and temperatures ranging from 1.1 to 2.3 GPa and 1300-1400 °C using a synthetic Ti-enriched Apollo ‘black glass’ composition in the CaO-FeO-MgO-Al₂O₃-TiO₂-SiO₂ system. Ilmenite-melt and armalcolite-melt partition coefficients (D) show highly incompatible values for the rare earth elements (REE) with the light REE more incompatible compared to the heavy REE ( 0.0020 ± 0.0010 to 0.069 ± 0.010 for ilmenite; 0.0048 ± 0.0023 to 0.041 ± 0.008 for armalcolite). D values for the high field strength elements vary from highly incompatible for Th, U and to a lesser extent W (for ilmenite: 0.0013 ± 0.0008, 0.0035 ± 0.0015 and 0.039 ± 0.005, and for armalcolite 0.008 ± 0.003, 0.0048 ± 0.0022 and 0.062 ± 0.03), to mildly incompatible for Nb, Ta, Zr, and Hf (e.g. 0.28 ± 0.05 and : 0.76 ± 0.07). Both minerals fractionate the high field strength elements with D_Ta/D_Nb and D_Hf/D_Zr between 1.3 and 1.6 for ilmenite and 1.3 and 1.4 for armalcolite. Armalcolite is slightly more efficient at fractionating Hf from W during lunar magma ocean crystallisation, with D_Hf/D_W = 12-13 compared to 6.7-7.5 for ilmenite. The transition metals vary from mildly incompatible to compatible, with the highest compatibilities for Cr in ilmenite (D ∼ 7.5) and V in armalcolite (D ∼ 8.1). D values show no clear variation with pressure in the small range covered.Crystal lattice strain modelling of D values for di-, tri- and tetravalent trace elements shows that in ilmenite, divalent elements prefer to substitute for Fe while armalcolite data suggest REE replacing Mg. Tetravalent cations appear to preferentially substitute for Ti in both minerals, with the exception of Th and U that likely substitute for the larger Fe or Mg cations. Crystal lattice strain modelling is also used to identify and correct for very small (∼0.3 wt.%) melt contamination of trace element concentration determinations in crystals.Our results are used to model the Lu-Hf-Ti concentrations of lunar high-Ti mare basalts. The combination of their subchondritic Lu/Hf ratios and high TiO₂ contents requires preferential dissolution of ilmenite or armalcolite from late-stage, lunar magma ocean cumulates into low-Ti partial melts of deeper pyroxene-rich cumulates.