The characteristics of seismic signals produced by lahars and pyroclastic flows: Volcán de Colima,México

The general characteristics of seismic signals produced by pyroclastic flows (generated by either the collapse of a lava dome or an eruptive column) and lahars at Volcán de Colima, México are discussed. The paper concentrates on the 2004–2006 activity associated with and following the extrusion of andesitic block-lava in October–November 2004. It is shown that the duration of the broad-band seismic records of pyroclastic flows lasts a few minutes while the duration of seismic records of lahars continues for tens of minutes or hours. The spectra of seismic records produced by pyroclastic flows are characterized by lower peak frequencies (around 3–4 Hz) than for lahars (around 6–8 Hz). This difference in the frequency content together with the difference in the duration of seismic signals allows early diagnostic of the events in real time.     

Coseismic and Aseismic Tilt Variations on Mount Etna

?—?In the last ten years (1990–1999), 21 discrete variations of continuous tilt signal have been recorded on Mount Etna, among which one episode was caused by the opening of the eruptive fracture. The remaining 20 anomalies can be classified into two categories: the first comprises 5 “instantaneous” tilt variations recorded in correspondence to the most energetic seismic events (M _L ?≥?3.3) localized on the high western part of the volcano; the second consists of 15 transient anomalies ranging from some hours to 1–2 days, observed at different times at the various tilt stations, with no correlation to seismic events or other evident volcanic episodes. The aseismic variations propagate through the volcanic edifice with a velocity between 4.5–6.0?km/day. Modeling studies suggest that the deformation is generated by a tensile source located 3–6?km SW from Etna volcano summit and 5–10?km depth.     

Controls on spatial and temporal variations in sand delivery to salmonid spawning riffles

Fine sediment infiltration into gravel interstices is known to be detrimental to incubating salmonid embryos. Infiltration into spawning riffles can show large spatial variations at the scale of a morphological unit and over time, with significant implications for embryo survival. Furthermore, some process‐based infiltration studies, and incubation‐to‐emergence models assume that fines are delivered to redds via suspension rather than bedload. This process‐based 12‐month study examined spatial patterns of predominantly sand infiltration into gravels in an upland trout stream, using infiltration baskets. An assessment of Rouse numbers for infiltrated sand indicated that it was transported predominantly as bedload at flow peaks. Clear temporal and spatial patterns existed, with highest rates of infiltration strongly associated with higher discharges (r² = 0.7, p < .05). Seasonal variations in the delivery of different grain sizes were also a feature, with enhanced contributions of 0.5–2 mm sediment during elevated winter flows and 0.125–0.5 mm sediment during spring and summer; the latter is potentially harmful to the later stages of embryo incubation. Clear spatial patterns were also evident across riffles, with highest rates of infiltration tending to occur in areas of lower relative roughness—the areas competent to transport sand for longer periods. Incubation‐to‐emergence models should take into consideration spatial patterns of fine sediment dynamics at the pool–riffle scale, to improve prediction.     

Under the surface: Pressure-induced planetary-scale waves,volcanic lightning,and gaseous clouds caused by the submarine eruption of Hunga Tonga-Hunga Ha'apai volcano

