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.     

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.     

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.     

Tagging-the-triangle algorithm for partitioning features with inconsistent boundaries

In spatial data sets, gaps or overlaps among features are frequently found in spatial tessellations due to the non-abutting edges with adjacent features. These non-abutting edges in loose tessellations are also called inconsistent boundaries or slivers; polygons containing at least one inconsistent boundary are called inconsistent polygons or sliver polygons. The existing algorithms to solve topological inconsistencies in sliver polygons suffer from one or more of three major issues, namely determination of tolerances, excessive CPU processing time for large data sets and loss of vertex history. In this article, we introduce a new algorithm that mitigates these three issues. Our algorithm efficiently searches the features with inconsistent polygons in a given spatial data set and logically partitions them among adjacent features. The proposed algorithm employs the constrained Delaunay triangulation technique to generate labelled triangles from which inconsistent polygons with gaps and overlaps are identified using label counts. These inconsistent polygons are then partitioned using the straight skeleton method. Moreover, each of these partitioned gaps or overlaps is distributed among the adjacent features to improve the topological consistency of the spatial data sets. We experimentally verified our algorithm using the real land cadastre data set. The comparison results show that the proposed algorithm is four times faster than the existing algorithm for data sets with 200,000 edges.     

Earthquakes as a potential contributing factor to climate change at multi-decadal scale

Climate changes at the multi-decadal scale are often associated with multi-decadal phase shifts of the dominant sea surface temperature （SST） pattern, such as the Pacific Decadal Oscillation （PDO）. The PDO may be associated with the North Pacific branch of the Thermohaline Circulation （THC）. Great earthquakes （M 〉8）, particularly along the route of the THC, might modulate the vertical mixing and bring deep, cold water to surface, contributing to multi-decadal changes in surface currents and the PDO. This may eventually lead to multi-decadal climate changes. We tested this hypothesis for the Pacific Ocean where great earthquakes have been frequently recorded. We found associations between the PDO and recurrent earthquakes along the route of the deep currents of the THC in the modern period since 1900, and relationships between hydroclimate change in Monsoonal Asia and historical earthquakes since 1300. However, it should be noted that this hypothesis is very preliminary and has many gaps that needs further evidences from more observational records and modeling studies.     

An improved processing sequence for uncorrelated Chirp sonar data

Chirp sonar systems can be used to obtain high resolution seismic reflection images of the sub-seafloor during marine surveys. The exact knowledge of the Chirp signature allows the use of deterministic algorithms to process the data, similarly to that applied to Vibroseis data on land. Here, it is described an innovative processing sequence to be applied to uncorrelated Chirp data, which can improve vertical and lateral resolution compared to conventional methods. It includes application of a Wiener filter to transform a frequency-modulated sweep into a minimum-phase pulse sequence. In this way, the data become causal and can undergo predictive deconvolution to reduce ringing and enhance vertical resolution. Afterwards, FX-deconvolution and Stolt migration can be applied to obtain an improved imaging of the subsurface. The result of this procedure is a seismic reflection image with higher resolution than traditional ones, which are normally represented using the envelope function of the signal. This technique can be particularly useful for engineering-geotechnical surveys and archaeological investigations that require a fine detail imaging of the uppermost meters of the sub-seafloor.     

Turbulent flux observations and modelling over a shallow lake and a wet grassland in the Nam Co basin,Tibetan Plateau

The Tibetan Plateau plays an important role in the global water cycle and is strongly influenced by climate change. While energy and matter fluxes have been more intensely studied over land surfaces, a large proportion of lakes have either been neglected or parameterised with simple bulk approaches. Therefore, turbulent fluxes were measured over wet grassland and a shallow lake with a single eddy-covariance complex at the shoreline in the Nam Co basin in summer 2009. Footprint analysis was used to split observations according to the underlying surface, and two sophisticated surface models were utilised to derive gap-free time series. Results were then compared with observations and simulations from a nearby eddy-covariance station over dry grassland, yielding pronounced differences. Observations and footprint integrated simulations compared well, even for situations with flux contributions including grassland and lake. The accessibility problem for EC measurements on lakes can be overcome by combining standard meteorological measurements at the shoreline with model simulations, only requiring representative estimates of lake surface temperature.     

Seasonal prediction of global sea level anomalies using an ocean–atmosphere dynamical model

Advanced warning of extreme sea level events is an invaluable tool for coastal communities, allowing the implementation of management policies and strategies to minimise loss of life and infrastructure damage. This study is an initial attempt to apply a dynamical coupled ocean–atmosphere model to the prediction of seasonal sea level anomalies (SLA) globally for up to 7 months in advance. We assess the ability of the Australian Bureau of Meteorology’s operational seasonal dynamical forecast system, the Predictive Ocean Atmosphere Model for Australia (POAMA), to predict seasonal SLA, using gridded satellite altimeter observation-based analyses over the period 1993–2010 and model reanalysis over 1981–2010. Hindcasts from POAMA are based on a 33-member ensemble of seasonal forecasts that are initialised once per month for the period 1981–2010. Our results show POAMA demonstrates high skill in the equatorial Pacific basin and consistently exhibits more skill globally than a forecast based on persistence. Model predictability estimates indicate there is scope for improvement in the higher latitudes and in the Atlantic and Southern Oceans. Most characteristics of the asymmetric SLA fields generated by El-Nino/La Nina events are well represented by POAMA, although the forecast amplitude weakens with increasing lead-time.     

Geophysical exploration on the subsurface geology of La Garrotxa monogenetic volcanic field (NE Iberian Peninsula)

We applied self-potential (SP) and electrical resistivity tomography (ERT) to the exploration of the uppermost part of the substrate geology and shallow structure of La Garrotxa monogenetic volcanic field, part of the European Neogene–Quaternary volcanic province. The aim of the study was to improve knowledge of the shallowest part of the feeding system of these monogenetic volcanoes and of its relationship with the subsurface geology. This study complements previous geophysical studies carried out at a less detailed scale and aimed at identifying deeper structures, and together will constitute the basis to establish volcanic susceptibility in La Garrotxa. SP study complemented previous smaller-scale studies and targeted key areas where ERT could be conducted. The main new results include the generation of resistivity models identifying dykes and faults associated with several monogenetic cones. The combined results confirm that shallow tectonics controlling the distribution of the foci of eruptive activity in this volcanic zone mainly correspond to NNW–SSE and accessorily by NNE–SSW Neogene extensional fissures and faults and concretely show the associated magmatic intrusions. These structures coincide with the deeper ones identified in previous studies, and show that previous Alpine tectonic structures played no apparent role in controlling the loci of this volcanism. Moreover, the results obtained show that the changes in eruption dynamics occurring at different vents located at relatively short distances in this volcanic area are controlled by shallow stratigraphical, structural and hydrogeological differences underneath these monogenetic volcanoes.