Influence of Boundary Condition Types on Unstable Density‐Dependent Flow

Boundary conditions are required to close the mathematical formulation of unstable density‐dependent flow systems. Proper implementation of boundary conditions, for both flow and transport equations, in numerical simulation are critical. In this paper, numerical simulations using the FEFLOW model are employed to study the influence of the different boundary conditions for unstable density‐dependent flow systems. A similar set up to the Elder problem is studied. It is well known that the numerical simulation results of the standard Elder problem are strongly dependent on spatial discretization. This work shows that for the cases where a solute mass flux boundary condition is employed instead of a specified concentration boundary condition at the solute source, the numerical simulation results do not vary between different convective solution modes (i.e., plume configurations) due to the spatial discretization. Also, the influence of various boundary condition types for nonsource boundaries was studied. It is shown that in addition to other factors such as spatial and temporal discretization, the forms of the solute transport equation such as divergent and convective forms as well as the type of boundary condition employed in the nonsource boundary conditions influence the convective solution mode in coarser meshes. On basis of the numerical experiments performed here, higher sensitivities regarding the numerical solution stability are observed for the Adams‐Bashford/Backward Trapezoidal time integration approach in comparison to the Euler‐Backward/Euler‐Forward time marching approach. The results of this study emphasize the significant consequences of boundary condition choice in the numerical modeling of unstable density‐dependent flow.     

Inhomogeneous reionization of the intergalactic medium regulated by radiative and stellar feedbacks

We study the inhomogeneous reionization in a critical density CDM universe resulting from stellar sources, including Population III objects. The spatial distribution of the sources is obtained from high-resolution numerical N -body simulations. We calculate the source properties, taking into account a self-consistent treatment of both radiative (i.e. ionizing and H₂-photodissociating photons) and stellar (i.e. SN explosions) feedbacks regulated by massive stars. This allows us to describe the topology of the ionized and dissociated regions at various cosmic epochs, and to derive the evolution of H, He and H₂ filling factors, soft UV background, cosmic star formation rate and the final fate of ionizing objects. The main results are: (i) galaxies reionize the intergalactic medium by z ≈10 (with some uncertainty related to the gas clumping factor), whereas H₂ is completely dissociated already by z ≈25; (ii) reionization is mostly caused by the relatively massive objects which collapse via H line cooling, while objects the formation of which relies on H₂ cooling alone are insufficient for this purpose; (iii) the diffuse soft UV background is the major source of radiative feedback effects for z ≤15; at higher z direct flux from neighbouring objects dominates; (iv) the match of the calculated cosmic star formation history with that observed at lower redshifts suggests that the conversion efficiency of baryons into stars is ≈1 per cent; (v) we find that a very large population of dark objects which failed to form stars is present by z ≈8. We discuss and compare our results with similar previous studies.     

Evidence for low‐temperature growth of fayalite and hedenbergite in MacAlpine Hills 88107, an ungrouped carbonaceous chondrite related to the CM‐CO clan

Abstract— The carbonaceous chondrite MacAlpine Hills (MAC) 88107 has bulk composition and mineralogy that are intermediate between those of CO and CM chondrites. This meteorite experienced minor alteration and a low degree of thermal metamorphism (petrologic type 3.1) and escaped post‐accretional brecciation. The alteration resulted in the formation of fayalite (Fa_90–100). Al‐free hedenbergite (～Fs₅₀Wo₅₀), phyllosilicates (saponite‐serpentine intergrowths), magnetite, and Ni‐bearing sulfides (pyrrhotite and pentlandite). Fayalite and hedenbergite typically occur as veins, which start at the opaque nodules in the chondrule peripheries, crosscut fine‐grained rims and either terminate at the boundaries with the neighboring fine‐grained rims or continue as layers between these rims. These observations suggest that fayalite and hedenbergite crystallized after accretion and compaction of the fine‐grained rims. Fayalite also overgrows isolated forsteritic (Fa_1–5) and fayalitic (Fa_20–40) olivine grains without any evidence for Fe‐Mg interdiffusion; it also replaces massive magnetite‐sulfide grains. The initial ⁵³Mn/⁵⁵Mn ratio of (1.58 ± 0.26) × 10^?6 in the MAC 88107 fayalite corresponds to an age difference between the formation of fayalite and refractory inclusions in Allende of either ～9 or 18 Ma, depending upon the value of the solar system initial abundance of ⁵³Mn used in age calculations. Formation of secondary fayalite and hedenbergite requires mobilization and transport of Ca, Si, and Fe either through a high‐temperature gaseous phase (Hua and Buseck, 1995) or low‐temperature aqueous solution (Krot et al., 1998a, b). The high‐temperature nebular model for the origin of fayalite (Hua and Buseck, 1995) fails to explain (a) formation of fayalite‐hedenbergite assemblages after accretion of fine‐grained rims that lack any evidence for high‐temperature processing; (b) extreme fractionation of refractory lithophile elements of similar volatility, Ca and Al, in hedenbergite; and (c) absence of Fe‐Mg interdiffusion along fayalite‐forsterite boundaries. We conclude that fayalite and hedenbergite in MAC 88107 formed during late‐stage, low‐temperature (approximately 150–200 °C) aqueous alteration. The data for MAC 88107 extend the evidence for an early onset of aqueous activity on chondrite parent bodies and reinforce the conclusion that liquid water played an important role in the chemical and mineralogical evolution of the first chondritic planetesimals.     