We present a narrative of the eruptive events culminating in the cataclysmic January 15, 2022 eruption of Hunga Tonga-Hunga Ha'apai Volcano by synthesizing diverse preliminary seismic, volcanological, sound wave, and lightning data available within the first few weeks after the eruption occurred. The first hour of eruptive activity produced fast-propagating tsunami waves, long-period seismic waves, loud audible sound waves, infrasonic waves, exceptionally intense volcanic lightning and an unsteady volcanic plume that transiently reached—at 58 ?km—the Earth's mesosphere. Energetic seismic signals were recorded worldwide and the globally stacked seismogram showed episodic seismic events within the most intense periods of phreatoplinian activity, and they correlated well with the infrasound pressure waveform recorded in Fiji. Gravity wave signals were strong enough to be observed over the entire planet in just the first few hours, with some circling the Earth multiple times subsequently. These large-amplitude, long-wavelength atmospheric disturbances come from the Earth's atmosphere being forced by the magmatic mixture of tephra, melt and gasses emitted by the unsteady but quasi-continuous eruption from 0402±1–1800 UTC on January 15, 2022. Atmospheric forcing lasted much longer than rupturing from large earthquakes recorded on modern instruments, producing a type of shock wave that originated from the interaction between compressed air and ambient (wavy) sea surface. This scenario differs from conventional ideas of earthquake slip, landslides, or caldera collapse-generated tsunami waves because of the enormous (～1000x) volumetric change due to the supercritical nature of volatiles associated with the hot, volatile-rich phreatoplinian plume. The time series of plume altitude can be translated to volumetric discharge and mass flow rate. For an eruption duration of ～12 ?h, the eruptive volume and mass are estimated at 1.9 ?km³ and ～2 900 ?Tg, respectively, corresponding to a VEI of 5–6 for this event. The high frequency and intensity of lightning was enhanced by the production of fine ash due to magma—seawater interaction with concomitant high charge per unit mass and the high pre-eruptive concentration of dissolved volatiles. Analysis of lightning flash frequencies provides a rapid metric for plume activity and eruption magnitude. Many aspects of this eruption await further investigation by multidisciplinary teams. It represents a unique opportunity for fundamental research regarding the complex, non-linear behavior of high energetic volcanic eruptions and attendant phenomena, with critical implications for hazard mitigation, volcano forecasting, and first-response efforts in future disasters.     

Estimates of eruption velocity and plume height from infrasonic recordings of the 2006 eruption of Augustine Volcano,Alaska

The 2006 eruption of Augustine Volcano, Alaska, began with an explosive phase comprising 13 discrete Vulcanian blasts. These events generated ash plumes reaching heights of 3–14 km. The eruption was recorded by a dense geophysical network including a pressure sensor located 3.2 km from the vent. Infrasonic signals recorded in association with the eruptions have maximum pressures ranging from 13–111 Pa. Eruption durations are estimated to range from 55–350 s. Neither of these parameters, however, correlates with eruption plume height. The pressure record, however, can be used to estimate the velocity and flux of material erupting from the vent, assuming that the sound is generated as a dipole source. Eruptive flux, in turn, is used to estimate plume height, assuming that the plume rises as a buoyant thermal. Plume heights estimated in this way correlate well with observations. Events that exhibit strongly impulsive waveforms are underestimated by the model, suggesting that flow may have been supersonic.     

New procedure for deriving multifrequential dynamic equivalent signals (LEMA_DES): a test-study based on Italian accelerometric records

The here proposed LEMA_DES (Levelled-Energy Multifrequential Analysis for Dynamic Equivalent Signals) procedure is a new approach for defining multifrequential dynamic equivalent signals from real accelerograms, to be applied for physic-analogue and numerical geotechnical modelling of induced seismic effects. In this approach, the resulting equivalent signals satisfy criteria of spectral, energetic and kinematic equivalence to the related real prototypes. The approach was tested to analyse the accelerometric records of the November 23rd, 1980 Irpinia (Italy) earthquake. Based on 48 selected records, correlations were studied between the characteristic parameters of both real and equivalent signals. These correlations demonstrate that the proposed approach guarantees: i) the energy equivalence of the derived signals, except for a half order of magnitude, and ii) the equivalence of the peak ground acceleration (PGA) values with relative errors below 105%. The computed relative error on the cumulative energy of the LEMA_DES signals (Δr V _eq %), which have spectral amplitudes at frequencies lower than 1 Hz, drops below 30%, while the same error increases above 2500%, in the same frequency range, for sinusoidal signals obtained according to traditional approaches. The PGAs of the LEMA_DES signals show a good fit with the PGA attenuation law proposed for the central-southern Apennines. Correlations between the Arias intensities and PGAs of the equivalent signals with respect to the actual ones demonstrate that their characteristic parameters: i) well represent the spatial variation in terms of energy and ground motion; ii) reproduce an analogue earthquake scenario with respect to the reference seismic event.     

Seismic characterization and monitoring of Fucino Basin (Central Italy)

