Evolution of the rates of mass wasting and fluvial sediment transfer from the epicentral area of the 1999, Mw 7.6 earthquake

The 1999 Chichi earthquake (Mw=7.6) triggered more than 20,000 landslides in the epicentral area in central west Taiwan, and subsequent typhoons have caused an even larger number of slope failures. As a result, the suspended sediment load of the epi- central Choshui River has increased dramatically. Measurements of suspended sedi- ment at a downstream gauging station indicate that the unit sediment concentration increased about six times due to the earthquake, and decreased exponentially due to flushing by subsequent typhoons. The e-folding time scale of the seismic perturbation of sediment transfer in the Choshui River is 3-5 years. Based on this estimate of the de- cay of the erosional response to the earthquake, a mass balance can be calculated for the earthquake, including co-seismic uplift and subsidence, post-seismic relaxation, and erosion. This mass balance shows that the Chi-Chi earthquake has acted to build ridge topography in the hanging wall of the fault, but in the far field, some destruc- tion of topography has occurred. However, our estimate of seismically-driven erosion may be incomplete. A detailed analysis of landsliding in the Chenyoulan tributary of the Choshui River indicates that most co-and post seismic landslide debris remains on hillslopes within the catchment. Recent typhoons have continued to cause high rates of landsliding high in the landscape, but rates of mass wasting near the stream net- work have decreased. The full geomorphic response to the Chi-Chi earthquake may be much larger, and more protracted than indicated by river gauging data

Exploring the origin of ice-filled craters in the north polar region of Mars

We investigate the origins of enigmatic ice-filled craters in the north polar region of Mars. We test several explanations for their origin, namely: (1) as polar cap remnants (2) accumulation independently of the polar cap, and (3) upwelling of subsurface water, analogous to either aufice or pingo formation on Earth. Each of these hypotheses has a significant impact on our understanding of Mars’ recent geological and climatic history and the behaviour of water and water ice at high latitudes. We used several lines of evidence to assess the most likely formation mechanism. We first performed a crater survey based on THEMIS visual data and MOLA elevation data to identify any craters that had domal central lumps which were different from normal central peaks. From this survey we identified 17 craters for further study. These include Louth, Korolev, Dokka and other unnamed craters. Using data from orbiting spectrometers; OMEGA on ESA’s Mars Express and CRISM on NASA’s Mars Reconnaissance Orbiter; we verified that the composition of the exposed central domes was predominantly water ice. We found the domes fell into three groups: (1) those completely covered by dunes, (2) those partially covered by dunes and (3) those with no dunes. We investigated the morphology and the relative position of the domes using MOLA elevation data. We found that the domes are always asymmetrically placed within the craters. However, this asymmetry could not easily be linked to wind directions as revealed by dune slip-faces [2]. The domes often have a moat-like structure and in some cases do not cover the entire crater floor, e.g. Louth Crater. From image data, we identified six craters which possessed internal stratigraphy, in the form of regularly spaced layers, and of these we have inspected three in detail. We found that the layers possess both strong sinuosity and high angle unconformities. We interpret the internal stratigraphy as representing a sequence of regular cyclic accumulations, which produced the layers, followed by asymmetric ablation and subsequent resumption of accumulation, to produce the unconformities. Hence, the present-day shape of the domes indicates that they are in a phase of ablation.We attribute the colour contrasts between layers to different levels of dust, or particulate content. This could form a source for the dunes, which are often located on the summits of these domes. We find that this sequence is best explained by a model of atmospheric condensation. Our measurements of internal layer spacing and observations of layer stratigraphy argues that these deposits are not linked directly to a former, more extensive polar cap

Cenozoic exhumation history of the Alborz Mountains, Iran: New constraints from low-temperature chronometry