Fine Resolution Palaeoecology Confirms Anthropogenic Impact During the Late Holocene in the Lower Hawkesbury Valley,NSW

This study reconstructs environmental conditions at Mill Creek, within the lower Hawkesbury-Nepean Valley, by the use of fine resolution palaeoecological studies of sediments. Archaeological surveys and historical accounts are used to investigate ecosystem response to known human activities. Research found that during prehistoric times (between 820 BP and the 1790s), the study area was well vegetated with dry sclerophyll communities on the valley sides, and a highly productive wetland community occupied the moist valley floor. Valley sides were generally stable, and little fire was evident in the landscape. Sediment was mainly deposited on the valley floor from the Hawkesbury-Nepean River during flood. After European settlement in the 1790s, sedimentation rates and charcoal preserved within sediments increased significantly, probably due to agricultural and clearing practices, both locally and upstream. Consequently, vegetation communities were altered. In 1967, the study area became part of Dharug National Park. Forestry and agricultural activities were excluded from the catchment, and the fire regime was reduced in frequency. Present-day vegetation communities appear to be becoming more like those that prevailed during prehistoric times. However, sedimentation rates remain elevated due to continued disturbance to the greater Hawkesbury River catchment. Recent poor water quality, coupled with a succession of floods on the Hawkesbury River, has increased wetland productivity levels, and their spatial extent, at Mill Creek.     

Improved understanding of spring and stream water responses in headwaters of the Indian Lesser Himalaya using stable isotopes,conductivity and temperature as tracers

Stable isotope data are presented for precipitation, spring and stream water in a headwater catchments in the Indian Lesser Himalaya. Isotopic contents of phreatic groundwater followed the local meteoric water line and showed minimal alteration by evaporation, suggesting fast recharge. Mean isotopic values for springs and the stream were close to the weighted annual mean for precipitation, indicating recharge was in synchrony with seasonal rainfall distribution. Precipitation exhibited isotopic declines of ?0.6‰ and ?0.2‰ δ¹⁸O per 100 m rise in elevation in July and August (monsoon), respectively. The time lag of one month between rainfall and spring discharge, combined with the isotopic lapse rate indicated a recharge elevation of 70–165 m above the spring outflow point, implying the water originated within the catchment. Time series of electrical conductivity and temperature of spring, seepage and stream waters confirmed the rapid recharge and limited storage capacity of the shallow aquifers.     

Automated Time Series Measurement of Microbial Concentrations in Groundwater-Derived Water Supplies

Fecal contamination by human and animal pathogens, including viruses, bacteria, and protozoa, is a potential human health hazard, especially with regards to drinking water. Pathogen occurrence in groundwater varies considerably in space and time, which can be difficult to characterize as sampling typically requires hundreds of liters of water to be passed through a filter. Here we describe the design and deployment of an automated sampler suited for hydrogeologically and chemically dynamic groundwater systems. Our design focused on a compact form to facilitate transport and quick deployment to municipal and domestic water supplies. We deployed a sampler to characterize water quality from a household well tapping a shallow fractured dolomite aquifer in northeast Wisconsin. The sampler was deployed from January to April 2017, and monitored temperature, nitrate, chloride, specific conductance, and fluorescent dissolved organic matter on a minute time step; water was directed to sequential microbial filters during three recharge periods that ranged from 5 to 20 days. Results from the automated sampler demonstrate the dynamic nature of the household water quality, especially with regard to microbial targets, which were shown to vary 1 to 2 orders of magnitude during a single sampling event. We believe assessments of pathogen occurrence and concentration, and related assessments of drinking well vulnerability, would be improved by the time-integrated characterization provided by this sampler.     

Offset‐variable density improves acoustic full‐waveform inversion: a shallow marine case study