The Fucino basin (Central Italy) is one of the largest intramountain alluvial plain in the Apennines range. It has a tectonic origin related to the presence of important systems of faults located in its northern and eastern edges. Some of these faults are still active and capable of generating strong seismic events. Site effects related to the soft soils filling the basin can be very important. In this paper we show the preliminary results of a seismic network installed in the Fucino area in order to collect information about site amplification effects and geometry of the basin. We analyze ambient seismic vibrations and recordings of about 150 local earthquakes mainly related to the seismic sequence of the April 6th 2009 Mw 6.3 L’Aquila event. Moreover the strongest events of L’Aquila sequence were analyzed at the three permanent strong-motion stations operating in the area. Using standard spectral techniques we investigate the variation of resonance frequencies within the basin. The ground motion recorded in the Fucino plain is mainly characterized by strong energy at low-frequencies (f < 1 Hz) affecting both horizontal and vertical components. This is particularly evident for stations deployed in correspondence of very thick deposits of sedimentary filling, where a significant increase of ground-motion amplitude and duration is likely caused by locally generated surface waves. The amplification at low-frequencies (<1 Hz) on the horizontal components can reach up a factor of 10 in comparison to nearby stiff sites. However, we found evidences of seismic amplification phenomena also for stiff sites surrounding the basin, including stations of the Italian strong motion network. The independent geological information and the shallow shear-velocity profiles available for the basin can be combined with resonance frequencies for deriving representative geological sections to be used as base for future numerical 2D–3D modeling of the basin.     

A novel technique for automatic seismic data processing using both integral and local feature of seismograms

A novel technique for automatic seismic data processing using both integral and local feature of seismograms was presented in this paper. Here, the term integral feature of seismograms refers to feature which may depict the shape of the whole seismograms. However, unlike some previous efforts which completely abandon the DIAL approach, i.e., signal detection, phase identification, association, and event localization, and seek to use envelope cross-correlation to detect seismic events directly, our technique keeps following the DIAL approach, but in addition to detect signals corresponding to individual seismic phases, it also detects continuous wave-trains and explores their feature for phase-type identification and signal association. More concrete ideas about how to define wave-trains and combine them with various detections, as well as how to measure and utilize their feature in the seismic data processing were expatiated in the paper. This approach has been applied to the routine data processing by us for years, and test results for a 16 days’ period using data from the Xinjiang seismic station network were presented. The automatic processing results have fairly low false and missed event rate simultaneously, showing that the new technique has good application prospects for improvement of the automatic seismic data processing.     

Sea-level changes and shelf break prograding sequences during the last 400 ka in the Aegean margins: Subsidence rates and palaeogeographic implications

The subsidence rates of the Aegean margins during the Middle-Upper Pleistocene were evaluated based on new and historical seismic profiling data. High-resolution seismic profiling (AirGun, Sparker and 3.5 kHz) have shown that (at least) four major oblique prograding sequences can be traced below the Aegean marginal slopes at increasing subbottom depths. These palaeo-shelf break glacial delta sediments have been developed during successive low sea-level stands (LST prograding sequences), suggesting continuous and gradual subsidence of the Aegean margins during the last 400 ka. Subsidence rates of the Aegean margins were calculated from the vertical displacement of successive topset-to-foreset transitions (palaeo-shelf break) of the LST prograding sediment sequences.The estimated subsidence rates that were calculated in the active boundaries of the Aegean microplate (North Aegean margins, Gulfs of Patras and Corinth) are high and range from 0.7 to 1.88 m ka^?1, while the lowest values (0.34–0.60 m ka^?1) are related to the low tectonic and seismic activity margins like the margin of Cyclades plateau. Lower subsidence rates (0.34–0.90 m ka^?1) were estimated for the period 146–18 ka BP (oxygen isotopic stages 6–2) and higher (1.46–1.88 m ka^?1) for the period from 425 to 250 ka BP (oxygen isotopic stages 12/10–8). A decrease of about 50% of the subduction rates in the Aegean margins was observed during the last 400 ka.During the isotopic stages 8, 10, 11 and 12, almost the 50–60% of the present Aegean Sea was land with extensive drainage systems and delta plains and large lakes in the central and North Aegean. Marine transgression in the North Aegean was rather occurred during the isotopic 9 interglacial period. The estimated palaeomorphology should imply fan delta development and sediment failures in the steep escarpments of the North Aegean margins and high sedimentation rates and turbidite sediment accumulation in the basins. It is deduced that the Black Sea was isolated from the Mediterranean during the Pleistocene prior oxygen isotopic stage 5.     

