A composite spessartite–appinite intrusion from Port-na-Blagh,County Donegal,Ireland

A sheet-like appinite from Breaghy Head, North Donegal, is closely associated with a spessartite intrusion. Field, mineralogical, and geochemical similarities indicate that the parent magmas were very similar, and the appinite magma may have been derived from the spessartite by small degrees of crystal fractionation. The appinite is divided into an amphibole-rich lower part and a felsic upper part; the mineralogical contrast between the two is best explained by two-pulse intrusion of a batch of spessartitic magma within which amphibole and pyroxene had already segregated due to crystal settling or flow sorting. Unlike other Caledonian appinites, the Breaghy Head intrusion is not one of a cluster, and appears to be an isolated development. Possible reasons for its location are discussed.     

The appraisal of surface orbital vibrators with buried geophone array for permanent reservoir monitoring

Time-lapse seismic is one of the main methods for monitoring changes in reservoir conditions caused by production or injection of fluids. One approach to time-lapse seismic is through permanent reservoir monitoring, whereby seismic sources and/or receivers are permanently deployed. Permanent reservoir monitoring can offer a more cost-effective and environmentally friendly solution than traditional campaign-based surveys that rely on temporarily deployed equipment while facilitating more frequent measurements. At the CO2CRC Otway Project, surface orbital vibrators were coupled to a buried geophone array to form a permanent reservoir monitoring system. These are fixed position seismic sources that provide both P and S waves using induction motor-driven eccentric masses. After an initial injection of CO₂ in February 2016, five months of continuous seismic data were acquired, and reflection imaging was used to assess the system performance. Analysis of the data showed the effects of weather variations on the near-surface conditions and the sweep signatures of surface orbital vibrators. Data processing flows of the continuous data was adapted from Vibroseis four-dimensional data processing flows. Ground roll proved a significant challenge to data processing. In addition, variations in the surface wave pattern were linked to major rainfall events. For the appraisal of surface orbital vibrators in imaging, a Vibroseis four-dimensional monitor survey data with similar geometry was also processed. Surface orbital vibrators are observed to be reliable sources with a potential to provide a repeatable signal, especially if the ground roll should fall outside the target window of interest. To guide future permanent reservoir monitoring applications, a repeatability analysis was performed for the various key data processing steps.     

Stable isotope and trace element status of subsistence-hunted bowhead and beluga whales in Alaska and gray whales in Chukotka

Tissues of bowhead, beluga, and gray whales were analyzed for Ag, Cd, Cu, Se, Zn, THg and MeHg (belugas only). Delta15N and delta13C in muscle were used to estimate trophic position and feeding habitat, respectively. Trace element concentrations in tissues were significantly different among whale species. Hepatic Ag was higher in belugas than bowheads and gray whales. Gray whales had lower Cd concentrations in liver and kidney than bowhead and belugas and a sigmoid correlation of Cd with length was noted for all whales. Renal and hepatic Se and THg were higher in belugas than in baleen whales. The hepatic molar ratio of Se:THg exceeded 1:1 in all species and was negatively correlated to body length. Hepatic and renal Zn in subsistence-harvested gray whales was lower than concentrations for stranded whales. Se:THg molar ratios and tissue concentrations of Zn may show promise as potential indicators of immune status and animal health.     

Patterns of inorganic phosphate uptake in Cassiopea xamachana: a bioindicator species

Nutrient levels in the nearshore waters of the Florida Keys have increased over the past few decades concomitant with a decline in the health of Florida's reef system. Phosphorus is a particular concern in the Florida Keys as it may be the limiting nutrient in nearshore waters. We demonstrate that the upside-down jellyfish, Cassiopea xamachana, decreases its rate of phosphate uptake following exposure to elevated levels of dissolved inorganic phosphate. We also show that this subsequent suppression of uptake rates persists for some time following exposure to elevated phosphates. Using these attributes, we experimentally investigated the use of C. xamachana as a bioindicator for dissolved inorganic phosphates in seawater. Our results show that these animals reveal comparative differences in environmental phosphates despite traditional testing methods yielding no detectable phosphates. We propose that C. xamachana is a bioindicator useful for integrating relevant information about phosphate availability in low nutrient environments.     

