Construction of Structural Geological Model Using Monte Carlo Simulation

Geotechnical and Geological Engineering - To optimize the prediction of structural geological conditions in the underground as of data collected at the surface, due to the usual great uncertainties...     

A cross-ecosystem comparison of spatial and temporal patterns of covariation in the recruitment of functionally analogous fish stocks

Temporal and spatial patterns of recruitment (R) and spawning stock biomass (S) variability were compared among functionally analogous species and similar feeding guilds from six marine ecosystems. Data were aggregated into four regions including the Gulf of Maine/Georges Bank, the Norwegian/Barents Seas, the eastern Bering Sea, and the Gulf of Alaska. Variability was characterized by calculating coefficients of variation and anomalies for three response variables: ln(R), ln(R/S), and stock–recruit model residuals. Patterns of synchrony and asynchrony in the response variables were examined among and between ecosystems, between- and within-ocean basins and among functionally analogous species groups using pair-wise correlation analysis corrected for within-time series autocorrelation, multivariate cross-correlation analyses and regime shift detectors. Time series trends in response variables showed consistent within basin similarities and consistent and coherent differences between the Atlantic and Pacific basin ecosystems. Regime shift detection algorithms identified two broad-scale regime shift time periods for the pelagic feeding guild (1972–1976 and 1999–2002) and possibly one for the benthic feeding guild (1999–2002). No spatial patterns in response variable coefficients of variation were observed. Results from multivariate cross-correlation analysis showed similar trends. The data suggest common external factors act in synchrony on stocks within ocean basins but temporal stock patterns, often of the same species or functional group, between basins change in opposition to each other. Basin-scale results (similar within but different between) suggest that the two geographically broad areas are connected by unknown mechanisms that, depending on the year, may influence the two basins in opposite ways. This work demonstrates that commonalities and synchronies in recruitment fluctuations can be found across geographically distant ecosystems but biophysical causes of the fluctuations remain difficult to identify.     

Glacial isostatic adjustment in Fennoscandia for a laterally heterogeneous earth

Glaciation and deglaciation in Fennoscandia during the last glacial cycles has significantly perturbed the Earth's equilibrium figure. Changes in the Earth's solid and geoidal surfaces due to external and internal mass redistributions are recorded in sequences of ancient coastlines, now either submerged or uplifted, and are still visible in observations of present‐day motions of the surface and glacially induced anomalies in the Earth's gravitational field. These observations become increasingly sophisticated with the availability of GPS measurements and new satellite gravity missions.
Observational evidence of the mass changes is widely used to constrain the radial viscosity structure of the Earth's mantle. However, lateral changes in earth model properties are usually not taken into account, as most global models of glacial isostatic adjustment assume radial symmetry for the earth model. This simplifying assumption contrasts with seismological evidence of significant lateral variations in the Earth's crust and upper mantle throughout the Fennoscandian region.
We compare predictions of glacial isostatic adjustment based on an ice model over the Fennoscandian region for the last glacial cycle for both radially symmetric and fully 3‐D earth models. Our results clearly reveal the importance of lateral variations in lithospheric thickness and asthenospheric viscosity for glacially induced model predictions. Relative sea‐level predictions can differ by up to 10–20 m, uplift rate predictions by 1–3 mm yr⁻¹ and free‐air gravity anomaly predictions by 2–4 mGal when a realistic 3‐D earth structure as proposed by seismic modelling is taken into account.     

Palaeovariations in the East-Asian Monsoon Regime Geochemically Recorded in Varved Sediments of Lake Sihailongwan (Northeast China, Jilin Province). Part 2: a 200-Year Record of Atmospheric Lead-210 Flux Variations and its Palaeoclimatic Implications

Palaeovariations of the atmospheric ²¹⁰Pb flux in Northeast China (Long Gang area, Jilin province), quantified by high resolution ²¹⁰Pb measurements on seasonally laminated sediments of Lake Sihailongwan are presented on decadal scale. The mean flux of unsupported ²¹⁰Pb between 1790 and 1970 for the centre of the lake basin is 517 Bq m⁻² yr⁻¹ with maximum deviations between −23% and +27% of this value. Flux rates above this average were found between 1783 and 1813, around 1836, and between 1860 and 1901. The mean ²¹⁰Pb^exc flux rates derived were clearly lower around 1821, 1908, 1930, and 1953. The ²¹⁰Pb^exc flux reached its minimum of 344 Bq m⁻² yr⁻¹ during the period 1977–1982. The atmospheric flux of unsupported ²¹⁰Pb is correlated with the precipitation frequency during the summer monsoon and shows coinciding variability with geochemical proxies that document the groundwater inflow into the lake. Al₂O₃-rich dust of remote provenance scavenged by wet-deposition in the rainy season is the major carrier of the atmospheric ²¹⁰Pb^exc flux.     

The onset of Pleistocene glaciation in the Barents Sea: implications for glacial isostatic adjustment