Performance of the Italian strong motion network during the 2009, L��Aquila seismic sequence (central Italy)

On April 6, 2009, the town of L’Aquila in the Abruzzo region (central Italy) was struck by a seismic event at 01:32 (UTC), of magnitude M_W = 6.3. The mainshock was followed by a long period of intense seismic activity and within seven days after the mainshock there were seven events of magnitude M_W ≥ 5 that occurred from April 6 to April 13. This long seismic sequence was characterized by a complex rupture mechanism that involved two major normal faults of the central Apennines: the Paganica and the Gorzano faults. The strong-motions of the mainshock were recorded by 64 stations of the Italian Strong-motion Network (RAN) operated by the National Civil Protection Department (DPC). Six stations of a local strong-motion array were working in NW L’Aquila suburb area. One of them, located at about 6 km from the Paganica fault surface tip-line, set up in trigger mode, recorded continuously for more than 20 min the mainshock and the aftershocks. Besides the mainshock, the RAN stations recorded in total 78 foreshocks and aftershocks of M_L ≥ 3.5, during the period from January to December 2009. The corresponding waveforms provide the most extensive digital strong ground motion data set ever recorded in Italy. Moreover, the 48 three-component observations of events of magnitude M_W ≥ 5, recorded at a distance less than 15 km from each of the major involved faults, provide a significant increasing of near-field records available for the Italian territory. Six days after the mainshock, the strong-motion dataset, referred to preliminary locations of the events with M_L ≥ 4.0, was made available on the DPC web site () and at the same time it was delivered to the ITACA database (). This dataset has been used by many authors in scientific papers and by engineers, geophysicists and geologists for professional technical works. In this paper, the present-day available strong-motion signals from the L’Aquila sequence and the performance of the Italian strong-motion network in terms of the number and quality of recorded data, the geometry and data transmission system are described. In addition the role of the temporary network that represents an extension of the permanent Italian strong-motion network, supporting the emergency response by civil protection authorities and improving the network coverage has been evaluated.     

Identification of blasting sources in the Dobrogea seismogenic region,Romania using seismo-acoustic signals

In order to discriminate between quarry blasts and earthquakes observed in the Dobrogea seismogenic region, a seismo-acoustic analysis was performed on 520 events listed in the updated Romanian seismic catalogue from January 2011 to December 2012. During this time interval, 104 seismo-acoustic events observed from a distance between 110 and 230 km and backazimuth interval of 110–160° from the IPLOR infrasound array were identified as explosions by associating with infrasonic signals. WinPMCC software for interactive analysis was applied to detect and characterize infrasonic signals in terms of backazimuth, speed and frequency content. The measured and expected values of both backazimuths and arrival times for the study events were compared in order to identify the sources of infrasound. Two predominant directions for seismo-acoustic sources’ aligning were observed, corresponding to the northern and central parts of Dobrogea, and these directions are further considered as references in the process of discriminating explosions from earthquakes. A predominance of high-frequency detections (above 1 Hz) is also observed in the infrasound data. The strong influence of seasonally dependent stratospheric winds on the IPLOR detection capability limits the efficiency of the discrimination procedure, as proposed by this study.     

Meso‐scale catchment sediment budgets: combining field surveys and modeling in the Dragonja catchment,southwest Slovenia

