Vegetation Dynamics in a Tidal Freshwater Wetland: A Long-Term Study at Differing Scales

Tidal freshwater wetlands are complex, species-rich ecosystems located at the interface between tidal estuaries and nontidal rivers. This study conducted on the Patuxent River estuary in Maryland was designed to assess vegetation dynamics over several decades to determine if there were directional changes in the dominant communities. Aerial photographs (1970, 1989, and 2007) documented broad-scale spatial changes in major plant communities. The coverage of areas dominated by Nuphar lutea and Phragmites australis expanded; mixed vegetation and scrub–shrub habitats were essentially unchanged; and Typha and Zizania aquatica communities fluctuated in coverage. Data collected between 1988 and 2010 from permanent plots and transects were used to examine fine-scale changes. Shifts in the importance of some species through time were observed, but there were no directional changes in community species composition. The lack of directional change as measured at a fine scale is characteristic of tidal freshwater wetlands in which variations in the abundance of individual species, especially annuals, are responsible for most short-term change in species composition. Changes in the composition of plant communities are interpreted as responses to variations in vertical accretion, stability of habitat types, invasive plant species, and herbivores. In the future, vegetation changes are likely to occur as a result of the intrusion of brackish water and increased flooding associated with global climate change and sea level rise. This long-term study establishes a baseline from which potential future changes to tidal freshwater wetlands can be better understood.     

Scaling radio astronomy signal correlation on heterogeneous supercomputers using variousdata distribution methodologies

Next generation radio telescopes will require orders of magnitude more computing power to provide a view of the universe with greater sensitivity. In the initial stages of the signal processing flow of a radio telescope, signal correlation is one of the largest challenges in terms of handling huge data throughput and intensive computations. We implemented a GPU cluster based software correlator with various data distribution models and give a systematic comparison based on testing results obtained using the Fornax supercomputer. By analyzing the scalability and throughput of each model, optimal approaches are identified across a wide range of problem sizes, covering the scale of next generation telescopes.     

Inversion of time-dependent nuclear well-logging data using neural networks

The purpose of this work was to investigate a new and fast inversion methodology for the prediction of subsurface formation properties such as porosity, salinity and oil saturation, using time‐dependent nuclear well logging data. Although the ultimate aim is to apply the technique to real‐field data, an initial investigation as described in this paper, was first required; this has been carried out using simulation results from the time‐dependent radiation transport problem within a borehole. Simulated neutron and γ‐ray fluxes at two sodium iodide (NaI) detectors, one near and one far from a pulsed neutron source emitting at ～14 MeV, were used for the investigation. A total of 67 energy groups from the BUGLE96 cross section library together with 567 property combinations were employed for the original flux response generation, achieved by solving numerically the time‐dependent Boltzmann radiation transport equation in its even parity form. Material property combinations (scenarios) and their correspondent teaching outputs (flux response at detectors) are used to train the Artificial Neural Networks (ANNs) and test data is used to assess the accuracy of the ANNs. The trained networks are then used to produce a surrogate model of the expensive, in terms of computational time and resources, forward model with which a simple inversion method is applied to calculate material properties from the time evolution of flux responses at the two detectors. The inversion technique uses a fast surrogate model comprising 8026 artificial neural networks, which consist of an input layer with three input units (neurons) for porosity, salinity and oil saturation; and two hidden layers and one output neuron representing the scalar photon or neutron flux prediction at the detector. This is the first time this technique has been applied to invert pulsed neutron logging tool information and the results produced are very promising. The next step in the procedure is to apply the methodology to real data.     

On Constraining Pilot Point Calibration with Regularization in PEST

Ground water model calibration has made great advances in recent years with practical tools such as PEST being instrumental for making the latest techniques available to practitioners. As models and calibration tools get more sophisticated, however, the power of these tools can be misapplied, resulting in poor parameter estimates and/or nonoptimally calibrated models that do not suit their intended purpose. Here, we focus on an increasingly common technique for calibrating highly parameterized numerical models—pilot point parameterization with Tikhonov regularization. Pilot points are a popular method for spatially parameterizing complex hydrogeologic systems; however, additional flexibility offered by pilot points can become problematic if not constrained by Tikhonov regularization. The objective of this work is to explain and illustrate the specific roles played by control variables in the PEST software for Tikhonov regularization applied to pilot points. A recent study encountered difficulties implementing this approach, but through examination of that analysis, insight into underlying sources of potential misapplication can be gained and some guidelines for overcoming them developed.     

