ALBERT

Description: Extreme extension across Seram and Ambon, eastern Indonesia: evidence for Banda slab rollback Solid Earth, 4, 277-314, 2013 Author(s): J. M. Pownall, R. Hall, and I. M. Watkinson 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.

Description: Changes in soil organic carbon and nitrogen capacities of Salix cheilophila Schneid along a revegetation chronosequence in semi-arid degraded sandy land of the Gonghe Basin, Tibet Plateau Solid Earth, 5, 1045-1054, 2014 Author(s): Y. Yu and Z. Q. Jia The Gonghe Basin is a sandified and desertified region of China, but the distribution of soil organic carbon (SOC) and total nitrogen (TN) along the cultivation chronosequence across this ecologically fragile region is not well understood. This study was carried out to understand the effects of restoration with Salix cheilophila for different periods of time (6, 11, 16, 21 years) to test whether it enhanced C and N storage. Soil samples, in four replications from seven depth increments (0–10, 10–20, 20–30, 30–50, 50–100, 100–150 and 150–200 cm), were collected in each stand. Soil bulk density, SOC, TN, aboveground biomass and root biomass were measured. Results indicated that changes occurred in both the upper and deeper soil layers with an increase in revegetation time. The 0–200 cm soil showed that the 6-year stand gained 3.89 Mg C ha −1 and 1.00 Mg N ha −1 , which accounted for 40.82% of the original SOC and 11.06% of the TN of the 0-year stand. The 11-year stand gained 7.82 Mg C ha −1 and 1.98 Mg N ha −1 in the 0–200 cm soil layers, accounting for 58.06% of the SOC and 19.80% of the TN of the 0-year stand. The 16-year stand gained 11.32 Mg C ha −1 and 3.30 Mg N ha −1 in the 0–200 cm soil layers, accounting for 66.71% of the SOC and 21.98% of the TN of the 0-year stand. The 21-year stand gained 13.05 Mg C ha −1 and 5.45 Mg N ha −1 from the same soil depth, accounting for 69.79% of the SOC and 40.47% of the TN compared with the 0-year stand. The extent of these changes depended on soil depth and plantation age. The results demonstrated that, as stand age increased, the storage of SOC and TN increased. These results further indicated that restoration with S. cheilophila has positive impacts on the Gonghe Basin and has increased the capacity of SOC sequestration and N storage. The shrub's role as carbon sink is compatible with system management and persistence. The findings are significant for assessing C and N sequestration accurately in semi-arid degraded high, cold sandy regions in the future.

Description: Preface: Environmental benefits of biochar Solid Earth, 5, 1301-1303, 2014 Author(s): J. Paz-Ferreiro, A. Méndez, A. M. Tarquis, A. Cerdà, and G. Gascó

Description: Exploring the potentials and limitations of the time-reversal imaging of finite seismic sources Solid Earth, 2, 95-105, 2011 Author(s): S. Kremers, A. Fichtner, G. B. Brietzke, H. Igel, C. Larmat, L. Huang, and M. Käser The characterisation of seismic sources with time-reversed wave fields is developing into a standard technique that has already been successful in numerous applications. While the time-reversal imaging of effective point sources is now well-understood, little work has been done to extend this technique to the study of finite rupture processes. This is despite the pronounced non-uniqueness in classic finite source inversions. The need to better constrain the details of finite rupture processes motivates the series of synthetic and real-data time reversal experiments described in this paper. We address questions concerning the quality of focussing in the source area, the localisation of the fault plane, the estimation of the slip distribution and the source complexity up to which time-reversal imaging can be applied successfully. The frequency band for the synthetic experiments is chosen such that it is comparable to the band usually employed for finite source inversion. Contrary to our expectations, we find that time-reversal imaging is useful only for effective point sources, where it yields good estimates of both the source location and the origin time. In the case of finite sources, however, the time-reversed field does not provide meaningful characterisations of the fault location and the rupture process. This result cannot be improved sufficiently with the help of different imaging fields, realistic modifications of the receiver geometry or weights applied to the time-reversed sources. The reasons for this failure are manifold. They include the choice of the frequency band, the incomplete recording of wave field information at the surface, the excitation of large-amplitude surface waves that deteriorate the depth resolution, the absence of a sink that should absorb energy radiated during the later stages of the rupture process, the invisibility of small slip and the neglect of prior information concerning the fault geometry and the inherent smoothness of seismologically inferred Earth models that prevents the beneficial occurrence of strong multiple-scattering. The condensed conclusion of our study is that the limitations of time-reversal imaging – at least in the frequency band considered here – start where the seismic source stops being effectively point-localised.

Description: Native American lithic procurement along the international border in the boot heel region of southwestern New Mexico Solid Earth, 2, 75-93, 2011 Author(s): K. E. Zeigler, P. Hogan, C. Hughes, and A. Kurota Multidisciplinary field projects can be very useful to a more fundamental understanding of the world around us, though these projects are not as common as they should be. In particular, the combination of archeology and geology combines our understanding of human behavior and human use of the landscape with an intimate knowledge of geologic processes and the materials available for human use in order to gain a broader understanding of human-Earth interaction. Here we present data from a cross-disciplinary project that uses a common dataset, archeological artifacts, to explore the anthropological and geologic implications of useage patterns. Archeological excavations and surveys conducted by the Office of Contract Archeology in 2007 along the route of the proposed international border fence reveal patterns of use of geologic materials by Archaic, Formative and Protohistoric Native Americans in the Boot Heel of southwestern New Mexico. Thousands of artifacts were recorded in multiple sites from Guadalupe Pass in the southern Peloncillo Mountains to the Carrizalillo Hills west of Columbus. We identified the lithologies of artifacts, ranging from projectile points to groundstones, and then constructed material movement maps based on either known procurement sites ("quarries") or outcrops identified as the closest source to a given site for each lithology. Not unexpectedly, the majority of the rock types utilized by native peoples are local siliceous volcanic materials. However, several artifacts constructed from obsidian were transported into the region from northern Mexico and eastern Arizona, indicating long-distance travel and/or trade routes. We also examine useage pattern difference between Archaic, Formative and Protohistoric sites. Additionally, a dramatic change in distribution of sources for geologic materials occurs between one pre-Spanish site and one post-Spanish site that are adjacent to one another.

Description: A re-evaluation of the Italian historical geomagnetic catalogue: implications for paleomagnetic dating at active Italian volcanoes Solid Earth, 2, 65-74, 2011 Author(s): F. D'Ajello Caracciolo, A. Pignatelli, F. Speranza, and A. Meloni Paleomagnetism is proving to represent one of the most powerful dating tools of volcanics emplaced in Italy during the last few centuries/millennia. This method requires that valuable proxies of the local geomagnetic field (paleo)secular variation ((P)SV) are available. To this end, we re-evaluate the whole Italian geomagnetic directional dataset, consisting of 833 and 696 declination and inclination measurements, respectively, carried out since 1640 AD at several localities. All directions were relocated via the virtual geomagnetic pole method to Stromboli (38.8° N, 15.2° E), the rough centre of the active Italian volcanoes. For declination-only measurements, missing inclinations were derived (always by pole method) by French data (for period 1670–1789), and by nearby Italian sites/years (for periods 1640–1657 and 1790–1962). Using post-1825 declination values, we obtain a 0.46 ± 0.19° yr −1 westward drift of the geomagnetic field for Italy. The original observation years were modified, considering such drift value, to derive at a drift-corrected relocated dataset. Both datasets were found to be in substantial agreement with directions derived from the field models by Jackson et al. (2000) and Pavon-Carrasco et al. (2009). However, the drift-corrected dataset minimizes the differences between the Italian data and both field models, and eliminates a persistent 1.6° shift of 1933–1962 declination values from Castellaccio with respect to other nearly coeval Italian data. The relocated datasets were used to calculate two post-1640 Italian SV curves, with mean directions calculated every 30 and 10 years before and after 1790, respectively. The curve comparison suggests that both available field models yield the best available SV curve to perform paleomagnetic dating of 1600–1800 AD Italian volcanics, while the Italian drift-corrected curve is probably preferable for the 19th century. For the 20th century, the global model by Jackson et al. (2000) yields more accurate inclination values, while the declinations from our drift-corrected curve seem to better represent the local field evolution, at least for the first half of the century.

Description: Open Plot Project: an open-source toolkit for 3-D structural data analysis Solid Earth, 2, 53-63, 2011 Author(s): S. Tavani, P. Arbues, M. Snidero, N. Carrera, and J. A. Muñoz In this work we present the Open Plot Project, an open-source software for structural data analysis, including a 3-D environment. The software includes many classical functionalities of structural data analysis tools, like stereoplot, contouring, tensorial regression, scatterplots, histograms and transect analysis. In addition, efficient filtering tools are present allowing the selection of data according to their attributes, including spatial distribution and orientation. This first alpha release represents a stand-alone toolkit for structural data analysis. The presence of a 3-D environment with digitalising tools allows the integration of structural data with information extracted from georeferenced images to produce structurally validated dip domains. This, coupled with many import/export facilities, allows easy incorporation of structural analyses in workflows for 3-D geological modelling. Accordingly, Open Plot Project also candidates as a structural add-on for 3-D geological modelling software. The software (for both Windows and Linux O.S.), the User Manual, a set of example movies (complementary to the User Manual), and the source code are provided as Supplement. We intend the publication of the source code to set the foundation for free, public software that, hopefully, the structural geologists' community will use, modify, and implement. The creation of additional public controls/tools is strongly encouraged.

Description: Interpretative modelling of a geological cross section from boreholes: sources of uncertainty and their quantification Solid Earth, 5, 1189-1203, 2014 Author(s): R. M. Lark, S. Thorpe, H. Kessler, and S. J. Mathers We conducted a designed experiment to quantify sources of uncertainty in geologists' interpretations of a geological cross section. A group of 28 geologists participated in the experiment. Each interpreted borehole record included up to three Palaeogene bedrock units, including the target unit for the experiment: the London Clay. The set of boreholes was divided into batches from which validation boreholes had been withheld; as a result, we obtained 129 point comparisons between the interpreted elevation of the base of the London Clay and its observed elevation in a borehole not used for that particular interpretation. Analysis of the results showed good general agreement between the observed and interpreted elevations, with no evidence of systematic bias. Between-site variation of the interpretation error was spatially correlated, and the variance appeared to be stationary. The between-geologist component of variance was smaller overall, and depended on the distance to the nearest borehole. There was also evidence that the between-geologist variance depends on the degree of experience of the individual. We used the statistical model of interpretation error to compute confidence intervals for any one interpretation of the base of the London Clay on the cross section, and to provide uncertainty measures for decision support in a hypothetical route-planning process. The statistical model could also be used to quantify error propagation in a full 3-D geological model produced from interpreted cross sections.

Description: Tunable diode laser measurements of hydrothermal/volcanic CO 2 and implications for the global CO 2 budget Solid Earth, 5, 1209-1221, 2014 Author(s): M. Pedone, A. Aiuppa, G. Giudice, F. Grassa, V. Francofonte, B. Bergsson, and E. Ilyinskaya Quantifying the CO 2 flux sustained by low-temperature fumarolic fields in hydrothermal/volcanic environments has remained a challenge, to date. Here, we explored the potential of a commercial infrared tunable laser unit for quantifying such fumarolic volcanic/hydrothermal CO 2 fluxes. Our field tests were conducted between April 2013 and March 2014 at Nea Kameni (Santorini, Greece), Hekla and Krýsuvík (Iceland) and Vulcano (Aeolian Islands, Italy). At these sites, the tunable laser was used to measure the path-integrated CO 2 mixing ratios along cross sections of the fumaroles' atmospheric plumes. By using a tomographic post-processing routine, we then obtained, for each manifestation, the contour maps of CO 2 mixing ratios in the plumes and, from their integration, the CO 2 fluxes. The calculated CO 2 fluxes range from low (5.7 ± 0.9 t d −1 ; Krýsuvík) to moderate (524 ± 108 t d −1 ; La Fossa crater, Vulcano). Overall, we suggest that the cumulative CO 2 contribution from weakly degassing volcanoes in the hydrothermal stage of activity may be significant at the global scale.

Description: A Mesoproterozoic continental flood rhyolite province, the Gawler Ranges, Australia: the end member example of the Large Igneous Province clan Solid Earth, 2, 25-33, 2011 Author(s): M. J. Pankhurst, B. F. Schaefer, P. G. Betts, N. Phillips, and M. Hand Rhyolite and dacite lavas of the Mesoproterozoic upper Gawler Range Volcanics (GRV) (〉30 000 km 3 preserved), South Australia, represent the remnants of one of the most voluminous felsic magmatic events preserved on Earth. Geophysical interpretation suggests eruption from a central cluster of feeder vents which supplied large-scale lobate flows 〉100 km in length. Pigeonite inversion thermometers indicate eruption temperatures of 950–1100 °C. The lavas are A-type in composition (e.g. high Ga/Al ratios) and characterised by elevated primary halogen concentrations (~1600 ppm fluorine, ~400 ppm chlorine). These depolymerised the magma such that temperature-composition-volatile non-Arrhenian melt viscosity modelling suggests they had viscosities of

Description: Metamorphic history and geodynamic significance of the Early Cretaceous Sabzevar granulites (Sabzevar structural zone, NE Iran) Solid Earth, 2, 219-243, 2011 Author(s): M. Nasrabady, F. Rossetti, T. Theye, and G. Vignaroli The Iranian ophiolites are part of the vast orogenic suture zones that mark the Alpine-Himalayan convergence zone. Few petrological and geochronological data are available from these ophiolitic domains, hampering a full assessment of the timing and regimes of subduction zone metamorphism and orogenic construction in the region. This paper describes texture, geochemistry, and the pressure-temperature path of the Early Cretaceous mafic granulites that occur within the Tertiary Sabzevar ophiolitic suture zone of NE Iran. Whole rock geochemistry indicates that the Sabzevar granulites are likely derived from a MORB-type precursor. They are thus considered as remnants of a dismembered dynamo-thermal sole formed during subduction of a back-arc basin (proto-Sabzevar Ocean) formed in the upper-plate of the Neotethyan slab. The metamorphic history of the granulites suggests an anticlockwise pressure-temperature loop compatible with burial in a hot subduction zone, followed by cooling during exhumation. Transition from a nascent to a mature stage of oceanic subduction is the geodynamic scenario proposed to accomplish for the reconstructed thermobaric evolution. When framed with the regional scenario, results of this study point to diachronous and independent tectonic evolutions of the different ophiolitic domains of central Iran, for which a growing disparity in the timing of metamorphic equilibration and of pressure-temperature paths can be expected to emerge with further investigations.

Description: Remobilization of silicic intrusion by mafic magmas during the 2010 Eyjafjallajökull eruption Solid Earth, 2, 271-281, 2011 Author(s): O. Sigmarsson, I. Vlastelic, R. Andreasen, I. Bindeman, J.-L. Devidal, S. Moune, J. K. Keiding, G. Larsen, A. Höskuldsson, and Th. Thordarson Injection of basaltic magmas into silicic crustal holding chambers and subsequent magma mingling or mixing is a process that has been recognised since the late seventies as resulting in explosive eruptions. Detailed reconstruction and assessment of the mixing process caused by such intrusion is now possible because of the exceptional time-sequence sample suite available from the tephra fallout of the 2010 summit eruption at Eyjafjallajökull volcano in South Iceland. Fallout from 14 to 19 April contains three glass types of basaltic, intermediate, and silicic compositions recording rapid magma mingling without homogenisation, involving evolved FeTi-basalt and silicic melt with composition identical to that produced by the 1821–1823 AD Eyjafjallajökull summit eruption. The time-dependent change in the magma composition suggests a binary mixing process with changing end-member compositions and proportions. Beginning of May, a new injection of primitive basalt was recorded by deep seismicity, appearance of Mg-rich olivine phenocrysts together with high sulphur dioxide output and presence of sulphide crystals. Thus, the composition of the basaltic injection became more magnesian and hotter with time provoking changes in the silicic mixing end-member from pre-existing melt to the solid carapace of the magma chamber. Finally, decreasing proportions of the mafic end-member with time in the erupted mixed-magma demonstrate that injections of Mg-rich basalt was the motor of the 2010 Eyjafjallajökull explosive eruption, and that its decreasing inflow terminated the eruption. Significant quantity of silicic magma is thus still present in the interior of the volcano. Our results show that detailed sampling during the entire eruption was essential for deciphering the complex magmatic processes at play, i.e. the dynamics of the magma mingling and mixing. Finally, the rapid compositional changes in the eruptive products suggest that magma mingling occurs on a timescale of a few hours to days whereas the interval between the first detected magma injection and eruption was several months.

