Search

The silicon isotopic composition ( δ Si 30 ) of the headwaters of the Ganges River, in the Himalaya, ranged from + 0.49 ± 0.01 ‰ to + 2.17 ± 0.04 ‰ at dissolved silicon (DSi) concentrations of 38 to 239 μM. Both the concentration and isotopic composition of DSi in the tributaries increased between the highest elevations to where the Ganges leaves the Himalayas at Rishikesh. The tributaries exhibit a linear correlation between δ Si 30 and DSi that may represent mixing between a low DSi, low δ Si 30 (e.g., 40 μM, + 0.5 ‰ ) component potentially reflecting fractionation during adsorption of a small fraction of silicon onto iron oxides and a high DSi, high δ Si 30 component (e.g., 240 μM, + 1.7 ‰ ) produced during higher intensity weathering with a greater proportional sequestration of weathered silicon into secondary minerals or biogenic silica. On the Ganges alluvial plain, in the Ganges and the Yamuna, Gomati, and their tributaries, DSi ranged from 122 to 218 μM while δ Si 30 ranged from + 1.03 ± 0.03 ‰ to + 2.46 ± 0.06 ‰ . Highest values of δ Si 30 occurred in the Gomati and its tributaries. In general, the lower DSi and higher δ Si 30 of DSi in these rivers suggests control of both by removal of DSi by secondary mineral formation and/or biogenic silica production. A simple 1-dimensional model with flow through a porous medium is introduced and provides a useful framework for understanding these results.

We reassess the applicability of the thermal plate cooling model to the subsidence of the North Pacific, Atlantic and North Indian Ocean Basins. We use a new numerical plate model in which the thermophysical parameters of the lithosphere vary with temperature according to the results of laboratory experiments, and the ridge temperature structure is consistent with the thickness of the oceanic crust. We first attempt to exclude thickened crust from our data set, and then to exclude swells and downwellings by masking regions of the data that remains that have significant gravity anomalies when there exists a clear regional correlation between intermediate-wavelength gravity and topography. We find that the average variation of depth with age is consistent with conventional half-space models until about 90 Myr. Thereafter, the departure from the half-space cooling curve is more rapid than predicted using simple conductive plate cooling models. The depth–age curves in the Pacific and Atlantic show ∼250 m of temporary shallowing between the ages of 90–130 Myr, a result consistent with the outcome of experiments on the initiation of small-scale boundary layer convection. The results do not change significantly if the estimated component of the gravity arising from plate cooling is subtracted prior to calculation of the correlation between gravity and topography. A 90-km-thick conductive plate is nevertheless a reasonable model for the average temperature structure of the oldest part of the Pacific ocean lithosphere. In the Pacific, the broad topographic undulations associated with the Line Island Swell, the Hawaiian Swell and surrounding basins have correlated gravity anomalies and an admittance of approximately 30 mGal km−1 and are likely to result from convective circulation in the upper mantle. In the Northeast Atlantic, the intermediate-wavelength admittance over the Cape Verde swell is similar; in the Northwest Atlantic over the Bermuda Swell it is slightly larger but not as well constrained.

We present evidence for the presence of complex anisotropy in the lowermost mantle from 3-D waveform modelling of observed core-diffracted shear waves that sample the southern edge of the African Large Low Shear Velocity Province (LLSVP). The anomalously strong amplitude of the SV component for the shear core-diffracted phase at large distances indicates the presence of anisotropy. We measure shear wave splitting parameters to determine which part of the elastic tensor is constrained by this particular data set. The modelling is performed using the spectral element method. The anisotropy is strong outside the LLSVP, weakens or rotates close to its boundary, and appears to be absent inside the LLSVP. The presence of the LLSVP margin may cause flow in the mantle to change direction. The occurrence of strong anisotropy in the region of fast seismic velocities is compatible with lattice-preferred orientation in post-perovskite due to accommodation of flow through dislocation creep.

Ongoing deglaciation in Iceland not only causes uplift at the surface but also increases magma production at depth due to decompression of the mantle. Here we study glacially induced decompression melting using 3-D models of glacial isostatic adjustment in Iceland since 1890. We find that the mean glacially induced pressure rate of change in the mantle increases melt production rates by 100–135%, or an additional 0.21–0.23 km3 of magma per year beneath Iceland. Approximately 50% of this melt is produced underneath central Iceland. The greatest volumetric increase is found directly beneath Iceland's largest ice cap, Vatnajökull, colocated with the most productive volcanoes. Our models of the effect of deglaciation on mantle melting predict a significantly larger volumetric response than previous models which only considered the effect of deglaciation of Vatnajökull, and only mantle melting directly below Vatnajökull. Although the ongoing deglaciation significantly increases the melt production rate, the increase in melt supply rate at the base of the lithosphere is delayed and depends on the melt ascent velocity through the mantle. Assuming that 25% of the melt reaches the surface, the upper limit on our deglaciation-induced melt estimates for central Iceland would be equivalent to an eruption the size of the 2010 Eyjafjallajökull summit eruption every seventh year.

Time-lapse, three-dimensional (3D) seismic surveys have imaged an accumulation of injected CO$_{2}$ adjacent to the Sleipner field in the North Sea basin. The changing pattern of reflectivity suggests that CO$_{2}$ accumulates within a series of interbedded sandstones and mudstones beneath a thick caprock of mudstone. Nine reflective horizons within the reservoir have been mapped on six surveys acquired between 1999 and 2008. These horizons have roughly elliptical planforms with eccentricities ranging between two and four. In the top half of the reservoir, horizon areas grow linearly with time. In the bottom half, horizon areas initially grow linearly for about eight years and then progressively shrink. The central portions of deeper reflective horizons dim with time. Amplitude analysis of horizons above, within, and below the reservoir show that this dimming is not solely caused by acoustic attenuation. Instead, it is partly attributable to CO$_{2}$ migration and/or CO$_{2}$ dissemination, which reduce the impedance contrast between sandstone and mudstone layers. Growth characteristics and permeability constraints suggest that each horizon grows by lateral spreading of a gravity current. This model is corroborated by the temporal pattern of horizon velocity pushdown beneath the reservoir. Horizon shrinkage may occur if the distal edge of a CO$_{2}$-filled layer penetrates the overlying mudstone, if the buoyant plume draws CO$_{2}$ upward, or if the effective permeability of deeper mudstone layers increases once interstitial brine has been expelled. Topographic control is evident at later times and produces elliptical planforms, especially toward the top of the reservoir. Our results show that quantitative mapping and analysis of time-lapse seismic surveys yield fluid dynamical insights which are testable, shedding light on the general problem of CO$_{2}$ sequestration.

Magma is transported in brittle rock through dikes and sills. This movement may be accompanied by the release of seismic energy that can be tracked from the Earth's surface. Locating dikes and deciphering their dynamics is therefore of prime importance in understanding and potentially forecasting volcanic eruptions. The Seismic Amplitude Ratio Analysis (SARA) method aims to track melt propagation using the amplitudes recorded across a seismic network without picking the arrival times of individual earthquake phases. This study validates this methodology by comparing SARA locations (filtered between 2 and 16 Hz) with the earthquake locations (same frequency band) recorded during the 2014–2015 Bár urn:x-wiley:jgrb:media:jgrb52508:jgrb52508-math-0003arbunga‐Holuhraun dike intrusion and eruption in Iceland. Integrating both approaches also provides the opportunity to investigate the spatiotemporal characteristics of magma migration during the dike intrusion and ensuing eruption. During the intrusion SARA locations correspond remarkably well to the locations of earthquakes. Several exceptions are, however, observed. (1) A low‐frequency signal was possibly associated with a subglacial eruption on 23 August. (2) A systematic retreat of the seismicity was also observed to the back of each active segment during stalled phases and was associated with a larger spatial extent of the seismic energy source. This behavior may be controlled by the dike's shape and/or by dike inflation. (3) During the eruption SARA locations consistently focused at the eruptive site. (4) Tremor‐rich signal close to ice cauldrons occurred on 3 September. This study demonstrates the power of the SARA methodology, provided robust site amplification; Quality Factors and seismic velocities are available.

