Phantom midge-based models for inferring past fish abundances

We sampled living and subfossil phantom midge (Diptera: Chaoboridae) larvae from surface sediments of 21 small lakes in Southern Sweden to examine the influence of fish and selected abiotic variables on the abundance and species composition of chaoborid assemblages. We expected total Chaoborus abundance to be inversely correlated with fish abundance and Chaoborus species most sensitive to fish predation to be found only in fishless lakes. We aimed to use the observed relationships to develop models to reconstruct past fish abundances from chaoborid remains and the abiotic environment. C. flavicans occurred in almost every lake, whereas subfossil C. obscuripes were found in the surface sediments of only one fishless lake. The density of living C. flavicans larvae correlated negatively with fish abundance, lake order and size. The concentration of C. flavicans subfossils was negatively associated with pH, lake size, water transparency and fish abundance. Regression models that included lake morphometry and landscape position as additional predictors of fish abundance performed better than models that used only Chaoborus predictors. The explained variance in fish abundance varied from 52 to 86%. Leave-one-out cross-validation indicated moderate performance of the two best models. These models explained 51 and 56% of the observed untransformed fish density and biomass, respectively. In addition, all Chaoborus models were unbiased in closely following the 1:1 reference line in plots of observed versus predicted values. These results are a promising step in developing midge-based paleolimnological reconstructions of past fish abundance, and the approach might be improved by including chironomid remains in the models.     

Direct Trace Element Determination in Oil and Gas Produced Waters with Inductively Coupled Plasma‐Optical Emission Spectrometry: Advantages of High‐Salinity Tolerance

Waters co‐produced during petroleum extraction are the largest waste streams from oil and gas development. Reuse or disposal of these waters is difficult due to their high salinities and the sheer volumes generated. Produced waters (PWs) may also contain valuable mineral commodities. While an understanding of produced water trace element composition is required for evaluating the associated resource and waste potential of these materials, measuring trace elements in brines is challenging due to the dilution requirements of typical methods. Alternatively, inductively coupled plasma‐optical emission spectrometry (ICP‐OES) has shown promise as being capable of direct measurements of trace elements within PWs with minimal dilution. Here, we evaluate direct ICP‐OES trace element quantification in PWs for seventeen trace elements (As, Al, Ba, Be, Cd, Cr, Co, Cu, Hg, Mo, Ni, Pb, Rb, Sb, U, V and Zn) within fifteen PWs from five U.S. continuous reservoirs. The total analytical uncertainties associated with the trace element levels determined using ICP‐OES were estimated to be better than ± 30% (2s) except for Rb, which could not be determined due to ionisation interferences. The ICP‐OES results are compared with trace element levels determined using inductively coupled plasma‐mass spectrometry from the same samples. Our results demonstrate the potential for direct analysis of high‐salinity waters using ICP‐OES with minimal dilution and provide trace element concentrations in waters from several important U.S. petroleum‐generating reservoirs where available data are sparse.     

Evolution of Heterogeneous Mixed Siliciclastic/Carbonate Aquifers Containing Metastable Sediments

Mixed carbonate and siliciclastic marine sediments commonly become freshwater aquifers in eastern coastal regions of the United States and many other global locations. As these deposits age, the carbonate fraction of the sediment is commonly removed by dissolution and the aquifer can become a solely siliciclastic system or contain zones or beds of pure quartz sand. During aquifer evolution, the sediment grain size characteristics, hydraulic conductivity, and porosity change. An investigation of these changes using mixed carbonate/siliciclastic sediment samples collected from a modern barrier island beach in southern Florida showed that the average mean grain diameter decreased with removal of the carbonate fraction, but the average hydraulic conductivity and porosity increased slightly, but not to statistical significance. This counterintuitive result occurs because of the change in the pore types from a combined shelter and intergranular pore system producing a dual porosity system in the mixed sediments to a single intergranular pore system in the siliciclastic sediment fraction. Within the mixed carbonate/siliciclastic sediment, in the pure carbonate fraction, large shell fractions form grain‐supported large pores, which become filled with sand‐sized quartz as the shell fragments decrease in size or as the sediment becomes compacted. The hydraulic conductivity increases because the shell fragments that were oriented perpendicular to flow caused an increase in the length of the flow path, or a larger scale tortuosity, compared with the flow through pure quartz sand.     

Simulation of Meteorological Fields Within and Above Urban and Rural Canopies with a Mesoscale Model

