Two-dimensional modelling of contaminant transport through saturated porous media using the radial point interpolation method (RPIM)

A new numerical tool is presented which models the two-dimensional contaminant transport through saturated porous media using a meshfree method called the radial point interpolation method (RPIM) with polynomial reproduction. In RPIM, an approximate solution is constructed entirely in terms of a set of nodes and no characterisation of the interrelationship of the nodes is needed. An advection-dispersion equation with sorption is considered to illustrate the applicability of the RPIM. The Galerkin weak form of the governing equation is formulated using two-dimensional meshfree shape functions constructed using thin plate spline radial basis functions. A computer program is developed for the implementation of the RPIM procedure. Three numerical examples are presented and the results are compared with those obtained from the analytical solution and finite element method. The experimental results are also used to validate the approach. The proposed RPIM has generated results with no oscillations and they are insensitive to Peclet constraints.     

Field scale heterogeneity of redox conditions in till-upscaling to a catchment nitrate model

Point scale studies in different settings of glacial geology show a large local variation of redox conditions. There is a need to develop an upscaling methodology for catchment scale models. This paper describes a study of field-scale heterogeneity of redox-interfaces in a till aquitard within an area of 600?×?600 m. The results showed significant variation of the depths to the redox-interface and thickness of the aquitard. Nitrate was present above the redox-interface but reduced to non-detectable levels a few metres below the interface. An upscaling approach for an area of 92 km² is proposed. Two models are proposed to predict the depth to the redox-interface in the aquitard and the resulting nitrate recharge concentrations to an underlying aquifer. The first model assumes that the depth to the redox-interface reflects the hydraulic head in the aquitard, and the second model assumes that the depth of the redox-interface is randomly distributed according to a log-normal probability distribution function. The upscaling approach using the random redox model estimated recharge concentrations comparable to the observed concentration in the underlying aquifer. The presented upscaling approach is applicable in distributed catchment models where sub-grid variability cannot be represented by the large grids.     

New basin modelling results from the Polish part of the Central European Basin system: implications for the Late Cretaceous–Early Paleogene structural inversion

Combined subsidence and thermal 1D modelling was performed on six well-sections located in the north-western Mid-Polish Trough/Swell in the eastern part of the Central European Basin system. The modelling allowed constraining quantitatively both the Mesozoic subsidence and the magnitude of the Late Cretaceous–Paleocene inversion and erosion. The latter most probably reached 2,400 m in the Mid-Polish Swell area. The modelled Upper Cretaceous thickness did not exceed 500 m, and probably corresponded to 200–300 m in the swell area as compared with more than 2,000 m in the adjacent non-inverted part of the basin. Such Upper Cretaceous thickness pattern implies early onset of inversion processes, probably in the Late Turonian or Coniacian. Our modelling, coupled with previous results of stratigraphic and seismic studies, demonstrates that the relatively low sedimentation rates in the inverted part of the basin during the Late Cretaceous were the net result of several discrete pulses of non-deposition and/or erosion that were progressively more pronounced towards the trough axis. The last phase of inversion started in the Late Maastrichtian and was responsible for the total amount of erosion, which removed also the reduced Upper Cretaceous deposits. According to our modelling results, a Late Cretaceous heat-flow regime which is similar to the present-day conditions (about 50 mW/m²) was responsible for the observed organic maturity of the Permian-Mesozoic rocks. This conclusion does not affect the possibility of Late Carboniferous–Permian and Late Permian–Early Triassic thermal events.     

Composition and pre-metamorphic geodynamic setting of the ultrahigh-pressure metabasic rocks from Dabie Shan,E-China