[1] The Alborz Mountains in the north of the Turkish-Iranian Plateau mainly developed in the Cenozoic as a consequence of the closure of Neo-Tethys and continental collision between Arabia and Eurasia. Cenozoic growth of the fold-and-thrust belt exploited an older Palaezoic-Mesozoic crustal fabric but the extent to which this governed the overall form of the mountain belt is unclear. To determine when and how the Alborz mountain belt has grown, apatite fission track (AFT) and (U-Th)/He (AHe) thermochronometry were performed on 46 bedrock samples collected along 8 transects across the range. AFT central ages range from 157 ± 24 Ma to 10 ± 1 Ma with most ages falling between 40 Ma and 10 Ma. AHe ages range from 17 ± 2 Ma to 6 ± 1 Ma. The data revealed enhanced exhumation ca. 35 ± 5 Ma, minor exhumation between ca. 30–20 Ma and an increase in exhumation thereafter. This pattern matches the tripartite Mid-Cenozoic stratigraphy of Central Iran, where Upper Oligocene-Lower Miocene carbonates are sandwiched between terrestrial clastic formations. The most intriguing thermochronometric signal found is a major acceleration of exhumation that initiated in the late Miocene to early Pliocene, recorded by the AHe data. There does not appear to be a direct tectonic cause for Pliocene intensified erosion, as convergence rates between Arabia and Eurasia have remained relatively constant. Enhanced exhumation at this time may reflect a climatically induced intensification of erosion during Caspian isolation and base level fall, or a regional tectonic re-organization of the Arabia-Eurasia collision, or both

Temporal response of mountain drainage basins in Taiwan to earthquake and typhoon perturbation.

In tectonically-active mountain belts, earthquake-triggered landslides deliver large amounts of sediment to rivers. In previous work, we have quantified the geomorphic impact of the 1999 Mw 7.6 Chi-Chi earthquake in Taiwan, which triggered >20,000 landslides and elevated suspended sediment loads in rivers by up to a factor of four. At the time, many coseismic landslides remained confined to hillslopes and, on the basis of four years of hydrometric data, we predicted that downslope transport of sediment would continue to occur during later storms. During the seven years since the Chi- Chi earthquake, several major typhoons storms have hit Taiwan (e.g., Typhoons Bilis, Toraji, Nari, Mindulle, Aere) and the Water Resources Agency of Taiwan has contin- ued to monitor water discharge and suspended sediment concentration. Here we use these new data to refine the spatial and temporal pattern of the decaying geomorphic response to the Chi-Chi earthquake in the face of several large typhoons. Our results indicate that the broad pattern of exponential decay in sediment concentration for a given river discharge (prevalent in winter seasons without typhoons) is punctuated by markedly elevated periods associated with typhoon storms. However, our analyses show that the change in unit sediment concentration (i.e., suspended sediment concen- tration for a unit water discharge) associated with each storm depends more strongly on the length of time elapsed since the earthquake than it does on the magnitude of the storm itself

Prediction of the area affected by earthquake-induced landsliding based on seismological parameters

We present an analytical, seismologically consistent expression for the surface area of the region within which most landslides triggered by an earthquake are located (landslide distribution area). This expression is based on scaling laws relating seismic moment, source depth, and focal mechanism with ground shaking and fault rupture length and assumes a globally constant threshold of acceleration for onset of systematic mass wasting. The seismological assumptions are identical to those recently used to propose a seismologically consistent expression for the total volume and area of landslides triggered by an earthquake. To test the accuracy of the model we gathered geophysical information and estimates of the landslide distribution area for 83 earthquakes. To reduce uncertainties and inconsistencies in the estimation of the landslide distribution area, we propose an objective definition based on the shortest distance from the seismic wave emission line containing 95 % of the total landslide area. Without any empirical calibration the model explains 56 % of the variance in our dataset, and predicts 35 to 49 out of 83 cases within a factor of 2, depending on how we account for uncertainties on the seismic source depth. For most cases with comprehensive landslide inventories we show that our prediction compares well with the smallest region around the fault containing 95 % of the total landslide area. Aspects ignored by the model that could explain the residuals include local variations of the threshold of acceleration and processes modulating the surface ground shaking, such as the distribution of seismic energy release on the fault plane, the dynamic stress drop, and rupture directivity. Nevertheless, its simplicity and first-order accuracy suggest that the model can yield plausible and useful estimates of the landslide distribution area in near-real time, with earthquake parameters issued by standard detection routines