In this paper, we present a methodology to construct a sediment budget for meso‐scale catchments. We combine extensive field surveys and expert knowledge of the catchment with a sediment delivery model. The meso‐scale Mediterranean drainage basin of the Dragonja (91 km²), southwest Slovenia, was chosen as case study area. During the field surveys, sheet wash was observed on sloping agricultural fields during numerous rainfall events, which was found to be the main source of sediment. With the sediment yield model WATEM/SEDEM the estimated net erosion on the hillslopes 4·1 t ha^–1 y^–1 (91% of inputs). The second source, bank erosion (4·2%; 0·25 t ha^–1 y^–1) was monitored during several years with erosion pins and photogrammetric techniques. The last source, channel incision, was derived from geomorphological mapping and lichenomery and provided 3·8% (0·17 t ha^–1 y^–1) of the sediment input. The river transports its suspended sediment mainly during high‐flow events (sampled with automated water samplers). About 27% (1·2 t ha^–1 y^–1) of the sediment delivered to the channel is deposited on floodplains and low terraces downstream (estimated with geomorphological mapping, coring and cesium‐137 measurements). The sediment transported as bedload disintegrates during transport to the outlet due to the softness of the bedrock material. As a result, the river carries no bedload when it reaches the sea. The results imply a build‐up of sediment in the valleys catchment. However, extreme flood events may flush large amounts of sediment stored in the lower parts of the system. Geomorphological evidence exists in the catchment that such high magnitude, low frequency events have happened in the past. Copyright © 2009 John Wiley & Sons, Ltd.     

Locating Small Aftershocks Using a Small-Aperture Temporary Seismic Array

In this paper, we developed a specialized method to locate small aftershocks using a small-aperture temporary seismic array. The array location technique uses the first P arrival times to determine the horizontal slowness vector of the incoming P wave, then combines it with S–P times to determine the event location. In order to reduce the influence of lateral velocity variation on the location determinations, we generated slowness corrections using events well-located by the permanent broadband network as calibration events, then we applied the corrections to the estimated slownesses. Applications of slowness corrections significantly improved event locations. This method can be a useful tool to locate events recorded by temporary fault-zone arrays in the near field but unlocated by the regional permanent seismic network. As a test, we first applied this method to 64 well-located aftershocks of the 1992 Landers, California, earthquake, recorded by both the Caltech/USGS Southern California Seismic Network and a small-aperture, temporary seismic array. The average horizontal and vertical separations between our locations and the well-determined catalogue locations are 1.35 and 1.75 km, respectively. We then applied this method to 132 unlocated aftershocks recorded only by the temporary seismic array. The locations show a clear tendency to follow the surface traces of the mainshock rupture.     

Detection of the Long Period Long Duration (LPLD) Events in Time-and Frequency-Domain

Long Period Long Duration (LPLD) signals are unusual seismic events that can be observed during hydraulic fracturing. These events are very similar in appearance to tectonic tremors sequences, which were first observed in subduction zones. Their nature is not well known. LPLD might be related to the productivity of the reservoir. Different methods of the LPLD events’ detection recorded during hydraulic fracturing are presented. The author applied two methods for LPLD detection–Butterworth filtering and Continuous Wavelet Transform (CWT). Additionally, a new approach to LPLD events detection–instantaneous seismic attributes–was used, common in a classical seismic interpretation but not in microseismic monitoring.     

Correlation of seismic activity and fumarole temperature at the Mt. Merapi volcano (Indonesia) in 2000

In this paper we present densely sampled fumarole temperature data, recorded continuously at a high-temperature fumarole of Mt. Merapi volcano (Indonesia). These temperature time series are correlated with continuous records of rainfall and seismic waveform data collected at the Indonesian–German multi-parameter monitoring network. The correlation analysis of fumarole temperature and precipitation data shows a clear influence of tropical rain events on fumarole temperature. In addition, there is some evidence that rainfall may influence seismicity rates, indicating interaction of meteoric water with the volcanic system. Knowledge about such interactions is important, as lava dome instabilities caused by heavy-precipitation events may result in pyroclastic flows. Apart from the strong external influences on fumarole temperature and seismicity rate, which may conceal smaller signals caused by volcanic degassing processes, the analysis of fumarole temperature and seismic data indicates a statistically significant correlation between a certain type of seismic activity and an increase in fumarole temperature. This certain type of seismic activity consists of a seismic cluster of several high-frequency transients and an ultra-long-period signal (<0.002 Hz), which are best observed using a broadband seismometer deployed at a distance of 600 m from the active lava dome. The corresponding change in fumarole temperature starts a few minutes after the ultra-long-period signal and simultaneously with the high-frequency seismic cluster. The change in fumarole temperature, an increase of 5 °C on average, resembles a smoothed step. Fifty-four occurrences of simultaneous high-frequency seismic cluster, ultra-long period signal and increase of fumarole temperature have been identified in the data set from August 2000 to January 2001. The observed signals appear to correspond to degassing processes in the summit region of Mt. Merapi.     

