Showing papers in "Solid Earth in 2013"

Kinematics of the South Atlantic rift

Abstract: . The South Atlantic rift basin evolved as a branch of a large Jurassic–Cretaceous intraplate rift zone between the African and South American plates during the final break-up of western Gondwana. While the relative motions between South America and Africa for post-break-up times are well resolved, many issues pertaining to the fit reconstruction and particularly the relation between kinematics and lithosphere dynamics during pre-break-up remain unclear in currently published plate models. We have compiled and assimilated data from these intraplated rifts and constructed a revised plate kinematic model for the pre-break-up evolution of the South Atlantic. Based on structural restoration of the conjugate South Atlantic margins and intracontinental rift basins in Africa and South America, we achieve a tight-fit reconstruction which eliminates the need for previously inferred large intracontinental shear zones, in particular in Patagonian South America. By quantitatively accounting for crustal deformation in the Central and West African Rift Zones, we have been able to indirectly construct the kinematic history of the pre-break-up evolution of the conjugate west African–Brazilian margins. Our model suggests a causal link between changes in extension direction and velocity during continental extension and the generation of marginal structures such as the enigmatic pre-salt sag basin and the Sao Paulo High. We model an initial E–W-directed extension between South America and Africa (fixed in present-day position) at very low extensional velocities from 140 Ma until late Hauterivian times (≈126 Ma) when rift activity along in the equatorial Atlantic domain started to increase significantly. During this initial ≈14 Myr-long stretching episode the pre-salt basin width on the conjugate Brazilian and west African margins is generated. An intermediate stage between ≈126 Ma and base Aptian is characterised by strain localisation, rapid lithospheric weakening in the equatorial Atlantic domain, resulting in both progressively increasing extensional velocities as well as a significant rotation of the extension direction to NE–SW. From base Aptian onwards diachronous lithospheric break-up occurred along the central South Atlantic rift, first in the Sergipe–Alagoas/Rio Muni margin segment in the northernmost South Atlantic. Final break-up between South America and Africa occurred in the conjugate Santos–Benguela margin segment at around 113 Ma and in the equatorial Atlantic domain between the Ghanaian Ridge and the Piaui-Ceara margin at 103 Ma. We conclude that such a multi-velocity, multi-directional rift history exerts primary control on the evolution of these conjugate passive-margin systems and can explain the first-order tectonic structures along the South Atlantic and possibly other passive margins.

Managing soil nitrate with cover crops and buffer strips in Sicilian vineyards

Abstract: . When soil nitrate levels are low, plants suffer nitrogen (N) deficiency but when the levels are excessive, soil nitrates can pollute surface and subsurface waters. Strategies to reduce the nitrate pollution are necessary to reach a sustainable use of resources such as soil, water and plant. Buffer strips and cover crops can contribute to the management of soil nitrates, but little is known of their effectiveness in semiarid vineyards plantations. The research was carried out in the south coast of Sicily (Italy) to evaluate nitrate trends in a vineyard managed both conventionally and using two different cover crops (Triticum durum and Vicia sativa cover crop). A 10 m-wide buffer strip was seeded with Lolium perenne at the bottom of the vineyard. Soil nitrate was measured monthly and nitrate movement was monitored by application of a 15N tracer to a narrow strip between the bottom of vineyard and the buffer and non-buffer strips. Lolium perenne biomass yield in the buffer strips and its isotopic nitrogen content were monitored. Vicia sativa cover crop management contributed with an excess of nitrogen, and the soil management determined the nitrogen content at the buffer areas. A 6 m buffer strip reduced the nitrate by 42% with and by 46% with a 9 m buffer strip. Thanks to catch crops, farmers can manage the N content and its distribution into the soil over the year, can reduced fertilizer wastage and reduce N pollution of surface and groundwater.

Spatial models for monitoring the spatio-temporal evolution of ashes after fire – a case study of a burnt grassland in Lithuania

Abstract: . Ash thickness is a key variable in the protection of soil against erosion agents after planned and unplanned fires. Ash thickness measurements were conducted along two transects (flat and sloping areas) following a grided experimental design. In order to interpolate data with accuracy and identify the techniques with the least bias, several interpolation methods were tested in the grided plot. Overall, the fire had a low severity. However, the fire significantly reduced the ground cover, especially on sloping areas, owing to the higher fire severity and/or less biomass previous to the fire. Ash thickness depended on fire severity and was thin where fire severity was higher and thicker in lower fire severity sites. The ash thickness decreased with time after the fire. Between 4 and 16 days after the fire, ash was transported by wind. The greatest reduction took place between 16 and 34 days after the fire as a result of rainfall, and was more efficient where fire severity was higher. Between 34 and 45 days after the fire, no significant differences in ash thickness were identified among ash colours and only traces of the ash layer remained. The omni-directional experimental variograms showed that variable structure did not change significantly with time. The ash spatial variability increased with time, particularly on the slope, as a result of water erosion.

Impacts of land use change in soil carbon and nitrogen in a Mediterranean agricultural area (Southern Spain)

Abstract: . The agricultural Mediterranean areas are dedicated to arable crops (AC), but in the last decades, a significant number of AC has led to a land use change (LUC) to olive grove (OG) and vineyards (V). A field study was conducted to determine the long-term effects (46 years) of LUC (AC by OG and V) and to determine soil organic carbon (SOC), total nitrogen (TN), C : N ratio and their stratification across the soil entire profile, in Montilla-Moriles denomination of origin (D.O.), in Calcic-Chromic Luvisols (LVcc/cr), an area under semiarid Mediterranean conditions. The experimental design consisted of studying the LUC on one farm between 1965 and 2011. Originally, only AC was farmed in 1965, but OG and V were farmed up to now (2011). This LUC principally affected the horizon thickness, texture, bulk density, pH, organic matter, organic carbon, total nitrogen and C : N ratio. The LUC had a negative impact in the soil, affecting the SOC and TN stocks. The conversion from AC to V and OG involved the loss of the SOC stock (52.7% and 64.9% to V and OG respectively) and the loss of the TN stock (42.6% and 38.1% to V and OG respectively). With respect to the stratification ratios (SRs), the effects were opposite; 46 years after LUC increased the SRs (in V and OG) of SOC, TN and C : N ratio.

