A parallel global-implicit 2-D solver for reactive transport problems in porous media based on a reduction scheme and its application to the MoMaS benchmark problem

In this article, an approach for the efficient numerical solution of multi-species reactive transport problems in porous media is described. The objective of this approach is to reformulate the given system of partial and ordinary differential equations (PDEs, ODEs) and algebraic equations (AEs), describing local equilibrium, in such a way that the couplings and nonlinearities are concentrated in a rather small number of equations, leading to the decoupling of some linear partial differential equations from the nonlinear system. Thus, the system is handled in the spirit of a global implicit approach (one step method) avoiding operator splitting techniques, solved by Newton’s method as the basic algorithmic ingredient. The reduction of the problem size helps to limit the large computational costs of numerical simulations of such problems. If the model contains equilibrium precipitation-dissolution reactions of minerals, then these are considered as complementarity conditions and rewritten as semismooth equations, and the whole nonlinear system is solved by the semismooth Newton method.     

The application of cluster analysis in geophysical data interpretation

A clustering algorithm that is based on density and is adaptive density-reachable is developed and presented for arbitrary data point distributions in some real-world applications, especially in geophysical data interpretation. Through comparisons of the new algorithm and other algorithms, it is shown that the new algorithm can reduce the dependency of domain knowledge and the sensitivity of abnormal data points, that it can improve the effectiveness of clustering results in which data are distributed in different shapes and different densities, and that it can get a better clustering efficiency. The application of the new clustering algorithm demonstrates that data mining techniques can be used in geophysical data interpretation and can get meaningful and useful results, and that the new clustering algorithm can be used in other real-world applications.     

Complete recycling of a magmatic arc: evidence from chemical and isotopic composition of Quaternary trench sediments in Chile (36°–40°S)

Marine Quaternary trench and slope sediments were sampled along the margin of the Southern Andes, Chile between 36° and 40°S. Major and trace element contents indicate only minor influence of weathering and transport fractionation. The whole rock composition of the sediments is similar to the average rock of the Cretaceous to Holocene magmatic arc of this section of the southern volcanic zone. Sr, Nd, and Pb isotope composition of the sediments also resembles closely the average composition of the magmatic arc. The contribution of compositionally distinct Palaeozoic crust, which makes up most of the volume of the forearc, is ~0–20% crustal Sr, Nd, and Pb according to the isotope record of the trench and slope sediments. Input of sediments from the continent into the subduction system was dominated by detritus from the magmatic arc at least for the last 20 My on the basis of the Oligocene to Holocene exhumation history of the margin.     

Integrated stratigraphy and 40Ar/39Ar chronology of the early to middle Miocene Upper Freshwater Molasse in western Bavaria (Germany)

A detailed integrated stratigraphic study (biostratigraphy and magnetostratigraphy) was carried out on five sections from the western part of the Bavarian Upper Freshwater Molasse of the North Alpine Foreland Basin (NAFB), greatly improving the chronostratigraphy of these sediments. The sections belong to the lithostratigraphic units Limnische Untere Serie (UL) and Fluviatile Untere Serie (UF) and contain 19 (mostly new) small-mammal bearing levels, significantly refining the local biostratigraphy. Radiometric ages obtained from glass shards from tuff horizons are used together with the biostratigraphic information for constructing and confirming the magnetostratigraphic correlation of the studied sections to the Astronomical Tuned Time Scale (ANTS04; Lourens et al. in Geologic Time Scale 2004, Cambridge University Press, 2004). This correlation implies that the UL lithostratigraphic unit corresponds to the latest Ottnangian and the Early Karpatian, whereas the UF corresponds to the Karpatian and the Early Badenian. This indicates that the Brackish- to Freshwater Molasse transition already occurred during the late Ottnangian. The pre-Riesian hiatus occurred in the latest Karpatian and lower Early Badenian in Eastern Bavaria and Bohemia and in the Late Karpatian and earliest Badenian in Western Bavaria. The geochemical and Ar–Ar data of volcanic ashes suggest that highly evolved silicic magmas from a single volcano or volcanic center, characterized by a uniform Nd isotopic composition, erupted repetitively over the course of at least 1.6 Myr. Three phases of eruptive activity were identified at 16.1 ± 0.2 Ma (Zahling-2), 15.6 ± 0.4 Ma (Krumbad), and 14.5 ± 0.2 Ma (Heilsberg, Hegau). The correlation of the local biostratigraphic zonation to the ANTS04 enables further the characterization of both the Ottnangian–Karpatian and Karpatian–Badenian boundaries in the NAFB by small-mammal biostratigraphy. According to these results the Ottnangian–Karpatian boundary is contemporaneous with the first appearance datum of Megacricetodon bavaricus (in the size of the type population) and the first common occurrence of Keramidomys thaleri, whereas Ligerimys florancei, Melissiodon dominans and Prodeinotherium aff. bavaricum have been already disappeared during the late Ottnangian. The Karpatian–Badenian boundary is characterized by a significant size increase of the large Megacricetodon lineage and possibly a (re-)immigration of Prodeinotherium bavaricum.     