Description: Pore formation during dehydration of a polycrystalline gypsum sample observed and quantified in a time-series synchrotron X-ray micro-tomography experiment Solid Earth, 3, 71-86, 2012 Author(s): F. Fusseis, C. Schrank, J. Liu, A. Karrech, S. Llana-Fúnez, X. Xiao, and K. Regenauer-Lieb We conducted an in-situ X-ray micro-computed tomography heating experiment at the Advanced Photon Source (USA) to dehydrate an unconfined 2.3 mm diameter cylinder of Volterra Gypsum. We used a purpose-built X-ray transparent furnace to heat the sample to 388 K for a total of 310 min to acquire a three-dimensional time-series tomography dataset comprising nine time steps. The voxel size of 2.2 μm 3 proved sufficient to pinpoint reaction initiation and the organization of drainage architecture in space and time. We observed that dehydration commences across a narrow front, which propagates from the margins to the centre of the sample in more than four hours. The advance of this front can be fitted with a square-root function, implying that the initiation of the reaction in the sample can be described as a diffusion process. Novel parallelized computer codes allow quantifying the geometry of the porosity and the drainage architecture from the very large tomographic datasets (2048 3 voxels) in unprecedented detail. We determined position, volume, shape and orientation of each resolvable pore and tracked these properties over the duration of the experiment. We found that the pore-size distribution follows a power law. Pores tend to be anisotropic but rarely crack-shaped and have a preferred orientation, likely controlled by a pre-existing fabric in the sample. With on-going dehydration, pores coalesce into a single interconnected pore cluster that is connected to the surface of the sample cylinder and provides an effective drainage pathway. Our observations can be summarized in a model in which gypsum is stabilized by thermal expansion stresses and locally increased pore fluid pressures until the dehydration front approaches to within about 100 μm. Then, the internal stresses are released and dehydration happens efficiently, resulting in new pore space. Pressure release, the production of pores and the advance of the front are coupled in a feedback loop.

Description: The regulation of the air: a hypothesis Solid Earth, 3, 87-96, 2012 Author(s): E. G. Nisbet, C. M. R. Fowler, and R. E. R. Nisbet We propose the hypothesis that natural selection, acting on the specificity or preference for CO 2 over O 2 of the enzyme rubisco (ribulose-1,5-bisphosphate carboxylase/oxygenase), has controlled the CO 2 :O 2 ratio of the atmosphere since the evolution of photosynthesis and has also sustained the Earth's greenhouse-set surface temperature. Rubisco works in partnership with the nitrogen-fixing enzyme nitrogenase to control atmospheric pressure. Together, these two enzymes control global surface temperature and indirectly the pH and oxygenation of the ocean. Thus, the co-evolution of these two enzymes may have produced clement conditions on the Earth's surface, allowing life to be sustained.

Description: Constraining fault interpretation through tomographic velocity gradients: application to northern Cascadia Solid Earth, 3, 53-61, 2012 Author(s): K. Ramachandran Spatial gradients of tomographic velocities are seldom used in interpretation of subsurface fault structures. This study shows that spatial velocity gradients can be used effectively in identifying subsurface discontinuities in the horizontal and vertical directions. Three-dimensional velocity models constructed through tomographic inversion of active source and/or earthquake traveltime data are generally built from an initial 1-D velocity model that varies only with depth. Regularized tomographic inversion algorithms impose constraints on the roughness of the model that help to stabilize the inversion process. Final velocity models obtained from regularized tomographic inversions have smooth three-dimensional structures that are required by the data. Final velocity models are usually analyzed and interpreted either as a perturbation velocity model or as an absolute velocity model. Compared to perturbation velocity model, absolute velocity models have an advantage of providing constraints on lithology. Both velocity models lack the ability to provide sharp constraints on subsurface faults. An interpretational approach utilizing spatial velocity gradients applied to northern Cascadia shows that subsurface faults that are not clearly interpretable from velocity model plots can be identified by sharp contrasts in velocity gradient plots. This interpretation resulted in inferring the locations of the Tacoma, Seattle, Southern Whidbey Island, and Darrington Devil's Mountain faults much more clearly. The Coast Range Boundary fault, previously hypothesized on the basis of sedimentological and tectonic observations, is inferred clearly from the gradient plots. Many of the fault locations imaged from gradient data correlate with earthquake hypocenters, indicating their seismogenic nature.

Description: Triplicated P-wave measurements for waveform tomography of the mantle transition zone Solid Earth, 3, 339-354, 2012 Author(s): S. C. Stähler, K. Sigloch, and T. Nissen-Meyer Triplicated body waves sample the mantle transition zone more extensively than any other wave type, and interact strongly with the discontinuities at 410 km and 660 km. Since the seismograms bear a strong imprint of these geodynamically interesting features, it is highly desirable to invert them for structure of the transition zone. This has rarely been attempted, due to a mismatch between the complex and band-limited data and the (ray-theoretical) modelling methods. Here we present a data processing and modelling strategy to harness such broadband seismograms for finite-frequency tomography. We include triplicated P-waves (epicentral distance range between 14 and 30°) across their entire broadband frequency range, for both deep and shallow sources. We show that is it possible to predict the complex sequence of arrivals in these seismograms, but only after a careful effort to estimate source time functions and other source parameters from data, variables that strongly influence the waveforms. Modelled and observed waveforms then yield decent cross-correlation fits, from which we measure finite-frequency traveltime anomalies. We discuss two such data sets, for North America and Europe, and conclude that their signal quality and azimuthal coverage should be adequate for tomographic inversion. In order to compute sensitivity kernels at the pertinent high body wave frequencies, we use fully numerical forward modelling of the seismic wavefield through a spherically symmetric Earth.

Description: Heat-flow and subsurface temperature history at the site of Saraya (eastern Senegal) Solid Earth, 3, 213-224, 2012 Author(s): F. Lucazeau and F. Rolandone New temperature measurements from eight boreholes in the West African Craton (WAC) reveal superficial perturbations down to 100 m below the alteration zone. These perturbations are both related to a recent increase in the surface air temperature (SAT) and to the site effects caused by fluid circulations and/or the lower conduction in the alterites. The ground surface temperature (GST), inverted from the boreholes temperatures, increased slowly in the past (~0.4 °C from 1700 to 1940) and then, more importantly, in recent years (~1.5 °C from 1940 to 2010). This recent trend is consistent with the increase of the SAT recorded at two nearby meteorological stations (Tambacounda and Kedougou), and more generally in the Sahel with a coeval rainfall decrease. Site effects are superimposed to the climatic effect and interpreted by advective (circulation of fluids) or conductive (lower conductivity of laterite and of high-porosity sand) perturbations. We used a 1-D finite differences thermal model and a Monte-Carlo procedure to find the best estimates of these site perturbations: all the eight borehole temperature logs can be interpreted with the same basal heat-flow and the same surface temperature history, but with some realistic changes of thermal conductivity and/or fluid velocity. The GST trend observed in Senegal can be confirmed by two previous borehole measurements made in 1983 in other locations of West Africa, the first one in an arid zone of northern Mali and the second one in a sub-humid zone in southern Mali. Finally, the background heat-flow is low (31±2 mW m −2 ), which makes this part of the WAC more similar with the observations in the southern part (33±8 mW m −2 ) rather than with those in the northern part and in the Pan-African domains where the surface heat-flow is 15–20 mW m −2 higher.

Description: Seismic imaging of sandbox experiments – laboratory hardware setup and first reflection seismic sections Solid Earth, 4, 93-104, 2013 Author(s): C. M. Krawczyk, M.-L. Buddensiek, O. Oncken, and N. Kukowski 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.

Description: Exhumation of (ultra-)high-pressure terranes: concepts and mechanisms Solid Earth, 4, 75-92, 2013 Author(s): C. J. Warren 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.

Description: Segmentation of the Izu-Bonin and Mariana slabs based on the analysis of the Benioff seismicity distribution and regional tomography results Solid Earth, 4, 59-73, 2013 Author(s): K. Jaxybulatov, I. Koulakov, and N. L. Dobretsov 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.

Description: Plate kinematics in the Cantabrian domain of the Pyrenean orogen Solid Earth, 3, 265-292, 2012 Author(s): S. Tavani The Cantabrian domain represents the western portion of the Pyrenean orogen, in the area where the Iberian continental lithosphere was subducted toward the north underneath the transitional to oceanic lithosphere of the Bay of Biscay. There, the about 100 km of orogenic convergence have been mostly accommodated in the northern portion of the orogen (i.e. the retro wedge) developed in the Bay of Biscay abyssal plain, while only crustal-scale folding with limited internal deformation occurred in the Cantabrian southern wedge (pro-wedge). Integrated meso- and macrostructural analyses and a reappraisal of available information from the transitional area between the Pyrenean and Cantabrian domains are presented in this work, allowing to set geometric and kinematic constraints on the entire Meso-Cenozoic history of the northern portion of the Iberian Plate, including subduction initiation and evolution in the western portion of the Pyrenean orogen. The structural record of the Late Jurassic to Early Cretaceous deformation stage, which was associated with rifting and seafloor spreading in the Bay of Biscay, indicates a ridge perpendicular (NNE-SSW oriented) extension, with no evidence of relevant strike-slip components during rifting. A Cenozoic NNW-SSE oriented shortening stage followed, related to the limited (about 100 km) north-directed subduction of the Iberian continental lithosphere underneath the transitional to oceanic lithosphere of the Bay of Biscay. Subduction led to the formation of the poorly-developed Cantabrian pro-wedge, which is laterally juxtaposed to the well-developed Pyrenean pro-wedge to the east. During this convergence stage, the structural framework in the Cantabrian pro-wedge, and particularly along its transition with the Pyrenean wedge to the east, was severely complicated by the reactivation of Paleozoic and Mesozoic inherited structures. Data presented in this work fully support the development of the Cantabrian Mountains as related to indentation and consequent thickening of the Bay of Biscay transitional lower crust during north-directed subduction of Iberian continental lithosphere. In essence, the Cantabrian pro-wedge is a lithospheric south-verging fault-propagation anticline developing above the subduction plane. The structural record in the area indicates that a lithospheric fault-propagation folding stage was predated, during the very early stages of orogenic shortening, by the development of a lithospheric-scale open syncline overlying the nucleation point of lithosphere sinking. Such a syncline is today partially preserved and represents one of the few natural examples of subduction initiation.

Description: Seismicity from February 2006 to September 2007 at the Rwenzori Mountains, East African Rift: earthquake distribution, magnitudes and source mechanisms Solid Earth, 3, 251-264, 2012 Author(s): M. Lindenfeld, G. Rümpker, A. Batte, and A. Schumann We have analysed the microseismic activity within the Rwenzori Mountains area in the western branch of the East African Rift. Seismogram recordings from a temporary array of up to 27 stations reveal approximately 800 events per month with local magnitudes ranging from –0.5 to 5.1. The earthquake distribution is highly heterogeneous. The majority of located events lie within faults zones to the east and west of the Rwenzoris with the highest seismic activity observed in the northeastern area, where the mountains are in contact with the rift shoulders. The hypocentral depth distribution exhibits a pronounced peak of seismic energy release at 15 km depth. The maximum extent of seismicity ranges from 20 to 32 km and correlates well with Moho depths that were derived from teleseismic receiver functions. We observe two general features: (i) beneath the rift shoulders, seismicity extends from the surface down to ca. 30 km depth; (ii) beneath the rift valley, seismicity is confined to depths greater than 10 km. From the observations there is no indication for a crustal root beneath the Rwenzori Mountains. The magnitude frequency distribution reveals a b -value of 1.1, which is consistent with the hypothesis that part of the seismicity is caused by magmatic processes within the crust. Fault plane solutions of 304 events were derived from P-polarities and SV/P amplitude ratios. More than 70 % of the source mechanisms exhibit pure or predominantly normal faulting. T-axis trends are highly uniform and oriented WNW–ESE, which is perpendicular to the rift axis and in good agreement with kinematic rift models. At the northernmost part of the region we observe a rotation of the T-axis trends to NEN–SWS, which may be indicative of a local perturbation of the regional stress field.

Description: Shallow water carbonate platforms (Late Aptian–Early Albian, Southern Apennines) in the context of supraregional to global changes: re-appraisal of palaeoecological events as reflectors of carbonate factory response Solid Earth, 3, 225-249, 2012 Author(s): A. Raspini This paper discusses the palaeoenvironmental significance of the "Orbitolina Level", the microbial carbonates and the Salpingoporella dinarica -rich deposits encased in the Aptian/Albian shallow water carbonate platform strata of Monte Tobenna and Monte Faito (Southern Italy). These facies show a peculiar field appearance due to their color and/or fossil content. In the shallow water carbonate strata, the Late Aptian "Orbitolina Level" was formed during a period of decreasing accommodation space. Microbial carbonates occur in different levels in the composite section. They reach their maximum thickness around the sequence boundaries just above the "Orbitolina Level" and close to the Aptian–Albian transition, and were not deposited during maximum flooding. S. dinarica -rich deposits occur in the lower part of the Monte Tobenna-Monte Faito composite section, in both restricted and more open lagoonal sediments. S. dinarica has its maximum abundance below the "Orbitolina Level" and disappears 11 m above this layer. On the basis of δ 13 C and δ 18 O values recorded at Tobenna-Faito, the succession has been correlated to global sea-level changes and to the main volcanic and climatic events during the Aptian. Deterioration of the inner lagoon environmental conditions was related to high trophic levels triggered by volcano-tectonic activity. Microbial carbonates were deposited especially in periods of third-order sea level lowering. In such a scenario, periods of increased precipitation during the Gargasian induced the mobilization of clay during flooding of the exposed platform due to high-frequency sea-level changes, with consequent terrigenous input to the lagoon. This and the high nutrient levels made the conditions unsuitable for the principle carbonate producers, and an opportunistic biota rich in orbitolinids ( Mesorbitolina texana and M. parva ) populated the platform. In the more open marine domain, the increased nutrient input enhanced the production of organic matter and locally led to the formation of black shales (e.g. the Niveau Fallot in the Vocontian Basin). It is argued that the concomitant low Mg/Ca molar ratio and high concentration of calcium in seawater could have favoured the development of the low-Mg calcite skeleton of the S. dinarica green algae. During third-order sea-level rise, no or minor microbial carbonates formed in the shallowlagoonal settings and S. dinarica disappeared. Carbonate neritic ecosystems were not influenced by the environmental changes inferred to have been induced by the mid-Cretaceous volcanism. The "Orbitolina Level", the microbial carbonates and the Salpingoporella dinarica -rich deposits in the studied Aptian/Albian shallow water carbonate strata are interpreted to be the response to environmental and oceanographic changes in shallow-water and deeper-marine ecosystems.

Description: High resolution reflection seismic profiling over the Tjellefonna fault in the Møre-Trøndelag Fault Complex, Norway Solid Earth, 3, 175-188, 2012 Author(s): E. Lundberg, C. Juhlin, and A. Nasuti The Møre-Trøndelag Fault Complex (MTFC) is one of the most prominent fault zones of Norway, both onshore and offshore. In spite of its importance, very little is known of the deeper structure of the individual fault segments comprising the fault complex. Most seismic lines have been recorded offshore or focused on deeper structures. This paper presents results from two reflection seismic profiles, located on each side of the Tingvollfjord, acquired over the Tjellefonna fault in the southeastern part of the MTFC. Possible kilometer scale vertical offsets, reflecting large scale northwest-dipping normal faulting, separating the high topography to the southeast from lower topography to the northwest have been proposed for the Tjellefonna fault or the Baeverdalen lineament. In this study, however, the Tjellefonna fault is interpreted to dip approximately 50–60° towards the southeast to depths of at least 1.3 km. Travel-time modeling of reflections associated with the fault was used to establish the geometry of the fault structure at depth, while detailed analysis of first P-wave arrivals in shot gathers, together with resistivity profiles, were used to define the near surface geometry of the fault zone. A continuation of the structure on the northeastern side of the Tingvollfjord is suggested by correlation of an in strike direction P-S converted reflection (generated by a fracture zone) seen on the reflection data from that side of the Tingvollfjord. The reflection seismic data correlate well with resistivity profiles and recently published near surface geophysical data. A highly reflective package forming a gentle antiform structure was also identified on both seismic profiles. This structure could be related to the folded amphibolite lenses seen on the surface or possibly by an important boundary within the gneissic basement rocks of the Western Gneiss Region. The fold hinge line of the structure is parallel with the Tjellefonna fault trace suggesting that the folding and faulting may have been related.

Description: Corrigendum to "Floating stones off El Hierro, Canary Islands: xenoliths of pre-island sedimentary origin in the early products of the October 2011 eruption" published in Solid Earth, 3, 97–110, 2012 Solid Earth, 3, 189-189, 2012 Author(s): V. R. Troll, A. Klügel, M.-A. Longpré, S. Burchardt, F. M. Deegan, J. C. Carracedo, S. Wiesmaier, U. Kueppers, B. Dahren, L. S. Blythe, T. H. Hansteen, C. Freda, D. A. Budd, E. M. Jolis, E. Jonsson, F. C. Meade, C. Harris, S. E. Berg, L. Mancini, M. Polacci, and K. Pedroza No abstract available.

Description: New developments in the analysis of column-collapse pyroclastic density currents through numerical simulations of multiphase flows Solid Earth, 3, 161-173, 2012 Author(s): S. Lepore and C. Scarpati A granular multiphase model has been used to evaluate the action of differently sized particles on the dynamics of fountains and associated pyroclastic density currents. The model takes into account the overall disequilibrium conditions between a gas phase and several solid phases, each characterized by its own physical properties. The dynamics of the granular flows (fountains and pyroclastic density currents) has been simulated by adopting a Reynolds-averaged Navier-Stokes model for describing the turbulence effects. Numerical simulations have been carried out by using different values for the eruptive column temperature at the vent, solid particle frictional concentration, turbulent kinetic energy, and dissipation. The results obtained provide evidence of the multiphase nature of the model and describe several disequilibrium effects. The low concentration (≤5 × 10 −4 ) zones lie in the upper part of the granular flow, above the fountain, and above the tail and body of pyroclastic density current as thermal plumes. The high concentration zones, on the contrary, lie in the fountain and at the base of the current. Hence, pyroclastic density currents are assimilated to granular flows constituted by a low concentration suspension flowing above a high concentration basal layer (boundary layer), from the proximal regions to the distal ones. Interactions among the solid particles in the boundary layer of the granular flow are controlled by collisions between particles, whereas the dispersal of particles in the suspension is determined by the dragging of the gas phase. The simulations describe well the dynamics of a tractive boundary layer leading to the formation of stratified facies during Strombolian to Plinian eruptions.