Exhumation of (ultra-)high-pressure terranes: concepts and mechanisms

Abstract: . The formation and exhumation of high and ultra-high-pressure, (U)HP, rocks of crustal origin appears to be ubiquitous during Phanerozoic plate subduction and continental collision events. Exhumation of (U)HP material has been shown in some orogens to have occurred only once, during a single short-lived event; in other cases exhumation appears to have occurred multiple discrete times or during a single, long-lived, protracted event. It is becoming increasingly clear that no single exhumation mechanism dominates in any particular tectonic environment, and the mechanism may change in time and space within the same subduction zone. Subduction zone style and internal force balance change in both time and space, responding to changes in width, steepness, composition of subducting material and velocity of subduction. In order for continental crust, which is relatively buoyant compared to the mantle even when metamorphosed to (U)HP assemblages, to be subducted to (U)HP conditions, it must remain attached to a stronger and denser substrate. Buoyancy and external tectonic forces drive exhumation, although the changing spatial and temporal dominance of different driving forces still remains unclear. Exhumation may involve whole-scale detachment of the terrane from the subducting slab followed by exhumation within a subduction channel (perhaps during continued subduction) or a reversal in motion of the entire plate (eduction) following the removal of a lower part of the subducting slab. Weakening mechanisms that may be responsible for the detachment of deeply subducted crust from its stronger, denser substrate include strain weakening, hydration, melting, grain size reduction and the development of foliation. These may act locally to form narrow high-strain shear zones separating stronger, less-strained crust or may act on the bulk of the subducted material, allowing whole-scale flow. Metamorphic reactions, metastability and the composition of the subducted crust all affect buoyancy and overall strength. Future research directions include identifying temporal and spatial changes in exhumation mechanisms within different tectonic environments, and determining the factors that influence those changes.

Seasonal changes in the soil hydrological and erosive response depending on aspect, vegetation type and soil water repellency in different Mediterranean microenvironments

Abstract: . Mediterranean areas are characterized by a strong spatial variability that makes the soil hydrological response highly complex. Moreover, Mediterranean climate has marked seasons that provoke dramatic changes on soil properties determining the runoff rates, such as soil water content or soil water repellency (SWR). Thus, soil hydrological and erosive response in Mediterranean areas can be highly time- as well as space-dependant. This study shows SWR, aspect and vegetation as factors of the soil hydrological and erosive response. Erosion plots were set up in the north- and the south-facing hillslope and rainfall, runoff, sediments and SWR were monitored. Soil water repellency showed a seasonal behaviour and it was presented in three out of four microenvironments after the summer, disappearing in the wet season. In general, runoff rate was higher in shrubs patches (0.47 ± 0.67 mm) than in inter-shrub soils (1.54 ± 2.14 mm), but it changed seasonally in different ways, depending on the aspect considered, decreasing in the north-facing hillslope and increasing in the south-facing one. The main factor determining the hydrological and erosive response was the rainfall intensity, regardless of the rainfall depth of the event. This response was modulated mainly by SWR in the north-facing hillslope and the vegetation pattern in the south-facing one.

Full-fit reconstruction of the Labrador Sea and Baffin Bay

Abstract: . Reconstructing the opening of the Labrador Sea and Baffin Bay between Greenland and North America remains controversial. Recent seismic data suggest that magnetic lineations along the margins of the Labrador Sea, originally interpreted as seafloor spreading anomalies, may lie within the crust of the continent–ocean transition. These data also suggest a more seaward extent of continental crust within the Greenland margin near Davis Strait than assumed in previous full-fit reconstructions. Our study focuses on reconstructing the full-fit configuration of Greenland and North America using an approach that considers continental deformation in a quantitative manner. We use gravity inversion to map crustal thickness across the conjugate margins, and assimilate observations from available seismic profiles and potential field data to constrain the likely extent of different crustal types. We derive end-member continental margin restorations following alternative interpretations of published seismic profiles. The boundaries between continental and oceanic crust (COB) are restored to their pre-stretching locations along small circle motion paths across the region of Cretaceous extension. Restored COBs are fitted quantitatively to compute alternative total-fit reconstructions. A preferred full-fit model is chosen based on the strongest compatibility with geological and geophysical data. Our preferred model suggests that (i) the COB lies oceanward of magnetic lineations interpreted as magnetic anomaly 31 (70 Ma) in the Labrador Sea, (ii) all previously identified magnetic lineations landward of anomaly 27 reflect intrusions into continental crust and (iii) the Ungava fault zone in Davis Strait acted as a leaky transform fault during rifting. This robust plate reconstruction reduces gaps and overlaps in Davis Strait and suggests that there is no need for alternative models proposed for reconstructions of this area including additional plate boundaries in North America or Greenland. Our favoured model implies that break-up and formation of continent–ocean transition (COT) first started in the southern Labrador Sea and Davis Strait around 88 Ma and then propagated north and southwards up to the onset of real seafloor spreading at 63 Ma in the Labrador Sea. In Baffin Bay, continental stretching lasted longer and actual break-up and seafloor spreading started around 61 Ma (chron 26).

Extreme extension across Seram and Ambon, eastern Indonesia: evidence for Banda slab rollback