Structural characteristics of the halite fabric type ‘Kristallbrocken’ from the Zechstein Basin with regard to its development

The Kristallbrocken are a characteristic centimetre- to decimetre-sized, laminated halite fabric type occurring in the Stassfurt Formation in the Zechstein Basin. Up to now, the nature of the Kristallbrocken, i.e. if they are relics of fine-grained, polycrystalline halite beds or clasts of ‘single crystal-layers’, as well as the deformation mechanisms of this halite type, were not clear from the literature. Drill core material from the salt deposit Teutschenthal at the southern rim of the Zechstein Basin now allowed investigating less intensely deformed samples for the first time. The deformational behaviour of these Kristallbrocken ranges from brittle to ductile, which is evidenced by fractured Kristallbrocken on the one hand and weakly bent or even folded Kristallbrocken on the other hand. Local X-ray texture measurements demonstrated that the Kristallbrocken are definitely single crystals and that they can be regarded as relics of formerly larger ‘single crystal-layers’ of up to several dm² in size. The folded Kristallbrocken clearly display by their single grain texture characteristics that their crystal lattice is bent, which was most likely enabled by a kind of flexural-shear folding, and did not develop after deformation from a fine-grained aggregate by recrystallisation. Due to their monocrystallinity, their originally large size, and the solid inclusions forming the internal lamination, the Kristallbrocken have clearly stronger rheological properties than the surrounding fine- to coarse-grained polycrystalline rock salt, and thus also deform by fracturing.     

Medium- and Long-term Recovery of Estuarine and Coastal Ecosystems: Patterns,Rates and Restoration Effectiveness

Many estuarine and coastal marine ecosystems have increasingly experienced degradation caused by multiple stressors. Anthropogenic pressures alter natural ecosystems and the ecosystems are not considered to have recovered unless secondary succession has returned the ecosystem to the pre-existing condition or state. However, depending upon the scales of time, space and intensity of anthropogenic disturbance, return along the historic trajectory of the ecosystem may: (1) follow natural restoration though secondary succession; (2) be re-directed through ecological restoration, or (3) be unattainable. In order to address the gaps in knowledge about restoration and recovery of estuarine and coastal ecosystems, this special feature includes the present overview and other contributions to provide a synthesis of our knowledge about recovery patterns, rates and restoration effectiveness. From the 51 examples collated in this contribution, we refine the recovery from the list of stressors into six recovery mechanisms: (1) recovery from sediment modification, which includes all aspects of dredging and disposal; (2) recovery by complete removal of stressors limiting natural ecosystem processes, which includes tidal marsh and inundation restoration; (3) recovery by speed of organic degradation, which includes oil discharge, fish farm wastes, sewage disposal, and paper mill waste; (4) recovery from persistent pollutants, which includes chemical discharges, such as TBT; (5) recovery from excessive biological removal, related to fisheries and (6) recovery from hydrological and morphological modification. Drawing upon experience both from these many examples and from an example of one comprehensive study, we show that although in some cases recovery can take <5 years, especially for the short-lived and high-turnover biological components, full recovery of coastal marine and estuarine ecosystems from over a century of degradation can take a minimum of 15–25 years for attainment of the original biotic composition and diversity may lag far beyond that period.     

