ESC Publications - Cambridge Univesity
Not a member yet
4214 research outputs found
Sort by
Active faulting within a megacity: the geometry and slip rate of the Pardisan thrust in central Tehran, Iran
Tehran, the capital city of Iran with a population of over 12 million, is one of the largest urban centres within the seismically active Alpine–Himalayan orogenic belt. Although several historic earthquakes have affected Tehran, their relation to individual faults is ambiguous for most. This ambiguity is partly due to a lack of knowledge about the locations, geometries and seismic potential of structures that have been obscured by dramatic urban growth over the past three decades, and which have covered most of the young geomorphic markers and natural exposures. Here we use aerial photographs from 1956, combined with an ∼1 m DEM derived from stereo Pleiades satellite imagery to investigate the geomorphology of a growing anticline above a thrust fault—the Pardisan thrust—within central Tehran. The topography across the ridge is consistent with a steep ramp extending from close to the surface to a depth of ∼2 km, where it presumably connects with a shallow-dipping detachment. No primary fault is visible at the surface, and it is possible that the faulting dissipates in the near surface as distributed shearing. We use optically stimulated luminescence to date remnants of uplifted and warped alluvial deposits that are offset vertically across the Pardisan fault, providing minimum uplift and slip-rates of at least 1 mm yr−1. Our study shows that the faults within the Tehran urban region have relatively rapid rates of slip, are important in the regional tectonics, and have a great impact on earthquake hazard assessment of the city and surrounding region
The 2008 Methoni earthquake sequence: the relationship between the earthquake cycle on the subduction interface and coastal uplift in SW Greece
Seismological, GPS and historical data suggest that most of the 40 mm yr−1 convergence at the Hellenic Subduction Zone is accommodated through aseismic creep, with earthquakes of MW ≲ 7 rupturing isolated locked patches of the subduction interface. The size and location of these locked patches are poorly constrained despite their importance for assessment of seismic hazard. We present continuous GPS time-series covering the 2008 MW 6.9 Methoni earthquake, the largest earthquake on the subduction interface since 1960. Post-seismic displacements from this earthquake at onshore GPS sites are comparable in magnitude with the coseismic displacements; elastic-dislocation modelling shows that they are consistent with afterslip on the subduction interface, suggesting that much of this part of the interface is able to slip aseismically and is not locked and accumulating elastic strain. In the Hellenic and other subduction zones, the relationship between earthquakes on the subduction interface and observed long-term coastal uplift is poorly understood. We use cGPS-measured coseismic offsets and seismological body-waveform modelling to constrain centroid locations and depths for the 2008 Methoni MW 6.9 and 2013 Crete MW 6.5 earthquakes, showing that the subduction interface reaches the base of the seismogenic layer SW of the coast of Greece. These earthquakes caused subsidence of the coast in regions where the presence of Pliocene–Quaternary marine terraces indicates recent uplift, so we conclude that deformation associated with the earthquake cycle on the subduction interface is not the dominant control on vertical motions of the coastline. It is likely that minor uplift on a short length scale (∼15 km) occurs in the footwalls of normal faults. We suggest, however, that most of the observed Plio-Quaternary coastal uplift in SW Greece is the result of thickening of the overriding crust of the Aegean by reverse faulting or distributed shortening in the accretionary wedge, by underplating of sediment of the Mediterranean seafloor, or a combination of these mechanisms
The impact of degassing on the oxidation state of basaltic magmas: A case study of Kīlauea volcano
Volcanic emissions link the oxidation state of the Earth's mantle to the composition of the atmosphere. Whether the oxidation state of an ascending magma follows a redox buffer – hence preserving mantle conditions – or deviates as a consequence of degassing remains under debate. Thus, further progress is required before erupted basalts can be used to infer the redox state of the upper mantle or the composition of their co-emitted gases to the atmosphere. Here we present the results of X-ray absorption near-edge structure (XANES) spectroscopy at the iron K-edge carried out for a series of melt inclusions and matrix glasses from ejecta associated with three eruptions of Kīlauea volcano (Hawai‘i). We show that the oxidation state of these melts is strongly correlated with their volatile content, particularly in respect of water and sulfur contents. We argue that sulfur degassing has played a major role in the observed reduction of iron in the melt, while the degassing of H2O and CO2 appears to have had a negligible effect on the melt oxidation state under the conditions investigated. Using gas–melt equilibrium degassing models, we relate the oxidation state of the melt to the composition of the gases emitted at Kīlauea. Our measurements and modelling yield a lower constraint on the oxygen fugacity of the mantle source beneath Kīlauea volcano, which we infer to be near the nickel nickel-oxide (NNO) buffer. Our findings should be widely applicable to other basaltic systems and we predict that the oxidation state of the mantle underneath most hotspot volcanoes is more oxidised than that of the associated lavas. We also suggest that whether the oxidation states of a basalt (in particular MORB) reflects that of its source, is primarily determined by the extent of sulfur degassing
Cooperative Jahn–Teller effect and the role of strain in the tetragonal-to-cubic phase transition in MgxCu1
Temperature and composition dependences of the I41/amd → [Fd\bar 3m] phase transition in the MgxCu1 − xCr2O4 spinel solid solution, due to the melting of the cooperative Jahn–Teller distortion, have been studied by means of single-crystal X-ray diffraction. Crystals with x = 0, 0.10, 0.18, 0.43, 0.46, 0.53, 1 were grown by flux decomposition methods. All crystals have been refined in the tetragonal I41/amd space group except for the Mg end-member, which has cubic symmetry. In MgxCu1 − xCr2O4 the progressive substitution of the Jahn–Teller, d9 Cu2+ cation with spherical and closed-shell Mg2+ has a substantial effect on the crystal structure, such that there is a gradual reduction of the splitting of a and c unit-cell parameters and flattening of the tetrahedra. Single-crystal diffraction data collected in situ up to T = 1173 K show that the tetragonal-to-cubic transition temperature decreases with increasing Mg content. The strength of the Cu—Cu interaction is, in effect, modulated by varying the Cu/Mg ratio. Structure refinements of diffraction data collected at different temperatures reveal that heating results in a gradual reduction in the tetrahedron compression, which remains significant until near the transition temperature, however, at which point the distortion of the tetrahedra rapidly vanishes. The spontaneous strain arising in the tetragonal phase is large, amounting to 10% shear strain, et, and ∼ 1% volume strain, Vs, in the copper chromite end-member at room temperature. Observed strain relationships are consistent with pseudoproper ferroelastic behaviour ([e_{\rm t}^2] ∝ Vs ∝ [q_{\rm JT}^2], where qJT is the order parameter). The I41/amd → [Fd\bar 3m] phase transition is first order in character for Cu-rich samples and then evolves towards second-order character. Although a third order term is permitted by symmetry in the Landau expansion, this behaviour appears to be more accurately represented by a 246 expansion with a change from negative to positive values of the fourth-order coefficient with progressive dilution of the Jahn–Teller cation
The validity of plagioclase-melt geothermometry for degassing-driven magma crystallisation.
Any quantitative interpretation of the formation conditions of igneous rocks requires methods for determining crystallisation temperature. Accurate application of such thermobarometers relies on the attainment of equilibrium in the system to be studied. This may be particularly difficult in silicic magmas, where diffusivities are low and crystallisation kinetics sluggish. Moreover, progressive degassing of volatileTrich magmas during ascent can result in continuous changes in effective undercooling, causing particular problems in achieving equilibrium between melt and crystals that grow in response to decompression. We consider these problems in the context of plagioclase-melt equilibria for magmas undergoing decompression and degassing-driven crystallisation, using two published thermometers. The two thermometers show similar trends with key parameters but absolute temperatures can vary by > 200 °C. Analysis of decompression experiments conducted at constant temperature shows systematic variations in calculated temperature and equilibrium constant with varying decompression rate and quench pressure. This indicates that an unrecognised lack of equilibration could result in significant temperature overestimates and potentially spurious results. This highlights the need to assess for equilibrium, and we discuss problems associated with some commonly used indicators of equilibration. Finally, retrospective analysis of published plagioclase-hosted melt inclusion suites from five subduction zone volcanoes shows systematic increases in calculated temperature and equilibrium constant with decreasing H₂O concentration. While this could represent the signature of latent heat of crystallisation, we suggest that such patterns should be treated with caution unless there is clear evidence of sustained equilibrium between plagioclase and melt during decompression
Flash flood events recorded by air temperature changes in caves: A case study in Covadura Cave (SE Spain)
On 28th September 2012, more than 150 mm rain fell in just two hours in some points of southeastern Spain, triggering intense flash floods that resulted in the death of ten people and widespread material damage. In the gypsum karst of Sorbas, rainfall intensity reached 33 mm/h. Air temperature monitoring in different levels of Covadura Cave, down to 85 m depth, enabled the effect of this extreme episode on the cave microclimate to be evaluated in real time. The cave air temperature increased by between 0.9 and 4.1 °C as a result of water flow into the cavity and intense mixing of air masses, in addition to the displacement of deeper air masses toward shallower levels produced by fast recharge of the surrounding karst aquifer. The lag between peak rainfall intensity and the highest cave air temperature was 5–6 h, indicating the response time of the karst to this rainfall event. No trends with depth were observed, suggesting that water not only flowed in through the main cave entrance but also through secondary accesses and fractures. Furthermore, the size of the cave passages and the intensity of air turbulence generated by waterfalls in the cave played an important role in producing these temperature differences. Even though the rainfall event lasted 10 h, cave air temperature did not return to pre-flash flood values until more than 20 days later. This indicates that, while waterflow through the cave might stop a few hours after the rainfall event, cave air temperature can be affected over a longer period. This can be explained by slow groundwater level decreasing of the surrounding karst aquifer and latent heat liberation produced by moisture condensation on the cave walls. Our results show how continuous monitoring of air temperature in caves can be a useful tool for evaluating the short-term effects of flash floods in subterranean karst systems
Not all feldspars are equal: a survey of ice nucleating properties across the feldspar group of minerals
Mineral dust particles from wind-blown soils are known to act as effective ice nucleating particles in the atmosphere and are thought to play an important role in the glaciation of mixed phase clouds. Recent work suggests that feldspars are the most efficient nucleators of the minerals commonly present in atmospheric mineral dust. However, the feldspar group of minerals is complex, encompassing a range of chemical compositions and crystal structures. To further investigate the ice-nucleating properties of the feldspar group we measured the ice nucleation activities of 15 characterized feldspar samples. We show that alkali feldspars, in particular the potassium feldspars, generally nucleate ice more efficiently than feldspars in the plagioclase series which contain significant amounts of calcium. We also find that there is variability in ice nucleating ability within these groups. While five out of six potassium-rich feldspars have a similar ice nucleating ability, one potassium rich feldspar sample and one sodium-rich feldspar sample were significantly more active. The hyper-active Na-feldspar was found to lose activity with time suspended in water with a decrease in mean freezing temperature of about 16 °C over 16 months; the mean freezing temperature of the hyper-active K-feldspar decreased by 2 °C over 16 months, whereas the "standard" K-feldspar did not change activity within the uncertainty of the experiment. These results, in combination with a review of the available literature data, are consistent with the previous findings that potassium feldspars are important components of arid or fertile soil dusts for ice nucleation. However, we also show that there is the possibility that some alkali feldspars may have enhanced ice nucleating abilities, which could have implications for prediction of ice nucleating particle concentrations in the atmosphere
Cenozoic Epeirogeny of the Indian Peninsula
Peninsular India is a cratonic region with asymmetric relief manifest by eastward tilting from the 1.5 km high Western Ghats escarpment toward the floodplains of eastward-draining rivers. Oceanic residual depth measurements on either side of India show that this west-east asymmetry is broader scale, occurring over distances of >2,000 km. Admittance analysis of free-air gravity and topography shows that the elastic thickness is 10 ±3 km, suggesting that regional uplift is not solely caused by flexural loading. To investigate how Indian physiography is generated, we have jointly inverted 530 river profiles to determine rock uplift rate as a function of space and time. Key erosional parameters are calibrated using independent geologic constraints (e.g. emergent marine deposits, elevated paleosurfaces, uplifted lignite deposits). Our results suggest that regional tilt grew at rates of up to 0.1 mm a– 1 between 25 Ma and the present day. Neogene uplift initiated in the south and propagated northward along the western margin. This calculated history is corroborated by low-temperature thermochronologic observations, by sedimentary flux of clastic deposits into the Krishna-Godavari delta, and by sequence stratigraphic architecture along adjacent rifted margins. Onset of regional uplift predates intensification of the Indian monsoon at 8 Ma, suggesting that rock uplift rather than climatic change is responsible for modern-day relief. A positive correlation between residual depth measurements and shear wave velocities beneath the lithosphere suggests that regional uplift is generated and maintained by temperature anomalies of ±100°C within a 200 ±25 km thick asthenospheric channel. This article is protected by copyright. All rights reserved
Organic membranes determine the pattern of the columnar prismatic layer of mollusc shells
The degree to which biological control is exercised compared to physical control of the organization of biogenic materials is a central theme in biomineralization. We show that the outlines of biogenic calcite domains with organic membranes are always of simple geometries, while without they are much more complex. Moreover, the mineral prisms enclosed within the organic membranes are frequently polycrystalline. In the prismatic layer of the mollusc shell, organic membranes display a dynamics in accordance with the von Neumann–Mullins and Lewis Laws for two-dimensional foam, emulsion and grain growth. Taken together with the facts that we found instances in which the crystals do not obey such laws, and that the same organic membrane pattern can be found even without the mineral infilling, our work indicates that it is the membranes, not the mineral prisms, that control the pattern, and the mineral enclosed within the organic membranes passively adjusts to the dynamics dictated by the latter
Inferences of mantle viscosity based on ice age data sets: Radial structure
We perform joint nonlinear inversions of glacial isostatic adjustment (GIA) data, including the following: postglacial decay times in Canada and Scandinavia, the Fennoscandian relaxation spectrum (FRS), late-Holocene differential sea level (DSL) highstands (based on recent compilations of Australian sea level histories), and the rate of change of the degree 2 zonal harmonic of the geopotential, J2. Resolving power analyses demonstrate the following: (1) the FRS constrains mean upper mantle viscosity to be ∼3 × 1020 Pa s, (2) postglacial decay time data require the average viscosity in the top ∼1500 km of the mantle to be 1021 Pa s, and (3) the J2 datum constrains mean lower mantle viscosity to be ∼5 × 1021 Pa s. To reconcile (2) and (3), viscosity must increase to 1022–1023 Pa s in the deep mantle. Our analysis highlights the importance of accurately correcting the J2 observation for modern glacier melting in order to robustly infer deep mantle viscosity. We also perform a large series of forward calculations to investigate the compatibility of the GIA data sets with a viscosity jump within the lower mantle, as suggested by geodynamic and seismic studies, and conclude that the GIA data may accommodate a sharp jump of 1–2 orders of magnitude in viscosity across a boundary placed in a depth range of 1000–1700 km but does not require such a feature. Finally, we find that no 1-D viscosity profile appears capable of simultaneously reconciling the DSL highstand data and suggest that this discord is likely due to laterally heterogeneous mantle viscosity, an issue we explore in a companion study