Abstract: . The island of Seram, which lies in the northern part of the 180°-curved Banda Arc, has previously been interpreted as a fold-and-thrust belt formed during arc-continent collision, which incorporates ophiolites intruded by granites thought to have been produced by anatexis within a metamorphic sole. However, new geological mapping and a re-examination of the field relations cause us to question this model. We instead propose that there is evidence for recent and rapid N–S extension that has caused the high-temperature exhumation of lherzolites beneath low-angle lithospheric detachment faults that induced high-temperature metamorphism and melting in overlying crustal rocks. These "Kobipoto Complex" migmatites include highly residual Al–Mg-rich garnet + cordierite + sillimanite + spinel + corundum granulites (exposed in the Kobipoto Mountains) which contain coexisting spinel + quartz, indicating that peak metamorphic temperatures likely approached 900 °C. Associated with these residual granulites are voluminous Mio-Pliocene granitic diatexites, or "cordierite granites", which crop out on Ambon, western Seram, and in the Kobipoto Mountains and incorporate abundant schlieren of spinel- and sillimanite-bearing residuum. Quaternary "ambonites" (cordierite + garnet dacites) emplaced on Ambon were also evidently sourced from the Kobipoto Complex migmatites as demonstrated by granulite-inherited xenoliths. Exhumation of the hot peridotites and granulite-facies Kobipoto Complex migmatites to shallower structural levels caused greenschist- to lower-amphibolite facies metapelites and amphibolites of the Tehoru Formation to be overprinted by sillimanite-grade metamorphism, migmatisation, and limited localised anatexis to form the Taunusa Complex. The extreme extension required to have driven Kobipoto Complex exhumation evidently occurred throughout Seram and along much of the northern Banda Arc. The lherzolites must have been juxtaposed against the crust at typical lithospheric mantle temperatures in order to account for such high-temperature metamorphism and therefore could not have been part of a cooled ophiolite. In central Seram, lenses of peridotites are incorporated with a major left-lateral strike-slip shear zone (the "Kawa Shear Zone"), demonstrating that strike-slip motions likely initiated shortly after the mantle had been partly exhumed by detachment faulting and that the main strike-slip faults may themselves be reactivated and steepened low-angle detachments. The geodynamic driver for mantle exhumation along the detachment faults and strike-slip faulting in central Seram is very likely the same; we interpret the extreme extension to be the result of eastward slab rollback into the Banda Embayment as outlined by the latest plate reconstructions for Banda Arc evolution.

A database of plagioclase crystal preferred orientations (CPO) and microstructures - implications for CPO origin, strength, symmetry and seismic anisotropy in gabbroic rocks

Abstract: . This study presents a unique database of 172 plagioclase Crystallographic Preferred Orientations (CPO) of variously deformed gabbroic rocks. The CPO characteristics as a function of the deformation regime (magmatic or crystal-plastic) are outlined and discussed. The studied samples are dominantly from slow- and fast-spread present-day ocean crust, as well as from the Oman ophiolite. Plagioclase is the dominant mineral phase in the studied samples. Plagioclase CPOs are grouped into three main categories: Axial-B, a strong point alignment of (010) with a girdle distribution of [100]; Axial-A, a strong point maximum concentration of [100] with parallel girdle distributions of (010) and (001); and P-type, point maxima of [100], (010), and (001). A majority of CPO patterns are Axial-B and P-type, in samples showing either magmatic or crystal-plastic deformation textures. Axial-A CPOs are less common; they represent 21% of the samples deformed by crystal-plastic flow. Although fabric strength (ODF J index) does not show any consistent variation as a function of the CPO patterns, there is a significant difference in the relationship between the ODF and pole figures J indices; the magmatic type microstructures have high (010) pole figures J indices, which increase linearly with ODF J index, whereas the high [100] pole figures J indices of plastically deformed samples vary in a more scattered manner with ODF J index. The multistage nature of plastic deformation superposed on a magmatic structure compared with magmatic flow, and the large number of possible slip-systems in plagioclase probably account for these differences. Calculated seismic properties (P wave and S wave velocities and anisotropies) of plagioclase aggregates show that anisotropy (up to 12% for P wave and 14% for S wave) tends to increase as a function of ODF J index. In comparison with the olivine 1998 CPO database, the magnitude of P wave anisotropy for a given J index is much less than olivine, whereas it is similar for S wave anisotropy. Despite a large variation of fabric patterns and geodynamic setting, seismic properties of plagioclase-rich rocks have similar magnitudes of anisotropy. There is a small difference in the aggregate elastic symmetry, with magmatic microstructures having higher orthorhombic and hexagonal components, whereas plastic deformation microstructures have a slightly higher monoclinic component, possibly correlated with predominant monoclinic simple shear flow in plastically deformed samples. Overall, plots for CPO strength (ODF J index), pole figure strength, CPO symmetry and seismic anisotropy show significant scattering. This could be related to sampling statistics, although our database is a factor of ten higher than the olivine database of 1998, or it could be related to the low symmetry (triclinic) structure of plagioclase resulting in the addition of degrees of freedom in the processes creating the CPOs.

Application of titanium-in-quartz thermobarometry to greenschist facies veins and recrystallized quartzites in the Hsüehshan range, Taiwan

Abstract: . The accuracy, reliability and best practises of Ti-in-quartz thermobarometry (TitaniQ) in greenschist facies rocks have not been established. To address these issues, we measured Ti concentrations in rutile-bearing samples of moderately deformed, partially recrystallized quartzite and vein quartz from the Hsuehshan range, Taiwan. The spread of Ti concentrations of recrystallized grains in quartzite correlates with recrystallized grain size. Recrystallized quartz (grain size ~100–200 μm) that formed during early deformation within the biotite stability field shows a marked increase in intermediate Ti-concentration grains (~1–10 ppm) relative to detrital porphyroclasts (Ti ~0.1–200 ppm). Fine recrystallized quartz (~5% of the samples by area, grain size ~10–20 μm) has a further restricted Ti concentration peaking at 0.8–2 ppm. This trend suggests equilibration of Ti in recrystallized quartz with a matrix phase during deformation and cooling. Unlike previously documented examples, Ti concentration in the quartzite is inversely correlated with blue cathodoluminescence. Deformation was associated with a minimum grain boundary diffusivity of Ti on the order of 10−22m2 s−1. Vein emplacement and quartzite recrystallization are independently shown to have occurred at 250–350 °C and 300–410 °C, respectively, with lithostatic pressure of 3–4 kbar (assuming a geothermal gradient of 25° km−1), and with hydrostatic fluid pressure. Estimates of the accuracy of TitaniQ at these conditions depend on whether lithostatic or fluid pressure is used in the TitaniQ calibration. Using lithostatic pressure and these temperatures, the Thomas et al. (2010) calibration yields Ti concentrations within error of concentrations measured by SIMS. If fluid pressure is instead used, predicted temperatures are ~30–40 °C too low. TitaniQ has potential to yield accurate PT information for vein emplacement and dynamic recrystallization of quartz at temperatures as low as ~250 °C, however clarification of the relevant pressure term and further tests in rutile-present rocks are warranted.

