Kalman filter-based algorithms for monitoring the ionosphere and plasmasphere with GPS in near-real time

Data collected from a GPS receiver located at low latitudes in the American sector are used to investigate the performance of the WinTEC algorithm [Anghel et al., 2008a, Kalman filter-based algorithm for near realtime monitoring of the ionosphere using dual frequency GPS data. GPS Solutions, accepted for publication; for different ionospheric modeling techniques: the single-shell linear, quadratic, and cubic approaches, and the multi-shell linear approach. Our results indicate that the quadratic and cubic approaches perform much better than the single-shell and multi-shell linear approaches in terms of post-fit residuals. The performance of the algorithm for the cubic approach is then further tested by comparing the vertical TEC predicted by WinTEC and USTEC [Spencer et al., 2004. Ionospheric data assimilation methods for geodetic applications. In: Proceedings of IEEE PLANS, Monterey, CA, 26–29 April, pp. 510–517] at five North American stations. In addition, since the GPS-derived total electron content (TEC) contains contributions from both ionospheric and plasmaspheric sections of the GPS ray paths, in an effort to improve the accuracy of the TEC retrievals, a new data assimilation module that uses background information from an empirical plasmaspheric model [Gallagher et al., 1988. An empirical model of the Earth's plasmasphere. Advances in Space Research 8, (8)15–(8)24] has been incorporated into the WinTEC algorithm. The new Kalman filter-based algorithm estimates both the ionospheric and plasmaspheric electron contents, the combined satellite and receiver biases, and the estimation error covariance matrix, in a single-site or network solution. To evaluate the effect of the plasmaspheric component on the estimated biases and total TEC and to assess the performance of the newly developed algorithm, we compare the WinTEC results, with and without the plasmaspheric term included, at three GPS receivers located at different latitudes in the American sector, during a solar minimum period characterized by quiet and moderate geomagnetic conditions. We also investigate the consistency of our plasmaspheric results by taking advantage of the specific donut-shaped geometry of the plasmasphere and applying the technique at 12 stations distributed roughly over four geomagnetic latitudes and three longitude sectors.     

Arctic microbial ecosystems and impacts of extreme warming during the International Polar Year

As a contribution to the International Polar Year program MERGE (Microbiological and Ecological Responses to Global Environmental change in polar regions), studies were conducted on the terrestrial and aquatic microbial ecosystems of northern Canada (details at: http://www.cen.ulaval.ca/merge/). The habitats included permafrost soils, saline coldwater springs, supraglacial lakes on ice shelves, epishelf lakes in fjords, deep meromictic lakes, and shallow lakes, ponds and streams. Microbiological samples from each habitat were analysed by HPLC pigment assays, light and fluorescence microscopy, and DNA sequencing. The results show a remarkably diverse microflora of viruses, Archaea (including ammonium oxidisers and methanotrophs), Bacteria (including filamentous sulfur-oxidisers in a saline spring and benthic mats of Cyanobacteria in many waterbodies), and protists (including microbial eukaryotes in snowbanks and ciliates in ice-dammed lakes). In summer 2008, we recorded extreme warming at Ward Hunt Island and vicinity, the northern limit of the Canadian high Arctic, with air temperatures up to 20.5 °C. This was accompanied by pronounced changes in microbial habitats: deepening of the permafrost active layer; loss of perennial lake ice and sea ice; loss of ice-dammed freshwater lakes; and 23% loss of total ice shelf area, including complete break-up and loss of the Markham Ice Shelf cryo-ecosystem. These observations underscore the vulnerability of Arctic microbial ecosystems to ongoing climate change.     

Evolution of overland flow connectivity in bare agricultural plots

Soil surface roughness not only delays overland flow generation but also strongly affects the spatial distribution and concentration of overland flow. Previous studies generally aimed at predicting the delay in overland flow generation by means of a single parameter characterizing soil roughness. However, little work has been done to find a link between soil roughness and overland flow dynamics. This is made difficult because soil roughness and hence overland flow characteristics evolve differently depending on whether diffuse or concentrated erosion dominates. The present study examined whether the concept of connectivity can be used to link roughness characteristics to overland flow dynamics. For this purpose, soil roughness of three 30‐m² tilled plots exposed to natural rainfall was monitored for two years. Soil micro‐topography was characterized by means of photogrammetry on a monthly basis. Soil roughness was characterized by the variogram, the surface stream network was characterized by network‐based indices and overland flow connectivity was characterized by Relative Surface Connection function (RSCf) functional connectivity indicator. Overland flow hydrographs were generated by means of a physically‐based overland flow model based on 1‐cm resolution digital elevation models. The development of eroded flow paths at the soil surface not only reduced the delay in overland flow generation but also resulted in a higher continuity of high flow velocity paths, an increase in erosive energy and a higher rate of increase of the overland flow hydrograph. Overland flow dynamics were found to be highly correlated to the RSCf characteristic points. By providing information regarding overland flow dynamics, the RSCf may thus serve as a quantitative link between soil roughness and overland flow generation in order to improve the overland flow hydrograph prediction. Copyright © 2016 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.     

Effects of streamflow reductions on aquatic macroinvertebrates: linking groundwater withdrawals and assemblage response in southern New Jersey streams,USA

Abstract

A comprehensive hydro-ecological investigation was conducted to determine the ecological response of increased groundwater withdrawals from the Kirkwood-Cohansey aquifer system, an important source of water supply in southern New Jersey, USA. Collocated observations were made of aquatic-macroinvertebrate assemblages and stream hydrologic attributes to develop flow–ecology response relations. A sub-regional transient groundwater flow model (MODFLOW) was used to simulate three plausible high-stress groundwater-withdrawal scenarios which resulted in stream baseflow reductions of approximately 0.12, 0.20, and 0.26 m³ s^-1. These reduction scenarios were used to construct flow-alteration ecological response models to evaluate aquatic-macroinvertebrate response to streamflow reduction. For example, flow-alteration ecological response models indicate that if groundwater withdrawals diminish mean annual streamflow from 1.1 to 0.6 m³ s^-1, the abundance of intolerant taxa could be reduced by as much as 20%. These flow-alteration ecological response modelling results could be used by resource professionals to evaluate alternative water management strategies to determine maximum basin withdrawal rates that meet ongoing human water demand while protecting biological integrity.

Editor D. Koutsoyiannis; Guest editor M. Acreman

Citation Kennen, J.G., Riskin, M.L., and Charles, E.G., 2014. Effects of streamflow reductions on aquatic macroinvertebrates: linking groundwater withdrawals and assemblage response in southern New Jersey streams, USA. Hydrological Sciences Journal, 59 (3–4), 545–561.     

Organic Contaminants in Portland Cements Used in Monitoring Well Construction

We report the results of two independent laboratory investigations to evaluate total and leachable concentrations of glycols, glycol ethers, phenol, and other compounds in representative Type I and Type I/II Portland cement products that United States Environmental Protection Agency (EPA), The American Society of Testing and Materials (ASTM) and others recommend as annular sealants in monitoring well completions. Water well drillers also use these cements in their well completions. The EPA has included some of these compounds for analysis in their National Hydraulic Fracturing Study to evaluate the effects of hydraulic fracturing on ground‐ and surface water resources. During any contaminant investigation, materials used in monitoring or water well drilling, completion, development, and sampling must be free of the chemicals being targeted by the regulatory agency. Three of five bulk cement products we tested contained part per million (ppm) (mg/kg) concentrations of diethylene glycol, ethylene glycol, tetraethylene glycol, and triethylene glycol; chemicals added as grinding aids during manufacture. Some cements also contained ppb (µg/kg) concentrations of benzoic acid, phenols, propylene glycol, and 2‐butoxyethanol. Leaching of cured cement samples in water produced ppm (mg/L) concentrations of glycols in the supernatant. These results show that cured cements in monitoring or water wells can contaminate groundwater samples with glycols and phenol. Our findings should help prevent future sample bias and false positives when testing for glycol compounds and phenol in groundwater samples from monitoring or water wells and highlight the need to test materials or products used in monitoring or water well drilling, completions, development, and sampling to avoid false positives when sampling and analyzing for less common analytical constituents.     

Statistical downscaling of historical monthly mean winds over a coastal region of complex terrain. II. Predicting wind components

A regression-based downscaling technique was applied to monthly mean surface wind observations from stations throughout western Canada as well as from buoys in the Northeast Pacific Ocean over the period 1979–2006. A predictor set was developed from principal component analysis of the three wind components at 500?hPa and mean sea-level pressure taken from the NCEP Reanalysis II. Building on the results of a companion paper, Curry et al. (Clim Dyn 2011, doi:10.1007/s00382-011-1173-3), the downscaling was applied to both wind speed and wind components, in an effort to evaluate the utility of each type of predictand. Cross-validated prediction skill varied strongly with season, with autumn and summer displaying the highest and lowest skill, respectively. In most cases wind components were predicted with better skill than wind speeds. The predictive ability of wind components was found to be strongly related to their orientation. Wind components with the best predictions were often oriented along topographically significant features such as constricted valleys, mountain ranges or ocean channels. This influence of directionality on predictive ability is most prominent during autumn and winter at inland sites with complex topography. Stations in regions with relatively flat terrain (where topographic steering is minimal) exhibit inter-station consistencies including region-wide seasonal shifts in the direction of the best predicted wind component. The conclusion that wind components can be skillfully predicted only over a limited range of directions at most stations limits the scope of statistically downscaled wind speed predictions. It seems likely that such limitations apply to other regions of complex terrain as well.     

Feldspar dissolution rates in the Topopah Spring Tuff, Yucca Mountain, Nevada

Two different field-based methods are used here to calculate feldspar dissolution rates in the Topopah Spring Tuff, the host rock for the proposed nuclear waste repository at Yucca Mountain, Nevada. The center of the tuff is a high silica rhyolite, consisting largely of alkali feldspar (60 wt%) and quartz polymorphs (35 wt%) that formed by devitrification of rhyolitic glass as the tuff cooled. First, the abundance of secondary aluminosilicates is used to estimate the cumulative amount of feldspar dissolution over the history of the tuff, and an ambient dissolution rate is calculated by using the estimated thermal history. Second, the feldspar dissolution rate is calculated by using measured Sr isotope compositions for the pore water and rock. Pore waters display systematic changes in Sr isotopic composition with depth that are caused by feldspar dissolution. The range in dissolution rates determined from secondary mineral abundances varies from 10⁻¹⁶ to 10⁻¹⁷ mol s⁻¹ kg tuff⁻¹ with the largest uncertainty being the effect of the early thermal history of the tuff. Dissolution rates based on pore water Sr isotopic data were calculated by treating percolation flux parametrically, and vary from 10⁻¹⁵ to 10⁻¹⁶ mol s⁻¹ kg tuff⁻¹ for percolation fluxes of 15 mm a⁻¹ and 1 mm a⁻¹, respectively. Reconciling the rates from the two methods requires that percolation fluxes at the sampled locations be a few mm a⁻¹ or less. The calculated feldspar dissolution rates are low relative to other measured field-based feldspar dissolution rates, possibly due to the age (12.8 Ma) of the unsaturated system at Yucca Mountain; because oxidizing and organic-poor conditions limit biological activity; and/or because elevated silica concentrations in the pore waters (50 mg L⁻¹) may inhibit feldspar dissolution.     

Early burial mud diapirism and its impact on stratigraphic architecture in the Carboniferous of the Shannon Basin,County Clare,Ireland

Delta fronts are often characterized by high rates of sediment supply that result in unstable slopes and a wide variety of soft‐sediment deformation, including the formation of overpressured and mobile muds that may flow plastically during early burial, potentially forming mud diapirs. The coastal cliffs of County Clare, western Ireland, expose Pennsylvanian (Namurian) delta‐front deposits of the Shannon Basin at large scale and in three dimensions. These deposits include decametre‐scale, internally chaotic mudstone masses that clearly impact the surrounding sedimentary strata. Evidence indicates that these were true mud (unlithified sediment) diapirs that pierced overlying strata. This study documents a well‐exposed ca 20 m tall mud diapir and its impact on the surrounding mouth‐bar deposits of the Tullig Cyclothem. A synsedimentary fault and associated rollover dome, evident from stratal thicknesses and the dip of the beds, define one edge of the diapir. These features are interpreted as recording the reactive rise of the mud diapir in response to extensional faulting along its margin. Above the diapir, heterolithic sandstones and siltstones contain evidence for the creation of localized accommodation, suggesting synsedimentary filling, tilting and erosion of a shallow sag basin accommodated by the progressive collapse of the diapir. Two other diapirs are investigated using three‐dimensional models built from ‘structure from motion’ drone imagery. Both diapirs are interpreted to have grown predominantly through passive rise (downbuilding). Stratal relationships for all three diapirs indicate that they were uncompacted and fluid‐rich mud beds that became mobilized through soft‐sediment deformation during early burial (i.e. <50 m, likely <10 m depth). Each diapir locally controlled the stratigraphic architecture in the shallow subsurface and potentially influenced local palaeocurrents on the delta. The mud diapirs studied herein are distinct from deeper ‘shale diapirs’ that have been inferred from seismic sections worldwide, now largely disputed.