Uncertainty in sewer sediment deposit modelling: Detailed vs simplified modelling approaches

The paper presents the results of a study in which the uncertainty levels associated with a detailed and a simplified/parsimonious sewer sediment modelling approach have been compared. The detailed approach used an Infoworks CS sewer network model combined with a user developed sediment transport code and the simplified approach used a conceptual sewer flow and quality model. The two approaches have been applied to a single case study sewer network and the simulation results compared. The case study was selected as moderate storm events had occurred during a 2 year rainfall and sewer flow monitoring period. Flooding had been observed and this was thought to be caused by significant solids accumulation in the sewer network. As a result sediment deposit measurements were carried out over a 6 month period. Model simulations were made of this period and predictions obtained of sediment deposit location and depth. The uncertainty analysis of both modelling approaches was carried out using Monte Carlo based computational methods. This was a limitation for the detailed approach with regards to computational time. Use of the simplified model was not constrained by this issue and so a more conventional assessment of the uncertainty was possible. The simplified approach, due to its structure, only provided a temporal estimate of uncertainty at the final section of the catchment. The detailed approach enabled an assessment of uncertainty at an individual pipe scale but only at the end of the simulation period. A comparison of the uncertainty estimations from both methods at the final section of the catchment and the end of the simulation period indicated comparable values of predicted uncertainty. Therefore a complementary use of both approaches would allow reasonably comparable estimations of levels of uncertainty at both a spatial and temporal scale. The use of such modelling approaches may provide a useful decision-making tool for sewer system management.     

Small-scale heterogeneities below the Lanzhou CTBTO seismic array,from seismic wave field fluctuations

Seismic wave field fluctuations below Lanzhou Comprehensive Test Ban Treaty Organization seismic array are determined and interpreted to describe the scattering characteristics of the crust beneath the northeastern margin of Qinghai–Tibet plateau in China. The frequency-dependent intensities of the mean and fluctuation wave fields are analyzed for 21 deep (206–632 km depth) teleseismic events. The observed wave field fluctuations in the frequency range 0.5–2.5 Hz can be explained by the scattering of the teleseismic P wave front at random media-type structures. The stable and narrow range of the fluctuation parameters from these teleseismic events indicates evidence for scattered seismic phases which are generated inside the crust. The reasonable structural models contain heterogeneities with 1–3% velocity fluctuations and 2.4–8.2 km correlation lengths in the crust with 53 km thickness.     

Testing the use of OSL from quartz grains for dating debris flows in Miyun,northeast Beijing,China

Extreme seasonal summer rain storms are common in the mountains to the north east of Beijing and these often result in mass movement of sediment slurries transported for up to a few km. These debris flows can be deadly and are very destructive to infrastructure and agriculture. This project tests the application of luminescence dating to determining the return frequency of such extreme events. The high sediment concentration and the very short flow duration gives very little opportunity for daylight resetting during transport and only a small fraction of the total mass is likely to be reset before transport begins. Here we examine the quartz single-grain dose distribution from a recent known-age (<25 years) debris flow from a small (∼3.9 km²) catchment ∼140 km north of Beijing and compare it with those from three samples from a sedimentary sequence containing the evidence of multiple flow events.Multi-grain quartz OSL signals are dominated by the fast component and <1% of the 150–200 μm grains give a detectable test dose (4.5 Gy) response. Single-grain beta dose recovery gave a ratio of 0.97 ± 0.06 (n = 30) with an over-dispersion of 23 ± 8% (CAM). Both the recent known age and the palaeo-distributions are highly dispersed with over-dispersions greater than 50%. The average weighted doses range between ∼3 mGy and ∼6.5 Gy, indicating that all deposits are no more than a few thousand years. Minimum age modelling give an age estimate for the youngest sample consistent with the known age, and minimum ages for the older palaeo-deposits suggest that there have been at least 3 major debris flows in this small catchment in the last 1000 years.     