Accurate simulation of air quality at neighbourhood scales (on order of 1-km horizontal grid spacing) requires detailed meteorological fields inside the roughness sub-layer (RSL). Since the assumptions of the roughness approach, used by most of the mesoscale models, are unsatisfactory at this scale, a detailed urban and rural canopy parameterisation, called DA-SM2-U, is developed inside the Penn State/NCAR Mesoscale Model (MM5) to simulate the meteorological fields within and above the urban and rural canopies. DA-SM2-U uses the drag-force approach to represent the dynamic and turbulent effects of the buildings and vegetation, and a modified version of the soil model SM2-U, called SM2-U(3D), to represent the thermodynamic effects of the canopy elements. The turbulence length scale is also modified inside the canopies. SM2-U(3D) assesses the sensible and latent heat fluxes from rural and urban surfaces in each of the computational layers inside the canopies by considering the shadowing effect, the radiative trapping by the street canyons, and the storage heat flux by the artificial surfaces. DA-SM2-U is tested during one simulated day above the city of Philadelphia, U.S.A. It is shown that DA-SM2-U is capable of simulating the important features observed in the urban and rural RSL, as seen in the vertical profiles of the shear stress, turbulent kinetic energy budget components, eddy diffusivity, potential air temperature, and specific humidity. Within the canopies, DA-SM2-U simulates the decrease of the wind speed inside the dense canopies, the skirting of the flow around the canopy blocks, warmer air inside the vegetation canopy than above open areas during the night and conversely during the day, and constantly warmer air inside the urban canopy. The comparison with measurements shows that the surface air temperature above rural and urban areas is better simulated by DA-SM2-U than by the `standard version' of MM5.     

Structure of the Ethiopian lithosphere: Xenolith evidence in the Main Ethiopian Rift

The lithospheric and sublithospheric processes associated with the transition from continental to oceanic magmatism during continental rifting are poorly understood, but may be investigated in the central Main Ethiopian Rift (MER) using Quaternary xenolith-bearing basalts. Explosive eruptions in the Debre Zeyit (Bishoftu) and Butajira regions, offset 20 km to the west of the contemporaneous main rift axis, host Al-augite, norite and lherzolite xenoliths, xenocrysts and megacrysts. Al-augite xenoliths and megacrysts derived from pressures up to 10 kb are the dominant inclusion in these recent basalts, which were generated as small degree partial melts of fertile peridotite between 15 and 25 kb. Neither the xenoliths nor the host basalts exhibit signs of carbonatitic or hydrous (amphibole + phlogopite) metasomatism, suggesting that infiltration of silicate melts resulting in pervasive Al-augite dyking/veining dominates the regional lithospheric mantle. Recent geophysical evidence has indicated that such veining/dyking is pervasive and segmented, supporting the connection of these Al-augite dykes/veins to the formation of a proto ridge axis. Al-augite xenoliths and megacrysts have been reported in other continental rift settings, suggesting that silicate melt metasomatism resulting in Al-augite dykes/veins is a fundamental processes attendant to continental rift development.     

Effects of hydrologic variability and remedial actions on first flush and metal loading from streams draining the Silverton caldera, 1992–2014

This study examined water quality in the upper Animas River watershed, a mined watershed that gained notoriety following the 2015 Gold King mine release of acid mine drainage to downstream communities. Water-quality data were used to evaluate trends in metal concentrations and loads over a two-decade period. Selected sites included three sites on tributary streams and one main-stem site on the Animas River downstream from the tributary confluences. During the study period, metal concentrations and loads varied seasonally and annually because of hydrologic variability and remedial actions designed to ameliorate the effects of acid mine drainage. Water-quality data were divided into two periods based on the timing of remedial activities in the watershed. The first period includes active water treatment, surface reclamation and installation of bulkheads in adits; the second period includes the decade following these activities. Water-quality data were used to estimate annual and monthly zinc loads using the Adjusted Maximum Likelihood Method (using LOADEST software) and U.S. Geological Survey streamflow data. This study presents one of the first applications of LOADEST focused on metal loads. Monthly flow-weighted concentrations were analysed using a Mann-Kendall trend test to determine the direction, magnitude, and significance of temporal trends in zinc loading in any given month and using t-test comparisons between the two periods. Zinc loads estimated for the Animas River below the tributaries indicate decreased zinc loading during the rising limb of the hydrograph in the second period, perhaps reflecting a reduction of snowmelt-derived zinc load following surface reclamation activities. In contrast, base-flow zinc loading increased at the main-stem site, perhaps because of the cessation of water treatment in tributary streams. Flow weighting of monthly load estimates yielded increased statistical significance and enabled more nuanced differentiation between the effects of hydrologic variability and remedial activities on zinc loading.     

Wide‐area estimates of evapotranspiration by red gum (Eucalyptus camaldulensis) and associated vegetation in the Murray–Darling River Basin,Australia

Floodplain red gum forests (Eucalyptus camaldulensis plus associated grasses, reeds and sedges) are sites of high biodiversity in otherwise arid regions of southeastern Australia. They depend on periodic floods from rivers, but dams and diversions have reduced flood frequencies and volumes, leading to deterioration of trees and associated biota. There is a need to determine their water requirements so environmental flows can be administered to maintain or restore the forests. Their water requirements include the frequency and extent of overbank flooding, which recharges the floodplain soils with water, as well as the actual amount of water consumed in evapotranspiration (ET). We estimated the flooding requirements and ET for a 38 134 ha area of red gum forest fed by the Murrumbidgee River in Yanga National Park, New South Wales. ET was estimated by three methods: sap flux sensors placed in individual trees; a remote sensing method based on the Enhanced Vegetation Index from MODIS satellite imagery and a water balance method based on differences between river flows into and out of the forest. The methods gave comparable estimates yet covered different spatial and temporal scales. We estimated flood frequency and volume requirements by comparing Normalized Difference Vegetation Index values from Landsat images with flood history from 1995 to 2014, which included both wet periods and dry periods. ET during wet years is about 50% of potential ET but is much less in dry years because of the trees' ability to control stomatal conductance. Based on our analyses plus other studies, red gum trees at this location require environmental flows of 2000 GL yr^?1 every other year, with peak flows of 20 000 ML d^?1, to produce flooding sufficient to keep them in good condition. However, only about 120–200 GL yr^?1 of river water is consumed in ET, with the remainder flowing out of the forest where it enters the Murray River system. Copyright © 2015 John Wiley & Sons, Ltd.     

Herschel SPIRE FTS telescope model correction

Emission from the Herschel telescope is the dominant source of radiation for the majority of SPIRE Fourier transform spectrometer (FTS) observations, despite the exceptionally low emissivity of the primary and secondary mirrors. Accurate modelling and removal of the telescope contribution is, therefore, an important and challenging aspect of FTS calibration and data reduction pipeline. A dust-contaminated telescope model with time invariant mirror emissivity was adopted before the Herschel launch. However, measured FTS spectra show a clear evolution of the telescope contribution over the mission and strong need for a correction to the standard telescope model in order to reduce residual background (of up to 7 Jy) in the final data products. Systematic changes in observations of dark sky, taken over the course of the mission, provide a measure of the evolution between observed telescope emission and the telescope model. These dark sky observations have been used to derive a time dependent correction to the telescope emissivity that reduces the systematic error in the continuum of the final FTS spectra to ～0.35 Jy.     

Herschel celestial calibration sources

Celestial standards play a major role in observational astrophysics. They are needed to characterise the performance of instruments and are paramount for photometric calibration. During the Herschel Calibration Asteroid Preparatory Programme approximately 50 asteroids have been established as far-IR/sub-mm/mm calibrators for Herschel. The selected asteroids fill the flux gap between the sub-mm/mm calibrators Mars, Uranus and Neptune, and the mid-IR bright calibration stars. All three Herschel instruments observed asteroids for various calibration purposes, including pointing tests, absolute flux calibration, relative spectral response function, observing mode validation, and cross-calibration aspects. Here we present newly established models for the four large and well characterized main-belt asteroids (1) Ceres, (2) Pallas, (4) Vesta, and (21) Lutetia which can be considered as new prime flux calibrators. The relevant object-specific properties (size, shape, spin-properties, albedo, thermal properties) are well established. The seasonal (distance to Sun, distance to observer, phase angle, aspect angle) and daily variations (rotation) are included in a new thermophysical model setup for these targets. The thermophysical model predictions agree within 5 % with the available (and independently calibrated) Herschel measurements. The four objects cover the flux regime from just below 1,000 Jy (Ceres at mid-IR N-/Q-band) down to fluxes below 0.1 Jy (Lutetia at the longest wavelengths). Based on the comparison with PACS, SPIRE and HIFI measurements and pre-Herschel experience, the validity of these new prime calibrators ranges from mid-infrared to about 700 μm, connecting nicely the absolute stellar reference system in the mid-IR with the planet-based calibration at sub-mm/mm wavelengths.     

Herschel SPIRE FTS relative spectral response calibration

Herschel/SPIRE Fourier transform spectrometer (FTS) observations contain emission from both the Herschel Telescope and the SPIRE Instrument itself, both of which are typically orders of magnitude greater than the emission from the astronomical source, and must be removed in order to recover the source spectrum. The effects of the Herschel Telescope and the SPIRE Instrument are removed during data reduction using relative spectral response calibration curves and emission models. We present the evolution of the methods used to derive the relative spectral response calibration curves for the SPIRE FTS. The relationship between the calibration curves and the ultimate sensitivity of calibrated SPIRE FTS data is discussed and the results from the derivation methods are compared. These comparisons show that the latest derivation methods result in calibration curves that impart a factor of between 2 and 100 less noise to the overall error budget, which results in calibrated spectra for individual observations whose noise is reduced by a factor of 2–3, with a gain in the overall spectral sensitivity of 23 % and 21 % for the two detector bands, respectively.