The mid-Atlantic plate boundary breaks up into a series of segments across Iceland. The South Iceland Seismic Zone (SISZ) is a complex transform zone where left-lateral E-W shear between the Reykjanes Peninsula Rift Zone and the Eastern Volcanic Zone is accommodated by bookshelf faulting along N-S lateral strike-slip faults. The SISZ is also a transient feature, migrating sideways in response to the southward propagation of the Eastern Volcanic Zone. Sequences of large earthquakes (M > 6) lasting from days to years and affecting most of the seismic zone have occurred repeatedly in historical time (last 1100 years), separated by intervals of relative quiescence lasting decades to more than a century. On May 29 2008, a Mw 6.1 earthquake struck the western part of the South Iceland Seismic Zone, followed within seconds by a slightly smaller event on a second fault ~5 km further west. Aftershocks, detected by a temporal array of 11 seismometers and three permanent Icelandic Meteorological Office stations were located using an automated Coalescence Microseismic Mapping technique. The epicenters delineate two major and several smaller N-S faults as well as an E-W zone of activity stretching further west into the Reykjanes Peninsula Rift Zone. Fault plane solutions show both right lateral and oblique strike slip mechanisms along the two major N-S faults. The aftershocks deepen from 3-5 km in the north to 8-9 km in the south, suggesting that the main faults dip southwards. The faulting is interpreted to be driven by the local stress due to transform motion between two parallel segments of the divergent plate boundary crossing Iceland.

We study the seismic response of layered basalts in the Faroe Islands using borehole data and\ud vertical seismic profiles from the Vestmanna borehole, combined with reflection and wideangle\ud seismic data recorded into arrays of both borehole and land multicomponent receivers.\ud Imaging through the basalt cover in the Faroe-Shetland Basin is a challenge for conventional\ud seismic surveys: scattering caused by the high reflectivity of the basalt as well as intra-basalt\ud multiples and high attenuation from the layered sequence make it difficult to image within\ud and beneath the basalts. This project allows us to correlate ultrasonic-scale velocity and\ud density measurements from the borehole together with ground-truthing from borehole logs\ud and core samples with the seismic-scale velocities and reflection images derived from VSP\ud and surface data. We find a good match of observed travel-times of borehole and wide-angle\ud P-wave data with those predicted from the borehole measurements, suggesting lateral\ud homogeneity over horizontal distances on the kilometre scale, and restricted transverse\ud anisotropy of the layered basalts. A pronounced intra-basalt reflector identified on the\ud multichannel surface seismic can be correlated with lithostratigraphic interpretation of the\ud borehole logs as caused by thick flows near the top of the Lower Basalt formation.

One-dimensional seismic reference models are known to be a nonunique solution to global seismic data, hampering an interpretation in terms of physical structure. Here we test the compatibility of the simplest hypothesis of a mantle convecting as a whole, with a constant pyrolitic composition with phase transitions, directly against the kind of seismic data that went into global seismic reference models, focusing on upper mantle structure down to 800 km depth. By randomly varying the elastic and anelastic parameters of the main mantle minerals within their uncertainty bounds, we generate a set of 100,000 adiabatic, pyrolitic models. A small number of these models ( 60) that is as satisfactory as the fit of preliminary reference Earth model (PREM) or AK135(-F). Although the accepted models have widely different combinations of the mineral parameters, there is a preference for a relatively high olivine shear modulus and its pressure and temperature derivatives, relatively low wadsleyite bulk and shear parameters, and relatively high ringwoodite bulk modulus derivatives. The resulting seismic profiles are very similar and, compared to PREM or AK135, have lower velocities above 400 km, larger jumps near "410," lower transition zone gradients, lower jumps around "660," and stronger gradients directly below. Such physical models that fit seismic data well enough to be useful as a seismic reference model can significantly facilitate physical interpretation of seismic structures. Copyright 2005 by the American Geophysical Union.

Migration of CO 2 through storage reservoirs can be monitored using time-lapse seismic reflection surveys. At the Sleipner Field, injected CO 2 is distributed throughout nine layers within the reservoir. These layers are too thin to be seismically resolvable by direct measurement of the separation between reflections from the top and bottom of each layer. Here, we develop and apply an inverse method for measuring thickness changes of the shallowest layer. Our approach combines differences in travel time down to \ud a specific reflection together with amplitude measurements to determine layer thicknesses from time-lapse surveys. A series of synthetic forward models were used to test the robustness of our inverse approach and to quantify uncer- \ud tainties. In the absence of ambient noise, this approach can unambiguously resolve layer thickness. If a realistic ambient noise distribution is included, layer thicknesses of 1–6 m are accurately retrieved with an uncertainty of ±0.5 m. We used this approach to generate a thickness map of the shallowest layer. The fidelity of this result was tested using measurements of layer thickness determined from the 2010 broadband seismic survey. The calculated volume of CO 2 within the shallowest layer increases at a rate that is quadratic in time, despite an approximately constant injection rate into the base of the reservoir. This result is consistent with a diminished growth rate of the areal extent of underlying layers. Finally, the relationship between caprock topography and layer thickness is explored and potential migration pathways that \ud charge this layer are identified.

Lithological variations in the mantle source regions under mid-ocean ridges and ocean islands have been proposed to play a key role in controlling melt generation and basalt composition. Here we combine compositional observations from Icelandic basalts and modelling of melting of a bi-lithologic peridotite-pyroxenite mantle to demonstrate that, while short-lengthscale major element variation is present in the mantle under Iceland, source heterogeneity does not make an important contribution to excess melt production. By identifying the major element characteristics of endmember Icelandic melts, we find enriched melts to be characterised by low SiO2 and CaO, but high FeO. We quantitatively compare endmember compositions to experimental partial melts generated from a range of lithologies, pressures and melt fractions. This comparison indicates that a single source composition cannot account for all the major element variation; depleted Icelandic melts can be produced by depleted peridotite melting, but the major element composition of enriched melts is best matched by melting of mantle sources that have been refertilised by the addition of up to 40% mid-ocean ridge basalt. The enriched source beneath Iceland is more fusible than the source of depleted melts, and as such will be over-represented in accumulated melts compared with its abundance in the source. Modelling of peridotite-pyroxenite melting, combined with our observational constraints on the composition of the Icelandic mantle, indicates that crustal thickness variations in the North Atlantic must be primarily due to mantle temperature and flow field variations.

An inverse model has been developed to determine the magnitude of denudation at\ud seabed and subsurface unconformities by using root-mean-square (RMS) stacking\ud velocity data derived from processing a set of seismic reflection profiles. This approach\ud provides superior spatial coverage in comparison to other methods, such as vitrinite\ud reflectance, apatite fission track, and sonic velocity modeling, which are restricted\ud to borehole locations. The model assumes exponential porosity decay with depth and a\ud standard velocity-porosity relationship in order to compute a synthetic RMS velocity\ud profile. Denudation values at two levels in the stratigraphy are then adjusted until the fit\ud between the model and the data is optimized. Successful modeling is dependent upon\ud independent estimates of the initial porosity of sediment since significant trade-off occurs\ud between initial porosity and denudation. Application to the west African shelf shows that\ud 0.5–1 km of denudation occurred along the entire margin, probably during late Neogene\ud times. The amount of denudation decreases oceanward and was probably triggered by\ud regional tilting associated with initiation and/or regeneration of continent-wide\ud mantle convective upwelling, which is thought to have affected much of subequatorial\ud Africa. A subsurface Oligocene unconformity represents as much as 2.5 km of denudation\ud and was probably produced by initiation of an oceanic current.