A workflow for building and calibrating 3-D geomechanical models &ndash a case study for a gas reservoir in the North German Basin

Abstract: . The optimal use of conventional and unconventional hydrocarbon reservoirs depends, amongst other things, on the local tectonic stress field. For example, wellbore stability, orientation of hydraulically induced fractures and – especially in fractured reservoirs – permeability anisotropies are controlled by the present-day in situ stresses. Faults and lithological changes can lead to stress perturbations and produce local stresses that can significantly deviate from the regional stress field. Geomechanical reservoir models aim for a robust, ideally "pre-drilling" prediction of the local variations in stress magnitude and orientation. This requires a numerical modelling approach that is capable to incorporate the specific geometry and mechanical properties of the subsurface reservoir. The workflow presented in this paper can be used to build 3-D geomechanical models based on the finite element (FE) method and ranging from field-scale models to smaller, detailed submodels of individual fault blocks. The approach is successfully applied to an intensively faulted gas reservoir in the North German Basin. The in situ stresses predicted by the geomechanical FE model were calibrated against stress data actually observed, e.g. borehole breakouts and extended leak-off tests. Such a validated model can provide insights into the stress perturbations in the inter-well space and undrilled parts of the reservoir. In addition, the tendency of the existing fault network to slip or dilate in the present-day stress regime can be addressed.

Grassland fire effect on soil organic carbon reservoirs in a semiarid environment

Abstract: . The aim of this work was to investigate the effect of an experimental fire used for grassland management on soil organic carbon (SOC) stocks. The study was carried out on Hyparrhenia hirta (L.) Stapf (Hh) grassland and Ampelodesmos mauritanicus (Desf.) T. Durand & Schinz (Am) grasslands located in the north of Sicily. Soil samples were collected at 0–5 cm before and after the experimental fire, and SOC was measured. During the grassland fire, soil surface temperature was monitored. Biomass of both grasses was analysed in order to determine dry weight and its chemical composition. The results showed that SOC varied significantly with vegetation type, while it is not affected in the short term by grassland fire. Am grassland stored more SOC compared with Hh grassland thanks to lower content in the biomass of the labile carbon pool. No significant difference was observed in SOC before and after fire, which could be caused by several factors: first, in both grassland types the measured soil temperature during fire was low due to thin litter layers; second, in a semiarid environment, a higher mineralization rate results in a lower soil carbon labile pool; and third, the SOC stored in the finest soil fractions, physically protected, is not affected by fire.

Application of soil quality indices to assess the status of agricultural soils irrigated with treated wastewaters

Abstract: . The supply of water is limited in some parts of the Mediterranean region, such as southeastern Spain. The use of treated wastewater for the irrigation of agricultural soils is an alternative to using better-quality water, especially in semi-arid regions. On the other hand, this practice can modify some soil properties, change their relationships and influence soil quality. In this work two soil quality indices were used to evaluate the effects of irrigation with treated wastewater in soils. The indices were developed studying different soil properties in undisturbed soils in SE Spain, and the relationships between soil parameters were established using multiple linear regressions. These indices represent the balance reached among properties in "steady state" soils. This study was carried out in four study sites from SE Spain irrigated with wastewater, including four study sites. The results showed slight changes in some soil properties as a consequence of irrigation with wastewater, the obtained levels not being dangerous for agricultural soils, and in some cases they could be considered as positive from an agronomical point of view. In one of the study sites, and as a consequence of the low quality wastewater used, a relevant increase in soil organic matter content was observed, as well as modifications in most of the soil properties. The application of soil quality indices indicated that all the soils of study sites are in a state of disequilibrium regarding the relationships between properties independent of the type of water used. However, there were no relevant differences in the soil quality indices between soils irrigated with wastewater with respect to their control sites for all except one of the sites, which corresponds to the site where low quality wastewater was used.

Post-processing scheme for modelling the lithospheric magnetic field

Abstract: . We investigated how the noise in satellite magnetic data affects magnetic lithospheric field models derived from these data in the special case where this noise is correlated along satellite orbit tracks. For this we describe the satellite data noise as a perturbation magnetic field scaled independently for each orbit, where the scaling factor is a random variable, normally distributed with zero mean. Under this assumption, we have been able to derive a model for errors in lithospheric models generated by the correlated satellite data noise. Unless the perturbation field is known, estimating the noise in the lithospheric field model is a non-linear inverse problem. We therefore proposed an iterative post-processing technique to estimate both the lithospheric field model and its associated noise model. The technique has been successfully applied to derive a lithospheric field model from CHAMP satellite data up to spherical harmonic degree 120. The model is in agreement with other existing models. The technique can, in principle, be extended to all sorts of potential field data with "along-track" correlated errors.

Segmentation of the Izu-Bonin and Mariana slabs based on the analysis of the Benioff seismicity distribution and regional tomography results

Abstract: . We present a new model of P and S velocity anomalies in the mantle down to a depth of 1300 km beneath the Izu-Bonin and Mariana (IBM) arcs. This model is derived based on tomographic inversion of global travel time data from the revised ISC catalogue. The results of inversion are thoroughly verified using a series of different tests. The obtained model is generally consistent with previous studies by different authors. We also present the distribution of relocated deep events projected to the vertical surface along the IBM arc system. Unexpectedly, the seismicity forms elongated vertical clusters instead of horizontal zones indicating phase transitions in the slab. We propose that these vertical seismicity zones mark zones of intense deformation and boundaries between semi-autonomous segments of the subducting plate. The P and S seismic tomography models consistently display the slab as prominent high-velocity anomalies coinciding with the distribution of deep seismicity. We can distinguish at least four segments which subduct differently. The northernmost segment of the Izu-Bonin arc has the gentlest angle of dipping which is explained by backward displacement of the trench. In the second segment, the trench stayed at the same location, and we observe the accumulation of the slab material in the transition zone and its further descending to the lower mantle. In the third segment, the trench is moving forward causing the steepening of the slab. Finally, for the Mariana segment, despite the backward displacement of the arc, the subducting slab is nearly vertical. Between the Izu-Bonin and Mariana arcs we clearly observe a gap which can be traced down to about 400 km in depth. Based on joint consideration of the tomography results and the seismicity distribution, we propose two different scenarios of the subduction evolution in the IBM zone during the recent time, depending on the reference frame of plate displacements. In the first case, we consider the movements in respect to the Philippine Plate, and explain the different styles of the subduction by the relative backward and forward migrations of the trench. In the second case, all the elements of the subduction system move westward in respect to the stable Asia. Different subduction styles are explained by the "anchoring" of selected segments of the slab, different physical properties of the subducting plate and the existence of buoyant rigid blocks related to sea mount and igneous provinces.