Description: The fate of fluids released from subducting slab in northern Cascadia Solid Earth, 3, 121-129, 2012 Author(s): K. Ramachandran and R. D. Hyndman Large amounts of water carried down in subduction zones are driven upward into the overlying forearc upper mantle and crust as increasing temperatures and pressure dehydrate the subducting crust. Through seismic tomography velocities we show (a) the overlying forearc mantle in northern Cascadia is hydrated to serpentinite, and (b) there is low Poisson's ratio at the base of the forearc lower crust that may represent silica deposited from the rising fluids. From the velocities observed in the forearc mantle, the volume of serpentinite estimated is ∼30 %. This mechanically weak hydrated forearc region has important consequences in limits to great earthquakes and to collision tectonics. An approximately 10 km thick lower crustal layer of low Poisson's ratio (σ = 0.22) in the forearc is estimated to represent a maximum addition of ∼14 % by volume of quartz (σ = 0.09). If this quartz is removed from rising silica-saturated fluids over long times, it represents a significant addition of silica to the continental crust and an important contributor to its average composition.

Description: Geomagnetic jerks characterization via spectral analysis Solid Earth, 3, 131-148, 2012 Author(s): B. Duka, A. De Santis, M. Mandea, A. Isac, and E. Qamili In this study we have applied spectral techniques to analyze geomagnetic field time-series provided by observatories, and compared the results with those obtained from analogous analyses of synthetic data estimated from models. Then, an algorithm is here proposed to detect the geomagnetic jerks in time-series, mainly occurring in the eastern component of the geomagnetic field. Applying such analysis to time-series generated from global models has allowed us to depict the most important space-time features of the geomagnetic jerks all over the globe, since the beginning of XXth century. Finally, the spherical harmonic power spectrum of the third derivative of the main geomagnetic field has been computed from 1960 to 2002.5, bringing new insights to understand the spatial evolution of these rapid changes of the geomagnetic field.

Description: Paleointensities on 8 ka obsidian from Mayor Island, New Zealand Solid Earth, 2, 259-270, 2011 Author(s): A. Ferk, R. Leonhardt, K.-U. Hess, and D. B. Dingwell The 8 ka BP (6050 BCE) pantelleritic obsidian flow on Mayor Island, Bay of Plenty, New Zealand, has been investigated using 30 samples from two sites. Due to a very high paramagnetic/ferromagnetic ratio, it was not possible to determine the remanence carriers. This is despite the fact that the samples were studied intensively at low, room, and high temperatures. We infer that a stable remanence within the samples is carried by single- or close to single-domain particles. Experiments to determine the anisotropy of thermoremanence tensor and the dependency on cooling rate were hampered due to alteration resulting from the repeated heating of the samples to temperatures just below the glass transition. Nonetheless, a well-defined mean paleointensity of 57.0 ± 1.0 μT, based on individual high quality paleointensity determinations, was obtained. This field value compares very well to a paleointensity of 58.1 ± 2.9 μT, which Tanaka et al. (2009) obtained for 5500 BCE at a site 100 km distant. Agreement with geomagnetic field models, however, is poor. Thus, gathering more high-quality paleointensity data for the Pacific region and for the southern hemisphere in general to better constrain global field models is very important.