Phytoplankton Interannual Variability at Cassino Beach, Southern Brazil (1992–2007), with Emphasis on the Surf Zone Diatom Asterionellopsis glacialis

A long-term study (monthly sampling, 1992 to 2007) was conducted in the surf zone of Cassino Beach, Southern Brazil, in order to detect possible natural and/or anthropogenic disturbances. Surface water temperature (6–29°C) was the only parameter with predictable seasonal variation; salinity (14–38) was inversely related to rainfall (3.1–485.2 mm month⁻¹) and low values followed extreme precipitation periods in 1997/1998 and 2002/2003 (El Ni?o years). Asterionellopsis glacialis and chlorophyll a presented high concentrations and peak frequency until 1998, when an intense mud deposition occurred with concomitant extreme rainfall. It affected the surf zone and beach, changing the hydrology and dissolved inorganic nutrient availability. Six phytoplankton species groups were recognized with distinct responses to this mud deposition. We conclude that large-scale climatic changes, like El Ni?o Southern Oscillation, in conjunction with human activities significantly altered the phytoplankton ecology of the highly dynamic Cassino Beach surf zone.     

Early warning of rainfall-induced shallow landslides and debris flows in the USA

The state of knowledge and resources available to issue alerts of precipitation-induced landslides vary across the USA. Federal and state agencies currently issue warnings of the potential for shallow, rapidly moving landslides and debris flows in a few areas along the Pacific coast and for areas affected by Atlantic hurricanes. However, these agencies generally lack resources needed to provide continuous support or to expand services to other areas. Precipitation thresholds that form the basis of landslide warning systems now exist for a few areas of the USA, but the threshold rainfall amounts and durations vary over three orders of magnitude nationwide and over an order of magnitude across small geographic areas such as a county. Antecedent moisture conditions also have a significant effect, particularly in areas that have distinct wet and dry seasons. Early warnings of shallow landslides that include specific information about affected areas, probability of landslide occurrence, and expected timing are technically feasible as illustrated by a case study from the Seattle, WA area. The four-level warning scheme (Null, Outlook, Watch, Warning) defined for Seattle is based on observed or predicted exceedance of a cumulative precipitation threshold and a rainfall intensity–duration threshold combined with real-time monitoring of soil moisture. Based on analysis of historical data, threshold performance varies according to precipitation characteristics, and threshold exceedance corresponds to a given probability of landslide occurrence. Experience in Seattle during December 2004 and January 2005 illustrates some of the challenges of providing landslide early warning on the USA West Coast.     

From whom Bells tolls: Reclassifying Bells among CR chondrites and implications for the formation conditions of CR parent bodies

The ungrouped carbonaceous chondrite (CC) Bells has long been considered petrographically similar to CM chondrites based on its matrix abundance and degree of aqueous alteration, but also shows significant isotopic affinities to CR chondrites. Its taxonomic status is thus important for clarifying the relationship of the CRHB (formerly “CR”) clan with other CCs. In this study, we measured the oxygen isotopic compositions of olivines in type I chondrules and isolated olivine grains in Bells. Bells olivines mostly have ∆¹⁷O > −4‰, similar to CR chondrites but unlike other CCs that are rich in refractory inclusions, in which chondrules are generally richer in ¹⁶O. Therefore, Bells is a CR chondrite (albeit an anomalous one), most similar to the rare, matrix-rich CRs like Al Rais. These chondrites (i) may not necessarily derive from the same primary parent body as mainstream CRs, (ii) bear witness to significant variations of the matrix/chondrule ratio within the CRHB clan, and (iii) may be a good analog for samples retrieved by the space mission OSIRIS-REx.