Redistribution of multi-phase particulate organic carbon in a marine shelf and canyon system during an exceptional river flood: Effects of Typhoon Morakot on the Gaoping River-Canyon system

This is the final published version of the article. It was originally published in Marine Geology (Sparkes RB, Lin I-T, Hovius N, Galy A, Liu JT, Xu X, Yang R, Marine Geology 2015, 363, 191–201, doi:10.1016/j.margeo.2015.02.013) http://dx.doi.org/10.1016/j.margeo.2015.02.013Volumetrically, turbidity currents are the prime suppliers of sediment to the deep sea, and conveyors of organic carbon from the terrestrial biosphere and submarine shelf into marine depositional basins. They result from complex processes of erosion, transport and deposition that can be difficult to study in detail. Here we present data from the Gaoping submarine canyon system, off SW Taiwan, which was perturbed in 2009 by the addition of flood deposits following Typhoon Morakot and sampled by gravity coring less than 2 months after the event. We use the different origins of organic carbon, distinguished by their carbon and nitrogen concentrations and δ13C and δ15N isotopic composition, to compare and contrast standard and extreme sedimentological conditions. Using well-constrained end-members, the results were de-convolved into inputs of metamorphic and sedimentary fossil organic carbon eroded within the Gaoping River basin, terrestrial non-fossil carbon and marine organic matter. In the upper Gaoping Canyon, sedimentation is dominated by the highly-localised hyperpycnal input of river washload and submarine sediment slumps, each associated with extensive flooding following Typhoon Morakot, whilst the shelf experienced deposition and reworking of hemi-pelagic marine sediments. A terrestrial signal is also found in the core-top of a fine-grained shelf sample over 20 km from the Gaoping Canyon, in a region normally dominated by marine carbon deposition, showing that Morakot was an unusually large flood event. Conversely, sediment from just above the canyon thalweg contains 0.23 wt.% depth-averaged marine organic carbon (37% of the TOC content) implying that terrestrial OC-dominated turbidites are tightly constrained within the canyon. Hyperpycnal processes can lead to the rapid and efficient transport of both terrestrial and submarine sediments to more permanent burial locations.RS was supported by an Engineering and Physical Sciences Research Council (EP/P502365/1 and EP/P504120/1) studentship. JTL was supported by grant number NSC95-2745-M-110-001 for the Fate of Terrestrial–Nonterrestrial Sediments in High Yield Particle–Export River–Sea Systems Program, which provided the cores in this study. We thank Peter Talling for his insightful and constructive comments on the manuscript and a further, anonymous reviewer for generous endorsement

The Exhumation of the Western Greater Caucasus; a Thermochronometric Study

This study provides 39 new thermochronometric analyses from the western part of the Greater Caucasus, a region in which existing data are extremely limited and of questionable quality. The new results are consistent with field studies that identify Triassic to Middle Jurassic (Cimmerian) and Oligo-Miocene (Alpine) orogenic erosional events. An inverse relationship between the fission track and depositional ages of Oligo-Miocene sedimentary samples also implies some degree of Eocene erosion of the Greater Caucasus and intermediate sediment storage. Cooling ages and field relationships within the core of the range, west of Mt Elbrus, require ~5 km of Permo-Triassic exhumation and restrict the overall amount of Cenozoic exhumation to ~2.5 km. Current exhumation rates are typically low, and do not support a Plio-Pleistocene increase in climate-driven denudation. High (~1 km Ma−1) rates of exhumation are restricted to the southern flank of the range in northwest Georgia. Despite a general lack of significant seismicity within the study region, this exhumation peak is close to the largest instrumentally recorded earthquake in the Caucasus (Ms = 7.0). This may suggest that exhumation is associated with the decoupling of the sedimentary succession from its crystalline basement in the southern part of the range and the inversion of the largely Jurassic fill of the Greater Caucasus basin. Rates of exhumation are compatible with this being driven by active shortening. Further sampling and analysis are required to provide a higher-resolution, low-temperature thermochronology of Alpine exhumation, to isolate the drivers for Palaeogene Dziruli Massif cooling and uplift, and to constrain better the extent of the current, localized phase of rapid exhumation