Relative depletion of niobium in some arc magmas and the continental crust: partitioning of K, Nb, La and Ce during melt/rock reaction in the upper mantle

Depletion of Nb relative to K and La is characteristic of lavas in subduction-related magmatic arcs, as distinct from mid-ocean ridge basalts. Nb depletion is also characteristic of the continental crust. This and other geochemical similarities between the continental crust and high-Mg# andesite magmas found in arcs suggests that the continental crust may have formed by accretion of andesites. Previous studies have shown that the major element characteristics of high-Mg# andesites may be produced by melt/rock reaction in the upper mantle. In this paper, new data on partitioning of K, Nb, La and Ce between garnet, orthopyroxene and clinopyroxene in mantle xenoliths, and on partitioning of Nb and La between orthopyroxene and liquid, show that garnet and orthopyroxene have Nb crystal/liquid distribution coefficients which are much larger than those of K and La. Similar fractionations of Nb from K and La are expected in spinel and olivine. For this reason, reactions between migrating melt and large masses of mantle peridotite can produce substantial depletion of Nb in derivative liquids. Modeling shows that reaction between ascending, mantle-derived melts and mantle peridotite is a viable mechanism for producing the trace element characteristics of high-Mg# andesite magmas and the continental crust.

Alternatively, small-degree melts of metabasalt and/or metasediment in the subducting slab may leave rutile in their residue, and will thus have large Nb depletions relative to K and La [1]. Slab melts are too rich in light rare earth elements and other incompatible elements, and too poor in compatible elements, to be parental to arc magmas. However, ascending slab melts may be modified by reaction with the mantle. Our new data permit modeling of the trace element effects of reaction between small-degree melts of the slab and mantle peridotite. Modeling shows that this type of reaction is also a viable mechanism for producing the trace element characteristics of high-Mg# andesites and the continental crust. These findings, in combination with previous results, suggest that melt/rock reaction in the upper mantle has been an important process in forming the continental crust and mantle lithosphere.     

Defining probability of saturation with indicator kriging on hard and soft data

In humid, well-vegetated areas, such as in the northeastern US, runoff is most commonly generated from relatively small portions of the landscape becoming completely saturated, however, little is known about the spatial and temporal behavior of these saturated regions. Indicator kriging provides a way to use traditional water table data to quantify probability of saturation to evaluate predicted spatial distributions of runoff generation risk, especially for the new generation of water quality models incorporating saturation excess runoff theory. When spatial measurements of a variable are transformed to binary indicators (i.e., 1 if above a given threshold value and 0 if below) and the resulting indicator semivariogram is modeled, indicator kriging produces the probability of the measured variable to exceed the threshold value. Indicator kriging gives quantified probability of saturation or, consistent with saturation excess runoff theory, runoff generation risk with depth to water table as the variable and the threshold set near the soil surface. The probability of saturation for a 120 m × 180 m hillslope based upon 43 measurements of depth to water table is investigated with indicator semivariograms for six storm events. The indicator semivariograms show high spatial structure in saturated regions with large antecedent rainfall conditions. The temporal structure of the data is used to generate interpolated (soft) data to supplement measured (hard) data. This improved the spatial structure of the indicator semivariograms for lower antecedent rainfall conditions. Probability of saturation was evaluated through indicator kriging incorporating soft data showing, based on this preliminary study, highly connected regions of saturation as expected for the wet season (April through May) in the Catskill Mountain region of New York State. Supplementation of hard data with soft data incorporates physical hydrology of the hillslope to capture significant patterns not available when using hard data alone for indicator kriging. With the need for water quality models incorporating appropriate runoff generation risk estimates on the rise, this manner of data will lay the groundwork for future model evaluation and development.     