In this paper the effect of a delayed onset of glaciation in the Barents Sea on glacial isostatic adjustment is investigated. The model calculations solve the sea-level equation governing the total mass redistributions associated with the last glaciation cycle on a spherically symmetric, linear, Maxwell viscoelastic earth for two different scenarios for the growth phase of the Barents Sea ice sheet. In the first ice model a linear growing history is used for the Barents Sea ice sheet, which closely relates its development to the build-up of other major Late Pleistocene ice sheets. In the second ice model the accumulation of the Barents Sea ice sheet is restricted to the last 6 ka prior to the last glacial maximum.
The calculations predict relative sea levels, present-day radial velocities, and gravity anomalies for the area formerly covered by the Weichselian ice sheet. The results show that observed relative sea levels in the Barents Sea are appropriate for distinguishing between the different glaciation histories. In particular, present-day observables such as the free-air gravity anomaly over the Barents Sea, and the present-day radial velocities are sensitive to changes in the glaciation history on this scale.
A palaeobathymetry derived from relative sea-level predictions before the last glacial maximum based on the second ice model essentially agrees with a palaeobathymetry derived by Lambeck (1995). The additional emerged areas provide centres for the build-up of an ice sheet and thus support the theory of Hald, Danielsen & Lorentzen (1990) and Mangerud et al. (1992) that the Barents Sea was an essentially marine environment shortly before the last glacial maximum.     

TiO2 exsolution from garnet by open-system precipitation: evidence from crystallographic and shape preferred orientation of rutile inclusions

We investigated rutile needles with a clear shape preferred orientation in garnet from (ultra) high-pressure metapelites from the Kimi Complex of the Greek Rhodope by electron microprobe, electron backscatter diffraction and TEM techniques. A definite though complex crystallographic orientation relationship between the garnet host and rutile was identified in that Rt[001] is either parallel to Grt<111> or describes cones with opening angle 27.6° around Grt<111>. Each Rt[001] small circle representing a cone on the pole figure displays six maxima in the density plots. This evidence together with microchemical observations in TEM, when compared to various possible mechanisms of formation, corroborates a precipitate origin. A review of exchange vectors for Ti substitution in garnet indicates that rutile formation from garnet cannot occur in a closed system. It requires that components are exchanged between the garnet interior and the rock matrix by solid-state diffusion, a process we refer to as “open-system precipitation” (OSP). The kinetically most feasible reaction of this type will dominate the overall process. The perhaps most efficient reaction involves internal oxidation of Fe²⁺ to Fe³⁺ and transfer from the dodecahedral to the octahedral site just vacated by $ {\text{Ti}}^{ 4+ }: 6\,{\text{M}}^{ 2+ }_{ 3} {\text{TiAl}}\left[ {{\text{AlSi}}_{ 2} } \right]{\text{O}}_{ 1 2} + 6\,{\text{M}}^{ 2+ }_{ 2, 5} {\text{TiAlSi}}_{ 3} {\text{O}}_{ 1 2} = 10\,{\text{M}}^{ 2+ }_{ 3.0} {\text{Al}}_{ 1. 8} {\text{Fe}}_{0. 2} {\text{Si}}_{ 3} {\text{O}}_{ 1 2} + {\text{M}}^{2+} + 2 {\text{e}}^{-} + 1 2\,{\text{TiO}}_{ 2} . $ OSP is likely to occur at conditions where the transition of natural systems to open-system behaviour becomes apparent, as in the granulite and high-temperature eclogite facies.     

Nuclear and incommensurate magnetic structure of NaFeGe2O6 between 5?K and 298?K and new data on multiferroic NaFeSi2O6

The compound NaFeGe₂O₆ was grown synthetically as polycrystalline powder and as large single crystals suitable for X-ray and neutron-diffraction experiments to clarify the low temperature evolution of secondary structural parameters and to determine the low temperature magnetic spins structure. NaFeGe₂O₆ is isotypic to the clinopyroxene-type compound aegirine and adopts the typical HT-C2/c clinopyroxene structure down to 2.5?K. The Na-bearing M2 polyhedra were identified to show the largest volume expansion between 2.5?K and room temperature, while the GeO₄ tetrahedra behave as stiff units. Magnetic susceptibility measurements show a broad maximum around 33?K, which marks the onset of low-dimensional magnetic ordering. Below 12?K NaFeGe₂O₆ transforms to an incommensurately modulated magnetic spin state, with k?=?[0.323, 1.0, 0.080] and a helical order of spins within the M1-chains of FeO₆ octahedra. This is determined by neutron-diffraction experiments on a single crystal. Comparison of NaFeGe₂O₆ with NaFeSi₂O₆ is given and it is shown that the magnetic ordering in the latter compound, aegirine, also is complex and is best described by two different spin states, a commensurate one with C2??/c?? symmetry and an incommensurate one, best being described by a spin density wave, oriented within the (1 0 1) plane.     

Modeling reactive transport in non-ideal aqueous–solid solution system

The numerical simulation of reactive mass transport processes in complex geochemical environments is an important tool for the performance assessment of future waste repositories. A new combination of the multi-component mass transport code GeoSys/RockFlow and the Gibbs Energy Minimization (GEM) equilibrium solver GEM-Selektor is used to calculate the accurate equilibrium of multiple non-ideal solid solutions which are important for the immobilization of radionuclides such as Ra. The coupled code is verified by a widely used benchmark of dissolution–precipitation in a calcite–dolomite system. A more complex application shown in this paper is the transport of Ra in the near-field of a nuclear waste repository. Depending on the initial inventories of Sr, Ba and sulfate, non-ideal sulfate and carbonate solid solutions can fix mobile Ra cations. Due to the complex geochemical interactions, the reactive transport simulations can describe the migration of Ra in a much more realistic way than using the traditional linear K_D approach only.