Subducted seafloor relief stops rupture in South American great earthquakes: Implications for rupture behaviour in the 2010 Maule, Chile earthquake

Great subduction earthquakes cause destructive surface deformation and ground shaking over hundreds of kilometres. Their rupture length is limited by the characteristic strength of the subduction plate interface, and by lateral variations in its mechanical properties. It has been proposed that subduction of topographic features such as ridges and seamounts can affect these properties and stop rupture propagation, but the required relief and physical mechanisms of topographic rupture limitation are not well understood. Here, we show that the rupture limits of thirteen historic great earthquakes along the South America-Nazca plate margin are strongly correlated with subducted topography with relief >1000m, including the Juan Fernandez Ridge. The northern limit of rupture in the Mw 8.8 Maule, Chile earthquake of 27 February 2010 is located where this ridge subducts. Analysis of intermediate-magnitude earthquakes shows that in most places, the subduction of high seafloor relief creates weak, aseismic zones at the plate interface, which prevent rupture propagation, but that the Juan Fernandez Ridge is associated with a locally strong plate interface. The maximum rupture length, and thus magnitude, of great subduction earthquakes is therefore determined by the size and lateral spacing of topographic features where they are present on the subducting plate. © 2010 Elsevier B.V

Automated analysis of carbon in powdered geological and environmental samples by Raman spectroscopy.

Raman spectroscopy can be used to assess the structure of naturally occurring carbonaceous materials (CM), which exist in a wide range of crystal structures. The sources of these geological and environmental materials include rocks, soils, river sediments, and marine sediment cores, all of which can contain carbonaceous material ranging from highly crystalline graphite to amorphous-like organic compounds. In order to fully characterize a geological sample and its intrinsic heterogeneity, several spectra must be collected and analyzed in a precise and repeatable manner. Here, we describe a suitable processing and analysis technique. We show that short-period ball-mill grinding does not introduce structural changes to semi-graphitized material and allows for easy collection of Raman spectra from the resulting powder. Two automated peak-fitting procedures are defined that allow for rapid processing of large datasets. For very disordered CM, Lorentzian profiles are fitted to five characteristic peaks, for highly graphitized material, three Voigt profiles are fitted. Peak area ratios and peak width measurements are used to classify each spectrum and allow easy comparison between samples. By applying this technique to samples collected in Taiwan after Typhoon Morakot, sources of carbon to offshore sediments have been identified. Carbon eroded from different areas of Taiwan can be seen mixed and deposited in the offshore flood sediments, and both graphite and amorphous-like carbon have been recycled from terrestrial to marine deposits. The practicality of this application illustrates the potential for this technique to be deployed to sediment-sourcing problems in a wide range of geological settings

An early diagnostic tool for diabetic peripheral neuropathy in rats

The skin's rewarming rate of diabetic patients is used as a diagnostic tool for early diagnosis of diabetic neuropathy. At present, the relationship between microvascular changes in the skin and diabetic neuropathy is unclear in streptozotocin (STZ) diabetic rats. The aim of this study was to investigate whether the skin rewarming rate in diabetic rats is related to microvascular changes and whether this is accompanied by changes observed in classical diagnostic methods for diabetic peripheral neuropathy. Computer-assisted infrared thermography was used to assess the rewarming rate after cold exposure on the plantar skin of STZ diabetic rats' hind paws. Peripheral neuropathy was determined by the density of intra-epidermal nerve fibers (IENFs), mechanical sensitivity, and electrophysiological recordings. Data were obtained in diabetic rats at four, six, and eight weeks after the induction of diabetes and in controls. Four weeks after the induction of diabetes, a delayed rewarming rate, decreased skin blood flow and decreased density of IENFs were observed. However, the mechanical hyposensitivity and decreased motor nerve conduction velocity (MNCV) developed 6 and 8 weeks after the induction of diabetes. Our study shows that the skin rewarming rate is related to microvascular changes in diabetic rats. Moreover, the skin rewarming rate is a non-invasive method that provides more information for an earlier diagnosis of peripheral neuropathy than the classical monofilament test and MNCV in STZ induced diabetic rats