The dynamics of laterally variable subductions: laboratory models applied to the Hellenides

Abstract: . We designed three-dimensional dynamically self-consistent laboratory models of subduction to analyse the relationships between overriding plate deformation and subduction dynamics in the upper mantle. We investigated the effects of the subduction of a lithosphere of laterally variable buoyancy on the temporal evolution of trench kinematics and shape, horizontal flow at the top of the asthenosphere, dynamic topography and deformation of the overriding plate. Two subducting units, which correspond to a negatively buoyant oceanic plate and positively buoyant continental one, are juxtaposed via a trench-perpendicular interface (analogue to a tear fault) that is either fully-coupled or shear-stress free. Differential rates of trench retreat, in excess of 6 cm yr−1 between the two units, trigger a more vigorous mantle flow above the oceanic slab unit than above the continental slab unit. The resulting asymmetrical sublithospheric flow shears the overriding plate in front of the tear fault, and deformation gradually switches from extension to transtension through time. The consistency between our models results and geological observations suggests that the Late Cenozoic deformation of the Aegean domain, including the formation of the North Aegean Trough and Central Hellenic Shear zone, results from the spatial variations in the buoyancy of the subducting lithosphere. In particular, the lateral changes of the subduction regime caused by the Early Pliocene subduction of the old oceanic Ionian plate redesigned mantle flow and excited an increasingly vigorous dextral shear underneath the overriding plate. The models suggest that it is the inception of the Kefalonia Fault that caused the transition between an extension dominated tectonic regime to transtension, in the North Aegean, Mainland Greece and Peloponnese. The subduction of the tear fault may also have helped the propagation of the North Anatolian Fault into the Aegean domain.

Short-wavelength undulatory extinction in quartz recording coseismic deformation in the middle crust – an experimental study

