The Contribution of Groundwater Discharge to Nutrient Exports from a Coastal Catchment: Post-Flood Seepage Increases Estuarine N/P Ratios

Four months of daily nutrient and radon (a natural groundwater tracer) observations at the outlet of a heavily drained coastal wetland illustrated how episodic floods and diffuse groundwater seepage influence the biogeochemistry of a sub-tropical estuary (Richmond River, New South Wales, Australia). Our observations downstream of the Tuckean Swamp (an acid sulphate soil floodplain) covered a dry stage, a flood triggered by a 213-mm rain event and a post-flood stage when surface water chemistry was dominated by groundwater discharge. Significant correlations were found between radon and ammonium and N/P ratios and between radon and dissolved organic nitrogen (DON) during the post-flood stage. While the flood lasted for 14?% of the time of the surface water time series, it accounted for 18?% of NH₄, 32?% of NO _x , 66?% of DON, 58?% of PO₄ and 55?% of dissolved organic phosphorus (DOP) catchment exports. Over the 4-month study period, groundwater fluxes of 35.0, 3.6, 36.3, 0.5 and 0.7?mmol?m^?2?day^?1 for NH₄, NO _x , DON, PO₄ and DOP, respectively, were estimated. The groundwater contribution to the total surface water catchment exports was nearly 100?% for ammonium, and <20?% for the other nutrients. Post-flood groundwater seepage shifted the system from a DON to a dissolved inorganic N-dominated system and doubled N/P ratios in surface waters. We hypothesise that the Richmond River Estuary N/P ratios may reflect a widespread trend of tidal rivers and estuaries becoming more groundwater-dominated and phosphorus-limited as coastal wetlands are drained for agriculture, grazing and development.     

Late Cretaceous ophiolite obduction and Paleocene India-Asia collision in the westernmost Himalaya

Abstract

Collision of the Kohistan island arc with Asia at ~100 Ma resulted in N-S compression within the Neo-Tethys at a spreading center north of the Indo-Pakistani craton. Subsequent India-Asia convergence converted the Neo-Tethyan spreading center into a short-lived subduction zone. The hanging wall of the subduction zone became the Waziristan, Khost and Jalalabad igneous complexes. During the Santonian- Campanian (late Cretaceous), thrusting of the NW IndoPakistani craton beneath Albian oceanic crust and a Cenomanian volcano-sedimentary complex, generated an ophiolite-radiolarite belt. Ophiolite obduction resulted in tectonic loading and flexural subsidence of the NW Indian margin and sub-CCD deposition of shelf-derived olistostromes and turbidites in the foredeep. Campanian-Maastriehtian calci- clastic and siliciclastic sediment gravity flows derived from both margins filled the foredeep as a huge allochthon of Triassic-Jurassic rise and slope strata was thrust ahead of the ophiolites onto the Indo-Pakistani craton. Shallow to intermediate marine strata covered the foredeep during the late Maastrichtian. As ophiolite obduction neared completion during the Maastrichtian, the majority of India-Asia convergence was accommodated along the southern margin of Asia. During the Paleocene, India was thrust beneath a second allochthon that included open marine middle Maastrichtian colored mélange which represents the Asian Makran-Indus-Tsangpo accretionary prism. Latérites that formed on the eroded ophiolites and structurally higher colored mélange during the Paleocene wei’e unconformably overlapped by upper Paleocene and Middle Eocene shallow marine limestone and shale that delineate distinct episodes of Paleocene collisional and Early Eocene post-collisional deformation.     

Discretizing the fracture-matrix interface to simulate solute transport

This article examines the required spatial discretization perpendicular to the fracture-matrix interface (FMI) for numerical simulation of solute transport in discretely fractured porous media. The discrete-fracture, finite-element model HydroGeoSphere ( Therrien et al. 2005 ) and a discrete-fracture implementation of MT3DMS ( Zheng 1990 ) were used to model solute transport in a single fracture, and the results were compared to the analytical solution of Tang et al. (1981) . To match analytical results on the relatively short timescales simulated in this study, very fine grid spacing perpendicular to the FMI of the scale of the fracture aperture is necessary if advection and/or dispersion in the fracture is high compared to diffusion in the matrix. The requirement of such extremely fine spatial discretization has not been previously reported in the literature. In cases of high matrix diffusion, matching the analytical results is achieved with larger grid spacing at the FMI. Cases where matrix diffusion is lower can employ a larger grid multiplier moving away from the FMI. The very fine spatial discretization identified in this study for cases of low matrix diffusion may limit the applicability of numerical discrete-fracture models in such cases.     

Elevation adjustments of paired natural levees during flooding of the Saskatchewan River

Natural levees control the exchange of water between an alluvial channel and its floodplain, but little is known about the spatial distribution and evolution of levee heights. The summer 2005 flood of the Saskatchewan River (Cumberland Marshes, east‐central Saskatchewan) inundated large areas of floodplain for up to seven weeks, forming prominent new deposits on natural levees along main‐stem channels. Measurements of flood‐deposit thickness and crest heights of 61 levee pairs show that the thickest deposits occur on the lower pre‐flood levee in 80% of the sites, though no clear relationship exists between deposit thickness and magnitude of height difference. Only 16% of the pairs displayed thicker deposits on the higher levee, half of which occurred at sites where relatively clear floodbasin waters re‐entered turbid channels during general flooding. Difference in crest elevation (ΔE) between paired levees is approximately log‐normally distributed, both before and after the flood, though with different mean values. Supplemental observations from tank experiments indicate that during near‐bankfull flows, temporally and spatially variable deposition and erosion occur on levees due to backwater effects associated with nearby channel bars and irregular rises of the channel bed forced by channel extension. During floods, preferential deposition in lows tends to even out crest heights. Copyright © 2009 John Wiley & Sons, Ltd.     