We have previously applied three‐dimensional acoustic, anisotropic, full‐waveform inversion to a shallow‐water, wide‐angle, ocean‐bottom‐cable dataset to obtain a high‐resolution velocity model. This velocity model produced an improved match between synthetic and field data, better flattening of common‐image gathers, a closer fit to well logs, and an improvement in the pre‐stack depth‐migrated image. Nevertheless, close examination reveals that there is a systematic mismatch between the observed and predicted data from this full‐waveform inversion model, with the predicted data being consistently delayed in time. We demonstrate that this mismatch cannot be produced by systematic errors in the starting model, by errors in the assumed source wavelet, by incomplete convergence, or by the use of an insufficiently fine finite‐difference mesh. Throughout these tests, the mismatch is remarkably robust with the significant exception that we do not see an analogous mismatch when inverting synthetic acoustic data. We suspect therefore that the mismatch arises because of inadequacies in the physics that are used during inversion. For ocean‐bottom‐cable data in shallow water at low frequency, apparent observed arrival times, in wide‐angle turning‐ray data, result from the characteristics of the detailed interference pattern between primary refractions, surface ghosts, and a large suite of wide‐angle multiple reflected and/or multiple refracted arrivals. In these circumstances, the dynamics of individual arrivals can strongly influence the apparent arrival times of the resultant compound waveforms. In acoustic full‐waveform inversion, we do not normally know the density of the seabed, and we do not properly account for finite shear velocity, finite attenuation, and fine‐scale anisotropy variation, all of which can influence the relative amplitudes of different interfering arrivals, which in their turn influence the apparent kinematics. Here, we demonstrate that the introduction of a non‐physical offset‐variable water density during acoustic full‐waveform inversion of this ocean‐bottom‐cable field dataset can compensate efficiently and heuristically for these inaccuracies. This approach improves the travel‐time match and consequently increases both the accuracy and resolution of the final velocity model that is obtained using purely acoustic full‐waveform inversion at minimal additional cost.     

Long-term monitoring using resident and caged mussels in Boston Harbor yield similar spatial and temporal trends in chemical contamination

Measurements of chemical contaminants in caged (transplanted) and resident mussel populations have become a routine tool for monitoring and assessing the status and trends of coastal water quality. However, few long-term data sets are available to assess the comparability and efficacy of these two monitoring approaches. Three long-term independent data sets exist for Boston Harbor: the National Mussel Watch program has analyzed resident blue mussels (Mytilus edulis) from the Boston Harbor/Massachusetts Bay region for over twenty years, the Massachusetts Water Resources Authority has annually deployed caged (transplanted) mussels (M. edulis) to assess bioaccumulation potential of sewage effluent discharged under its NPDES permit for over fourteen years, and the GulfWatch program has analyzed resident blue mussel populations for over twelve years. Together, these data provide consistent and comparable information on temporal and spatial changes in chemical contamination in Boston Harbor as steps were taken to reduce contaminant loading. The data also demonstrate the complementary nature of resident and caged (transplanted) mussels for assessing contaminant trends even when the basic approaches and sampling frequency differ. These fifteen-year data sets demonstrate contaminant concentrations in mussels from Boston Harbor are similar and with few exceptions have significantly decreased since the early 1990s. The observed trends also demonstrate broad scale improvements to the quality of Boston Harbor and expand understanding of the response of coastal systems to interventions that reduce the load of chemicals to the ocean.     

Episodic riverine influence on surface DIC in the coastal Gulf of Maine

Anomalously high precipitation and river discharge during the spring of 2005 caused considerable freshening and depletion of dissolved inorganic carbon (DIC) in surface waters along the coastal Gulf of Maine. Surface pCO₂ and total alkalinity (TA) were monitored by repeated underway sampling of a cross-shelf transect in the western Gulf of Maine (GOM) during 2004–05 to examine how riverine fluxes, mixing, and subsequent biological activity exert control on surface DIC in this region. Most of the variability in surface DIC concentration was attributable to mixing of low DIC river water with higher DIC, saline GOM waters, but net biological uptake of DIC was significant especially during the spring and summer seasons. The extent and persistence of the coastal freshwater intrusion exerted considerable influence on net carbon dynamics. Integrated over the 10-m surface layer of our study region (∼5 × 10⁴ km²), net biological DIC uptake was 0.48 × 10⁸ mol C during April–July of 2004 compared to 1.33 × 10⁸ mol C during April–July of 2005. We found the temporal signature and magnitude of DIC cycling to be different in adjacent plume-influenced and non-plume regions. Extreme events such as the freshwater anomaly observed in 2005 will affect mean estimates of coastal carbon fluxes, thus budgets based on short time series of observations may be skewed and should be viewed with caution.