The structure of the Earth's inner core is not well known between depths of ∼100–200 km beneath the inner core boundary. This is a result of the PKP core phase triplication and the existence of strong precursors to PKP phases, which hinder the measurement of inner core compressional PKIKP waves at epicentral distances between roughly 143 and 148°. Consequently, interpretation of the detailed structure of deeper regions also remains difficult. To overcome these issues we stack seismograms in slowness and time, separating the PKP and PKIKP phases which arrive simultaneously but with different slowness. We apply this method to study the inner core's Western hemisphere beneath South and Central America using paths travelling in the quasi-polar direction between 140 and 150° epicentral distance, which enables us to measure PKiKP–PKIKP differential traveltimes up to greater epicentral distance than has previously been done. The resulting PKiKP–PKIKP differential traveltime residuals increase with epicentral distance, which indicates a marked increase in seismic velocity for polar paths at depths greater than 100 km compared to reference model AK135. Assuming a homogeneous outer core, these findings can be explained by either (i) inner core heterogeneity due to an increase in isotropic velocity or (ii) increase in anisotropy over the studied depth range. Although this study only samples a small region of the inner core and the current data cannot distinguish between the two alternatives, we prefer the latter interpretation in the light of previous work.

The active magmatic processes beneath volcanoes in continental rifts is poorly understood. For example, until recently in the East African rift (EAR), the majority of the young volcanoes were thought to be inactive. More recent studies have shown that numerous volcanoes in the EAR are seismically active and deforming rapidly. However, an unambiguous sign of actively degassing magma hosted in shallow magma bodies has eluded most investigators. Here we present detailed analysis of the first low-frequency (LF) earthquake swarms to be observed in the Main Ethiopian Rift. The earthquakes locate to beneath Tullu Moye volcano and are directly related to the presence of a shallow magma body with a high fluid content. Using spectral modelling we show that the LF earthquakes appear to have low stress-drops (1–50 kPa) which we interpret in terms of low rupture velocities and high pore-fluid pressure. Careful relocation of the LF earthquakes place them approximately 4 km below the surface within one of two possible clusters. However, analysis of the correlation between earthquake waveforms show that each swarm contains a range of earthquake families and as such a diversity of earthquake source mechanisms. To explain these observations, we propose the seismicity is induced by H2O/CO2 fluid pulses from the shallow magma body into a highly fractured region. Fluid pulses cause high pore fluid pressures, which also cause the low rupture velocities.

[1] Despite its importance, the spatial and temporal pattern of dynamic topography generated by mantle convective circulation is poorly known. We present accurate estimates of dynamic topography from oceanic basins and continental margins surrounding Australia. Our starting point is measurement of residual depth anomalies on the oldest oceanic floor adjacent to the continental shelf. These anomalies were determined from a combined dataset of ~200 seismic reflection and wide-angle images of well-sedimented oceanic crust. They have amplitudes of between −1 km and +0.5 km, and their spatial variation is broadly consistent with long-wavelength free-air gravity and shallow seismic tomographic anomalies. Along the Northwest Shelf, a regional depth anomaly of −300 to −700 m intersects the adjacent continental shelf. The temporal evolution of this anomaly was determined by analyzing the stratigraphic architecture of an extensive carbonate platform, which fringes the shelf and records a dramatic switch from progradation to aggradation during Neogene times. Three-dimensional seismic mapping calibrated by boreholes was used to calculate water-loaded subsidence histories at rollover points of clinoforms along the shelf. At 9 ± 3 Ma, the rate of subsidence increases from 5 to up 75 m Myr−1, generating a subsidence anomaly of −300 to −700 m. The amplitude of this anomaly varies along the shelf and cannot be generated by glacio-eustatic sea-level variation. Instead, we propose that a regional subsidence episode, which affects both the proximal shelf and the distal oceanic basin, was generated by convective drawdown. By combining our results with other published estimates of uplift and subsidence, a map of Australia, which shows the spatial and temporal pattern of dynamic topography is presented. Most, but not all, of Australia's epeirogeny can be attributed to rapid northward motion of the Australian plate over a pre-existing pattern of convective circulation.

In situ cosmogenic 10Be and 26Al measurements on sedimentary quartzites, together with analysis of abandoned river terraces, can be used to show how drainage responded to the lateral propagation of a late Quaternary anticline forming above a blind reverse fault in Central Otago, New Zealand. A close link between the progression of cosmogenic ages and the tectonic geomorphology allows us to confirm that uplift and propagation rates on the anticline are in the ranges 0.08–0.12 mm yr−1 and 0.8–1.5 mm yr−1, respectively, over the last 550,000 years. The agreement between the isotopic and geomorphological evidence in turn requires that minimum 10Be exposure ages as great as 660 ka are not in steady state with respect to erosion. This is a geochemical result of significance, as it opens the possibility of a more regional analysis of landscape evolution in this region, using the same rocks. On an adjacent anticline, thought to be even older on geomorphological arguments, 10Be concentrations give minimum exposure ages in the range 750–1400 ka. These extremely old minimum ages are rare worldwide, as they are usually limited to smaller values by erosion. They are attributable to the very resistant nature of the quartz-rich boulders in which the measurements were made, which in places are shown to have erosion rates lower than 0.5 mm kyr−1.

[1] Longitudinal river profiles, where elevation of a river bed is plotted as a function of distance along the river bed, contain information about uplift rate. When a region adjacent to a reference level (e.g., sea level) is uplifted, a rapid change in gradient occurs near the river mouth. The erosional process causes this change in gradient to migrate upstream. Thus a river profile is effectively a ‘tape recording’ of the uplift rate history, provided that the erosional process can be adequately parameterized. Here, we use a non-linear equation to relate the shape of a river profile, z(x), to uplift rate history, U(t). If erosion is assumed to be dominated by knickpoint retreat, an inverse model can be formulated and used to calculate uplift rate histories. Our model builds upon standard stream profile analysis, which focuses on the relationship between profile slope and drainage area. We have applied this analytical approach to river profiles from the Bie Dome, Angola. Calculated uplift rate histories agree with independent geologic estimates.

The extent to which persistent, rather than transient, fissures (wide planar voids) can exist along upper crustal faults is important in assessing fault permeability to mineral and hydrocarbon-bearing fluids. Variscan (late Carboniferous) faults cutting Dinantian (Lower Carboniferous) limestones on the Gower peninsula, South Wales, host clear evidence for fissures up to several metres wide. Evidence includes dendritic hematite growth and elongate calcite growth into open voids, spar ball and cockade breccia formation, laminated sediment infill and void-collapse breccias. Detailed mapping reveals cross-cutting geometries and brecciation of earlier fissure fills, showing that fissures were formed during, rather than after, active faulting. Fissures therefore probably formed by geometric mismatch between displaced fault walls, rather than by solution widening along inactive faults.