Abstract: . Deformation experiments are carried out on natural vein quartz in a modified Griggs-type solid medium apparatus to explore the preservation potential of microfabrics created by crystal-plastic deformation at high stress, overprinted during subsequent creep at lower stress. A corresponding stress history is expected for the upper plastosphere, where fault slip during an earthquake causes quasi-instantaneous loading to high stress, followed by stress relaxation. The question is whether evidence of crystal-plastic deformation at high stress, hence an indicator of past seismic activity, can still be identified in the microstructure after overprint by creep at lower stresses. First, quartz samples are deformed at a temperature of 400 °C and constant strain rate of 10−4 s−1 ("kick"), and then held at 900 to 1000 °C at residual stress ("creep"). In quartz exclusively subject to high-stress deformation, lamellar domains of slightly differing crystallographic orientation (misorientation angle

Quantification of magma ascent rate through rockfall monitoring at the growing/collapsing lava dome of Volcán de Colima, Mexico

Abstract: . The most recent eruptive phase of Volcan de Colima, Mexico, started in 1998 and was characterized by dome growth with a variable effusion rate, interrupted intermittently by explosive eruptions. Between November 2009 and June 2011, activity at the dome was mostly limited to a lobe on the western side where it had previously started overflowing the crater rim, leading to the generation of rockfall events. As a consequence of this, no significant increase in dome volume was perceivable and the rate of magma ascent, a crucial parameter for volcano monitoring and hazard assessment could no longer be quantified via measurements of the dome's dimensions. Here, we present alternative approaches to quantify the magma ascent rate. We estimate the volume of individual rockfalls through the detailed analysis of sets of photographs (before and after individual rockfall events). The relationship between volume and infrared images of the freshly exposed dome surface and the seismic signals related to the rockfall events were then investigated. Larger rockfall events exhibited a correlation between its previously estimated volume and the surface temperature of the freshly exposed dome surface, as well as the mean temperature of rockfall mass distributed over the slope. We showed that for larger events, the volume of the rockfall correlates with the maximum temperature of the newly exposed lava dome as well as a proxy for seismic energy. It was therefore possible to calibrate the seismic signals using the volumes estimated from photographs and the count of rockfalls over a certain period was used to estimate the magma extrusion flux for the period investigated. Over the course of the measurement period, significant changes were observed in number of rockfalls, rockfall volume and hence averaged extrusion rate. The extrusion rate was not constant: it increased from 0.008 ± 0.003 to 0.02 ± 0.007 m3 s−1 during 2010 and dropped down to 0.008 ± 0.003 m3 s−1 again in March 2011. In June 2011, magma extrusion had come to a halt. The methodology presented represents a reliable tool to constrain the growth rate of domes that are repeatedly affected by partial collapses. There is a good correlation between thermal and seismic energies and rockfall volume. Thus it is possible to calibrate the seismic records associated with the rockfalls (a continuous monitoring tool) to improve volcano monitoring at volcanoes with active dome growth.

3-D reflection seismic imaging of the Hontomín structure in the Basque–Cantabrian Basin (Spain)

Abstract: . The Basque–Cantabrian Basin of the northern Iberia Peninsula constitutes a unique example of a major deformation system, featuring a dome structure developed by extensional tectonics followed by compressional reactivation. The occurrence of natural resources in the area and the possibility of establishing a geological storage site for carbon dioxide motivated the acquisition of a 3-D seismic reflection survey in 2010, centered on the Jurassic Hontomin dome. The objectives of this survey were to obtain a geological model of the overall structure and to establish a baseline model for a possible geological CO2 storage site. The 36 km2 survey included approximately 5000 mixed (Vibroseis and explosives) source points recorded with a 25 m inline source and receiver spacing. The target reservoir is a saline aquifer, at approximately 1450 m depth, encased and sealed by carbonate formations. Acquisition and processing parameters were influenced by the rough topography and relatively complex geology. A strong near-surface velocity inversion is evident in the data, affecting the quality of the data. The resulting 3-D image provides constraints on the key features of the geologic model. The Hontomin structure is interpreted to consist of an approximately 107 m2 large elongated dome with two major (W–E and NW–SE) striking faults bounding it. Preliminary capacity estimates indicate that about 1.2 Gt of CO2 can be stored in the target reservoir.

Experimental deformation and recrystallization of olivine – processes and timescales of damage healing during postseismic relaxation at mantle depths

Abstract: . Experiments comprising sequences of deformation (at 300 or 600 °C) and annealing at varying temperature (700 to 1100 °C), time (up to 144 h) and stress (up to 1.5 GPa) were carried out in a Griggs-type apparatus on natural olivine-rich peridotite samples to simulate deformation and recrystallization processes in deep shear zones that reach mantle depth as continuations of seismically active faults. The resulting olivine microfabrics were analysed by polarization and electron microscopy (SEM/EBSD, TEM). Core-and-mantle-like microstructures are the predominant result of our experiments simulating rapid stress relaxation (without or with minor creep) after a high-stress deformation event: porphyroclasts (> 100 μm) are surrounded by new grains comprising fragments and recrystallized grains with a wide range in size (2 to 40 μm). Areas with small grains (l 10 μm) trace former high-strain zones generated during initial high-stress deformation even after annealing at a temperature of 1100 °C for 70 h. A weak crystallographic preferred orientation (CPO) of new olivine grains is related to the orientation of the original host crystals but appears unrelated to the strain field. Based on these findings, we propose that olivine microstructures in natural shear-zone peridotites with a large range in grain size, localized fine-grained zones, and a weak CPO not related to the strain field are diagnostic for a sequence of high-stress deformation followed by recrystallization at low stresses, as to be expected in areas of seismic activity. We extended the classic Avrami-kinetics equation by accounting for time-dependent growth kinetics and constrained the involved parameters relying on our results and previous studies devoted to the kinetics of defect processes in olivine. Extrapolation to natural conditions suggests that the observed characteristic microstructure may develop within as little as tens of years and less than ten thousands of years. These recrystallization microstructures have a great diagnostic potential for past seismic activity because they are expected to be stable over geological timescales, since driving forces for further modification are not sufficient to erase the characteristic heterogeneities.

Reprocessed height time series for GPS stations

Abstract: . Precise weekly positions of 403 Global Positioning System (GPS) stations located worldwide are obtained by reprocessing GPS data of these stations for the time span from 4 January 1998 until 29 December 2007. The processing algorithms and models used as well as the solution and results obtained are presented. Vertical velocities of 266 GPS stations having a tracking history longer than 2.5 yr are computed; 107 of them are GPS stations located at tide gauges (TIGA observing stations). The vertical velocities calculated in this study are compared with the estimates from the co-located tide gauges and other GPS solutions. The formal errors of the estimated vertical velocities are 0.01–0.80 mm yr−1. The vertical velocities of our solution agree within 1 mm yr−1 with those of the recent solutions (ULR5 and ULR3) of the Universite de La Rochelle for about 67–75 per cent of the common stations. Examples of typical behaviour of station height changes are given and interpreted. The derived height time series and vertical motions of continuous GPS at tide gauges stations can be used for correcting the vertical land motion in tide gauge records of sea level changes.

An objective rationale for the choice of regularisation parameter with application to global multiple-frequency S -wave tomography

Abstract: . In a linear ill-posed inverse problem, the regularisation parameter (damping) controls the balance between minimising both the residual data misfit and the model norm. Poor knowledge of data uncertainties often makes the selection of damping rather arbitrary. To go beyond that subjectivity, an objective rationale for the choice of damping is presented, which is based on the coherency of delay-time estimates in different frequency bands. Our method is tailored to the problem of global multiple-frequency tomography (MFT), using a data set of 287 078 S-wave delay times measured in five frequency bands (10, 15, 22, 34, and 51 s central periods). Whereas for each ray path the delay-time estimates should vary coherently from one period to the other, the noise most likely is not coherent. Thus, the lack of coherency of the information in different frequency bands is exploited, using an analogy with the cross-validation method, to identify models dominated by noise. In addition, a sharp change of behaviour of the model e∞-norm, as the damping becomes lower than a threshold value, is interpreted as the signature of data noise starting to significantly pollute at least one model component. Models with damping larger than this threshold are diagnosed as being constructed with poor data exploitation. Finally, a preferred model is selected from the remaining range of permitted model solutions. This choice is quasi-objective in terms of model interpretation, as the selected model shows a high degree of similarity with almost all other permitted models (correlation superior to 98% up to spherical harmonic degree 80). The obtained tomographic model is displayed in the mid lower-mantle (660–1910 km depth), and is shown to be compatible with three other recent global shear-velocity models. A wider application of the presented rationale should permit us to converge towards more objective seismic imaging of Earth's mantle.

Seismic imaging of sandbox experiments – laboratory hardware setup and first reflection seismic sections

Abstract: . With the study and technical development introduced here, we combine analogue sandbox simulation techniques with seismic physical modelling of sandbox models. For that purpose, we designed and developed a new mini-seismic facility for laboratory use, comprising a seismic tank, a PC-driven control unit, a positioning system, and piezoelectric transducers used here for the first time in an array mode. To assess the possibilities and limits of seismic imaging of small-scale structures in sandbox models, different geometry setups were tested in the first 2-D experiments that also tested the proper functioning of the device and studied the seismo-elastic properties of the granular media used. Simple two-layer models of different materials and layer thicknesses as well as a more complex model comprising channels and shear zones were tested using different acquisition geometries and signal properties. We suggest using well sorted and well rounded grains with little surface roughness (glass beads). Source receiver-offsets less than 14 cm for imaging structures as small as 2.0–1.5 mm size have proven feasible. This is the best compromise between wide beam and high energy output, and is applicable with a consistent waveform. Resolution of the interfaces of layers of granular materials depends on the interface preparation rather than on the material itself. Flat grading of interfaces and powder coverage yields the clearest interface reflections. Finally, sandbox seismic sections provide images of high quality showing constant thickness layers as well as predefined channel structures and indications of the fault traces from shear zones. Since these were artificially introduced in our test models, they can be regarded as zones of disturbance rather than tectonic shear zones characterized by decompaction. The multiple-offset surveying introduced here, improves the quality with respect to S / N ratio and source signature even more; the maximum depth penetration in glass-bead layers thereby amounts to 5 cm. Thus, the presented mini-seismic device is already able to resolve structures within simple models of saturated porous media, so that multiple-offset seismic imaging of shallow sandbox models, that are structurally evolving, is generally feasible.

Paleosols in the Transantarctic Mountains: indicators of environmental change

Abstract: . The Transantarctic Mountains (TAMs), a 3500 km long chain that subdivides East Antarctica from West Antarctica, are important for reconstructing the tectonic, glacial, and climatic history of Antarctica. With an ice-free area of 24 200 km2 (50% of the total in Antarctica), the TAMs contain an unusually high proportion of paleosols, including relict and buried soils. The unconsolidated paleosols range from late Quaternary to Miocene in age, the semi-consolidated paleosols are of early Miocene to Oligocene age, and the consolidated paleosols are of Paleozoic age. Paleosols on unconsolidated deposits are emphasized in this study. Examples are given from the McMurdo Dry Valleys (78° S) and two outlet glaciers in the central and southern TAMS, including the Hatherton–Darwin Glacier region (80° S) and the Beardmore Glacier region (85°30' S). Relict soils constitute 73% of all of the soils examined; 10% of the soils featured burials. About 26% of the soils examined are from the last glaciation (

Two-dimensional numerical investigations on the termination of bilinear flow in fractures

Abstract: . Bilinear flow occurs when fluid is drained from a permeable matrix by producing it through an enclosed fracture of finite conductivity intersecting a well along its axis. The terminology reflects the combination of two approximately linear flow regimes: one in the matrix with flow essentially perpendicular to the fracture, and one along the fracture itself associated with the non-negligible pressure drop in it. We investigated the characteristics, in particular the termination, of bilinear flow by numerical modeling allowing for an examination of the entire flow field without prescribing the flow geometry in the matrix. Fracture storage capacity was neglected relying on previous findings that bilinear flow is associated with a quasi-steady flow in the fracture. Numerical results were generalized by dimensionless presentation. Definition of a dimensionless time that, other than in previous approaches, does not use geometrical parameters of the fracture permitted identifying the dimensionless well pressure for the infinitely long fracture as the master curve for type curves of all fractures with finite length from the beginning of bilinear flow up to fully developed radial flow. In log–log scale the master curve's logarithmic derivative initially follows a 1/4-slope straight line (characteristic for bilinear flow) and gradually bends into a horizontal line (characteristic for radial flow) for long times. During the bilinear flow period, isobars normalized to well pressure propagate with the fourth and second root of time in fracture and matrix, respectively. The width-to-length ratio of the pressure field increases proportional to the fourth root of time during the bilinear period, and starts to deviate from this relation close to the deviation of well pressure and its derivative from their fourth-root-of-time relations. At this time, isobars are already significantly inclined with respect to the fracture. The type curves of finite fractures all deviate counterclockwise from the master curve instead of clockwise or counterclockwise from the 1/4-slope straight line as previously proposed. The counterclockwise deviation from the master curve was identified as the arrival of a normalized isobar reflected at the fracture tip 16 times earlier. Nevertheless, two distinct regimes were found in regard to pressure at the fracture tip when bilinear flow ends. For dimensionless fracture conductivities TD 10, the pressure at the fracture tip has reached substantial fractions of the associated change in well pressure when the flow field transforms towards intermittent formation linear flow at times that scale inversely with the fourth power of dimensionless fracture conductivity. Our results suggest that semi-log plots of normalized well pressure provide a means for the determination of hydraulic parameters of fracture and matrix after shorter test duration than for conventional analysis.

Indications for different types of brittle failure due to active coal mining using waveform similarities of induced seismic events

Abstract: . Longwall mining activity in the Ruhr coal mining district leads to mining-induced seismicity. For detailed studies the seismicity of a single longwall panel beneath the town of Hamm-Herringen in the eastern Ruhr area was monitored between June 2006 and July 2007 with a dense temporary network of 15 seismic stations. More than 7000 seismic events with magnitudes between –1.7 ≤ ML ≤ 2.0 were detected and localized in this period. Most of the events occurred in the vicinity of the moving longwall face. In order to find possible differences in the brittle failure types of these events an association of the events to distinct clusters is performed based on their waveform characteristics. This task is carried out using a new clustering algorithm utilizing a network similarity matrix which is created by combining all available 3-component single station similarity matrices. The resultant network matrix is then sorted with respect to the similarity of its rows leading to a sorted matrix immediately indicating the clustering of the event catalogue. Finally, clusters of similar events are extracted by visual inspection. This approach results in the identification of several large clusters which are distinct with respect to their spatial and temporal characteristics as well as their frequency magnitude distributions. Comparable clusters are also found with a conventional single linkage approach, however, the new routine seems to be able to associate more events to specific clusters without merging the clusters. The nine largest observed clusters can be tentatively divided into three different groups that indicate different types of brittle failure. The first group consists of the two largest clusters which constitute more than half of all recorded events. Results of a relative relocation using cross-correlation data suggest that these events are confined to the extent of the mined out longwall and cluster close to the edges of the active longwall at the depth of active mining. These events occur in lockstep with the longwall advance and exhibit a high b value of the Gutenberg–Richter relation (GR) of about 1.5 to 2.5 and consist of small magnitude events. Thus, these events represent the immediate energy release adjacent to the mined out area. The second group consists of clusters located either slightly above or below the depth of active mining and occurring at the current position of the longwall face within the confines of the longwall. They consist of generally stronger events and do not follow GR. This activity might be linked to the failure of more competent layers above and below the mined out seam resulting in larger magnitude events. Finally, one cluster represents seismic activity with a rather low b value below 1 and events located partly towards the north of the longwall which are delayed with respect to the advance of the longwall face. These events are interpreted as brittle failure on pre-existing tectonic structures reactivated by the mining activity.

Regional wave propagation using the discontinuous Galerkin method

Abstract: . We present an application of the discontinuous Galerkin (DG) method to regional wave propagation. The method makes use of unstructured tetrahedral meshes, combined with a time integration scheme solving the arbitrary high-order derivative (ADER) Riemann problem. This ADER-DG method is high-order accurate in space and time, beneficial for reliable simulations of high-frequency wavefields over long propagation distances. Due to the ease with which tetrahedral grids can be adapted to complex geometries, undulating topography of the Earth's surface and interior interfaces can be readily implemented in the computational domain. The ADER-DG method is benchmarked for the accurate radiation of elastic waves excited by an explosive and a shear dislocation source. We compare real data measurements with synthetics of the 2009 L'Aquila event (central Italy). We take advantage of the geometrical flexibility of the approach to generate a European model composed of the 3-D EPcrust model, combined with the depth-dependent ak135 velocity model in the upper mantle. The results confirm the applicability of the ADER-DG method for regional scale earthquake simulations, which provides an alternative to existing methodologies.

Strain localisation in mechanically layered rocks beneath detachment zones: insights from numerical modelling

Abstract: . We have designed a series of fully dynamic numerical simulations aimed at assessing how the orientation of mechanical layering in rocks controls the orientation of shear bands and the depth of penetration of strain in the footwall of detachment zones. Two parametric studies are presented. In the first one, the influence of stratification orientation on the occurrence and mode of strain localisation is tested by varying initial dip of inherited layering in the footwall with regard to the orientation of simple shear applied at the rigid boundary simulating a rigid hanging wall, all scaling and rheological parameter kept constant. It appears that when Mohr–Coulomb plasticity is being used, shear bands are found to localise only when the layering is being stretched. This corresponds to early deformational stages for inital layering dipping in the same direction as the shear is applied, and to later stages for intial layering dipping towards the opposite direction of shear. In all the cases, localisation of the strain after only γ=1 requires plastic yielding to be activated in the strong layer. The second parametric study shows that results are length-scale independent and that orientation of shear bands is not sensitive to the viscosity contrast or the strain rate. However, decreasing or increasing strain rate is shown to reduce the capacity of the shear zone to localise strain. In the later case, the strain pattern resembles a mylonitic band but the rheology is shown to be effectively linear. Based on the results, a conceptual model for strain localisation under detachment faults is presented. In the early stages, strain localisation occurs at slow rates by viscous shear instabilities but as the layered media is exhumed, the temperature drops and the strong layers start yielding plastically, forming shear bands and localising strain at the top of the shear zone. Once strain localisation has occured, the deformation in the shear band becomes extremely penetrative but the strength cannot drop since the shear zone has a finite thickness.

Can vesicle size distributions assess eruption intensity during volcanic activity

Abstract: . We studied three-dimensional (3-D) vesicle size distributions by X-ray microtomography in scoria collected during the relatively quiescent Phase II of the April–May 2010 eruption at Eyjafjallajokull volcano, Iceland. Our goal was to compare cumulative vesicle size distributions (VSDs) measured in these samples with those found in Stromboli volcano, Italy. Stromboli was chosen because its VSDs are well-characterized and show a correlation with eruption intensity: typical Strombolian activity produces VSDs with power-law exponents near 1, whereas larger and more energetic vulcanian-type explosions and Plinian eruptions produce VSDs with power-law exponents near 1.5. The first hypothesis to be tested was whether or not the samples studied in this work would contain VSDs similar to normal Strombolian products, display higher power-law exponents, or be described by exponential functions. Before making this comparison, we tested a second hypothesis, which was that the magma–water interactions in the Eyjafjallajokull eruption might have a significant effect on the VSDs. We performed 1 bar bubble-growth experiments in which the samples were inundated with water and compared them to similar control experiments without water inundation. No significant differences between the VSDs of the two sets of experiments were found, and the second hypothesis is not supported by the experimental evidence. The Phase II Eyjafjallajokull VSDs are described by power-law exponents of ~0.8, typical of normal Strombolian eruptions, and support the first hypothesis. The comparable VSDs and behavior of Phase II of the Eyjafjallajokull 2010 eruption to Stromboli are interpreted to be a reflection of similar conduit systems in both volcanoes that are being constantly fed by the ascent of mingled/mixed magma from depth. Such behavior implies that continued activity during Phase II of the Eyjafjallajokull eruption could be expected and would have been predicted, had our VSDs been measured in real time during the eruption. However, the products studied show no peculiar feature that could herald the renewed eruption intensity observed in the following Phase III of the eruption.

Crustal 3-D geometry of the Kristineberg area (Sweden) with implications on VMS deposits

Abstract: . Structural analysis of the Palaeoproterozoic volcanogenic massive sulfide (VMS) hosting Kristineberg area, Sweden, constrained by existing magnetotelluric (MT) and seismic reflection data, reveals that the complex geometry characterized by non-cylindrical antiformal structures is due to transpression along the termination of a major high-strain zone. Similar orientations of the host rock deformation fabrics and the VMS ore lenses indicate that the present-day geometry of the complex VMS deposits in the Kristineberg area may be attributed to tectonic transposition. The tectonic transposition was dominantly controlled by reverse shearing and related upright to overturned folding, with increasing contribution of strike-slip shearing and sub-horizontal flow towards greater crustal depths. Furthermore, the northerly dip of the previously recognized subsurface crustal reflector within the Kristineberg area is attributed to formation of crustal compartments with opposite polarities within the scale of the whole Skellefte district. The resulting structural framework of the main geological units is visualized in a 3-D model which is available as a 3-D PDF document through the publication website.