As part of the Yangtze plate, segments of the Dabie Shan terrane of Central China underwent ultra-high pressure metamorphism during Triassic subduction. We studied the geochemistry of the abundant eclogites to evaluate the nature of the protoliths and their geodynamic setting. Although some previous geochemical work exists, the analyses and interpretation herein are based on a new subdivision of the ultra-high pressure sequence into basement and cover units (Changpu and Ganghe Unit), revealing new and important results. In addition, eclogites of the so-called HP Unit south of the UHP units were studied. Whereas the large ion lithophile elements indicate postmagmatic, metasomatic changes of some samples, the high-field strength elements and the rare earth elements display original magmatic trends. The geochemical characteristics of the eclogites of the ultra-high pressure areas display a strong dependence on their “structural” and geographic position. The eclogites of the basement and the Changpu Unit indicate melt intrusion and extrusion in a continental rift system, i.e. during extensional tectonics. In contrast, the Ganghe Unit is characterized by a pronounced chemical homogeneity. The composition of the eclogites indicates generation from a mantle source highly influenced by slab-derived fluids. Those of the HP Unit show similar characteristics. Magmatism of the Ganghe and HP Unit probably occurred in a continental arc setting. A similar age for both units, geographically and/or tectonically separated, is possible. The geodynamic interpretation based on the geochemistry of the four units points to a Neoproterozoic scenario in which the protoliths of the HP and the cover units could have been of similar age and deposited in one evolving geological system. A rift-related larger-scale basin might have formed, e.g. a continental back-arc basin behind a magmatic arc after or simultaneous to sedimentation and magmatism in the magmatic arc. Alternatively, magmatism occurred in independent geodynamic settings, distinct in time and space. The units were juxtaposed during exhumation, after subduction to varying depths.     

Environmental Factors Contributing to the Development and Demise of a Toxic Dinoflagellate (<Emphasis Type="Italic">Karlodinium veneficum</Emphasis>) Bloom in a Shallow,Eutrophic, Lagoonal Estuary

A dense bloom of the ichthyotoxic dinoflagellate Karlodinium veneficum was discovered in the Neuse River Estuary, North Carolina, on 19 October 2006 and was associated with four subsequent fish kills. Microscopic, photopigment, DNA, and toxicological techniques confirmed bloom identity and toxicity. High-resolution spatio-temporal data from ship-board and fixed automated sampling stations provided a unique opportunity to investigate the environmental conditions that initiated, maintained, and terminated the K. veneficum bloom. Bloom initiation and growth were favored by high nutrient availability and reduced dispersal during the period of declining riverine discharge after Tropical Storm Ernesto. K. veneficum out-competed other co-occurring dinoflagellates, perhaps because of the production of karlotoxins that are known to act as grazing deterrents and to facilitate mixotrophic feeding. Once the bloom was established, small-scale hydrodynamic processes, coupled with vertical migration, concentrated cells along a frontal convergence to high densities (>200,000 cells per milliliter). By 26 October 2006, wind mixing and possible nutrient stress disrupted the bloom. Release of cell-bound toxins during the bloom collapse likely accounted for the associated fish kill events where fish were reported as exhibiting typical symptoms of karlotoxin poisoning. The dynamics of this bloom underscore the tight control of harmful algal blooms by meteorological forcing, hydrology, and sediment nutrient input in this shallow lagoonal estuary.     

Below- and Aboveground Spartina alterniflora Production in a Louisiana Salt Marsh

The monthly variations of below- and aboveground biomass of Spartina alterniflora were documented for a south Louisiana salt marsh from March 2004 to March 2005, and in March 2006 and 2007. The annual production rate above- and belowground was 1821 and 11,676 g m^?2, respectively (Smalley method), and the annual production rate per biomass belowground was 10.7 g dry weight^?1, which are highs along the latitudinal distributions of the plant’s range. The average root + rhizome/shoot ratio (R&R/S) was 2.6:1, which is lower than the R&R/S ratios of 4 to 5.1 reported for Spartina sp. marshes in the northeastern US. The belowground biomass increased from July to September and fluctuated between October and November, after which it declined until February when the growing season began. The belowground biomass was dominated by rhizomes, which declined precipitously in spring and then rose to a seasonal high in the month before declining again as the late summer rise in inflorescence began. Over half of the root biomass in a 30-cm soil profile was in the upper 10 cm, and in the 10- to 20-cm profile for rhizomes. The maximum March biomass above- and belowground was four to five times that of the minimum biomass over the four sampling years. The net standing stock (NSS) of N and P in live biomass aboveground compared to that in the belowground biomass was about 1.7 times higher and equal, respectively, but the NSS of N and P for the live + dead biomass was about six times higher belowground. The average nitrogen/phosphorous molar ratios of 16:1 aboveground is in agreement with the often tested N limitation of biomass accumulation aboveground, whereas the 37:1 belowground ratio suggests that there is an influence of P on R&R foraging for P belowground. Some implications for management and restoration are, in part, that salt marshes should be evaluated and examined using information on the plant’s physiology and production both below- and aboveground.     