Stream geochemical response to reductions in acid deposition in headwater streams: Chronic versus episodic acidification recovery

In the decades since dramatic reductions in acid emissions and subsequent deposition in North America and Europe, the vast majority of research on aquatic chemical recovery has focused on trends in acid–base concentrations during baseflow conditions. Missing from such assessments is consideration of higher flow periods, when flow paths and chemical concentrations change and episodic acidification may occur. Stream samples collected weekly and bi‐hourly during high‐flow storm events from 1992 to 2015 at three watersheds representing a gradient in response to acid deposition in the south‐eastern United States were used to evaluate temporal trends in acid anions, base cations, acid neutralizing capacity (ANC), and pH for each of five distinct flow exceedance regimes, as well as trends in concentration–discharge (C–Q) relationships. For the most well‐buffered watershed (mean ANC = 220 μeq L^?1), ANC and pH increased at a similar rate for baseflow and storm flow conditions. Decreases in sulfate controlled ANC trends at higher flows, whereas smaller sulfate decreases in combination with base cation increases resulted in similar ANC improvements for low‐flow periods. For the most acid sensitive watershed (mean ANC = 10 μeq L^?1), no trends in ANC or pH were observed for the lowest flows (>90% flow exceedance), whereas reductions in sulfate resulted in significant increases in ANC during higher flow conditions. At all sites, greater rates of sulfate decline during high‐flow, as compared with low‐flow, conditions are likely a result of a reduced capacity of near surface soils to adsorb and retain sulfur in these non‐glaciated watersheds. Overall, consistent increases in pH (~0.01–0.02 pH units year^?1) during higher flow conditions (<10% flow exceedance) in contrast to the variable trends observed during lower flows (>50% flow exceedance) illustrate that episodic acidification is recovering at an equal or greater rate than chronic acidification in these watersheds.     

Fast data assimilation using a nonlinear Kalman filter and a model surrogate: An application to the Columbia River estuary

A fast and adjoint-free nonlinear data assimilation (DA) system was developed to simulate 3D baroclinic circulation in estuaries, leveraging two recently developed technologies: (1) a nonlinear model surrogate that executes forward simulation three orders of magnitude faster than a forward numerical circulation code and (2) a nonlinear extension to the reduced-dimension Kalman filter that estimates the state of the model surrogate. The noise sources in the Kalman filter were calibrated using empirical cross-validation and accounted for errors in model and model forcing.The DA system was applied to assimilate in situ measurements of water levels, salinities, and temperatures in simulations of the Columbia River estuary. To validate the DA results, we used a combination of cross-validation studies, process-oriented studies, and tests of statistical and dynamical consistency. The validation studies showed that DA improved the representation of several important processes in the estuary, including nonlinear tidal propagation, salinity intrusion, estuarine residual circulation, heat balance, and response of the estuary to coastal winds.     

Observations Of A Progressive Warm-Fresh Oceanic Front and A Cool-Moist Atmospheric Front At Duck, North Carolina

This paper presents a unique set of observations of nearlycoincident and progressive oceanic and marine atmospheric boundary-layer (MABL) fronts in a coastal zone. The event was observed during the afternoon of 12 May 1996 at the United States Army Corps of Engineers, Coastal Engineering Research Center, Field Research Facility pier at Duck, North Carolina. The oceanic front was warm and fresh. Current variabilityaccompanied the oceanic front. A marked MABL front preceded the oceanic front by several minutes and had characteristics of a sea-breezefront. This MABL front separated warmer and dryer pre-frontal air from cooler and moister post-frontal air. Wind direction and wind speedvariability accompanied the MABL front. The sea-surface roughness signatures of both fronts were detected by an X-band pulsed Doppler radar. Supporting data are used to identify each front detected by the radar and to calculate each front's velocity. In an attempt to explain the sea-surface roughness variations associated with each front, the radar data are compared to corresponding variations in wind speed, wind direction, and air-sea temperature difference.