Description: An open ocean record of the Toarcian oceanic anoxic event Solid Earth, 2, 245-257, 2011 Author(s): D. R. Gröcke, R. S. Hori, J. Trabucho-Alexandre, D. B. Kemp, and L. Schwark Oceanic anoxic events were time intervals in the Mesozoic characterized by widespread distribution of marine organic matter-rich sediments (black shales) and significant perturbations in the global carbon cycle. These perturbations are globally recorded in sediments as carbon isotope excursions irrespective of lithology and depositional environment. During the early Toarcian, black shales were deposited on the epi- and pericontinental shelves of Pangaea, and these sedimentary rocks are associated with a pronounced (ca. 7 ‰) negative (organic) carbon isotope excursion (CIE) which is thought to be the result of a major perturbation in the global carbon cycle. For this reason, the lower Toarcian is thought to represent an oceanic anoxic event (the T-OAE). If the T-OAE was indeed a global event, an isotopic expression of this event should be found beyond the epi- and pericontinental Pangaean localities. To address this issue, the carbon isotope composition of organic matter (δ 13 C org of lower Toarcian organic matter-rich cherts from Japan, deposited in the open Panthalassa Ocean, was analysed. The results show the presence of a major (〉6 ‰) negative excursion in δ 13 C org that, based on radiolarian biostratigraphy, is a correlative of the lower Toarcian negative CIE known from Pangaean epi- and pericontinental strata. A smaller negative excursion in δ 13 C org (ca. 2 ‰) is recognized lower in the studied succession. This excursion may, within the current biostratigraphic resolution, represent the excursion recorded in European epicontinental successions close to the Pliensbachian/Toarcian boundary. These results from the open ocean realm suggest, in conjunction with other previously published datasets, that these Early Jurassic carbon cycle perturbations affected the active global reservoirs of the exchangeable carbon cycle (deep marine, shallow marine, atmospheric).

Description: Effect of glacial-interglacial sea-level changes on the displacement and stress field in the forearc and along the plate interface of subduction zones Solid Earth, 3, 63-70, 2012 Author(s): T. Li and A. Hampel Combined seismological, space-geodetic and numerical studies have shown that the seismicity at subduction zones may be modulated by tides and glacier fluctuations on timescales of 1–100 a, because these changes in loads on Earth's surface are able to alter the stress field in the upper plate and along the plate interface. Here we use a two-dimensional finite-element model of a subduction zone to investigate how glacial-interglacial sea-level changes affect the forearc region and the plate interface. The model results show that a sea-level fall by 125 m over 100 ka causes up to 0.7 m of vertical displacement, with the maximum uplift occurring between the trench and the coast. The uplift signal induced by the sea-level fall decreases to zero ~20 km landward of the coastline. A subsequent sea-level rise by 125 m over 20 ka causes subsidence, which is again most pronounced in the submarine part of the forearc. The sea-level changes cause horizontal displacements of up to 0.12 m, which are directed seaward during sea-level fall and landward during sea-level rise. With respect to the stress field, the sea-level changes lead to variations in the vertical stress and the shear stress of up to 1.23 MPa and 0.4 MPa, respectively. The shear stress variations are highest beneath the coast, i.e. in the area where the sea-level changes cause the strongest flexure. The resulting Coulomb stress changes on the plate interface are of the order of 0.2–0.5 MPa and indicate that earthquakes are promoted during sea-level fall and delayed during sea-level rise. Our findings imply that eustatic sea-level changes during glacial-interglacial periods may have induced displacements and stress changes that were large enough to affect the seismic cycle of subduction thrusts.

Description: Geophysical characterisation of two segments of the Møre-Trøndelag Fault Complex, Mid Norway Solid Earth, 2, 125-134, 2011 Author(s): A. Nasuti, C. Pascal, J. Ebbing, and J. F. Tønnesen The Møre-Trøndelag Fault Complex (MTFC) has controlled the tectonic evolution of Mid Norway and its shelf for the past 400 Myr through repeated reactivations during Palaeozoic, Mesozoic and perhaps Cenozoic times, the very last phase of reactivation involving normal to oblique-slip faulting. Despite its pronounced signature in the landscape, its deep structure has largely remained unresolved until now. We focused on two specific segments of the MTFC (i.e. the Tjellefonna and Bæverdalen faults) and acquired multiple geophysical datasets (i.e. gravity, magnetic, resistivity and shallow refraction profiles). A 100–200 m-wide zone of gouge and/or brecciated bedrock steeply dipping to the south is interpreted as being the Tjellefonna fault sensu stricto. The fault appears to be flanked by two additional but minor damage zones. A secondary normal fault also steeply dipping to the south but involving indurated breccias was detected ~1 km farther north. The Bæverdalen fault, ~12 km farther north, is interpreted as a ~700 m-wide and highly deformed zone involving fault gouge, breccias and lenses of intact bedrock. As such, it is probably the most important fault segment in the studied area and accommodated most of the strain during presumably Late Jurassic normal faulting. Our geophysical data are indicative of a Bæverdalen fault dipping steeply towards the south, in agreement with the average orientation of the local tectonic grain. Our findings suggest that the influence of Mesozoic normal faulting along the MTFC on landscape development is more complex than previously thought.

Description: New zircon data supporting models of short-lived igneous activity at 1.89 Ga in the western Skellefte District, central Fennoscandian Shield Solid Earth, 2, 205-217, 2011 Author(s): P. Skyttä, T. Hermansson, J. Andersson, M. Whitehouse, and P. Weihed New U-Th-Pb zircon data (SIMS) from three intrusive phases of the Palaeoproterozoic Viterliden intrusion in the western Skellefte District, central Fennoscandian Shield, dates igneous emplacement in a narrow time interval at about 1.89 Ga. A locally occurring quartz-plagioclase porphyritic tonalite, here dated at 1889 ± 3 Ma, is considered the youngest of the intrusive units, based on the new age data and field evidence. This supports an existing interpretation of its fault-controlled emplacement after intrusion of the dominating hornblende-tonalite units, in this study dated at 1892 ± 3 Ma. The Viterliden magmatism was synchronous with the oldest units of the Jörn type early-orogenic intrusions in the eastern part of the district (1.89–1.88 Ga; cf. Gonzàles Roldán, 2010). A U-Pb zircon age for a felsic metavolcanic rock from the hanging-wall to the Kristineberg VMS deposit, immediately south of the Viterliden intrusion, is constrained at 1883 ± 6 Ma in this study. It provides a minimum age for the Kristineberg ore deposit and suggests contemporaneous igneous/volcanic activity throughout the Skellefte District. Furthermore, it supports the view that the Skellefte Group defines a laterally continuous belt throughout this "ore district". Tentative correlation of the 1889 ± 3 Ma quartz-plagioclase porphyritic tonalite with the Kristineberg "mine porphyry" suggests that these units are coeval at about 1.89 Ga. Based on the new age determinations, the Viterliden intrusion may equally well have intruded into or locally acted as a basement for the ore-hosting Skellefte Group volcanic rocks.

Description: The model of self-generated seismo-electromagnetic oscillations of the LAI system Solid Earth, 2, 17-23, 2011 Author(s): M. K. Kachakhidze, Z. A. Kereselidze, and N. K. Kachakhidze Very low frequency (VLF) electromagnetic radiation (in diapason 1 kHz–1 MHz) in the atmosphere, generated during an earthquake preparation period, may be connected with the linear size characterising the expected earthquake focus. In order to argue this hypothesis, a very simple quasi-electrostatic model is used: the local VLF radiation may represent the self-generated (own) electromagnetic oscillations of interactive seismoactive segments of the lithosphere-atmosphere system. This model qualitatively explains the well-known precursor effects of earthquakes. In addition, using this model after diagnosing existing data makes it principally possible to forecast an expected earthquake with certain precision. As a physical basis of the working hypothesis is the atmospheric effect of polarization charges occurring in the surface layer of the Earth, it is possible to test the following constructed model in the Earth's crust, where the reason for polarization charge generation may be different from piezo-electric mechanism, e.g., some other mechanism.

Description: Optical method for measuring bed topography and flow depth in an experimental flume Solid Earth, 2, 143-154, 2011 Author(s): A. Limare, M. Tal, M. D. Reitz, E. Lajeunesse, and F. Métivier We describe an optical method known as moiré for acquiring quasi-simultaneous measurements of bed topography and flow depth in laboratory experiments. The moiré method is based on projecting a fringe pattern (grating) on the bed and analyzing the deformation of the pattern caused by the topography with respect to a reference plane. The height of the object is encoded in the phase of the pattern and can be retrieved either through Fourier transform or phase shifting algorithms. The methodology enables image-based non-contact measurements over a continuous surface at very high spatial and temporal resolutions. We use a commercial software package of a moiré method called Light3D to map bed topography and flow depth in an experimental braided channel and demonstrate how the method can be used to characterize a full range of statistics not previously possible.

Description: First observational evidence for the CO 2 -driven origin of Stromboli's major explosions Solid Earth, 2, 135-142, 2011 Author(s): A. Aiuppa, M. Burton, P. Allard, T. Caltabiano, G. Giudice, S. Gurrieri, M. Liuzzo, and G. Salerno We report on the first detection of CO 2 flux precursors of the till now unforecastable "major" explosions that intermittently occur at Strombolivolcano (Italy). An automated survey of the crater plume emissions in the period 2006–2010, during which 12 such explosions happened, demonstrated that these events are systematically preceded by a brief phase of increasing CO 2 /SO 2 weight ratio (up to 〉40) and CO 2 flux (〉1300 t d −1 ) with respect to the time-averaged values of 3.7 and ~500 t d −1 typical for standard Stromboli's activity. These signals are best explained by the accumulation of CO 2 -rich gas at a discontinuity of the plumbing system (decreasing CO 2 emission at the surface), followed by increasing gas leakage prior to the explosion. Our observations thus supports the recent model of Allard (2010) for a CO 2 -rich gas trigger of recurrent major explosions at Stromboli, and demonstrates the possibility to forecast these events in advance from geochemical precursors. These observations and conclusions have clear implications for monitoring strategies at other open-vent basaltic volcanoes worldwide.

Description: 3-D thermo-mechanical laboratory modeling of plate-tectonics: modeling scheme, technique and first experiments Solid Earth, 2, 35-51, 2011 Author(s): D. Boutelier and O. Oncken We present an experimental apparatus for 3-D thermo-mechanical analogue modeling of plate tectonic processes such as oceanic and continental subductions, arc-continent or continental collisions. The model lithosphere, made of temperature-sensitive elasto-plastic analogue materials with strain softening, is submitted to a constant temperature gradient causing a strength reduction with depth in each layer. The surface temperature is imposed using infrared emitters, which allows maintaining an unobstructed view of the model surface and the use of a high resolution optical strain monitoring technique (Particle Imaging Velocimetry). Subduction experiments illustrate how the stress conditions on the interplate zone can be estimated using a force sensor attached to the back of the upper plate and adjusted via the density and strength of the subducting lithosphere or the lubrication of the plate boundary. The first experimental results reveal the potential of the experimental set-up to investigate the three-dimensional solid-mechanics interactions of lithospheric plates in multiple natural situations.

Description: Reflection seismic studies over the end-glacial Burträsk fault, Skellefteå, Sweden Solid Earth, 2, 9-16, 2011 Author(s): C. Juhlin and B. Lund Reflection seismic data were acquired along a ca. 22 km long profile over the end-glacial Burträsk fault with a nominal receiver and source spacing of 20 m. A steeply dipping reflection can be correlated to the Burträsk fault, indicating that the fault dips at about 55° to the southeast near the surface. The reflection from the fault is rather poorly imaged, probably due to a lateral offset in the fault of about 1 km at this location and the crookedness of the seismic profile in the vicinity of the fault. A more pronounced steeply dipping reflection is observed about 4 km southeast of the Burträsk fault. Based on its correlation with a topographic low at the surface this reflection is interpreted to originate from a fracture zone. There are no signs of large displacements along this zone as the glacial ice receded, but earthquakes could be associated with it today. Other reflections on the processed seismic section may originate from changes in lithological variations in the supra-crustal rocks or from intrusions of more mafic rock. Constraints on the fault geometry provided by the reflection seismic data will help determine what stresses were required to activate the fault when the major rupture along it occurred ca. 9500 years ago.

Description: Some improvements in subbasalt imaging using pre-stack depth migration Solid Earth, 2, 1-7, 2011 Author(s): I. Flecha, R. Carbonell, R. W. Hobbs, and H. Zeyen Subbasalt imaging can be improved by carefully applying pre-stack depth migration. Pre-stack depth migration requires a detailed velocity model and an accurate traveltime calculation. Ray tracing methods are fast but, often fail in calculating traveltimes in complex models, specially, when they feature high velocity contrasts. Finitte difference solutions of the eikonal are more stable and can produce a traveltime field for the whole model avoiding shadow zones. A synthetic test was carried out to check the performance of a new pre-stack depth migration algorithm in a model that features a high velocity layer surrounded by lower velocities. The results reasonably reproduce the original model. The same scheme was used to process long-offset reflection data from the Faroe Shelf where conventional techniques (stack) were insufficient to assess the structure under a basalt layer. Pre-stack depth migration produced an improved image which recovered the main features in the stacked section and allowed to identify some subbasalt coherent events.

Description: 3-D reflection seismic imaging of the Hontomín structure in the Basque–Cantabrian Basin (Spain) Solid Earth, 4, 481-496, 2013 Author(s): J. Alcalde, D. Martí, C. Juhlin, A. Malehmir, D. Sopher, E. Saura, I. Marzán, P. Ayarza, A. Calahorrano, A. Pérez-Estaún, and R. Carbonell 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 Hontomín 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 CO 2 storage site. The 36 km 2 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 Hontomín structure is interpreted to consist of an approximately 10 7 m 2 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 CO 2 can be stored in the target reservoir.

Description: Two-dimensional numerical investigations on the termination of bilinear flow in fractures Solid Earth, 4, 331-345, 2013 Author(s): A. E. Ortiz R., R. Jung, and J. Renner 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 T D 〈 1, a significant pressure increase is not observed at the fracture tip until bilinear flow is succeeded by radial flow at a fixed dimensionless time. For T D 〉 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.

Description: An objective rationale for the choice of regularisation parameter with application to global multiple-frequency S -wave tomography Solid Earth, 4, 357-371, 2013 Author(s): C. Zaroli, M. Sambridge, J.-J. Lévêque, E. Debayle, and G. Nolet 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 ℓ ∞ -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.

Description: A workflow for building and calibrating 3-D geomechanical models &ndash a case study for a gas reservoir in the North German Basin Solid Earth, 4, 347-355, 2013 Author(s): K. Fischer and A. Henk 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.

Description: The microstructural record of porphyroclasts and matrix of partly serpentinized peridotite mylonites – from brittle and crystal-plastic deformation to dissolution–precipitation creep Solid Earth, 4, 315-330, 2013 Author(s): J. Bial and C. A. Trepmann We present microfabrics in high-pressure, metamorphic, partly serpentinized peridotite mylonites from the Voltri Massif, in which porphyroclasts and matrix record independent deformation events. The microfabrics are analysed using polarization microscopy and electron microscopy (SEM/EBSD, EMP). The mylonites contain diopside and olivine porphyroclasts originating from the mantle protolith embedded in a fine-grained matrix consisting mainly of antigorite and minor olivine and pyroxene. The porphyroclasts record brittle and crystal-plastic deformation of the peridotite at upper-mantle conditions and differential stresses of a few hundred MPa. After the peridotites became serpentinized, deformation occurred mainly by dissolution–precipitation creep resulting in a pronounced foliation of the antigorite matrix, crenulation cleavages and newly precipitated olivine and pyroxene from the pore fluid at sites of dilation, i.e. in strain shadows next to porphyroclasts and folded fine-grained antigorite layers. Antigorite reveals a pronounced associated shape preferred orientation (SPO) and crystallographic preferred orientation (CPO) with the basal (001) cleavage plane oriented in the foliation plane. In monomineralic antigorite aggregates at sites of stress concentration around porphyroclasts, a characteristically reduced grain size and deflecting CPO as well as sutured grain boundaries indicate also some contribution of crystal-plastic deformation and grain-boundary migration of antigorite. In contrast, the absence of any intragranular deformation features in newly precipitated olivine in strain shadows reveals that stresses were not sufficiently high to allow for significant dislocation creep of olivine at conditions at which antigorite is stable. The porphyroclast microstructures are not associated with the microstructures of the mylonitic matrix, but are inherited from an independent earlier deformation. The porphyroclasts record a high-stress deformation of the peridotite with dislocation creep of olivine in the upper mantle probably related to rifting processes, whereas the serpentinite matrix records dominantly dissolution–precipitation creep and low stresses during subduction and exhumation.

Description: Can vesicle size distributions assess eruption intensity during volcanic activity? Solid Earth, 4, 373-380, 2013 Author(s): A. LaRue, D. R. Baker, M. Polacci, P. Allard, and N. Sodini 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 Eyjafjallajökull 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 Eyjafjallajökull 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 Eyjafjallajökull 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 Eyjafjallajökull 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 Eyjafjallajökull 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.

Description: New insights on the occurrence of peperites and sedimentary deposits within the silicic volcanic sequences of the Paraná Magmatic Province, Brazil Solid Earth, 5, 121-130, 2014 Author(s): A. C. F. Luchetti, A. J. R. Nardy, F. B. Machado, J. E. O. Madeira, and J. M. Arnosio The PMP (Paraná Magmatic Province) is characterized by lava flows of the Early Cretaceous Serra Geral Formation which covers about 75% of the Paraná Basin (southern and southeastern Brazil), composed of a thick (up to 1600 m) volcanic sequence formed by a succession of petrographically and geochemically distinct units of basic and silicic composition. The whole package must have been emplaced during approximately 3 million years of nearly uninterrupted activity. A few aeolian sandstone layers, indicating arid environmental conditions (the Botucatu Formation), are interlayered in the lower basalts. Above the basalts, the Palmas and Chapecó Members are composed of silicic volcanic rocks (quartz latites, dacites, rhyodacites and rhyolites) and basalts. This paper presents new evidence of sedimentation episodes separating silicic volcanic events, expressed by the occurrence of sedimentary deposits. Interaction between the volcanic bodies and the coeval unconsolidated sediments formed peperites. The sediments were observed between basaltic lava flows and silicic rocks or interlayered in the Palmas-type rocks, between the Chapecó-type rocks and overlying basaltic flows, between silicic bodies of the Palmas and Chapecó types, and interlayered within Palmas-type units. The observed structures indicate that the sediments were still wet and unconsolidated, or weakly consolidated, at the time of volcanism, which, coupled with the sediment features, reflect environmental conditions that are different from those characterizing the Botucatu arid conditions.

Description: Thermal shock and splash effects on burned gypseous soils from the Ebro Basin (NE Spain) Solid Earth, 5, 131-140, 2014 Author(s): J. León, M. Seeger, D. Badía, P. Peters, and M. T. Echeverría Fire is a natural factor of landscape evolution in Mediterranean ecosystems. The middle Ebro Valley has extreme aridity, which results in a low plant cover and high soil erodibility, especially on gypseous substrates. The aim of this research is to analyze the effects of moderate heating on physical and chemical soil properties, mineralogical composition and susceptibility to splash erosion. Topsoil samples (15 cm depth) were taken in the Remolinos mountain slopes (Ebro Valley, NE Spain) from two soil types: Leptic Gypsisol (LP) in a convex slope and Haplic Gypsisol (GY) in a concave slope. To assess the heating effects on the mineralogy we burned the soils at 105 and 205 °C in an oven and to assess the splash effects we used a rainfall simulator under laboratory conditions using undisturbed topsoil subsamples (0–5 cm depth of Ah horizon). LP soil has lower soil organic matter (SOM) and soil aggregate stability (SAS) and higher gypsum content than GY soil. Gypsum and dolomite are the main minerals (〉80%) in the LP soil, while gypsum, dolomite, calcite and quartz have similar proportions in GY soil. Clay minerals (kaolinite and illite) are scarce in both soils. Heating at 105 °C has no effect on soil mineralogy. However, heating to 205 °C transforms gypsum to bassanite, increases significantly the soil salinity (EC) in both soil units (LP and GY) and decreases pH only in GY soil. Despite differences in the content of organic matter and structural stability, both soils show no significant differences (P 〈 0.01) in the splash erosion rates. The size of pores is reduced by heating, as derived from variations in soil water retention capacity.

Description: Seismic visibility of a deep subduction channel – insights from numerical simulation of high-frequency seismic waves emitted from intermediate depth earthquakes Solid Earth, 5, 141-159, 2014 Author(s): W. Friederich, L. Lambrecht, B. Stöckhert, S. Wassmann, and C. Moos Return flow in a deep subduction channel (DSC) has been proposed to explain rapid exhumation of high pressure–low temperature metamorphic rocks, entirely based on the fossil rock record. Supported by thermo-mechanical models, the DSC is envisioned as a thin layer on top of the subducted plate reaching down to minimum depths of about 150 km. We perform numerical simulations of high-frequency seismic wave propagation (1–5 Hz) to explore potential seismological evidence for the in situ existence of a DSC. Motivated by field observations, for modeling purposes we assume a simple block-in-matrix (BIM) structure with eclogitic blocks floating in a serpentinite matrix. Homogenization calculations for BIM structures demonstrate that effective seismic velocities in such composites are lower than in the surrounding oceanic crust and mantle, with nearly constant values along the entire length of the DSC. Synthetic seismograms for receivers at the surface computed for intermediate depth earthquakes in the subducted oceanic crust for models with and without DSC turn out to be markedly influenced by its presence or absence. While for both models P and S waveforms are dominated by delayed high-amplitude guided waves, models with DSC exhibit a very different pattern of seismic arrivals compared to models without DSC. The main reason for the difference is the greater length and width of the low-velocity channel when a DSC is present. Seismic velocity heterogeneity within the DSC or oceanic crust is of minor importance. The characteristic patterns allow for definition of typical signatures by which models with and without DSC may be discriminated. The signatures stably recur in slightly modified form for earthquakes at different depths inside subducted oceanic crust. Available seismological data from intermediate depth earthquakes recorded in the forearc of the Hellenic subduction zone exhibit similar multi-arrival waveforms as observed in the synthetic seismograms for models with DSC. According to our results, observation of intermediate depth earthquakes along a profile across the forearc may allow to test the hypothesis of a DSC and to identify situations where such processes could be active today.

Description: The permeability and elastic moduli of tuff from Campi Flegrei, Italy: implications for ground deformation modelling Solid Earth, 5, 25-44, 2014 Author(s): M. J. Heap, P. Baud, P. G. Meredith, S. Vinciguerra, and T. Reuschlé The accuracy of ground deformation modelling at active volcanoes is a principal requirement in volcanic hazard mitigation. However, the reliability of such models relies on the accuracy of the rock physical property (permeability and elastic moduli) input parameters. Unfortunately, laboratory-derived values on representative rocks are usually rare. To this end we have performed a systematic laboratory study on the influence of pressure and temperature on the permeability and elastic moduli of samples from the two most widespread lithified pyroclastic deposits at the Campi Flegrei volcanic district, Italy. Our data show that the water permeability of Neapolitan Yellow Tuff and a tuff from the Campanian Ignimbrite differ by about 1.5 orders of magnitude. As pressure (depth) increases beyond the critical point for inelastic pore collapse (at an effective pressure of 10–15 MPa, or a depth of about 750 m), permeability and porosity decrease significantly, and ultrasonic wave velocities and dynamic elastic moduli increase significantly. Increasing the thermal stressing temperature increases the permeability and decreases the ultrasonic wave velocities and dynamic elastic moduli of the Neapolitan Yellow Tuff; whereas the tuff from the Campanian Ignimbrite remains unaffected. This difference is due to the presence of thermally unstable zeolites within the Neapolitan Yellow Tuff. For both rocks we also find, under the same pressure conditions, that the dynamic (calculated from ultrasonic wave velocities) and static (calculated from triaxial stress-strain data) elastic moduli differ significantly. The choice of elastic moduli in ground deformation modelling is therefore an important consideration. While we urge that these new laboratory data should be considered in routine ground deformation modelling, we highlight the challenges for ground deformation modelling based on the heterogeneous nature (vertically and laterally) of the rocks that comprise the caldera at Campi Flegrei.

Description: Seasonal changes in the soil hydrological and erosive response depending on aspect, vegetation type and soil water repellency in different Mediterranean microenvironments Solid Earth, 4, 497-509, 2013 Author(s): M. A. Gabarrón-Galeote, J. F. Martínez-Murillo, M. A. Quesada, and J. D. Ruiz-Sinoga 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.

Description: Review of some significant claimed irregularities in Scandinavian postglacial uplift on timescales of tens to thousands of years – earthquakes in Denmark? Solid Earth, 5, 109-118, 2014 Author(s): S. Gregersen and P. H. Voss The postglacial uplift and surrounding subsidence in Scandinavia is well described as close to regular, and the phenomenon is similar on timescales of tens, hundreds and thousands of years studied via geodesy, seismology and geology. Searches for irregularities in the form of earthquakes claimed in the scientific literature have disclosed many earthquakes right after the Ice Age, 9000 yr ago, and some later cases for further evaluation. In a previous report, the present authors have mentioned doubts about the validity of some of the most significant claimed irregularities. In the present paper, a review is made of these significant claimed irregularities in the south-western flank of the Scandinavian postglacial uplift and subsidence via literature studies of geodetic and geological claims of earthquakes as well as discussions in the field. Geodetic observations exist for all of Scandinavia that describe the phenomenon on a scale of 10s–100s of years. Earthquake observations in seismology are of relevance in the same timescales. Geological studies of dated shore lines describe the postglacial vertical earth-surface motion in a quite different timescale of 100s–1000s of years. There is a need for integration of these observations geographically. This is happening in the various timescales in the DynaQlim project. The review finds the claims unlikely to be earthquakes concerning the following: (1) geodynamical motion in the Copenhagen area, (2) a palaeo-earthquake in Læsø and (3) the recently proposed water level discrepancy in the southern part of Denmark. The assessment is less certain, but falls to improbable when concerning (4) proposed palaeo-earthquakes by Hallandsåsen in south-western Sweden.

Description: Use of phytoremediation and biochar to remediate heavy metal polluted soils: a review Solid Earth, 5, 65-75, 2014 Author(s): J. Paz-Ferreiro, H. Lu, S. Fu, A. Méndez, and G. Gascó Anthropogenic activities are resulting in an increase of the use and extraction of heavy metals. Heavy metals cannot be degraded and hence accumulate in the environment, having the potential to contaminate the food chain. This pollution threatens soil quality, plant survival and human health. The remediation of heavy metals deserves attention, but it is impaired by the cost of these processes. Phytoremediation and biochar are two sound environmental technologies which could be at the forefront to mitigate soil pollution. This review provides an overview of the state of the art of the scientific research on phytoremediation and biochar application to remediate heavy-metal-contaminated soils. Research to date has attempted only in a limited number of occasions to combine both techniques, however we discuss the potential advantages of combining both, and the potential mechanisms involved in the interaction between phytoremediators and biochar. We identified specific research needs to ensure a sustainable use of phytoremediation and biochar as remediation tools.

Description: Editorial Note A case of plagiarism: "Modelling of the wave fields by the modification of the matrix method in anisotropic media" published in Solid Earth Discuss., 6, 1–19, 2014 Solid Earth, 5, 161-161, 2014 Author(s): F. Storti No abstract available.

Description: Corrigendum to "3-D reflection seismic imaging of the Hontomín structure in the Basque–Cantabrian Basin (Spain)" published in Solid Earth, 4, 481–496, 2013 Solid Earth, 5, 181-181, 2014 Author(s): J. Alcalde, D. Martí, C. Juhlin, A. Malehmir, D. Sopher, E. Saura, I. Marzán, P. Ayarza, A. Calahorrano, A. Pérez-Estaún, and R. Carbonell No abstract available.

Description: Lithospheric-scale structures in New Guinea and their control on the location of gold and copper deposits Solid Earth, 5, 163-179, 2014 Author(s): L. T. White, M. P. Morse, and G. S. Lister The locations of major gold and copper deposits on the island of New Guinea are considered by many to be controlled by a series of transfer faults that strike N–S to NE–SW, perpendicular to the long axis of the island. The premise is that these faults dilate perpendicular to the regional stress field, forming conduits for metalliferous gases and fluids to drop out of solution. However, the data on which this idea was first proposed were often not presented or, when the data were presented, were of poor quality or low resolution. We therefore present a review of the existing structural interpretations and compare these with several recently published geophysical data sets to determine if the mineralization controlling transfer faults could be observed. These data were used to produce a new lineament map of New Guinea. A comparison of the lineaments with the location of major gold and copper deposits indicates there is a link between the arc-normal structures and mineralization. However, it is only those deposits that are less than 4.5 million years old that could be associated with these structures. Gravity and seismic tomography data indicate that some of these structures could penetrate deep levels of the lithosphere, providing some support to the earlier idea that the arc-normal structures act as conduits for the younger mineral deposits of New Guinea. The gravity data can also be used to infer the location of igneous intrusions at depth, which could have brought metal-bearing fluids and gases closer to the Earth's surface. These regions might be of interest for future exploration campaigns, particularly those areas that are crosscut by deep, vertical faults. However, new exploration models are needed to explain the location of the deposits that are older than 5 Ma.

Description: The ring-shaped thermal field of Stefanos crater, Nisyros Island: a conceptual model Solid Earth, 5, 183-198, 2014 Author(s): M. Pantaleo and T. R. Walter Fumarole fields related to hydrothermal processes release the heat of the underground through permeable pathways. Thermal changes, therefore, are likely to depend also on the size and permeability variation of these pathways. There may be different explanations for the observed permeability changes, such as fault control, lithology, weathering/alteration, heterogeneous sediment accumulation/erosion and physical changes of the fluids (e.g., temperature and viscosity). A common difficulty, however, in surface temperature field studies at active volcanoes is that the parameters controlling the ascending routes of fluids are poorly constrained in general. Here we analyze the crater of Stefanos, Nisyros (Greece), and highlight complexities in the spatial pattern of the fumarole field related to permeability conditions. We combine high-resolution infrared mosaics and grain-size analysis of soils, aiming to elaborate parameters controlling the appearance of the fumarole field. We find a ring-shaped thermal field located within the explosion crater, which we interpret to reflect near-surface contrasts of the soil granulometry and volcanotectonic history at depth. We develop a conceptual model of how the ring-shaped thermal field formed at the Stefanos crater and similarly at other volcanic edifices, highlighting the importance of local permeability contrast that may increase or decrease the thermal fluid flux.

Description: DInSAR Coseismic Deformation of the May 2011 M w 5.1 Lorca Earthquake (southeastern Spain) Solid Earth, 3, 111-119, 2012 Author(s): T. Frontera, A. Concha, P. Blanco, A. Echeverria, X. Goula, R. Arbiol, G. Khazaradze, F. Pérez, and E. Suriñach The coseismic superficial deformation at the region of Lorca (Murcia, southeastern Spain) due to the M w 5.1 earthquake on 11 May 2011 was characterized by a multidisciplinary team, integrating information from DInSAR, GPS and numerical modelling techniques. Despite the moderate magnitude of the event, quantitative information was obtained from the interferometric study of a pair of TerraSAR-X images. The DinSAR results defined the trace of the fault plane and evidenced uplift of the hanging wall block in agreement with the estimated deformation obtained through an elastic rupture dislocation numerical model. Meanwhile for the footwall block, interferometric results showed that tectonic deformation is masked by an important subsidence related to groundwater extraction previously identified at the area of study. Horizontal crustal deformation rates and velocity vectors, obtained from GPS stations existent at the area, were also coherent with the tectonic setting of the southern margin of the Iberian Peninsula and with the focal mechanism calculated for the Lorca event. The analysis of a continuous GPS site in Lorca showed good agreement with the horizontal N–S direction component relative to the numerical model and tectonics of the region. This is the first time at this seismic active area that a multi-technique analysis has been performed immediately after the occurrence of a seismic event, comparing the existing deformation data with a theoretical numerical model based on estimated seismic rupture dislocation.

Description: Floating stones off El Hierro, Canary Islands: xenoliths of pre-island sedimentary origin in the early products of the October 2011 eruption Solid Earth, 3, 97-110, 2012 Author(s): V. R. Troll, A. Klügel, M.-A. Longpré, S. Burchardt, F. M. Deegan, J. C. Carracedo, S. Wiesmaier, U. Kueppers, B. Dahren, L. S. Blythe, T. H. Hansteen, C. Freda, D. A. Budd, E. M. Jolis, E. Jonsson, F. C. Meade, C. Harris, S. E. Berg, L. Mancini, M. Polacci, and K. Pedroza A submarine eruption started off the south coast of El Hierro, Canary Islands, on 10 October 2011 and continues at the time of this writing (February 2012). In the first days of the event, peculiar eruption products were found floating on the sea surface, drifting for long distances from the eruption site. These specimens, which have in the meantime been termed "restingolites" (after the close-by village of La Restinga), appeared as black volcanic "bombs" that exhibit cores of white and porous pumice-like material. Since their brief appearance, the nature and origin of these "floating stones" has been vigorously debated among researchers, with important implications for the interpretation of the hazard potential of the ongoing eruption. The "restingolites" have been proposed to be either (i) juvenile high-silica magma (e.g. rhyolite), (ii) remelted magmatic material (trachyte), (iii) altered volcanic rock, or (iv) reheated hyaloclastites or zeolite from the submarine slopes of El Hierro. Here, we provide evidence that supports yet a different conclusion. We have analysed the textures and compositions of representative "restingolites" and compared the results to previous work on similar rocks found in the Canary Islands. Based on their high-silica content, the lack of igneous trace element signatures, the presence of remnant quartz crystals, jasper fragments and carbonate as well as wollastonite (derived from thermal overprint of carbonate) and their relatively high oxygen isotope values, we conclude that "restingolites" are in fact xenoliths from pre-island sedimentary layers that were picked up and heated by the ascending magma, causing them to partially melt and vesiculate. As they are closely resembling pumice in appearance, but are xenolithic in origin, we refer to these rocks as "xeno-pumice". The El Hierro xeno-pumices hence represent messengers from depth that help us to understand the interaction between ascending magma and crustal lithologies beneath the Canary Islands as well as in similar Atlantic islands that rest on sediment-covered ocean crust (e.g. Cape Verdes, Azores). The occurrence of "restingolites" indicates that crustal recycling is a relevant process in ocean islands, too, but does not herald the arrival of potentially explosive high-silica magma in the active plumbing system beneath El Hierro.

Description: Re-evaluation of the Mentelle Basin, a polyphase rifted margin basin, offshore southwest Australia: new insights from integrated regional seismic datasets Solid Earth, 2, 107-123, 2011 Author(s): D. Maloney, C. Sargent, N. G. Direen, R. W. Hobbs, and D. R. Gröcke Vintage 2-D (two-dimensional) seismic reflection surveys from the sparsely explored Mentelle Basin (western Australian margin) have been reprocessed and integrated with a recent high-quality 2-D seismic survey and stratigraphic borehole data. Interpretation of these data sets allows the internal geometry of the Mentelle Basin fill and depositional history to be reanalysed and new insights into its formation revealed. Basin stratigraphy can be subdivided into several seismically defined megasequences separated by major unconformities related to both breakup between India-Madagascar and Australia-Antarctica in the Valanginian-Late Hauterivian and tectonically-driven switches in deposition through the Albian. Resting on the Valanginian-Late Hauterivian breakup unconformity are several kilometre-scale mounded structures that formed during Late Jurassic to Early Cretaceous extension. These have previously been interpreted as volcanic edifices although direct evidence of volcanic feeder systems is lacking. An alternative interpretation is that these features may be carbonate build-ups. The latter interpretation carries significant climatic ramifications since carbonate build-ups would have formed at high palaeolatitude, ~60° S. Soon after breakup, initial subsidence resulted in a shallow marine environment and deposition of Barremian-Aptian silty-sandy mudstones. As subsidence continued, thick successions of Albian ferruginous black clays were deposited. Internally, seismic megasequences composed of successions of black clays show previously unresolved unconformities, onlapping and downlapping packages, which reflect a complex depositional, rifting and subsidence history at odds with their previous interpretation as open marine sediments. Southwestwards migration of the Kerguelen hotspot led to thermal contraction and subsidence to the present day water depth (~3000 m). This was accompanied by Turonian-Santonian deposition of massive chalk beds, which are unconformably overlain by pelagic Palaeocene-Holocene sediments. This substantial unconformity is related to the diachronous breakup and onset of slow spreading between Australia and Antarctica, which may have led to the reactivation and inversion of basement faults and was followed by rapid seafloor spreading from the Middle Eocene to the present.

Description: Dynamical geochemistry of the mantle Solid Earth, 2, 159-189, 2011 Author(s): G. F. Davies The reconciliation of mantle chemistry with the structure of the mantle inferred from geophysics and dynamical modelling has been a long-standing problem. This paper reviews three main aspects. First, extensions and refinements of dynamical modelling and theory of mantle processing over the past decade. Second, a recent reconsideration of the implications of mantle heterogeneity for melting, melt migration, mantle differentiation and mantle segregation. Third, a recent proposed shift in the primitive chemical baseline of the mantle inferred from observations of non-chondritic 142 Nd in the Earth. It seems most issues can now be resolved, except the level of heating required to maintain the mantle's thermal evolution. A reconciliation of refractory trace elements and their isotopes with the dynamical mantle, proposed and given preliminary quantification by Hofmann, White and Christensen, has been strengthened by work over the past decade. The apparent age of lead isotopes and the broad refractory-element differences among and between ocean island basalts (OIBs) and mid-ocean ridge basalts (MORBs) can now be quantitatively accounted for with some assurance. The association of the least radiogenic helium with relatively depleted sources and their location in the mantle have been enigmatic. The least radiogenic helium samples have recently been recognised as matching the proposed non-chondritic primitive mantle. It has also been proposed recently that noble gases reside in a so-called hybrid pyroxenite assemblage that is the result of melt from fusible pods reacting with surrounding refractory peridotite and refreezing. Hybrid pyroxenite that is off-axis may not remelt and erupt at MORs, so its volatile constituents would recirculate within the mantle. Hybrid pyroxenite is likely to be denser than average mantle, and thus some would tend to settle in the D" zone at the base of the mantle, along with some old subducted oceanic crust. Residence times in D" are longer, so the hybrid pyroxenite there would be less degassed. Plumes would sample both the degassed, enriched old oceanic crust and the gassy, less enriched hybrid pyroxenite and deliver them to OIBs. These findings can account quantitatively for the main He, Ne and Ar isotopic observations. It has been commonly inferred that the MORB source is strongly depleted of incompatible elements. However it has recently been argued that conventional estimates of the MORB source composition fail to take full account of mantle heterogeneity, and in particular focus on an ill-defined "depleted" mantle component while neglecting less common enriched components. Previous estimates have also been tied to the composition of peridotites, but these probably do not reflect the full complement of incompatible elements in the heterogeneous mantle. New estimates that account for enriched mantle components suggest the MORB source complement of incompatibles could be as much as 50–100 % larger than previous estimates. A major difficulty has been the inference that mass balances of incompatible trace elements could only be satisfied if there is a deep enriched layer in the mantle, but the Earth's topography precludes such a layer. The difficulty might be resolved if either the Earth is depleted relative to chondritic or the MORB source is less depleted than previous estimates. Together these factors can certainly resolve the mass balance difficulties.

Description: Positive geothermal anomalies in oceanic crust of Cretaceous age offshore Kamchatka Solid Earth, 2, 191-198, 2011 Author(s): G. Delisle Heat flow measurements were carried out in 2009 offshore Kamchatka during the German-Russian joint-expedition KALMAR. An area with elevated heat flow in oceanic crust of Cretaceous age – detected ~30 yr ago in the course of several Russian heat flow surveys – was revisited. One previous interpretation postulated anomalous lithospheric conditions or a connection between a postulated mantle plume at great depth (〉200 km) as the source for the observed high heat flow. However, the positive heat flow anomaly – as our bathymetric data show – is closely associated with the fragmentation of the western flank of the Meiji Seamount into a horst and graben structure initiated during descent of the oceanic crust into the subduction zone offshore Kamchatka. This paper offers an alternative interpretation, which connects high heat flow primarily with natural convection of fluids in the fragmented rock mass and, as a potential additional factor, high rates of erosion, for which evidence is available from our collected bathymetric image. Given high erosion rates, warm rock material at depth rises to nearer the sea floor, where it cools and causes temporary elevated heat flow.

Description: Domains of Archean mantle lithosphere deciphered by seismic anisotropy – inferences from the LAPNET array in northern Fennoscandia Solid Earth, 2, 303-313, 2011 Author(s): J. Plomerová, L. Vecsey, V. Babuška, and LAPNET Working Group An international LAPNET array (2007–2009, http://www.oulu.fi/sgo-oty/lapnet ) of the POLENET/LAPNET sub-project of the POLENET-IPY consortium, related to seismic and geodetic studies in the Arctic regions, consisted of about 60 broadband seismic stations located on the territory of northern Finland and adjacent parts of Sweden, Norway and Russia. We analyze relative P-wave travel-time deviations evaluated for a subset of 90 teleseismic events recorded by the LAPNET array and show examples of lateral variations of shear-wave splitting to demonstrate variability of fabrics of the Archean mantle lithosphere. The initial results clearly demonstrate the Archean mantle lithosphere consists of domains with consistent fabrics reflecting fossil anisotropic structures. 3-D self-consistent anisotropic models with inclined symmetry axes accommodate two independent sets of body-wave anisotropic observations. Individual domains are delimited by boundaries (sutures), where the anisotropic parameters change. The results obtained from the LAPNET array fill a gap in structural studies of the upper mantle beneath northern Fennoscandia.

Description: Erosion rates deduced from seasonal mass balance along the upper Urumqi River in Tianshan Solid Earth, 2, 283-301, 2011 Author(s): Y. Liu, F. Métivier, J. Gaillardet, B. Ye, P. Meunier, C. Narteau, E. Lajeunesse, T. Han, and L. Malverti We report measurements performed during two complete flow seasons on the Urumqi River, a proglacial mountain stream in the northeastern flank of the Tianshan, an active mountain range in Central Asia. This survey of flow dynamics and sediment transport (dissolved, suspended and bed loads), together with a 25-year record of daily discharge, enables the assessment of secular denudation rates on this high mountain catchment of Central Asia. Our results show that chemical weathering accounts for more than one-third of the total denudation rate. Sediment transported as bed load cannot be neglected in the balance, given that sand and gravel transport accounts for one third of the solid load of the river. Overall, the mean denudation rates are low, averaging 46 t × km −2 × yr −1 (17–18 m Myr −1 ). We furthermore analyse the hydrologic record to show that the long-term sediment budget is not dominated by extreme and rare events but by the total amount of rainfall or annual runoff. The rates we obtain are in agreement with rates obtained from the mass balance reconstruction of the Plio-Quaternary gravely deposits of the foreland but signicantly lower than the rates recently obtained from cosmogenic dating of the Kuitun River sands, west of the Urumqi River. We show that the resolution of this incompatibility may have an important consquence for our understanding of the interplay between erosion and tectonics in the semi-humid ranges of Central Asia.

Description: Influence of the Ringwoodite-Perovskite transition on mantle convection in spherical geometry as a function of Clapeyron slope and Rayleigh number Solid Earth, 2, 315-326, 2011 Author(s): M. Wolstencroft and J. H. Davies We investigate the influence on mantle convection of the negative Clapeyron slope ringwoodite to perovskite and ferro-periclase mantle phase transition, which is correlated with the seismic discontinuity at 660 km depth. In particular, we focus on understanding the influence of the magnitude of the Clapeyron slope (as measured by the Phase Buoyancy parameter, P ) and the vigour of convection (as measured by the Rayleigh number, Ra ) on mantle convection. We have undertaken 76 simulations of isoviscous mantle convection in spherical geometry, varying Ra and P . Three domains of behaviour were found: layered convection for high Ra and more negative P , whole mantle convection for low Ra and less negative P , and transitional behaviour in an intervening domain. The boundary between the layered and transitional domain was fit by a curve P = α Ra β where α = −1.05, and β = −0.1, and the fit for the boundary between the transitional and whole mantle convection domain was α = −4.8, and β = −0.25. These two curves converge at Ra ≈ 2.5 × 10 4 (well below Earth mantle vigour) and P ≈ −0.38. Extrapolating to high Ra , which is likely earlier in Earth history, this work suggests a large transitional domain. It is therefore likely that convection in the Archean would have been influenced by this phase change, with Earth being at least in the transitional domain, if not the layered domain.

Description: Bio-chemostratigraphy of the Barremian-Aptian shallow-water carbonates of the southern Apennines (Italy): pinpointing the OAE1a in a Tethyan carbonate platform Solid Earth, 3, 1-28, 2012 Author(s): M. Di Lucia, A. Trecalli, M. Mutti, and M. Parente Low biostratigraphic resolution and lack of chronostratigraphic calibration hinder precise correlations between platform carbonates and coeval deep-water successions. These are the main obstacle when studying the record of Mesozoic oceanic anoxic events in carbonate platforms. In this paper carbon and strontium isotope stratigraphy are used to produce the first chronostratigraphic calibration of the Barremian-Aptian biostratigraphy of the Apenninic carbonate platform of southern Italy. According to this calibration, the segment of decreasing δ 13 C values, leading to the negative peak that is generally taken as the onset of the Selli event, starts a few metres above the last occurrence of Palorbitolina lenticularis and Voloshinoides murgensis . The following rise of δ 13 C values, corresponding to the interval of enhanced accumulation of organic matter in deep-water sections, ends just below the first acme of Salpingoporella dinarica , which roughly corresponds to the segment of peak δ 13 C values. The whole carbon isotope excursion associated with the oceanic anoxic event 1a is bracketed in the Apenninic carbonate platform between the last occurrence of Voloshinoides murgensis and the "Orbitolina level", characterized by the association of Mesorbitolina parva and Mesorbitolina texana . Since these bioevents have been widely recognized beyond the Apenninic platform, the calibration presented in this paper can be used to pinpoint the interval corresponding to the Early Aptian oceanic anoxic event in other carbonate platforms of central and southern Tethys. This calibration will be particularly useful to interpret the record of the Selli event in carbonate platform sections for which a reliable carbon isotope stratigraphy is not available.

Description: Phanerozoic environments of black shale deposition and the Wilson Cycle Solid Earth, 3, 29-42, 2012 Author(s): J. Trabucho-Alexandre, W. W. Hay, and P. L. de Boer The spatial and temporal distribution of black shales is related to the development of environments in which they accumulate and to a propitious combination of environmental variables. In recent years, much has been done to improve our understanding of the mechanisms behind the temporal distribution of black shales in the Phanerozoic and of the environmental variables that result in their deposition. However, the interpretation of ancient black shale depositional environments is dominated by an oversimplistic set of three depositional models that do not capture their complexity and dynamics. These three models, the restricted circulation, the (open) ocean oxygen minimum and the continental shelf models, are an oversimplification of the variety of black shale depositional environments that arise and coexist throughout the course of a basin's Wilson Cycle, i.e. the dynamic sequence of events and stages that characterise the evolution of an ocean basin, from the opening continental rift to the closing orogeny. We examine the spatial distribution of black shales in the context of the Wilson Cycle using examples from the Phanerozoic. It is shown that the geographical distribution of environments of black shale deposition and the position of black shales in the basin infill sequence strongly depend on basin evolution, which controls the development of sedimentary environments where black shales may be deposited. The nature of the black shales that are deposited, i.e. lithology and type of organic matter, also depends on basin evolution and palaeogeography. We propose that in studies of black shales more attention should be given to the sedimentary processes that have led to their formation and to the interpretation of their sedimentary environments.

Description: Tomography of the 2011 Iwaki earthquake (M 7.0) and Fukushima nuclear power plant area Solid Earth, 3, 43-51, 2012 Author(s): P. Tong, D. Zhao, and D. Yang High-resolution tomographic images of the crust and upper mantle in and around the area of the 2011 Iwaki earthquake (M 7.0) and the Fukushima nuclear power plant are determined by inverting a large number of high-quality arrival times with both the finite-frequency and ray tomography methods. The Iwaki earthquake and its aftershocks mainly occurred in a boundary zone with strong variations in seismic velocity and Poisson's ratio. Prominent low-velocity and high Poisson's ratio zones are revealed under the Iwaki source area and the Fukushima nuclear power plant, which may reflect fluids released from the dehydration of the subducting Pacific slab under Northeast Japan. The 2011 Tohoku-oki earthquake (Mw 9.0) caused static stress transfer in the overriding Okhotsk plate, resulting in the seismicity in the Iwaki source area that significantly increased immediately following the Tohoku-oki mainshock. Our results suggest that the Iwaki earthquake was triggered by the ascending fluids from the Pacific slab dehydration and the stress variation induced by the Tohoku-oki mainshock. The similar structures under the Iwaki source area and the Fukushima nuclear power plant suggest that the security of the nuclear power plant site should be strengthened to withstand potential large earthquakes in the future.

Description: Asymmetry of high-velocity lower crust on the South Atlantic rifted margins and implications for the interplay of magmatism and tectonics in continental breakup Solid Earth, 5, 1011-1026, 2014 Author(s): K. Becker, D. Franke, R. Trumbull, M. Schnabel, I. Heyde, B. Schreckenberger, H. Koopmann, K. Bauer, W. Jokat, and C. M. Krawczyk High-velocity lower crust (HVLC) and seaward-dipping reflector (SDR) sequences are typical features of volcanic rifted margins. However, the nature and origin of HVLC is under discussion. Here we provide a comprehensive analysis of deep crustal structures in the southern segment of the South Atlantic and an assessment of HVLC along the margins. Two new seismic refraction lines off South America fill a gap in the data coverage and together with five existing velocity models allow for a detailed investigation of the lower crustal properties on both margins. An important finding is the major asymmetry in volumes of HVLC on the conjugate margins. The seismic refraction lines across the South African margin reveal cross-sectional areas of HVLC 4 times larger than at the South American margin, a finding that is opposite to the asymmetric distribution of the flood basalts in the Paraná–Etendeka Large Igneous Province. Also, the position of the HVLC with respect to the SDR sequences varies consistently along both margins. Close to the Falkland–Agulhas Fracture Zone in the south, a small body of HVLC is not accompanied by SDRs. In the central portion of both margins, the HVLC is below the inner SDR wedges while in the northern area, closer to the Rio Grande Rise-Walvis Ridge, large volumes of HVLC extend far seaward of the inner SDRs. This challenges the concept of a simple extrusive/intrusive relationship between SDR sequences and HVLC, and it provides evidence for formation of the HVLC at different times during the rifting and breakup process. We suggest that the drastically different HVLC volumes are caused by asymmetric rifting in a simple-shear-dominated extension.

Description: Focal mechanism and depth of the 1956 Amorgos twin earthquakes from waveform matching of analogue seismograms Solid Earth, 5, 1027-1044, 2014 Author(s): A. Brüstle, W. Friederich, T. Meier, and C. Gross Historic analogue seismograms of the large 1956 Amorgos twin earthquakes which occurred in the volcanic arc of the Hellenic subduction zone (HSZ) were collected, digitized and reanalyzed to obtain refined estimates of their depth and focal mechanism. In total, 80 records of the events from 29 European stations were collected and, if possible, digitized. In addition, bulletins were searched for instrument parameters required to calculate transfer functions for instrument correction. A grid search based on matching the digitized historic waveforms to complete synthetic seismograms was then carried out to infer optimal estimates for depth and focal mechanism. Owing to incomplete or unreliable information on instrument parameters and frequently occurring technical problems during recording, such as writing needles jumping off mechanical recording systems, much less seismograms than collected proved suitable for waveform matching. For the first earthquake, only seven seismograms from three different stations at Stuttgart (STU), Göttingen (GTT) and Copenhagen (COP) could be used. Nevertheless, the waveform matching grid search yields two stable misfit minima for source depths of 25 and 50 km. Compatible fault plane solutions are either of normal faulting or thrusting type. A separate analysis of 42 impulsive first-motion polarities taken from the International Seismological Summary (ISS bulletin) excludes the thrusting mechanism and clearly favors a normal faulting solution with at least one of the potential fault planes striking in SW–NE direction. This finding is consistent with the local structure and microseismic activity of the Santorini–Amorgos graben. Since crustal thickness in the Amorgos area is generally less than 30 km, a source depth of 25 km appears to be more realistic. The second earthquake exhibits a conspicuously high ratio of body wave to surface wave amplitudes suggesting an intermediate-depth event located in the Hellenic Wadati–Benioff zone. This hypothesis is supported by a focal mechanism analysis based on first-motion polarities, which indicates a mechanism very different from that of the first event. A waveform matching grid search done to support the intermediate-depth hypothesis proved not to be fruitful because the body wave phases are overlain by strong surface wave coda of the first event inhibiting a waveform match. However, body to surface wave amplitude ratios of a modern intermediate-depth event with an epicenter close to the island of Milos observed at stations of the German Regional Seismic Network (GRSN) exhibit a pattern similar to the one observed for the second event with high values in a frequency band between 0.05 Hz and 0.3 Hz. In contrast, a shallow event with an epicenter in western Crete and nearly identical source mechanism and magnitude, shows very low ratios of body and surface wave amplitude up to 0.17 Hz and higher ratios only beyond that frequency. Based on this comparison with a modern event, we estimate the source depth of the second event to be greater than 100 km. The proximity in time and space of the two events suggests a triggering of the second, potentially deep event by the shallow first one.

Description: Coffee husk mulch on soil erosion and runoff: experiences under rainfall simulation experiment Solid Earth, 5, 851-862, 2014 Author(s): H. Moreno-Ramón, S. J. Quizembe, and S. Ibáñez-Asensio The high erosion rates found in the agriculture land make valuable the use of mulches to control the soil and water losses. Coffee husk ( Coffea canephora var. robusta ) can be one of those mulches. This paper evaluates how to apply the mulch in order to obtain the best effectiveness. An experimental factorial design 4 × 3 × 2 with two replicates was designed in a greenhouse with a total number of 48 cases. All the samples were deposited in trays of 0.51 m 2 and applied a simulated rain of 122 mm h −1 during 21 min. The factors examined were the following: four soil classes; three treatments – buried (B), surface (S) and non-residue (C) – and the presence (WC) or absence (WOC) of the soil surface crusting. The coffee husk residue (S and B treatments) reduced runoff by 10.2 and 46% respectively, soil losses by 78.3 and 88.7% and sediment concentration by 77 and 84.4%. The infiltration rate increased on average by 104 and 167%, and time to runoff by 1.58 and 2.07 min respectively. Coffee husk is an efficient mulch to reduce the soil and water losses, although it could not completely cushion the influence of crust.

Description: Exploring the shallow structure of the San Ramón thrust fault in Santiago, Chile (~33.5° S), using active seismic and electric methods Solid Earth, 5, 837-849, 2014 Author(s): D. Díaz, A. Maksymowicz, G. Vargas, E. Vera, E. Contreras-Reyes, and S. Rebolledo The crustal-scale west-vergent San Ramón thrust fault system, which lies at the foot of the main Andean Cordillera in central Chile, is a geologically active structure with manifestations of late Quaternary complex surface rupture on fault segments along the eastern border of the city of Santiago. From the comparison of geophysical and geological observations, we assessed the subsurface structural pattern that affects the sedimentary cover and rock-substratum topography across fault scarps, which is critical for evaluating structural models and associated seismic hazard along the related faults. We performed seismic profiles with an average length of 250 m, using an array of 24 geophones (Geode), with 25 shots per profile, to produce high-resolution seismic tomography to aid in interpreting impedance changes associated with the deformed sedimentary cover. The recorded travel-time refractions and reflections were jointly inverted by using a 2-D tomographic approach, which resulted in variations across the scarp axis in both the velocities and the reflections that are interpreted as the sedimentary cover-rock substratum topography. Seismic anisotropy observed from tomographic profiles is consistent with sediment deformation triggered by west-vergent thrust tectonics along the fault. Electrical soundings crossing two fault scarps were used to construct subsurface resistivity tomographic profiles, which reveal systematic differences between lower resistivity values in the hanging wall with respect to the footwall of the geological structure, and clearly show well-defined east-dipping resistivity boundaries. These boundaries can be interpreted in terms of structurally driven fluid content change between the hanging wall and the footwall of the San Ramón fault. The overall results are consistent with a west-vergent thrust structure dipping ~55° E in the subsurface beneath the piedmont sediments, with local complexities likely associated with variations in fault surface rupture propagation, fault splays and fault segment transfer zones.

Description: Evidence of magma activation beneath the Harrat Lunayyir basaltic field (Saudi Arabia) from attenuation tomography Solid Earth, 5, 873-882, 2014 Author(s): I. Koulakov, S. El Khrepy, N. Al-Arifi, I. Sychev, and P. Kuznetsov We present a seismic attenuation model for the crust beneath the Cenozoic basaltic field of Harrat Lunayyir (western Saudi Arabia), where a strong seismic swarm occurred in 2009. The tomography inversion uses the envelope shape of the S wave seismograms from over 300 strong events ( M 〈 3.5). The resulting attenuation structures appear to be consistent with the distribution of seismic velocities. The obtained 3-D attenuation model distinguishes the low-attenuation zones down to 5 km depth corresponding to the rigid basaltic cover. At greater depths, we detect a high-attenuation anomaly coinciding with the main seismicity cluster. We propose that this zone corresponds to the upper part of the conduit area ascending from deeper magma sources. According to the distributions of local events, fluids and melts from this conduit appear to reach a depth of ∼2 km, but were not able to reach the surface and cause the eruption in 2009.

Description: Velocity structure and the role of fluids in the West Bohemia Seismic Zone Solid Earth, 5, 863-872, 2014 Author(s): C. Alexandrakis, M. Calò, F. Bouchaala, and V. Vavryčuk In this study, we apply the double-difference tomography to investigate the detailed 3-D structure within and around the Nový Kostel Seismic Zone, an area in the Czech Republic known for frequent occurrences of earthquake swarms. We use data from the 2008 swarm since it has already been analysed in terms of earthquake focal mechanisms, principal faults, tectonic stress and foci migration. We selected about 500 microearthquakes recorded at 22 local seismic stations of the West Bohemia seismic monitoring network (WEBNET). Applying double-difference tomography, combined with weighted average model (WAM) post-processing to correct for parameter dependence effects, we produce and interpret 3-D models of the Vp-to-Vs ratio (Vp/Vs) in and around the focal zone. The modelled Vp/Vs ratio shows several distinct structures, namely an area of high Vp/Vs ratio correlating with the foci of the microearthquakes, and a layer of low values directly above it. These structures may reflect changes in lithology and/or fluid concentration. The overlaying low Vp/Vs ratio layer coincides with the base of the Fichtelgebirge (Smrčiny) granitic intrusion. It is possible that the base of the layer acts as a fluid trap and an upper limit to the seismicity, resulting in observed periodic swarms.

Description: Transport processes at quartz–water interfaces: constraints from hydrothermal grooving experiments Solid Earth, 5, 883-899, 2014 Author(s): K. Klevakina, J. Renner, N. Doltsinis, and W. Adeagbo We performed hydrothermal annealing experiments on quartzite samples at temperatures of 392 to 568 °C and fluid pressures of 63 to 399 MPa for up to 120 h, during which hydrothermal grooves developed on the free surfaces of the samples. An analysis of surface topology and groove characteristics with an atomic force microscope revealed a range of surface features associated with the simultaneous and successive operation of several processes partly depending on crystal orientation during the various stages of an experiment. Initially, dissolution at the quartzite-sample surface occurs to saturate the fluid in the capsule with SiO 2 . Subsequently, grooving controlled by diffusion processes takes place parallel to dissolution and precipitation due to local differences in solubility. Finally, quench products develop on grain surfaces during the termination of experiments. The average groove-root angle amounts to about 160°, varying systematically with misorientation between neighboring grains and depending slightly on temperature and run duration. The grooving is thermally activated, i.e., groove depth ranging from 5 nm to several micrometers for the entire suite of experiments generally increases with temperature and/or run time. We use Mullins' classical theories to constrain kinetic parameters for the transport processes controlling the grooving. In the light of previous measurements of various diffusion coefficients in the system SiO 2 –H 2 O, interface diffusion of Si is identified as the most plausible rate-controlling process. Grooving could potentially proceed faster by diffusion through the liquid if the fluid were not convecting in the capsule. Characteristic times of healing of microfractures in hydrous environments constrained from these kinetic parameters are consistent with the order of magnitude of timescales over which postseismic healing occurs in situ according to geophysical surveys and recurrence intervals of earthquakes.

Description: Simulation of seismic waves at the earth's crust (brittle–ductile transition) based on the Burgers model Solid Earth, 5, 1001-1010, 2014 Author(s): J. M. Carcione, F. Poletto, B. Farina, and A. Craglietto The earth's crust presents two dissimilar rheological behaviors depending on the in situ stress-temperature conditions. The upper, cooler part is brittle, while deeper zones are ductile. Seismic waves may reveal the presence of the transition but a proper characterization is required. We first obtain a stress–strain relation, including the effects of shear seismic attenuation and ductility due to shear deformations and plastic flow. The anelastic behavior is based on the Burgers mechanical model to describe the effects of seismic attenuation and steady-state creep flow. The shear Lamé constant of the brittle and ductile media depends on the in situ stress and temperature through the shear viscosity, which is obtained by the Arrhenius equation and the octahedral stress criterion. The P and S wave velocities decrease as depth and temperature increase due to the geothermal gradient, an effect which is more pronounced for shear waves. We then obtain the P−S and SH equations of motion recast in the velocity-stress formulation, including memory variables to avoid the computation of time convolutions. The equations correspond to isotropic anelastic and inhomogeneous media and are solved by a direct grid method based on the Runge–Kutta time stepping technique and the Fourier pseudospectral method. The algorithm is tested with success against known analytical solutions for different shear viscosities. A realistic example illustrates the computation of surface and reverse-VSP synthetic seismograms in the presence of an abrupt brittle–ductile transition.

Description: 3-D geomechanical–numerical model of the contemporary crustal stress state in the Alberta Basin (Canada) Solid Earth, 5, 1123-1149, 2014 Author(s): K. Reiter and O. Heidbach In the context of examining the potential usage of safe and sustainable geothermal energy in the Alberta Basin, whether in deep sediments or crystalline rock, the understanding of the in situ stress state is crucial. It is a key challenge to estimate the 3-D stress state at an arbitrarily chosen point in the crust, based on sparsely distributed in situ stress data. To address this challenge, we present a large-scale 3-D geomechanical–numerical model (700 km × 1200 km × 80 km) from a large portion of the Alberta Basin, to provide a 3-D continuous quantification of the contemporary stress orientations and stress magnitudes. To calibrate the model, we use a large database of in situ stress orientation (321 S Hmax ) as well as stress magnitude data (981 S V , 1720 S hmin and 2 (+11) S Hmax ) from the Alberta Basin. To find the best-fit model, we vary the material properties and primarily the displacement boundary conditions of the model. This study focusses in detail on the statistical calibration procedure, because of the large amount of available data, the diversity of data types, and the importance of the order of data tests. The best-fit model provides the total 3-D stress tensor for nearly the whole Alberta Basin, and allows estimation of stress orientation and stress magnitudes in advance of any well. First-order implications for the well design and configuration of enhanced geothermal systems are revealed. Systematic deviations of the modelled stress from the in situ data are found for stress orientations in the Peace River and the Bow Island Arch as well as for leak-off test magnitudes.

Description: Wave-equation-based travel-time seismic tomography – Part 2: Application to the 1992 Landers earthquake ( M w 7.3) area Solid Earth, 5, 1169-1188, 2014 Author(s): P. Tong, D. Zhao, D. Yang, X. Yang, J. Chen, and Q. Liu High-resolution 3-D P and S wave crustal velocity and Poisson's ratio models of the 1992 Landers earthquake ( M w 7.3) area are determined iteratively by a wave-equation-based travel-time seismic tomography (WETST) technique. The details of data selection, synthetic arrival-time determination, and trade-off analysis of damping and smoothing parameters are presented to show the performance of this new tomographic inversion method. A total of 78 523 P wave and 46 999 S wave high-quality arrival-time data from 2041 local earthquakes recorded by 275 stations during the period of 1992–2013 are used to obtain the final tomographic models, which cost around 10 000 CPU hours. Checkerboard resolution tests are conducted to verify the reliability of inversion results for the chosen seismic data and the wave-equation-based travel-time seismic tomography method. Significant structural heterogeneities are revealed in the crust of the 1992 Landers earthquake area which may be closely related to the local seismic activities. Strong variations of velocity and Poisson's ratio exist in the source regions of the Landers and three other nearby strong earthquakes. Most seismicity occurs in areas with high-velocity and low Poisson's ratio, which may be associated with the seismogenic layer. Pronounced low-velocity anomalies revealed in the lower crust along the Elsinore, the San Jacinto, and the San Andreas faults may reflect the existence of fluids in the lower crust. The recovery of these strong heterogeneous structures is facilitated by the use of full wave equation solvers and WETST and verifies their ability in generating high-resolution tomographic models.

Description: Wave-equation-based travel-time seismic tomography – Part 1: Method Solid Earth, 5, 1151-1168, 2014 Author(s): P. Tong, D. Zhao, D. Yang, X. Yang, J. Chen, and Q. Liu In this paper, we propose a wave-equation-based travel-time seismic tomography method with a detailed description of its step-by-step process. First, a linear relationship between the travel-time residual Δ t = T obs –T syn and the relative velocity perturbation δ c( x )/c( x ) connected by a finite-frequency travel-time sensitivity kernel K ( x ) is theoretically derived using the adjoint method. To accurately calculate the travel-time residual Δ t , two automatic arrival-time picking techniques including the envelop energy ratio method and the combined ray and cross-correlation method are then developed to compute the arrival times T syn for synthetic seismograms. The arrival times T obs of observed seismograms are usually determined by manual hand picking in real applications. Travel-time sensitivity kernel K ( x ) is constructed by convolving a~forward wavefield u ( t , x ) with an adjoint wavefield q ( t , x ). The calculations of synthetic seismograms and sensitivity kernels rely on forward modeling. To make it computationally feasible for tomographic problems involving a large number of seismic records, the forward problem is solved in the two-dimensional (2-D) vertical plane passing through the source and the receiver by a high-order central difference method. The final model is parameterized on 3-D regular grid (inversion) nodes with variable spacings, while model values on each 2-D forward modeling node are linearly interpolated by the values at its eight surrounding 3-D inversion grid nodes. Finally, the tomographic inverse problem is formulated as a regularized optimization problem, which can be iteratively solved by either the LSQR solver or a~nonlinear conjugate-gradient method. To provide some insights into future 3-D tomographic inversions, Fréchet kernels for different seismic phases are also demonstrated in this study.

Description: Preface: Soil processes in cold-climate environments Solid Earth, 5, 1205-1208, 2014 Author(s): M. Oliva, P. Pereira, J. Bockheim, and A. Navas

Description: Thermal conditions during deformation of partially molten crust from TitaniQ geothermometry: rheological implications for the anatectic domain of the Araçuaí belt, eastern Brazil Solid Earth, 5, 1223-1242, 2014 Author(s): G. C. G. Cavalcante, A. Vauchez, C. Merlet, M. Egydio-Silva, M. H. Bezerra de Holanda, and B. Boyer During the Neoproterozoic orogeny, the middle crust of the Araçuaí belt underwent widespread partial melting. At the regional scale, this anatectic domain is characterized by a progressive rotation of the flow direction from south to north, suggesting a 3-D deformation of the anatectic middle crust. To better determine whether melt volumes present in the anatectic middle crust of the Araçuaí orogen were large enough to allow a combination of gravity-driven and convergence-driven deformation, we used the titanium-in-quartz (TitaniQ) geothermometer to estimate the crystallization temperatures of quartz grains in the anatectic rocks. When possible, we compared these estimates with thermobarometric estimates from traditional exchange geothermobarometers applied to neighboring migmatitic kinzigites. TitaniQ temperatures range from 750 to 900 °C, suggesting that quartz starts crystallizing at minimum temperatures of ≥ 800 °C. These results, combined with the bulk-rock chemical composition of diatexites, allows the estimation of a minimum of ~ 30% melt and a corresponding viscosity of ~ 10 9 –10 10 Pa s. Such a minimum melt content and low viscosity are in agreement with interconnected melt networks observed in the field. Considering that these characteristics are homogeneous over a wide area, this supports the finding that the strength of the middle crust was severely weakened by extensive partial melting, making it prone to gravity-driven flow and lateral extrusion.

Description: Relative tectonic activity classification in the Kermanshah area, western Iran Solid Earth, 5, 1277-1291, 2014 Author(s): M. Arian and Z. Aram Due to the closing of the subduction zone and the collision of the Arabian and Eurasian plates, the High Zagros region has always been affected by a wide range of tectonic variations. In this research, the Gharasu River basin, which is located in the Kermanshah area, was selected as the study area, six geomorphic indices were calculated, and the results of each one were divided into three classes. Then, using the indices, relative tectonic activity was calculated, and their values were classified and analysed in four groups. Regions were identified as very high, high, moderate and low. In analysing the results and combining them with field observations and regional geology, the results were often associated and justified with field evidence. The highest value is located on the Dokeral anticline in the crush zone in Zagros. Most of the areas with high and moderate values of Index of Active Tectonics (IAT) are also located in the crush zone in Zagros. Crushing in this zone is because of the main fault mechanism of the Zagros region. The result in this paper confirms previous research in this region. At the eastern end of the study area, the value of Iat is high, which could be the result of the Sarab and Koh-e Sefid fault mechanisms.

Description: Picroilmenites in Yakutian kimberlites: variations and genetic models Solid Earth, 5, 915-938, 2014 Author(s): I. V. Ashchepkov, N. V. Alymova, A. M. Logvinova, N. V. Vladykin, S. S. Kuligin, S. I. Mityukhin, H. Downes, Yu. B. Stegnitsky, S. A. Prokopiev, R. F. Salikhov, V. S. Palessky, and O. S. Khmel'nikova Major and trace element variations in picroilmenites from Late Devonian kimberlite pipes in Siberia reveal similarities within the region in general, but show individual features for ilmenites from different fields and pipes. Empirical ilmenite thermobarometry (Ashchepkov et al., 2010), as well as common methods of mantle thermobarometry and trace element geochemical modeling, shows long compositional trends for the ilmenites. These are a result of complex processes of polybaric fractionation of protokimberlite melts, accompanied by the interaction with mantle wall rocks and dissolution of previous wall rock and metasomatic associations. Evolution of the parental magmas for the picroilmenites was determined for the three distinct phases of kimberlite activity from Yubileynaya and nearby Aprelskaya pipes, showing heating and an increase of Fe# (Fe# = Fe / (Fe + Mg) a.u.) of mantle peridotite minerals from stage to stage and splitting of the magmatic system in the final stages. High-pressure (5.5–7.0 GPa) Cr-bearing Mg-rich ilmenites (group 1) reflect the conditions of high-temperature metasomatic rocks at the base of the mantle lithosphere. Trace element patterns are enriched to 0.1–10/relative to primitive mantle (PM) and have flattened, spoon-like or S- or W-shaped rare earth element (REE) patterns with Pb 〉 1. These result from melting and crystallization in melt-feeding channels in the base of the lithosphere, where high-temperature dunites, harzburgites and pyroxenites were formed. Cr-poor ilmenite megacrysts (group 2) trace the high-temperature path of protokimberlites developed as result of fractional crystallization and wall rock assimilation during the creation of the feeder systems prior to the main kimberlite eruption. Inflections in ilmenite compositional trends probably reflect the mantle layering and pulsing melt intrusion during melt migration within the channels. Group 2 ilmenites have inclined REE enriched patterns (10–100)/PM with La / Yb n ~ 10–25, similar to those derived from kimberlites, with high-field-strength elements (HFSE) peaks (typical megacrysts). A series of similar patterns results from polybaric Assimilation + fractional crystallization (AFC) crystallization of protokimberlite melts which also precipitated sulfides (Pb 〈 1) and mixed with partial melts from garnet peridotites. Relatively low-Ti ilmenites with high-Cr content (group 3) probably crystallized in the metasomatic front under the rising protokimberlite source and represent the product of crystallization of segregated partial melts from metasomatic rocks. Cr-rich ilmenites are typical of veins and veinlets in peridotites crystallized from highly contaminated magma intruded into wall rocks in different levels within the mantle columns. Ilmenites which have the highest trace element contents (1000/PM) have REE patterns similar to those of perovskites. Low Cr contents suggest relatively closed system fractionation which occurred from the base of the lithosphere up to the garnet–spinel transition, according to monomineral thermobarometry for Mir and Dachnaya pipes. Restricted trends were detected for ilmenites from Udachnaya and most other pipes from the Daldyn–Alakit fields and other regions (Nakyn, Upper Muna and Prianabarie), where ilmenite trends extend from the base of the lithosphere mainly up to 4.0 GPa. Interaction of the megacryst forming melts with the mantle lithosphere caused heating and HFSE metasomatism prior to kimberlite eruption.

Description: Biochar increases plant-available water in a sandy loam soil under an aerobic rice crop system Solid Earth, 5, 939-952, 2014 Author(s): M. T. de Melo Carvalho, A. de Holanda Nunes Maia, B. E. Madari, L. Bastiaans, P. A. J. van Oort, A. B. Heinemann, M. A. Soler da Silva, F. A. Petter, B. H. Marimon Jr., and H. Meinke The main objective of this study was to assess the impact of biochar rate (0, 8, 16 and 32 Mg ha −1 ) on the water retention capacity (WRC) of a sandy loam Dystric Plinthosol. The applied biochar was a by-product of slow pyrolysis (∼450 °C) of eucalyptus wood, milled to pass through a 2000 μm sieve that resulted in a material with an intrinsic porosity ≤10 μm and a specific surface area of ∼3.2 m 2 g −1 . The biochar was incorporated into the top 15 cm of the soil under an aerobic rice system. Our study focused on both the effects on WRC and rice yields 2 and 3 years after its application. Undisturbed soil samples were collected from 16 plots in two soil layers (5–10 and 15–20 cm). Soil water retention curves were modelled using a nonlinear mixed model which appropriately accounts for uncertainties inherent of spatial variability and repeated measurements taken within a specific soil sample. We found an increase in plant-available water in the upper soil layer proportional to the rate of biochar, with about 0.8% for each Mg ha −1 biochar amendment 2 and 3 years after its application. The impact of biochar on soil WRC was most likely related to an effect in overall porosity of the sandy loam soil, which was evident from an increase in saturated soil moisture and macro porosity with 0.5 and 1.6% for each Mg ha −1 of biochar applied, respectively. The increment in soil WRC did not translate into an increase in rice yield, essentially because in both seasons the amount of rainfall during the critical period for rice production exceeded 650 mm. The use of biochar as a soil amendment can be a worthy strategy to guarantee yield stability under short-term water-limited conditions. Our findings raise the importance of assessing the feasibility of very high application rates of biochar and the inclusion of a detailed analysis of its physical and chemical properties as part of future investigations.

Description: Biochar can be used to capture essential nutrients from dairy wastewater and improve soil physico-chemical properties Solid Earth, 5, 953-962, 2014 Author(s): T. A. Ghezzehei, D. V. Sarkhot, and A. A. Berhe Recently, the potential for biochar use to recapture excess nutrients from dairy wastewater has been a focus of a growing number of studies. It is suggested that biochar produced from locally available excess biomass can be important in reducing release of excess nutrient elements from agricultural runoff, improving soil productivity, and long-term carbon (C) sequestration. Here we present a review of a new approach that is showing promise for the use of biochar for nutrient capture. Using batch sorption experiments, it has been shown that biochar can adsorb up to 20–43% of ammonium and 19–65% of the phosphate in flushed dairy manure in 24 h. These results suggest a potential of biochar for recovering essential nutrients from dairy wastewater and improving soil fertility if the enriched biochar is returned to soil. Based on the sorption capacity of 2.86 and 0.23 mg ammonium and phosphate, respectively, per gram of biochar and 10–50% utilization of available excess biomass, in the state of California (US) alone, 11 440 to 57 200 tonnes of ammonium-N and 920–4600 tonnes of phosphate can be captured from dairy waste each year while at the same time disposing up to 8–40 million tons of excess biomass.

Description: Evaluating the importance of surface soil contributions to reservoir sediment in alpine environments: a combined modelling and fingerprinting approach in the Posets-Maladeta Natural Park Solid Earth, 5, 963-978, 2014 Author(s): L. Palazón, L. Gaspar, B. Latorre, W. H. Blake, and A. Navas Soil in alpine environments plays a key role in the development of ecosystem services and in order to maintain and preserve this important resource, information is required on processes that lead to soil erosion. Similar to other mountain alpine environments, the Benasque catchment is characterised by temperatures below freezing that can last from November to April, intense rainfall events, typically in spring and autumn, and rugged topography which makes assessment of erosion challenging. Indirect approaches to soil erosion assessment, such as combined model approaches, offer an opportunity to evaluate soil erosion in such areas. In this study (i) the SWAT (Soil and Water Assessment Tool) hydrological and erosion model and (ii) sediment fingerprinting procedures were used in parallel to assess the viability of a combined modelling and tracing approach to evaluate soil erosion processes in the area of the Posets-Maladeta Natural Park (central Spanish Pyrenees). Soil erosion rates and sediment contribution of potential sediment sources defined by soil type (Kastanozems/Phaeozems; Fluvisols and Cambisols) were assessed. The SWAT model suggested that, with the highest specific sediment yields, Cambisols are the main source of sediment in the Benasque catchment and Phaeozems and Fluvisols were identified as the lowest sediment contributors. Spring and winter model runs gave the highest and lowest specific sediment yield, respectively. In contrast, sediment fingerprinting analysis identified Fluvisols, which dominate the riparian zone, as the main sediment source at the time of sampling. This indicates the importance of connectivity as well as potential differences in the source dynamic of material in storage versus that transported efficiently from the system at times of high flow. The combined approach enabled us to better understand soil erosion processes in the Benasque alpine catchment, wherein SWAT identified areas of potential high sediment yield in large flood events but sediment fingerprinting identified areas that, due to high connectivity, contribute more to channel-stored sediment deposits.

Description: Degradation of buried ice and permafrost in the Veleta cirque (Sierra Nevada, Spain) from 2006 to 2013 as a response to recent climate trends Solid Earth, 5, 979-993, 2014 Author(s): A. Gómez-Ortiz, M. Oliva, F. Salvador-Franch, M. Salvà-Catarineu, D. Palacios, J. J. de Sanjosé-Blasco, L. M. Tanarro-García, J. Galindo-Zaldívar, and C. Sanz de Galdeano The Veleta cirque is located at the foot of the Veleta peak, one of the highest summits of the Sierra Nevada National Park (southern Spain). This cirque was the source of a glacier valley during the Quaternary cold periods. During the Little Ice Age it sheltered a small glacier, the most southerly in Europe, about which we have possessed written records since the 17th century. This glacier still had ice residues until the mid-20th century. This ice is no longer visible, but a residue persists along with discontinuous permafrost trapped under strata of rock blocks that make up an incipient rock glacier. From 2006 to 2013, this rock glacier was monitored by measurement of the temperature of the active layer, the degree of snow cover on the ground, movements of the body of the rock glacier and geophysical prospection inside it. The results show that the relict ice and trapped permafrost have been steadily declining. The processes that explain this degradation occur in chain, starting from the external radiation that affects the ground in summer, which is when the temperatures are higher. In effect, when this radiation steadily melts the snow on the ground, the thermal expansive wave advances into the heart of the active layer, reaching the ceiling of the frozen mass, which it then degrades and melts. In this entire linked process, the circulation of meltwaters fulfils a highly significant function, as they act as heat transmitters. The complementary nature of these processes explains the subsidence and continuous changes in the entire clastic pack and the melting of the frozen ceiling on which it rests. This happens in summer in just a few weeks. All these events, in particular the geomorphological ones, take place on the Sierra Nevada peaks within certain climate conditions that are at present unfavourable to the maintenance of snow on the ground in summer. These conditions could be related to recent variations in the climate, starting in the mid-19th century and most markedly since the second half of the 20th century. The work and results highlight the climate sensitivity of the peaks of the Sierra Nevada to the effect of climate change and its impact on the dynamics of ecosystems, which is a benchmark for evaluating the current evolution of landscapes of Mediterranean high mountains.

Description: Future accreted terranes: a compilation of island arcs, oceanic plateaus, submarine ridges, seamounts, and continental fragments Solid Earth, 5, 1243-1275, 2014 Author(s): J. L. Tetreault and S. J. H. Buiter Allochthonous accreted terranes are exotic geologic units that originated from anomalous crustal regions on a subducting oceanic plate and were transferred to the overriding plate by accretionary processes during subduction. The geographical regions that eventually become accreted allochthonous terranes include island arcs, oceanic plateaus, submarine ridges, seamounts, continental fragments, and microcontinents. These future allochthonous terranes (FATs) contribute to continental crustal growth, subduction dynamics, and crustal recycling in the mantle. We present a review of modern FATs and their accreted counterparts based on available geological, seismic, and gravity studies and discuss their crustal structure, geological origin, and bulk crustal density. Island arcs have an average crustal thickness of 26 km, average bulk crustal density of 2.79 g cm −3 , and three distinct crustal units overlying a crust–mantle transition zone. Oceanic plateaus and submarine ridges have an average crustal thickness of 21 km and average bulk crustal density of 2.84 g cm −3 . Continental fragments presently on the ocean floor have an average crustal thickness of 25 km and bulk crustal density of 2.81 g cm −3 . Accreted allochthonous terranes can be compared to these crustal compilations to better understand which units of crust are accreted or subducted. In general, most accreted terranes are thin crustal units sheared off of FATs and added onto the accretionary prism, with thicknesses on the order of hundreds of meters to a few kilometers. However, many island arcs, oceanic plateaus, and submarine ridges were sheared off in the subduction interface and underplated onto the overlying continent. Other times we find evidence of terrane–continent collision leaving behind accreted terranes 25–40 km thick. We posit that rheologically weak crustal layers or shear zones that were formed when the FATs were produced can be activated as detachments during subduction, allowing parts of the FAT crust to accrete and others to subduct. In many modern FATs on the ocean floor, a sub-crustal layer of high seismic velocities, interpreted as ultramafic material, could serve as a detachment or delaminate during subduction.

Description: Polycyclic aromatic hydrocarbons in post-fire soils of drained peatlands in western Meshchera (Moscow region, Russia) Solid Earth, 5, 1305-1317, 2014 Author(s): A. Tsibart, A. Gennadiev, T. Koshovskii, and A. Watts Polycyclic aromatic hydrocarbons (PAHs) are priority pollutants that arrive in the environment from numerous anthropogenic and natural sources, but the data on their natural sources including wildfires remain insufficient. The level of contamination and the composition of PAHs in soils of the areas affected by wildfires were studied in this work. The study was conducted in the Moscow region (Russia) in areas occupied by drained peatland and strongly damaged by fires in 2002, 2010 and 2012. The features of PAH accumulation and the profile distributions in histosols and histic podzols after the fires of different times were analyzed. It was shown that new soil horizons formed after the fires – Cpir, Hpir and incipient O horizons – and that these horizons differ in PAH accumulation rate. Maximal total concentrations of 14 PAHs were detected in charred peat horizons Hpir (up to 330 ng g −1 ) and in post-fire incipient O horizons (up to 180 ng g −1 ), but the high-molecular-weight PAHs (benz(ghi)perylene, benz(a)pyrene, benz(k)fluoranthene) were revealed only in charry peat horizons. The trends of higher PAH concentrations were found in cases when smoldering combustion resulted in rather thick residual peat horizons. In cases of almost complete pyrogenic destruction of He horizons, total PAH concentrations were no more than 50 ng g −1 . Also, PAH accumulation in upper horizons of soils near the sites of the latest fires was observed.

Description: Characterization and interaction of driving factors in karst rocky desertification: a case study from Changshun, China Solid Earth, 5, 1329-1340, 2014 Author(s): E. Q. Xu and H. Q. Zhang As the most severe ecological issue in southwest China, karst rocky desertification (KRD) has both threatened and constrained regional sustainable development. Comprehensively understanding the relationship between the evolution of KRD and relevant driving data would provide more information to combat KRD in such complex karst environments. Past studies have been limited in quantifying the relative importance of driving factors influencing fine-scale KRD evolution, and have also lacked insight into their interactive impacts. To address these issues, we have used geographical information system techniques and a geographical detector model to explore the spatial consistency of driving factors and their interactions in relation to the evolution of KRD. Changshun County in China was selected as a representative area for the study. Nine relevant driving factors, including both natural and anthropogenic factors, were studied in regard to their relationships with KRD transformation between 2000 and 2010. Our results demonstrate the relative importance of driving data in influencing the improvement and deterioration of KRD. Lithology, soil type and road influence are identified as the leading factors. Interestingly, to our study at least, there is no significant difference between the impacts of natural and anthropogenic factors influencing KRD improvement, and even natural factors have a higher impact on KRD deterioration. Factors were found to enhance the influence of each other for KRD transformation. In particular, the results show a non-linearly enhanced effect between driving factors, which significantly aggravates KRD. New information found in our study helps to effectively control and restore areas afflicted by KRD.

Description: Microscale strain partitioning? Differential quartz crystallographic fabric development in Phyllite, Hindu Kush, Northwestern Pakistan Solid Earth, 5, 1319-1327, 2014 Author(s): K. P. Larson, J. L. Lamming, and S. Faisal Spatially referenced quartz c axis fabrics demonstrate the preservation of multiple, distinct fabrics in a specimen collected from northwestern Pakistan. The overall fabric yielded by the specimen is dominated by a single population of quartz grains, while the fabric signatures of two other unique, spatially distinct populations are overwhelmed. It is these minor fabrics, however, that provide information on temperature of deformation (403 ± 50 °C), differential stress (8.6 + 2.6/−1.5 MPa to 15.0 +3.8/−2.5 MPa), strain rate (10 −16 s −1 to 10 −15 s −1 ), and strain partitioning recorded by the specimen.

Description: Features of the Earth surface deformations in the Kamchatka peninsula and their relation to geoacoustic emission Solid Earth, 5, 1293-1300, 2014 Author(s): I. A. Larionov, Y. V. Marapulets, and B. M. Shevtsov The paper presents the results of investigations of deformation processes in the near-surface sedimentary rocks, which have been carried out in a seismically active region of the Kamchatka peninsula since 2007. The peculiarity of the experiments on registration of geodeformations is the application of a laser strainmeter–interferometer constructed according to the Michelson interferometer scheme. Besides rock deformations, geoacoustic emission in the frequency range from several hertz to the first tens of kilohertz is under investigation. Piezoceramic hydrophones installed in artificial water reservoirs are applied. It is shown that periods of primary rock compression and tension with a duration of up to several months are distinguished in the geodeformation process at the observation site. During the direction change in the deformations, when the geodeformation process rate grows, an increase in geoacoustic radiation is observed.

Description: Using open sidewalls for modelling self-consistent lithosphere subduction dynamics Solid Earth, 3, 313-326, 2012 Author(s): M. V. Chertova, T. Geenen, A. van den Berg, and W. Spakman Subduction modelling in regional model domains, in 2-D or 3-D, is commonly performed using closed (impermeable) vertical boundaries. Here we investigate the merits of using open boundaries for 2-D modelling of lithosphere subduction. Our experiments are focused on using open and closed (free slip) sidewalls while comparing results for two model aspect ratios of 3:1 and 6:1. Slab buoyancy driven subduction with open boundaries and free plates immediately develops into strong rollback with high trench retreat velocities and predominantly laminar asthenospheric flow. In contrast, free-slip sidewalls prove highly restrictive on subduction rollback evolution, unless the lithosphere plates are allowed to move away from the sidewalls. This initiates return flows pushing both plates toward the subduction zone speeding up subduction. Increasing the aspect ratio to 6:1 does not change the overall flow pattern when using open sidewalls but only the flow magnitude. In contrast, for free-slip boundaries, the slab evolution does change with respect to the 3:1 aspect ratio model and slab evolution does not resemble the evolution obtained with open boundaries using 6:1 aspect ratio. For models with open side boundaries, we could develop a flow-speed scaling based on energy dissipation arguments to convert between flow fields of different model aspect ratios. We have also investigated incorporating the effect of far-field generated lithosphere stress in our open boundary models. By applying realistic normal stress conditions to the strong part of the overriding plate at the sidewalls, we can transfer intraplate stress to influence subduction dynamics varying from slab roll-back, stationary subduction, to advancing subduction. The relative independence of the flow field on model aspect ratio allows for a smaller modelling domain. Open boundaries allow for subduction to evolve freely and avoid the adverse effects (e.g. forced return flows) of free-slip boundaries. We conclude that open boundaries in combination with intraplate stress conditions are to be preferred for modelling subduction evolution (rollback, stationary or advancing) using regional model domains.

Description: The 11 May 2011 earthquake at Lorca (SE Spain) viewed in a structural-tectonic context Solid Earth, 2, 199-204, 2011 Author(s): R. L. M. Vissers and B. M. L. Meijninger The Lorca earthquake of 11 May 2011 in the Betic Cordillera of SE Spain occurred almost exactly on the Alhama de Murcia fault, a marked fault that forms part of a NE-SW trending belt of faults and thrusts. The fault belt is reminiscent of a strike-slip corridor, but recent structural studies have provided clear evidence for reverse motions on these faults. Focal mechanisms of the main earthquake, but also of a foreshock, are strikingly consistent with structural observations on the Alhama de Murcia fault. This strengthens the conclusion that, rather than a strike-slip fault, the fault is at present a contractional fault with an oblique reverse sense of motion, presumably in response to the NW-directed motion of Africa with respect to Europe.

Description: The role of karst in engineering and environmental geosciences Solid Earth, 2, 155-158, 2011 Author(s): H. C. Ho Karst is a unique landform developed by soluble rock. It usually relates to the groundwater drainage system, and provides important water resources. Current researches indicate that karst is closely related to the Earth system and environmental protection, and it can also create potential natural hazards such as sinkhole flooding and land subsidence in urban area. Its relationship with hydrogeology has also been an important factor for studying water pollution and nutrient cycles in engineering geosciences and agricultural geology.

Description: Strength and permeability recovery of tuffisite-bearing andesite Solid Earth, 3, 191-198, 2012 Author(s): S. Kolzenburg, M. J. Heap, Y. Lavallée, J. K. Russell, P. G. Meredith, and D. B. Dingwell Tuffisites, the products of subsurface fragmentation, transport and deposition, are common in explosive volcanic environments. Their study provides direct insight to the mechanical processes operating within volcanic conduits. Here we document the influence of the presence of coherent tuffisite veins on the physical properties of andesitic rocks. We find that (1) compressive strength is unaffected by the presence and/or orientation of tuffisites, (2) permeability doubles when tuffisites are oriented favorably (at 45° to the fluid flow direction), and (3) ultrasonic wave velocities show a continuous increase with depth, independent of vein presence and orientation. Although the influence of tuffisites on andesitic rock properties determined here is modest, we emphasize that the material tested represents the post-eruptive state of tuffisite. Thus, these results likely delineate the upper and lower boundaries of strength vs. permeability and porosity, respectively. Our evidence suggests that, via compaction and lithification, tuffisites may restore the strength of the volcanic host-rocks to that of their pre-tuffisite values.

Description: Thermal structure and intermediate-depth seismicity in the Tohoku-Hokkaido subduction zones Solid Earth, 3, 355-364, 2012 Author(s): P. E. van Keken, S. Kita, and J. Nakajima The cause of intermediate-depth (〉40 km) seismicity in subduction zones is not well understood. The viability of proposed mechanisms, which include dehydration embrittlement, shear instabilities and the presence of fluids in general, depends significantly on local conditions, including pressure, temperature and composition. The well-instrumented and well-studied subduction zone below Northern Japan (Tohoku and Hokkaido) provides an excellent testing ground to study the conditions under which intermediate-depth seismicity occurs. This study combines new finite element models that predict the dynamics and thermal structure of the Japan subduction system with a high-precision hypocenter data base. The upper plane of seismicity is principally contained in the crustal portion of the subducting slab and appears to thin and deepen within the crust at depths 〉80 km. The disappearance of seismicity overlaps in most of the region with the predicted phase change of blueschist to hydrous eclogite, which forms a major dehydration front in the crust. The correlation between the thermally predicted blueschist-out boundary and the disappearance of seismicity breaks down in the transition from the northern Japan to Kurile arc below western Hokkaido. Adjusted models that take into account the seismically imaged modified upper mantle structure in this region fail to adequately recover the correlation that is seen below Tohoku and eastern Hokkaido. We conclude that the thermal structure below Western Hokkaido is significantly affected by time-dependent, 3-D dynamics of the slab. This study generally supports the role of fluids in the generation of intermediate-depth seismicity.

Description: Morphology and surface features of olivine in kimberlite: implications for ascent processes Solid Earth, 5, 313-326, 2014 Author(s): T. J. Jones, J. K. Russell, L. A. Porritt, and R. J. Brown Most kimberlite rocks contain large proportions of ellipsoidal-shaped xenocrystic olivine grains that are derived mainly from disaggregation of peridotite. Here, we describe the shapes, sizes and surfaces of olivine grains recovered from kimberlite lavas erupted from the Quaternary Igwisi Hills volcano, Tanzania. The Igwisi Hills kimberlitic olivine grains are compared to phenocrystic olivine, liberated from picritic lavas, and mantle olivine, liberated from a fresh peridotite xenolith. Image analysis, scanning electron microscopy imagery and laser microscopy reveal significant differences in the morphologies and surface features of the three crystal populations. The kimberlitic olivine grains form smooth, rounded to ellipsoidal shapes and have rough flaky micro-surfaces that are populated by impact pits. Mantle olivine grains are characterised by flaked surfaces and indented shapes consistent with growth as a crystal aggregate. Phenocrystic olivine exhibit faceted, smooth-surfaced crystal faces. We suggest that the unique shape and surface properties of the Igwisi Hills kimberlitic olivine grains are products of the transport processes attending kimberlite ascent from mantle source to surface. We infer that the unique shapes and surfaces of kimberlitic olivine grains result from three distinct mechanical processes attending their rapid transport through the thick cratonic mantle lithosphere: (1) penetrative flaking from micro-tensile failure induced by rapid decompression; (2) sustained abrasion and attrition driven by particle–particle collisions between grains within a turbulent, volatile-rich flow regime; and (3) higher-energy particle–particle collisions producing impact cavities superimposed on decompression structures. The combination of these processes during the rapid ascent of kimberlite magmas is responsible for the distinctive ellipsoidal shape of olivine xenocrysts found in kimberlites worldwide.

Description: Comparing a thermo-mechanical Weichselian Ice Sheet reconstruction to reconstructions based on the sea level equation: aspects of ice configurations and glacial isostatic adjustment Solid Earth, 5, 371-388, 2014 Author(s): P. Schmidt, B. Lund, J-O. Näslund, and J. Fastook In this study we compare a recent reconstruction of the Weichselian Ice Sheet as simulated by the University of Maine ice sheet model (UMISM) to two reconstructions commonly used in glacial isostatic adjustment (GIA) modelling: ICE-5G and ANU (Australian National University, also known as RSES). The UMISM reconstruction is carried out on a regional scale based on thermo-mechanical modelling, whereas ANU and ICE-5G are global models based on the sea level equation. The three models of the Weichselian Ice Sheet are compared directly in terms of ice volume, extent and thickness, as well as in terms of predicted glacial isostatic adjustment in Fennoscandia. The three reconstructions display significant differences. Whereas UMISM and ANU includes phases of pronounced advance and retreat prior to the last glacial maximum (LGM), the thickness and areal extent of the ICE-5G ice sheet is more or less constant up until the LGM. During the post-LGM deglaciation phase ANU and ICE-5G melt relatively uniformly over the entire ice sheet in contrast to UMISM, which melts preferentially from the edges, thus reflecting the fundamental difference in the reconstruction scheme. We find that all three reconstructions fit the present-day uplift rates over Fennoscandia equally well, albeit with different optimal earth model parameters. Given identical earth models, ICE-5G predicts the fastest present-day uplift rates, and ANU the slowest. Moreover, only for ANU can a unique best-fit model be determined. For UMISM and ICE-5G there is a range of earth models that can reproduce the present-day uplift rates equally well. This is understood from the higher present-day uplift rates predicted by ICE-5G and UMISM, which result in bifurcations in the best-fit upper- and lower-mantle viscosities. We study the areal distributions of present-day residual surface velocities in Fennoscandia and show that all three reconstructions generally over-predict velocities in southwestern Fennoscandia and that there are large differences in the fit to the observational data in Finland and northernmost Sweden and Norway. These difference may provide input to further enhancements of the ice sheet reconstructions.

Description: On the complexity of surface ruptures during normal faulting earthquakes: excerpts from the 6 April 2009 L'Aquila (central Italy) earthquake ( M w 6.3) Solid Earth, 5, 389-408, 2014 Author(s): L. Bonini, D. Di Bucci, G. Toscani, S. Seno, and G. Valensise Over the past few years the assessment of the earthquake potential of large continental faults has increasingly relied on field investigations. State-of-the-art seismic hazard models are progressively complementing the information derived from earthquake catalogs with geological observations of active faulting. Using these observations, however, requires full understanding of the relationships between seismogenic slip at depth and surface deformation, such that the evidence indicating the presence of a large, potentially seismogenic fault can be singled out effectively and unambiguously. We used observations and models of the 6 April 2009, M w 6.3, L'Aquila, normal faulting earthquake to explore the relationships between the activity of a large fault at seismogenic depth and its surface evidence. This very well-documented earthquake is representative of mid-size yet damaging earthquakes that are frequent around the Mediterranean basin, and was chosen as a paradigm of the nature of the associated geological evidence, along with observational difficulties and ambiguities. Thanks to the available high-resolution geologic, geodetic and seismological data aided by analog modeling, we reconstructed the full geometry of the seismogenic source in relation to surface and sub-surface faults. We maintain that the earthquake was caused by seismogenic slip in the range 3–10 km depth, and that the slip distribution was strongly controlled by inherited discontinuities. We also contend that faulting was expressed at the surface by pseudo-primary breaks resulting from coseismic crustal bending and by sympathetic slip on secondary faults. Based on our results we propose a scheme of normal fault hierarchization through which all surface occurrences related to faulting at various depths can be interpreted in the framework of a single, mechanically coherent model. We stress that appreciating such complexity is crucial to avoiding severe over- or under-estimation of the local seismogenic potential.

Description: AxiSEM: broadband 3-D seismic wavefields in axisymmetric media Solid Earth, 5, 425-445, 2014 Author(s): T. Nissen-Meyer, M. van Driel, S. C. Stähler, K. Hosseini, S. Hempel, L. Auer, A. Colombi, and A. Fournier We present a methodology to compute 3-D global seismic wavefields for realistic earthquake sources in visco-elastic anisotropic media, covering applications across the observable seismic frequency band with moderate computational resources. This is accommodated by mandating axisymmetric background models that allow for a multipole expansion such that only a 2-D computational domain is needed, whereas the azimuthal third dimension is computed analytically on the fly. This dimensional collapse opens doors for storing space–time wavefields on disk that can be used to compute Fréchet sensitivity kernels for waveform tomography. We use the corresponding publicly available AxiSEM ( www.axisem.info ) open-source spectral-element code, demonstrate its excellent scalability on supercomputers, a diverse range of applications ranging from normal modes to small-scale lowermost mantle structures, tomographic models, and comparison with observed data, and discuss further avenues to pursue with this methodology.