Zircon M257 ‐ a Homogeneous Natural Reference Material for the Ion Microprobe U‐Pb Analysis of Zircon

We introduce and propose zircon M257 as a future reference material for the determination of zircon U‐Pb ages by means of secondary ion mass spectrometry. This light brownish, flawless, cut gemstone specimen from Sri Lanka weighed 5.14 g (25.7 carats). Zircon M257 has TIMS‐determined, mean isotopic ratios (2s uncertainties) of 0.09100 ± 0.00003 for ²⁰⁶pb/²³⁸U and 0.7392 ± 0.0003 for ²⁰⁷pb/²³⁵U. Its ²⁰⁶pb/²³⁸U age is 561.3 ± 0.3 Ma (unweighted mean, uncertainty quoted at the 95% confidence level); the U‐Pb system is concordant within uncertainty of decay constants. Zircon M257 contains ～ 840 μg g^?1 U (Th/U ～ 0.27). The material exhibits remarkably low heterogeneity, with a virtual absence of any internal textures even in cathodoluminescence images. The uniform, moderate degree of radiation damage (estimated from the expansion of unit‐cell parameters, broadening of Raman spectral parameters and density) corresponds well, within the “Sri Lankan trends”, with actinide concentrations, U‐Pb age, and the calculated alpha fluence of 1.66 × 10¹⁸ g^?1. This, and a (U+Th)/He age of 419 ± 9 Ma (2s), enables us to exclude any unusual thermal history or heat treatment, which could potentially have affected the retention of radiogenic Pb. The oxygen isotope ratio of this zircon is 13.9%o VSMOW suggesting a metamorphic genesis in a marble or calc‐silicate skarn.     

Provenance of the Bosnian Flysch

Sandwiched between the Adriatic Carbonate Platform and the Dinaride Ophiolite Zone, the Bosnian Flysch forms a c. 3000 m thick, intensely folded stack of Upper Jurassic to Cretaceous mixed carbonate and siliciclastic sediments in the Dinarides. New petrographic, heavy mineral, zircon U/Pb and fission-track data as well as biostratigraphic evidence allow us to reconstruct the palaeogeology of the source areas of the Bosnian Flysch basin in late Mesozoic times. Middle Jurassic intraoceanic subduction of the Neotethys was shortly followed by exhumation of the overriding oceanic plate. Trench sedimentation was controlled by a dual sediment supply from the sub-ophiolitic high-grade metamorphic soles and from the distal continental margin of the Adriatic plate. Following obduction onto Adria, from the Jurassic–Cretaceous transition onwards a vast clastic wedge (Vranduk Formation) was developed in front of the leading edge, fed by continental basement units of Adria that experienced Early Cretaceous synsedimentary cooling, by the overlying ophiolitic thrust sheets and by redeposited elements of coeval Urgonian facies reefs grown on the thrust wedge complex. Following mid-Cretaceous deformation and thermal overprint of the Vranduk Formation, the depozone migrated further towards SW and received increasing amounts of redeposited carbonate detritus released from the Adriatic Carbonate Platform margin (Ugar Formation). Subordinate siliciclastic source components indicate changing source rocks on the upper plate, with ophiolites becoming subordinate. The zone of the continental basement previously affected by the Late Jurassic–Early Cretaceous thermal imprint has been removed; instead, the basement mostly supplied detritus with a wide range of pre-Jurassic cooling ages. However, a c. 80 Ma, largely synsedimentary cooling event is also recorded by the Ugar Formation, that contrasts the predominantly Early Cretaceous cooling of the Adriatic basement and suggests, at least locally, a fast exhumation.