We present a suite of programs that implement decades-old algorithms for computation of seismic surface wave reflection and transmission coefficients at a welded contact between two laterally homogeneous quarter-spaces. For Love as well as Rayleigh waves, the algorithms are shown to be capable of modelling multiple mode conversions at a lateral discontinuity, which was not shown in the original publications or in the subsequent literature. Only normal incidence at a lateral boundary is considered so there is no Love–Rayleigh coupling, but incidence of any mode and coupling to any (other) mode can be handled. The code is written in Python and makes use of SciPy's Simpson's rule integrator and NumPy's linear algebra solver for its core functionality. Transmission-side results from this code are found to be in good agreement with those from finite-difference simulations. In today's research environment of extensive computing power, the coded algorithms are arguably redundant but SWRT can be used as a valuable testing tool for the ever evolving numerical solvers of seismic wave propagation. SWRT is available via GitHub (https://github.com/arjundatta23/SWRT.git).

We used seismic body waves, radar interferometry and field investigation to examine the source processes of the destructive earthquake of 2005 February 22 near Zarand, in south–central Iran. The earthquake ruptured an intramountain reverse fault, striking E–W and dipping north at ∼60° to a depth of about 10 km. It produced a series of coseismic scarps with up to 1 m vertical displacement over a total distance of ∼13 km, continuous for 7 km. The line of the coseismic ruptures followed a known geological fault of unknown, but probably pre-Late Cenozoic, age and involved bedding-plane slip where the scarps were continuous at the surface. However, any signs of earlier coseismic ruptures along this fault had been obliterated by the time of the 2005 earthquake, probably by land sliding and weathering, so that the fault could not reasonably have been identified as active beforehand. The 2005 fault is at an oblique angle to the range-bounding Kuh Banan strike-slip fault, and may represent a splay from that fault, related to its southern termination. Other intramountain reverse faulting earthquakes have occurred in Iran, but this is the first to have produced a clear, mapped surface rupture, and to have been studied with InSAR. Faults of this type represent a serious seismic hazard in Iran and are difficult to assess, because their geomorphological expression is much less clear than the range-bounding reverse faults, which are more common and have been better studied.

The Mw 6.1 2016 Parina earthquake led to extension of the south Peruvian Andes along a normal fault with evidence of Holocene slip. We use interferometric synthetic aperture radar, seismology, and field mapping to determine a source model for this event and show that extension at Parina is oriented NE-SW, which is parallel to the shortening direction in the adjacent sub-Andean lowlands. In addition, we use earthquake source models and GPS data to demonstrate that shortening within the sub-Andes is parallel to topographic gradients. Both observations imply that forces resulting from spatial variations in gravitational potential energy are important in controlling the geometry of the deformation in the Andes. We calculate the horizontal forces per unit length acting between the Andes and South America due to these potential energy contrasts to be 4–8 ×1012 N/m along strike of the mountain range. Normal faulting at Parina implies that the Andes in south Peru have reached the maximum elevation that can be supported by the forces transmitted across the adjacent foreland, which requires that the foreland faults have an effective coefficient of friction ≤ 0.2. Additionally, the onset of extension in parts of the central Andes following orogen-wide compression in the late Miocene suggests that there has been a change in the force balance within the mountains. We propose that shortening on weak detachment faults within the Andean foreland since ∼5–9 Ma reduced the shear tractions acting along the base of the upper crust in the eastern Andes, leading to extension in the highest parts of the range. ©2018. American Geophysical Union. All Rights Reserved.

Buoyant plumes of CO2 spreading through water-saturated permeable rock, bounded by layers of lower permeability, tend to spread laterally. As they advance, they may gradually leak through fractures or discontinuities in the lower permeability boundary, leading to a gradual waning of the original plume and dispersal of CO2 higher in the formation. With a finite release of CO2, the trailing edge of the plume recedes with time, and capillary forces tend to trap a fraction of this CO2 within the pore space. This also leads to a gradual waning of the plume with time and limits the mass of CO2 which can leak through the boundary and rise higher into the formation. We explore the balance between these two effects and calculate some of the controls on the fraction of the CO2 plume which becomes trapped within the original layer of rock.

The effects of dispersion acting on gravity currents propagating through porous media are considered theoretically and experimentally. We exploit the large aspect ratio of these currents to formulate a depth-averaged model of the evolution of the mass and buoyancy. Dispersion, acting predominantly at the interface between the current and the ambient, is velocity dependent and acts to entrain fluid into the gravity current, in direct analogy to turbulent mixing. Here, we show that when the gravity current is fed by a constant buoyancy and mass flux the buoyancy of the current is self-similar and recovers the classical solution for gravity currents in porous media. In contrast, the profile and the depth-averaged concentration of the current evolve in a non-self-similar manner. The total volume of the current increases with time as $t^{1/3}$ due to this dispersive entrainment. We test our theoretical predictions using a suite of laboratory experiments in which the evolution of the concentration within the current was mapped using a dye-attenuation technique. These experimental results show good agreement with the early-time limits of our theoretical model, and in particular accurately predict the evolution of the depth-averaged concentration profile. These results suggest that mixing within porous media may be modelled using an effective dispersive entrainment, the magnitude of which may be set by the underlying structure of the porous medium.

The coherence between Bouguer gravity anomalies and topography is widely used to estimate the value of $$T_{e}$$ , the effective elastic thickness of the lithosphere. In areas where there is little topography but substantial free air gravity anomalies there is often little coherence between the free air anomalies and topography. In such regions the Bouguer coherence method generally gives estimates of $$T_{e}$$ of 90 km or more. A detailed analysis shows that, under these conditions, the value of the Bouguer coherence $$\gamma ^2_b$$ is entirely controlled by the ratio of the power spectra of the free air gravity anomalies and the uncompensated topography, and contains no information about the value of $$T_{e}$$ . What is worse, under these circumstances the variation of $$\gamma ^2_b$$ with wavelength closely resembles that expected for large values of $$T_{e}$$ . These results show that neither the Bouguer coherence method nor the admittance method can produce meaningful estimates of $$T_{e}$$ when the free air gravity anomalies are incoherent with the topography.

Although agriculturally accelerated soil erosion is implicated in the unsustainable environmental degradation of mountain environments, such as in the Himalaya, the effects of land use can be challenging to quantify in many mountain settings because of the high and variable natural background rates of erosion. In this study, we present new long-term denudation rates, derived from cosmogenic 10Be analysis of quartz in river sediment from the Likhu Khola, a small agricultural river basin in the Middle Hills of central Nepal. Calculated long-term denudation rates, which reflect background natural erosion processes over 1000+ years prior to agricultural intensification, are similar to present-day sediment yields and to soil loss rates from terraces that are well maintained. Similarity in short- and long-term catchment-wide erosion rates for the Likhu is consistent with data from elsewhere in the Nepal Middle Hills but contrasts with the very large increases in short-term erosion rates seen in agricultural catchments in other steep mountain settings. Our results suggest that the large sediment fluxes exported from the Likhu and other Middle Hills rivers in the Himalaya are derived in large part from natural processes, rather than from soil erosion as a result of agricultural activity. Catchment-scale erosional fluxes may be similar over short and long timescales if both are dominated by mass wasting sources such as gullies, landslides, and debris flows (e.g., as is evident in the landslide-dominated Khudi Khola of the Nepal High Himalaya, based on compiled data). As a consequence, simple comparison of catchment-scale fluxes will not necessarily pinpoint land use effects on soils where these are only a small part of the total erosion budget, unless rates of mass wasting are also considered. Estimates of the mass wasting contribution to erosion in the Likhu imply catchment-averaged soil production rates on the order of ~ 0.25–0.35 mm yr−1, though rates of mass wasting are poorly constrained. The deficit between our best estimates for soil production rates and measurements of soil loss rates supports conclusions from previous studies that terraced agriculture in the Likhu may not be associated with a large systematic soil deficit, at least when terraces are well maintained, but that poorly managed terraces, forest, and scrubland may lead to rapid depletion of soil resources.

Seismic noise interferometry allows the continuous and real-time measurement of relative seismic velocity through a volcanic edifice. Because seismic velocity is sensitive to the pressurization state of the system, this method is an exciting new monitoring tool at active volcanoes. Despite the potential of this tool, no studies have yet comprehensively compared velocity to other geophysical observables on a short-term time scale at a volcano over a significant length of time. We use volcanic tremor (~0.3 to 1.0 Hz) at Kīlauea as a passive source for interferometry to measure relative velocity changes with time. By cross-correlating the vertical component of day-long seismic records between ~230 station pairs, we extract coherent and temporally consistent coda wave signals with time lags of up to 120 s. Our resulting time series of relative velocity shows a remarkable correlation between relative velocity and the radial tilt record measured at Kīlauea summit, consistently correlating on a time scale of days to weeks for almost the entire study period (June 2011 to November 2015). As the summit continually deforms in deflation-inflation events, the velocity decreases and increases, respectively. Modeling of strain at Kīlauea suggests that, during inflation of the shallow magma reservoir (1 to 2 km below the surface), most of the edifice is dominated by compression—hence closing cracks and producing faster velocities—and vice versa. The excellent correlation between relative velocity and deformation in this study provides an opportunity to understand better the mechanisms causing seismic velocity changes at volcanoes, and therefore realize the potential of passive interferometry as a monitoring tool.

Hu, X-F, Jeans, C.V. and Dickson, J.A.D. 2012. Geochemical and stable isotope patterns of calcite cementation in the Upper Cretaceous Chalk, UK: Direct evidence from calcite-filled vugs in brachiopods. Acta GeologicaPolonica, 62 (2), 143-172. Warszawa. The history of research into the cementation of the Upper Cretaceous Chalk of the UK is reviewed. Calcitefilled vugs within the shell cavities of terebratulid brachiopods from the Cenomanian Chalk of eastern England have been investigated by cathodoluminesence imaging, staining, electron microprobe and stable isotope analysis. This has provided the first detailed analysis of the geochemistry of the Chalk’s cement. Two cement series, suboxic and anoxic, are recognized. Both start with a Mg-rich calcite with positive δ 13 C values considered to have been precipitated under oxic conditions influenced by aerobic ammonification. The suboxic series is characterized by positive δ 13 C values that became increasingly so as cementation progressed, reaching values of 3.5‰. Manganese is the dominant trace element in the earlier cement, iron in the later cement. Mnand Fe-reducing microbes influenced cement precipitation and the trace element and δ 13 C patterns. The anoxic series is characterized by δ 13 C values that became increasingly negative as cementation progressed, reaching values of -6.5‰. Trace elements are dominated by iron and manganese. Sulphate-reducing microbes influenced cement precipitation and the trace element and δ 13 C patterns. Both cement series are related closely to lithofacies and early lithification pre-dating the regional hardening of the Chalk. The suboxic series occurs in chalk which was continuously deposited and contained hematite pigment and limited organic matter. The anoxic series was associated with slow to nil deposition and hardground development in chalks that originally contained hematite pigment but no longer do so, and an enhanced supply of organic matter.

The source parameters and slip distribution of the 2013 May 11 Mw 6.1 Minab earthquake are studied using seismology, geodesy and field observations. We observe left-lateral strike-slip motion on a fault striking ENE–WSW; approximately perpendicular to previously studied faults in the Minab–Zendan–Palami fault zone. The fault that ruptured in 2013 is one of a series of ∼E–W striking left-lateral faults visible in the geology and geomorphology. These accommodate a velocity field equivalent to right-lateral shear on ∼N–S striking planes by clockwise rotations about vertical axes. The presence of these faults can reconcile differences in estimates of fault slip rates in the western Makran from GPS and Quaternary dating. The longitudinal range of shear in the western Makran is likely to be controlled by the distance over which the underthrusting Arabian lithosphere deepens in the transition from continent–continent collision in the Zagros to oceanic subduction in the Makran.

Pagination differs from hardbound copy deposited in Cambridge University Library., Reactions between CO2-charged brines and reservoir minerals might either enhance the long-term storage of CO2 in geological reservoirs or facilitate leakage by corroding cap rocks and fault seals. Modelling the progress of such reactions is frustrated by uncertainties in the absolute mineral surface reaction rates and the significance of other rate limiting steps in natural systems. This study uses the chemical evolution of groundwater from the Jurassic Navajo Sandstone, part of a leaking natural accumulation of CO2 at Green River, Utah, in the Colorado Plateau, USA, to place constraints on the rates and potential controlling mechanisms of the mineral-fluid reactions,under elevated CO2 pressures, in a natural system. The progress of individual reactions, inferred from changes in groundwater chemistry is modelled using mass balance techniques. The mineral reactions are close to stoichiometric with plagioclase and K-feldspar dissolution largely balanced by precipitation of clay minerals and carbonate. Mineral modes, in conjunction with published surface area measurements and flow rates estimated from hydraulic head measurements, are then used to quantify the kinetics of feldspar dissolution. Maximum estimated dissolution rates for plagioclase and K-feldspar are 2x10-14 and 4x10-16 mol·m-2·s-1, respectively. Fluid ion-activity products are close to equilibrium (e.g. DGr for plagioclase between -2 and -10 kJ/mol) and lie in the region in which mineral surface reaction rates show a strong dependence on DGr. Local variation in DGr is attributed to the injection and disassociation of CO2 which initially depresses silicate mineral saturation in the fluid, promoting feldspar dissolution. With progressive flow through the aquifer, feldspar hydrolysis reactions consume H+ and liberate solutes to solution which increase mineral saturation in the fluid and rates slow as a consequence. The measured plagioclase dissolution rates at low DGr would be compatible with far-from-quilibrium rates of ~1x10-13 mol·m-2·s-1 as observed in some experimental studies. This suggests that the discrepancy between field and laboratory reaction rates may in part be explained by the differences in the thermodynamic state of natural and experimental fluids, with field-scale reactions occurring close to equilibrium whereas most laboratory experiments are run far-from-equilibrium. Surface carbonate deposits and cementation within the footwall of the local fault systems record multiple injections of CO2 into the Navajo Aquifer and leakage of CO2 from the site over ca. 400,000 years. The d18O, d13C and 87Sr/86Sr of these deposits record rapid rates of CO2 leakage (up to ~1000 tonnes/a) following injection of CO2, but rates differ by an order of magnitude between each fault, due to differences in the fault architecture. Elevated pCO2 enhances rates of feldspar dissolution in the host aquifer and carbonate precipitation in fracture conduits. Silicate mineral dissolution rates decline and carbonate precipitation rates increase as pH and the CO2 charge dissipate. The Sr/Ca of calcite cements record average precipitation rates of ~2x10-6 mol/m2/s, comparable to laboratory derived calcite precipitation rates in fluids with elevated Mn/Ca and Fe/Ca, at cc of ~1 to 3. This suggests that far-from-equilibrium carbonate precipitation, which blocks fracture conduits and causes the leaking system to self-seal, driven by CO2 degassing in the shallow subsurface, can be accurately modeled with laboratory derived rates. Sandstones altered in CO2 leakage conduits exhibit extensive dissolution of hematite grain coatings and are chemically bleached as a result. Measurements of Eh-pH conditions in the modern fluid, and modeling of paleo-Eh-pH conditions using calcite Fe and Mn concentrations, suggests that the CO2-charged groundwaters are reducing, due to their low dissolved O2 content and that pH suppression due to high pCO2 is capable of dissolving and transporting large concentrations of metals. Exhumed paleo-CO2 reservoirs along the crest of the Green River anticline have been identified using volatile hosting fluid inclusions. Paleo-CO2-charged fluids mobilized hydrocarbons and CH4 from deeper formations, enhancing the reductive dissolution of hematite, which produced spectacular km-scale bleached patterns in these sediment., This work was funded by a post-graduate research grant from Shell Global Solutions International awarded to Niko Kampman.

The competition between wind and buoyancy forces acting on a naturally ventilated\ud building with one upwind and one downwind opening and a distributed heat source\ud can lead to multiple steady states. This study considers the impact of a second\ud downwind opening on the existence of these multiple steady states. It is found that if\ud this new opening exceeds a critical area, then the multiple steady states are eliminated,\ud and the transition from buoyancy to wind dominated flow is continuous. This critical\ud area is shown to be a function of the relative heights of the three windows. The\ud predictions of the model are in accord with new laboratory experiments, and the\ud implications for building design are discussed.

SUMMARY We use improved focal mechanisms and centroid depth estimates of earthquakes, combined with GPS velocities, to examine the tectonics of the Hellenic subduction zone, and in particular the processes occurring at both ends of the Hellenic Arc. Nubia-Aegean convergence is accommodated by shallowly dipping thrust-faulting along the subduction–zone interface, as well as by steeper splay faults in the overriding material. From a comparison of observed and expected seismic moment release over the last 100 yr, combined with existing knowledge of the longer-term documented historical record, we confirm earlier suggestions that most (80 per cent) of this convergence is accommodated aseismically, that is, that the subduction zone is uncoupled. This conclusion is robust, even allowing for rare very large earthquakes on splay faults, such as that of AD 365, and also allowing for the contribution of small earthquakes. The downgoing Nubian plate deforms by arc-parallel contraction at all depths, from 200 km seaward of Crete to at least 100 km within the subducting slab. Extensional (T) axes of earthquakes are aligned downdip within the descending slab suggesting that, even if the aseismic prolongation of the slab has reached the 670 km mantle discontinuity, it does not transmit stresses to shallower depths. Shallow thrust-faulting earthquakes on the subduction interface show a divergence of slip vectors round the arc, and GPS measurements show that this is accommodated mainly by E–W extension on normal faults in the overriding Aegean material. The eastern end of the subduction zone, south of Rhodes, displays distributed deformation in the overriding material, including a mixture of strike-slip and splay-thrust faulting, and probably involves rotations about a vertical axes. Here slip on the interface itself is by thrust faulting with slip vectors oblique to the arc but parallel to the overall Nubia-Aegean convergence: there is no evidence for slip-partitioning in the traditional sense. In the west, the subduction zone terminates in a distributed zone of parallel NE–SW strike-slip faults, of which the most prominent is the Kefalonia Transform Fault (KTF). A flexural gravity anomaly confirms that the deep bathymetric escarpment of the KTF is a lateral ramp, formed as the Ionian islands are emplaced SW onto the Apulian lithosphere, and enhanced by minor thrust faulting with slip vectors perpendicular to the scarp. Distributed parallel strike-slip faults both SW and NE of mainland central Greece terminate in E–W graben in central Greece, which accommodate the overall NE–SW shear by clockwise block rotation. Central Greece therefore acts as a relay zone between the strike-slip faulting of the NE Aegean and the Ionian Islands–western Peloponnese.

\ud \ud The thermochemical structure of lithospheric and asthenospheric mantle exert primary controls on surface topography and volcanic activity. Volcanic rock compositions and mantle seismic velocities provide indirect observations of this structure. Here, we compile and analyze a global database of the distribution and composition of Neogene-Quaternary intraplate volcanic rocks. By integrating this database with seismic tomographic models, we show that intraplate volcanism is concentrated in regions characterized by slow upper mantle shear-wave velocities and by thin lithosphere (i.e.

Factors such as network geometry, network size and phase-picking accuracy have significant\ud effects on the precision of seismic hypocentre locations. In turn, the precision of the hypocentral locations\ud dictates the degree to which morphological details within seismic swarms may be resolved. The Icelandic\ud national seismic network (SIL) is designed to monitor seismic activity across large expanses of Iceland in realtime\ud using automated earthquake detection and location software. Here we examine the performance of the\ud SIL network relative to a much denser, local network of seismometers deployed around the Askja volcano in\ud the Northern Volcanic Zone. A subset of earthquakes from the 2007–2008 dyke intrusion beneath Mt. Upptyppingar\ud is used to compare single- and multi-event hypocentral locations. Specifically, we highlight 288, high\ud signal-to-noise ratio events that occurred during an intensive sequence of earthquakes from 6–24 July 2007,\ud when the temporary Askja network was active. A careful refinement of phase onsets recorded by our wellconfigured,\ud dense network of receivers reveals hypocentres clustered tightly on a planar structure, interpreted\ud as a dyke dipping at 49. The root-mean-square (RMS) misfit to the plane (114 m) is only slightly greater than\ud the uncertainties in relative locations of the earthquakes themselves, and constitutes a three-fold reduction in\ud RMS misfit over SIL relative locations. The improved precision, facilitated predominantly by a more favourable\ud network size and configuration, permits a more detailed analysis of the intrusion.

The Eastern Layered Intrusion of the Rum Layered Suite\ud comprises paired peridotite and allivalite (troctolite and gabbro)layers forming 16 macro-rhythmic units. Whereas the majority of these macro-units are believed to have formed by a process of crystal^liquid differentiation involving successive accumulation of crystals from multiple picritic replenishments of the chamber, the Unit 9 peridotite is interpreted as a layer-parallel picrite intrusion. Closely correlated with this discontinuous peridotite body\ud is a distinctive feature generally known as theWavy Horizon, which divides the overlying allivalite into a lower troctolite and an upper gabbro along a well-defined undulating contact. We propose that the Wavy Horizon is a metasomatic feature formed consequent to the removal of clinopyroxene from an original gabbroic mush.\ud Foundering of the mush into the picritic sill resulted in the replacement of the original interstitial liquid by one saturated only in olivine ( plagioclase). Progressive through-flow of this liquid resulted in the stripping out of clinopyroxene from the lower parts of the allivalite.We interpret theWavy Horizon as a reaction front, representing the point at which the invading liquid became saturated in clinopyroxene.The distinctive pyroxene-enriched zone immediately above theWavy Horizon could have formed when\ud mixing of the interstitial liquids on either side of the reaction front formed a supercooled liquid oversaturated in pyroxene, as a result of the curvature of the olivine plagioclase clinopyroxene cotectic. The presence of many such approximately layer-parallel features, defined by differences in pyroxene content, in the Eastern Layered Intrusion of Rum suggests that such an infiltration^\ud reaction process was not unique to Unit 9.

A preliminary study of the aftershocks of three earthquakes that occurred near to Corinth (Greece) in 1981 is combined with observations of the morphology and faulting to understand the evolution of the Eastern Gulf of Corinth. The well located aftershocks form a zone 60km long and 20km wide. They do not lie on the main fault planes and are mostly located between the north-dipping faulting on which the first two earthquakes occurred and the south-dipping faulting associated with the third event. A cluster of aftershocks also lies in the footwall of the eastern end of the south-dipping fault of the third event.

The southern East African Rift has an unusually large seismogenic thickness (35–40 km), which is responsible for wide tilted basins and extremely long faults with the potential for M7-8 normal-faulting earthquakes. From 6–19 December 2009, a shallow earthquake sequence (four of Mw > 5.5) hit the Karonga region of northern Lake Malawi. The location is 50 km west of the rift-bounding Livingstone Fault, within the hanging-wall. We used seismology and InSAR to obtain source parameters and combined this with information on rift structure from geomorphology and seismic profiles. The deformation is consistent with rupture of a shallow, west-dipping fault, with no evidence for the involvement of magmatic fluids. Although the Livingstone Fault dominates local geomorphology, the Karonga earthquakes demonstrate that the hanging-wall block is actively breaking up, reflecting temporal and spatial migration of activity or the release of stresses within it.

We present a compilation of 2131 high-fidelity mechanisms and centroid depths of moderate-magnitude earthquakes derived using synthetic body-waveform modeling (the Global Waveform-Modelled Earthquake Catalog v1.0—gWFM), which can be visualized and downloaded online (see Data and Resources). In this article, we describe the methods used to construct the gWFM and present a comparison between the earthquake depths and focal mechanisms in the gWFM with those derived by the International Seismological Centre, Global Centroid Moment Tensor (Global CMT) project, and the U.S. Geological Survey (USGS) W-phase, as well as 60 events studied using geodesy. We find that 20%–30% of the earthquakes in routine global catalogs have depths that differ by more than 10 km from those in the gWFM. Shallow-crustal earthquakes of Mw 5–6 are typically the worst located in depth by routine catalogs. Over 90% of the earthquakes in the gWFM are within ±15° in strike, ±5° in dip, and ±15° in rake of the Global CMT and USGS W-phase best double-couple moment tensor solutions. However, the mechanisms of shallow Mw 5–6 earthquakes in the routine catalogs can be inaccurate, due to the well-known insensitivity of long-period surface waves to the vertical dip-slip components of the moment tensor. The gWFM is an archive of well-constrained earthquake source parameters, though it will continue to update as new earthquake mechanisms and depths are published, thereby remaining an up-to-date research tool.

The Earth's inner core displays transverse velocity anisotropy with cylindrical symmetry, which causes the anomalous zonal splitting of inner core sensitive normal modes. In this paper, we extend existing theory for calculating normal mode splitting from models of cylindrical velocity anisotropy to include models of anelastic attenuation anisotropy. Furthermore, we derive the equations that can be used to rewrite the attenuation anisotropy parameters natural to normal mode considerations in terms of body wave attenuation anisotropy for compressional and shear waves.

Along mid-ocean ridges the extending crust is segmented1 on length scales of 10–1,000 km. Where rift segments are offset from one another, motion between segments is accommodated by transform faults that are oriented orthogonally to the main rift axis. Where segments overlap, non-transform offsets with a variety of geometries2 accommodate shear motions. Here we use micro-seismic data to analyse the geometries of faults at two overlapping rift segments exposed on land in north Iceland. Between the rift segments, we identify a series of faults that are aligned sub-parallel to the orientation of the main rift. These faults slip through left-lateral strike-slip motion. Yet, movement between the overlapping rift segments is through right-lateral motion. Together, these motions induce a clockwise rotation of the faults and intervening crustal blocks in a motion that is consistent with a bookshelf-faulting mechanism, named after its resemblance to a tilting row of books on a shelf3. The faults probably reactivated existing crustal weaknesses, such as dyke intrusions, that were originally oriented parallel to the main rift and have since rotated about 15° clockwise. Reactivation of pre-existing, rift-parallel weaknesses contrasts with typical mid-ocean ridge transform faults and is an important illustration of a non-transform offset accommodating shear motion between overlapping rift segments.

The apparent success of inverse modeling of continent‐wide drainage inventories is perplexing. An ability to obtain reasonable fits between observed and calculated longitudinal river profiles implies that drainage networks behave simply and predictably at length scales of O(102–103) km and time scales of O(100–102) Ma. This behavior suggests that rivers respond in an organized way to large‐scale tectonic forcing. On the other hand, stream power laws are empirical approximations since fluvial processes are complex, nonlinear, and probably susceptible to disparate temporal and spatial shocks. To bridge the gap between these different perceptions of landscape evolution, we present and interpret a suite of power spectra for African river profiles that traverse different climatic zones, lithologic boundaries, and biotic distributions. At wavelengths urn:x-wiley:jgrf:media:jgrf20973:jgrf20973-math-0001102 km, power spectra have slopes of −2, consistent with red noise, demonstrating that profiles are self‐similar at these length scales. At wavelengths urn:x-wiley:jgrf:media:jgrf20973:jgrf20973-math-0002102 km, there is a crossover transition to slopes of −1, consistent with pink noise, for which power scales according to the inverse of wavenumber. Onset of this transition suggests that spatially correlated noise, perhaps generated by instabilities in water flow and by lithologic heterogeneities, becomes more prevalent at wavelengths shorter than ∼100 km. At longer wavelengths, this noise gradually diminishes and self‐similar scaling emerges. Our analysis is consistent with the concept that complexities of river profile development can be characterized by an adaptation of the Langevin equation, by which simple advective models of erosion are driven by a combination of external forcing and noise.

We have studied the seismicity of the western margin of the South Caspian Basin (SCB) and the neighbouring Talesh fold and thrust belt. We have used the hypocentroidal decomposition multiple-event location technique to obtain accurate location of events recorded during 2 yr of observation. Data from a temporary seismic network in northwest Iran and other national and regional networks were combined to make an accurate assessment of seismicity in the region. Significant offshore seismicity is observed in a 50-km wide margin of the SCB. East of the Talesh Fault along the Caspian coastline, the depth of seismicity varies from 20 to 47 km. This pattern extends inland about 20–25 km west of the North Talesh Fault. This pattern of seismicity indicates that the basement slab of the South Caspian is undergoing intense seismic deformation as it is underthrusting beneath the northern Talesh, whereas the sedimentary cover deforms aseismically. The seismicity, depths, and previous focal mechanisms of the larger offshore events are consistent with low-angle underthrusting of the South Caspian floor. Within the Talesh, seismicity is mostly concentrated in the northern and southern structural arcs of the range, where deformation is more intense and complicated. Shallow crustal seismicity in the eastern flank of the Talesh is much less intense than in the western flank, where it signifies the deformation of the upper continental crust. One major observation is the lack of any significant N–S alignment of shallow epicentres inside the central Talesh to match the observed right-lateral shear deformation there. This suggests that shear deformation inside the Talesh may have a distributed nature, rather than being concentrated on a single thorough-going fault zone, as the Talesh moves northward relative to the South Caspian. We have determined a new moment tensor solution in the southwestern Talesh, with a dominant N–S trending right-lateral motion, the only solution so far confirming along-strike shear deformation in the Talesh.

We investigate fault friction from dynamic modeling of fault slip prior to and following the Mw 6.0 earthquake which ruptured the Parkfield segment of the San Andreas Fault in 2004. The dynamic modeling assumes a purely rate-strengthening friction law, with a logarithmic dependency on sliding rate: μ=μ *+a-blnvv*. The initial state of stress is explicitly taken into account, and afterslip is triggered by the stress change induced by the earthquake source model given a priori. We consider different initial stress states and two coseismic models, and invert for the other model parameters using a nonlinear inversion scheme. The model parameters include the reference friction μ*, the friction rate dependency characterized by the quantity a-b, assumed to be either uniform or depth dependent. The model parameters are determined from fitting the transient postseismic geodetic signal measured at continuous GPS stations. Our study provides a view of frictional properties at the kilometers scale over the 0-15 km depth illuminated by the coseismic stress change induced by the Parkfield earthquake. The reference friction is estimated to be between 0.1 and 0.5. With independent a priori constraints on the amplitude of differential stress, the range of possible values narrows down to 0.1-0.17. The friction rate coefficient a-b is estimated to be � 10- 3 - 10- 2 with a hint that it increases upward from about 1-3 � 10-3 at 3-7 km depth to about 4-7 � 10-3 at 0-1 km depth. It is remarkable that our results are consistent with frictional properties measured on rock samples recovered from the fault zone thanks to the San Andreas Fault Observatory at Depth experiment. ©2013. American Geophysical Union. All Rights Reserved.

We present BurnMan, an open-source mineral physics toolbox to determine elastic properties for specified compositions in the lower mantle by solving an Equation of State (EoS). The toolbox, written in Python, can be used to evaluate seismic velocities of new mineral physics data or geodynamic models, and as the forward model in inversions for mantle composition. The user can define the composition from a list of minerals provided for the lower mantle or easily include their own. BurnMan provides choices in methodology, both for the EoS and for the multiphase averaging scheme. The results can be visually or quantitatively compared to observed seismic models. Example user scripts show how to go through these steps. This paper includes several examples realized with BurnMan: First, we benchmark the computations to check for correctness. Second, we exemplify two pitfalls in EoS modeling: using a different EoS than the one used to derive the mineral physical parameters or using an incorrect averaging scheme. Both pitfalls have led to incorrect conclusions on lower mantle composition and temperature in the literature. We further illustrate that fitting elastic velocities separately or jointly leads to different Mg/Si ratios for the lower mantle. However, we find that, within mineral physical uncertainties, a pyrolitic composition can match PREM very well. Finally, we find that uncertainties on specific input parameters result in a considerable amount of variation in both magnitude and gradient of the seismic velocities.

The Mentawai segment of the Sumatra subduction zone is locked and likely to produce a large earthquake in the near future. A part of this locked zone ruptured on 12 September 2007 producing twin earthquakes of Mw = 8.5 and 7.9. Recently, a third earthquake of Mw = 7.8 occurred on the 25th October 2010, SW of Pagai Island, Sumatra. The earthquake generated an unexpected very large tsunami on Pagai Islands with run‐up height of up to 8 m. Here we present seismic reflection and bathymetry images from the 2010 epicentral region acquired before the earthquake. We find that the frontal thrust is the main active fault in this region and might have ruptured up to the seafloor at 6 km water depth uplifting the water column and producing a large tsunami. Furthermore, finite fault models indicate that this earthquake ruptured the frontal section of the subduction zone, which is generally believed to slip aseismically and be incapable of producing large earthquakes. The presence of aftershocks near the subduction front further confirms that the frontal section of the subduction zone is not aseismic. If the rest of the Mentawai locked zone ruptures the frontal section of the subduction zone during a megathrust, then the resulting tsunami in the Indian Ocean could be devastating.

This dissertation uses a variety of seismological and mathematical tools to investigate the mechanisms\ud of deep earthquakes, the properties of fault planes and rupture at high pressure within the Earth,\ud and the local conditions where deep earthquakes occur. This work is split into two parts, focusing\ud respectively on the properties and conditions of deep faulting.

Convergence in the eastern Mediterranean of oceanic Nubia with Anatolia and the Aegean is complex and poorly understood. Large volumes of sediment obscure the shallow structure of the subduction zone, and since much of the convergence is accommodated aseismically, there are limited earthquake data to constrain its kinematics. We present new source models for recent earthquakes, combining these with field observations, published GPS velocities and reflection-seismic data to investigate faulting in three areas: the Florence Rise, SW Turkey and the Pliny and Strabo Trenches. The depths and locations of earthquakes reveal the geometry of the subducting Nubian plate NE of the Florence Rise, a bathymetric high that is probably formed by deformation of sediment at the surface projection of the Anatolia-Nubia subduction interface. In SW Turkey, the presence of a strike-slip shear zone has often been inferred despite an absence of strike-slip earthquakes. We show that the GPS-derived strain-rate field is consistent with extension on the orthogonal systems of normal faults observed in the region and that strike-slip faulting is not required to explain observed GPS velocities. Further SW, the Pliny and Strabo Trenches are also often interpreted as strike-slip shear zones, but almost all nearby earthquakes have either reverse-faulting or normal-faulting focal mechanisms. Oblique convergence across the trenches may be accommodated either by a partitioned system of strike-slip and reverse faults or by oblique slip on the Aegean-Nubia subduction interface. The observed late-Quaternary vertical motions of coastlines close to the subduction zone are influenced by the interplay between: (1) thickening of the material overriding the subduction interface associated with convergence, which promotes coastal uplift; and (2) subsidence due to extension and associated crustal thinning. Long-wavelength gravity data suggest that some of the observed topographic contrasts in the eastern Mediterranean are supported by mantle convection. However, whether the convection is time dependent and whether its pattern moves relative to Nubia are uncertain, and its contribution to present-day rates of vertical coastal motions is therefore hard to constrain. The observed extension of the overriding material in the subduction system is probably partly related to buoyancy forces arising from topographic contrasts between the Aegean, Anatolia and the Mediterranean sea floor, but the reasons for regional variations are less clear.

Topographic domes that are located far from active plate boundaries are often characterised by rapid, youthful\ud uplift, contemporaneous mafic volcanism, radial drainage patterns, and positive long-wavelength gravity anomalies. There is increasing evidence that they are underlain by anomalously low sub-plate seismic velocities. Despite their\ud well-known geomorphological expression, the origin of these epeirogenic features remains enigmatic and is much debated. Here, we investigate potential mechanisms of uplift by combining disparate observations from the Borborema\ud and Angolan plateaux that straddle the Brazilian and southwest African margins, respectively. Oceanic residual depth\ud measurements, drainage analysis, stratigraphic architecture, emergent marine terraces and basement denudation are\ud used to constrain their regional uplift histories. In both cases, the bulk of topographic growth occurred within the\ud last 30 Ma in the absence of significant tectonic deformation. We derive present-day mantle temperature and lithospheric thickness from Neogene to recent volcanic trace element compositions and upper mantle shear wave velocities. Volcanic geochemistry in northeast Brazil is compatible with decompression melting of warm asthenosphere and potentially a minor contribution from metasomatised lithospheric mantle. In Angola, melting of metasomatised lithosphere is most likely triggered by injection of small-degree asthenospheric-derived melts. We find no evidence for an\ud asthenospheric thermal anomaly > 50 ◦ C above ambient beneath either region. Present-day lithospheric thickness is ∼ 100 km beneath Angola and could be as thin as 60 km in the Borborema Province. For Angola, thermobarometry on\ud mantle xenocrysts from Cretaceous kimberlites is used to estimate palaeogeothermal gradients. Results indicate a pre-\ud existing gradient in lithospheric thickness between the edge of the Congo craton and the centre of the Angolan dome\ud at ∼ 120 Ma. This gradient likely steepened as a result of additional Neogene thinning by ∼ 30 ± 10 km beneath the\ud centre of the dome. We conclude that the mechanism for Neogene epeirogenic uplift of the Borborema and Angolan\ud domes is introduction of a small positive temperature anomaly into the asthenosphere that causes thermomechanical\ud thinning of the overlying lithospheric mantle.

We agree that there is now convincing evidence which favours the idea that a region encompassing the British Isles is dynamically supported by convective circulation beneath the lithospheric plate. For example, the coincidence of a long wavelength (>800 km) free-air gravity anomaly with anomalously slow velocities beneath the lithospheric plate suggests that regional elevation of Northern Britain is maintained by a tongue of hot asthenosphere protruding away from Iceland (Jones et al. 2002a; Arrowsmith et al. 2005). These observations are corroborated by the crustal thickness measurements of Davis et al. (2012). The Icelandic plume initiated 62 Ma and there is excellent evidence in the Irminger and Icelandic basins that this plume has waxed and waned with different periodicities (Jones et al. 2002b; Poore et al. 2009). This transient activity has manifest itself in geochemical and earthquake records along the Reykjanes Ridge (Poore et al. 2011; Parnell-Turner et al. 2013), in the Neogene record of deep-water overflow across the Greenland-Scotland Ridge (Poore et al. 2006), and along the fringing continental margins where spectacular buried ephemeral landscapes have been encountered (Shaw-Champion et al. 2008; Rudge et al. 2008; Hartley et al. 2011). Around the British Isles, there is considerable onshore and offshore evidence favouring a history of transient and/or permanent vertical motions which may reflect Palaeogene and Neogene plume activity (White & Lovell 1997; Jones et al. 2001; Jones et al. 2002a; Al-Kindi et al. 2003; Jones & White 2003; Mackay et al. 2005). This regional pattern of epeirogeny is locally complicated by faulting and folding, which is often difficult to date with precision.

We study the longitudinal spreading of a passive tracer by a two-dimensional pressure-driven flow through a composite layer of porous rock which is bounded above and below by impermeable seal rock. We focus on the flow across the interface between two neighbouring zones of the rock. First, we show that, with isotropic permeability, if the interface between the two zones is tilted relative to the boundaries, then this results in a difference in travel times across the formation which in turns leads to a net shear flow. We explore the strength of this shear as a function of (a) the permeability ratio across the interface, and (b) the interface angle. Second, we show that if one zone of the rock is cross-bedded, then with uniform flow, the pressure gradient is directed at an angle to the boundary. As a result, there is a transition zone across the interface, which again leads to a net shear, even if the interface is normal to the boundaries of the layer. We explore the competition between these effects, showing how they may combine constructively to produce a larger shear, or may negate one another, reducing or reversing the sign of the shear, depending on the angle of the interface, the degree of anisotropy and the change in effective downstream permeability across the interface. We discuss some of the implications of this shear for modelling flow in such composite rocks.

In this paper we examine the role of thermal mass in buffering the interior temperature of a naturally ventilated building from the diurnal fluctuations in the environment. First, we show that the effective thermal mass which is in good thermal contact with the air is limited by the diffusion distance into the thermal mass over one diurnal temperature cycle. We also show that this effective thermal mass may be modelled as an isothermal mass. Temperature fluctuations in the effective thermal mass are attenuated and phase-shifted from those of the interior air, and therefore heat is exchanged with the interior air. The evolution of the interior air temperature is then controlled by the relative magnitudes of (i) the time for the heat exchange between the effective thermal mass and the air; (ii) the time for the natural ventilation to replace the air in the space with air from the environment; and (iii) the period of the diurnal oscillations of the environment. Through analysis and numerical solution of the governing equations, we characterize a number of different limiting cases. If the ventilation rate is very small, then the thermal mass buffers the interior air temperature from fluctuations in the environment, creating a near-isothermal interior. If the ventilation rate increases, so that there are many air changes over the course of a day, but if there is little heat exchange between the thermal mass and interior air, then the interior air temperature locks on to the environment temperature. If there is rapid thermal equilibration of the thermal mass and interior air, and a high ventilation rate, then both the thermal mass and the interior air temperatures lock on to the environment temperature. However, in many buildings, the more usual case is that in which the time for thermal equilibration is comparable to the period of diurnal fluctuations, and in which ventilation rates are moderate. In this case, the fluctuations of the temperature of the thermal mass lag those of the interior air, which in turn lag those of the environment. We consider the implications of these results for the use of thermal mass in naturally ventilated buildings.