A Survey of Techniques for Predicting Earthquake Ground Motions for Engineering Purposes

Over the past four or five decades many advances have been made in earthquake ground-motion prediction and a variety of procedures have been proposed. Some of these procedures are based on explicit physical models of the earthquake source, travel-path and recording site while others lack a strong physical basis and seek only to replicate observations. In addition, there are a number of hybrid methods that seek to combine benefits of different approaches. The various techniques proposed have their adherents and some of them are extensively used to estimate ground motions for engineering design purposes and in seismic hazard research. These methods all have their own advantages and limitations that are not often discussed by their proponents. The purposes of this article are to: summarise existing methods and the most important references, provide a family tree showing the connections between different methods and, most importantly, to discuss the advantages and disadvantages of each method.

Magnetic perturbations to the acoustic modes of roAp stars

We study the effect of a large-scale surface magnetic field on the non-radial acoustic modes of roAp stars. Special attention is given to the use of a variational principle which is used for determining the shifts in the frequencies with relative ease, enabling us to avoid having to calculate the perturbed eigenfunctions. With knowledge of the frequency shifts we then estimate the eigenfunctions in a simpler, albeit approximate way. The results indicate frequency shifts of the order of few μHz, which depend on the order, degree and azimuthal order of the mode. The loss of energy through Alfvén waves is also estimated from the imaginary parts of the frequency shifts. The results indicate that the loss is particularly high near specific frequencies. This might indicate the presence of a selection effect, which could make some modes more likely to be excited than others. However, our results do not explain why the modes observed appear always to be aligned with the axis of the magnetic field. Finally, the estimated perturbed eigenfunctions contain strong components of spherical harmonics that differ from those of the original unperturbed modes.     

Mg2+ removal in the system Mg2+—amorphous SiO2—H2O by adsorption and Mg-hydroxysilicate precipitation

Two chemical processes can remove Mg²⁺ from suspensions containing amorphous silica (am-SiO₂) at low temperatures: adsorption and precipitation of a Mg-hydroxysilicate resembling sepiolite. Mg²⁺ removal from am-SiO₂ suspensions was investigated, and the relative role of the two removal processes evaluated, as a function of: pH, ionic strength, Mg²⁺ concentration, and temperature.The extent of Mg²⁺ adsorption onto am-SiO₂ decreases with increasing NaCl concentration due to displacement of Mg²⁺ by Na⁺. At NaCl concentrations of 0.05 M and above, adsorption occurs only at pH values above 8.5, where rapid dissolution of am-SiO₂ gives rise to high concentrations of dissolved silica, resulting in supersaturation with respect to sepiolite. Removal of Mg²⁺, at concentrations of 40 to 650 μM, from am-SiO₂ suspensions in 0.70 M NaCl at 25 °C occurs at pH 9.0 and above. Experiments show that under these conditions adsorption and Mg-hydroxysilicate precipitation remove Mg²⁺ at similar rates. For 0.05 M Mg²⁺, at 0.70 M ionic strength and 25 °C, measurable Mg²⁺ removal occurs down to ca. pH 7.5 but is primarily due to Mg-hydroxysilicate precipitation. For the same solution conditions at 5°C, Mg²⁺ removal occurs above pH 8.0 and is primarily due to adsorption.Assuming that increasing pressure does not greatly enhance adsorption, Mg²⁺ adsorption onto am-SiO₂ is an insignificant process in sea water. The surface charge of pristine am-SiO₂ in sea water is primarily controlled by interactions with Na⁺. The principal reaction between Mg²⁺ and am-SiO₂ in marine sediments is sepiolite precipitation.The age distribution of sepiolite in siliceous pelagic sediments is influenced by temperatures of bottom waters and by geothermal gradients.     

Trace element distribution in mineral separates of the Allende inclusions and their genetic implications

Concentrations of the REE, Sc, Co, Fe, Zn, Ir, Na and Cr were determined by instrumental neutron activation and mass spectrometric isotope dilution analysis for mineral separates of the coarseand fine-grained types (group I and II of Martin and Mason's classification) of the Allende inclusions.These data, combined with data on mineral/liquid partition coefficients, oxygen isotope distributions and diffusion calculations, suggest the following: (1) Minerals in the coarse-grained inclusions (group I) crystallized in a closed system with respect to refractory elements. On the other hand, differences in oxygen isotope distributions among minerals preclude a totally molten stage in the history of the inclusion. Group I inclusions were formed by rapid condensation (either to liquid or solid) in a supercooled solar nebula; extrasolar pyroxene and spinel dust were included but not melted in the condensing inclusions, thus preserving their extrasolar oxygen isotope composition. REE were distributed by diffusion during the subsequent heating at subsolidus temperatures; because oxygen diffuses much more slowly at these temperatures, the oxygen isotope anomalies were preserved. (2) The fine-grained (group II) inclusions were also formed by condensation from a super-cooled nebular gas; however, REE-rich clinopyroxene and spinel were formed early and REE-poor sodalite and nepheline were formed later and mechanically mixed with clinopyroxene and spinel to form the inclusions. The REE patterns of the bulk inclusions and the mineral separates are fractionated, indicating that REE abundances in the gaseous phase were already fractionated at the time of condensation of the minerals. (3) Pre-existing Mg isotope anomalies in the coarse-grained inclusions must have been erased during the heating stage thus resetting the ²⁶Al-²⁶Mg chronometer.