Thermobaric structure of the Kokchetav ultrahigh-pressure–high-pressure massif deduced from a north–south transect in the Kulet and Saldat–Kol regions, northern Kazakhstan

Abstract To investigate the regional thermobaric structure of the diamondiferous Kokchetav ultrahigh‐pressure and high‐pressure (UHP–HP) massif and adjacent units, eclogite and other metabasites in the Kulet and Saldat–Kol regions, northern Kazakhstan, were examined. The UHP–HP massif is subdivided into four units, bounded by subhorizontal faults. Unit I is situated at the lowest level of the massif and consists of garnet–amphibolite and acidic gneiss with minor pelitic schist and orthogneiss. Unit II, which structurally overlies Unit I, is composed mainly of pelitic schist and gneiss, and whiteschist locally with abundant eclogite blocks. The primary minerals observed in Kulet and Saldat–Kol eclogites are omphacite, sodic augite, garnet, quartz, rutile and minor barroisite, hornblende, zoisite, clinozoisite and phengite. Rare kyanite occurs as inclusions in garnet. Coesite inclusions occur in garnet porphyroblasts in whiteschist from Kulet, which are closely associated with eclogite masses. Unit III consists of alternating orthogneiss and amphibolite with local eclogite masses. The structurally highest unit, Unit IV, is composed of quartzitic schist with minor pelitic, calcareous, and basic schist intercalations. Mineral assemblages and compositions, and occurrences of polymorphs of SiO₂ (quartz or coesite) in metabasites and associated rocks in the Kulet and Saldat–Kol regions indicate that the metamorphic grades correspond to epidote–amphibolite, through high‐pressure amphibolite and quartz–eclogite, to coesite–eclogite facies conditions. Based on estimations by several geothermobarometers, eclogite from Unit II yielded the highest peak pressure and temperature conditions in the UHP–HP massif, with metamorphic pressure and temperature decreasing towards the upper and lower structural units. The observed thermobaric structure is subhorizontal. The UHP–HP massif is overlain by a weakly metamorphosed unit to the north and is underlain by the low‐pressure Daulet Suite to the south; boundaries are subhorizontal faults. There is a distinct pressure gap across these boundaries. These suggest that the highest grade unit, Unit II, has been selectively extruded from the greatest depths within the UHP–HP unit during the exhumation process, and that all of the UHP–HP unit has been tectonically intruded and juxtaposed into the adjacent lower grade units at shallower depths of about 10 km.     

Profiles of Fecal Porphyrins in River Otters Following the Exxon Valdez Oil Spill

Median levels of Coproporphyrin III (Copro III) in fecal samples of river otters (Lontra canadensis) collected from an oiled area in Prince William Sound, Alaska, USA, during 1990 were significantly higher than in samples collected from the same oiled area during 1996 (p=0.011, one way analysis of variance), a nonoiled reference area in Prince William Sound during 1996 (p=0.002) and a reference area in southeast Alaska during 1998 (p=0.004). An overall test of significance that combined probabilities from the statistical analysis of this porphyrin study with those from other biomarker studies revealed a significant difference in physiological response of river otters between oiled and nonoiled areas of the Sound for 1990 (p < 0.01). We demonstrated that changes in levels of fecal porphyrins may serve as a biomarker that may contribute to a health assessment of wild river otters.     

Low- and intermediate-mass close binary evolution and the initial–final mass relation

Using Eggleton's stellar evolution code, we carry out 150 runs of Population I binary evolution calculations with the initial primary mass between 1 and 8 M_⊙, the initial mass ratio     between 1.1 and 4, and the onset of Roche lobe overflow (RLOF) at an early, middle or late Hertzsprung-gap stage. We assume that RLOF is conservative in the calculations, and find that the remnant mass of the primary may change by more than 40 per cent over the range of initial mass ratio or orbital period, for a given primary mass. This is contrary to the often-held belief that the remnant mass depends only on the progenitor mass if mass transfer begins in the Hertzsprung gap. We fit a formula, with an error less than 3.6 per cent, for the remnant (white dwarf) mass as a function of the initial mass M _1i of the primary, the initial mass ratio q _i and the radius of the primary at the onset of RLOF. We also find that a carbon–oxygen white dwarf with mass as low as 0.33 M_⊙ may be formed if the initial mass of the primary is around 2.5 M_⊙.