Perturbation of fluvial sediment fluxes following the 2008 Wenchuan earthquake

Quantifying the removal of co‐seismic landslide material after a large‐magnitude earthquake is central to our understanding of geomorphic recovery from seismic events and the topographic evolution of tectonically active mountain ranges. In order to gain more insight into the fluvial erosion response to co‐seismic landslides, we focus on the sediment fluxes of rivers flowing through the rupture zone of the 2008 M_w 7.9 Wenchuan earthquake in the Longmen Shan of the eastern Tibetan Plateau. Over the post‐seismic period of 2008–2013, we annually collected river sediment samples (0.25–1 mm) at 19 locations and measured the concentration of cosmogenic ¹⁰Be in quartz. When compared with published pre‐earthquake data, the ¹⁰Be concentrations declined dramatically after the earthquake at all sampling sites, but with significant spatial differences in the amplitude of this decrease, and were starting to increase toward pre‐earthquake level in several basins over the 5‐year survey. Our analysis shows that the amplitude of ¹⁰Be decrease is controlled by the amount of landslides directly connected to the river network. Calculations based on ¹⁰Be mixing budgets indicate that the sediment flux of the 0.25–1 mm size fraction increased up to sixfold following the Wenchuan earthquake. Our results also suggest that fluvial erosion became supply limited shortly after the earthquake, and predict that it could take a few years to several decades for fluvial sediment fluxes to go back to pre‐earthquake characteristics, depending on catchment properties. We also estimate that it will take at least decades and possibly up to thousands of years to remove the co‐seismic landslide materials from the catchments in the Longmen Shan. Copyright © 2017 John Wiley & Sons, Ltd.     

Examining sediment infill dynamics at Naracoorte cave megafauna sites using multiple luminescence dating signals

Relatively little is known about the long-term sediment accumulation dynamics of Naracoorte Cave Complex (NCC) solution pipe cavities, and many of the megafauna-bearing infill deposits at this globally significant Australian Pleistocene fossil locality remain partially dated or lack any numerical age control. In this study, we assess the suitability of three different luminescence dating signals for improving existing chronologies at six Late and Middle Pleistocene NCC sites (n = 22 samples), and we undertake multi-site examinations of NCC sediment infill dynamics spanning the last 550 thousand years (ka). Modern analogue samples collected from above and beneath two active cave entrances confirm that single-grain OSL, single-grain TT-OSL and multi-grain pIR-IRSL signals can be reset down to insignificant residual levels (<10⁻¹–10⁰ Gy) when compared with the natural dose ranges of interest for most NCC palaeontological applications. Replicate luminescence dating comparisons performed at six NCC fossil sites (n = 15 samples) reveal consistent ages for twenty-eight out of thirty-one paired OSL–TT-OSL, OSL–pIR-IRSL and TT-OSL–pIR-IRSL datasets. Nineteen of the twenty Middle to Late Pleistocene samples analysed from the NCC sites produce homogeneous OSL, TT-OSL and pIR-IRSL D_e datasets suggesting that the NCC solution pipe deposits considered here are generally not affected by syn-depositional mixing complications that can take place within closed karst cavities (i.e., the remobilisation of unbleached grains from pre-existing cave floor sediments, and their subsequent translocation and incorporation into deposits within deeper parts of the cavity along with externally bleached, allochthonous grain populations). Detailed examination of solution pipe dynamics at Smoke Tortoise Cave (SMT) reveals a complex accumulation history focused on the marine isotope stage (MIS) 9 and MIS 7 interglacial complexes, as well as the MIS 8e interstadial. The SMT case study highlights that NCC solution pipes are not simply associated with short-lived opening and sediment accumulation events, but may involve multiple, discontinuous deposition episodes and reactivation events. An initial multi-site examination of all published NCC infill chronologies (n = 70) appears to suggest statistically significant, preferential solution pipe development during the relatively wet parts of interglacial or interstadial cycles. The non-uniform infill age distribution implies that NCC solution pipe dynamics may have exerted taphonomic biases on fossil accumulation, which should be taken into consideration when reconstructing long-term palaeoecological histories from NCC solution pipe cavities.