Crustal structure of the Reykjanes Ridge near 62°N,on the basis of seismic refraction and gravity data

Explosion deep seismic sounding data sections of high quality had been obtained with RV Meteor in the Reykjanes Iceland Seismic Project (RRISP77 [Angenheister, G., Gebrande, H., Miller, H., Goldflam, P., Weigel, W., Jacoby, W.R., Pálmason, G., Björnsson, S., Einarsson, P., Pavlenkova, N.I., Zverev, S., Litvinenko, I.V., Loncarecic, B., Solomon, S., 1980. Reykjanes Ridge Iceland Seismic Experiment (RRISP 77). J. Geophys. 47, 228–238]) which close an information gap near 62°N. Preliminary results were presented by Weigel [Weigel, W., 1980. Aufbau des Reykjanes Rückens nach refraktionsseismischen Messungen. In: Weigel, W. (Ed.), Reykjanes Rücken, Island, Norwegischer Kontinentalrand. Abschlusskolloquium, Hamburg zur Meteor-Expedition, vol. 45. DFG, Bonn, pp. 53–61], and here we report on the data and results of interpretation. Clear refracted phases to 90 km distance permit crustal and uppermost mantle structure to be modelled by ray tracing. The apparent P-wave velocities are around 4.5, 6–6.5, 7–7.6 and 8.2–8.7 km/s, but no wide-angle reflections have been clearly seen. Accompanying sparker reflection data reveal thin sediment ponds in the axial zone and up to 400 m thick sediments at 10 Ma crustal age. Ray tracing reveals the following model below the sediments: (1) a distinct, 1–2 km thick upper crust (layer 2A) with Vp increasing with age (to 10 Ma) from <3.4 to 4.9 km/s and with a vertical gradient of 0.1–0.2 km/s/km, (2) a lower crust or layer 3 beginning at depths of 2 (axis) to 4 km (10 Ma age) below sea level with 6.1–6.8 km/s and similar vertical gradients as above, (3) the lower crust bottoms at 5.2–9.5 km depth below sea level (0–10 Ma) with a marked discontinuity, underneath which (4) Vp rises from about 7.5–7.8 km/s (0–10 Ma) with a positive vertical gradient of, again, 0.1–0.2 km/s/km such that 8 km/s would be reached at 12 km and deeper near the axis. Our preferred interpretation is that the mantle begins at the distinct discontinuity (“Moho”), but a deeper “Moho” of Vp ≈ 8 km/s cannot be excluded. From Iceland southward to 60°N several experiments show a decrease of crustal thickness from 14 to 8 km. Velocity trends with age across the ridge reflect cooling and filling of cracks, and thickness trends probably suggest volcanic productivity variations as previously suggested.Gravity inversion concentrates on a profile across the ridge with the above seismic a priori information; with 0.2–0.5 km depth uncertainty it leads to a good fit (±2.5 mGal where seismic data exist). Best fitting densities are (in kg/m³) for sediments, 2180; upper crust, 2450–2570; lower crust, 2850–2940; mantle lithosphere, 3215–3240 with a deficit for an asthenospheric wedge of no more than −100 kg/m³. The morphological ridges and troughs superimposed on the SE ridge flank are partly correlated, partly anti-correlated with the Bouguer anomaly and suggest that variable crustal density variations accompany the morphology variations.     

1-g Experimental investigation of bi-layer soil response and kinematic pile bending

The effect of soil inhomogeneity and material nonlinearity on kinematic soil–pile interaction and ensuing bending under the passage of vertically propagating seismic shear waves in layered soil, is investigated by means of 1-g shaking table tests and nonlinear numerical simulations. To this end, a suite of scale model tests on a group of five piles embedded in two-layers of sand in a laminar container at the shaking table facility in BLADE Laboratory at University of Bristol, are reported. Results from white noise and sine dwell tests were obtained and interpreted by means of one-dimensional lumped parameter models, suitable for inhomogeneous soil, encompassing material nonlinearity. A frequency range from 0.1 Hz to 100 Hz and 5 Hz to 35 Hz for white noise and sine dwell tests, respectively, and an input acceleration range from 0.015 g to 0.1 g, were employed. The paper elucidates that soil nonlinearity and inhomogeneity strongly affect both site response and kinematic pile bending, so that accurate nonlinear analyses are often necessary to predict the dynamic response of pile foundations.     

Kinematics of the eastern part of the North Anatolian Fault Zone

The North Anatolian Fault Zone (NAFZ), which marks the boundary between Anatolia and the Eurasian plate, is one of the world's most seismically active structures. Although the eastern part of NAFZ has high seismic hazard, there is a lack of geodetic information about the present tectonics of this region. Even though many scientists would like to study this area, geographical and logistical problems make performing scientific research difficult. In order to investigate contemporary neotectonic deformation on the eastern NAFZ and in its neighborhood, a relatively dense Global Positioning System (GPS) monitoring network was established in 2003. Geodetic observations were performed in three GPS campaigns in an area of 350 km × 200 km with 12-month intervals. In addition, 14 new GPS stations were measured far from the deforming area. Since this region includes the intersection of the NAFZ and the East Anatolian Fault Zone (EAFZ), deformation is complex and estimating seismic hazard is difficult. One important segment is the Yedisu segment and it has not broken since the 1784 earthquake. After the 1992 Erzincan and 2003 Pulumur earthquakes, the Coulomb stress loading on the Yedisu segment of the NAFZ has increased significantly, emphasizing the need to monitor this region. We computed the horizontal velocity field with respect to Eurasia and strain rates field as well. GPS-derived velocities relative to Eurasia are in the range of 16–24 mm/year, which are consistent with the regional tectonics. The principal strain rates were derived from the velocity field. Results show that strain is accumulating between the NAFZ and EAFZ along small secondary fault branches such as the Ovacik Fault (OF).     

Sensitivity of ground-motion scenarios to earthquake source parameters in the Tehran metropolitan area,Iran

The northern Tehran fault (NTF) is a principal active fault of the Alborz mountain belt in the northern Iran. The fault is located north of the highly populated Metropolitan Area of Tehran. Historical records and paleoseismological studies have shown that the NTF poses a high seismic risk for the Tehran region and the surrounding cities (e.g. Karaj). A series of ground-motion simulations are carried out using a hybrid kinematic-stochastic model to calculate broadband (0.1–20 Hz) ground-motion time histories for deterministic earthquake scenarios (M7.2) on the NTF. We will describe the source characteristics of the target event to develop a list of scenario earthquakes that are probably similar to a large earthquake on the NTF. The effect of varying different rupture parameters such as rupture velocity and rise time on the resulting broadband strong motions has been investigated to evaluate the range of uncertainty in seismic scenarios. The most significant parameters in terms of ground-shaking level are the rise time and the value of the rupture velocity. For the worst-case scenario, the maximum expected horizontal acceleration, and velocity at rock sites in Tehran range between 128 and 1315 cm/s/s and 11–191 cm/s, respectively. For the lowest scenario, the corresponding values range between 102 and 776 cm/s/s and 12 to 81 cm/s. Nonlinear soil effects may change these results but are not accounted for in this study. The largest variability of ground motion is observed in neighborhood of asperity and also in the direction of rupture propagation. The calculated standard deviation of all ground-motion scenarios is less than 30% of the mean. The capability of the simulation method to synthesize expected ground motions and the appropriateness of the key parameters used in the simulations are confirmed by comparing the synthetic peak ground motions (PGA, PGV and response spectra) with empirical ground-motion prediction equations.     

Magnetic susceptibility and Cs-137 inventory variability as influenced by land use change and slope positions in a hilly,semiarid region of west-central Iran

The objective of this study was to explore the slope position and land use change effects on the variability in magnetic susceptibility and ¹³⁷Cs inventory as the soil redistribution indicators in a hilly semiarid calcareous area in Iran. The selected study area is located in a hilly region with pasture and cultivation land use of Fereydunshahr, Isfahan Province in west-central Iran. In the two mentioned dominant ecosystems, four slope positions including summit, shoulder, backslope and footslope were identified and in each land use and slope position, three cores were selected to collect 72 soil samples from three depths (0–10, 10–20, 20–30 cm) in an area of 15 × 15 cm. Additional 28 soil samples were collected from the reference site for soil loss and deposition calculations by using the Cs-137 measurement. The results of the study with the use of the Cs-137 technique showed that the average soil loss in the pasture land (46.4 t ha^− 1 yr^− 1) was significantly (p < 0.05) lower than the average soil loss in the cultivated land (80.4 t ha^− 1 yr^− 1). The highest soil loss in both land uses was obtained in the shoulder position, 60.1 and 84.4 t ha^− 1 yr^− 1, respectively, for the pasture and cultivated lands. Moreover, the highest rates of soil deposition was observed in a footslope position in both land uses and they were 34 and 32.4 t ha^− 1 yr^− 1 for the pasture and cultivated lands, respectively. Magnetic susceptibility was significantly (p < 0.05) greater in pasture (χ_lf = 41.51 × 10^− 8 m³/kg) than in the cultivated land (χ_lf = 34.90 × 10^− 8 m³/kg). The pasture land with a lower soil loss rate, indicated significantly higher magnetic susceptibility in all landform positions as compared to that in the cultivated land. The results of the correlation analysis showed that among the studied soil physico-chemical properties, χ_lf (r = 0.83, p < 0.01) in the pasture land had the highest correlation with the Cs-137 inventory. Throughout the non-linear regression analysis, χ_lf was introduced for relating soil parameters and the cesium inventory explained 68% and 79% of the total variability of ¹³⁷Cs in the pasture and cultivated lands, respectively. The results implied that the variability in the magnetic susceptibility within the hillslope is consistent with the variation of the Cs-inventory; and the results thus demonstrate the slope and land use effects on soil redistribution.     

Stability of massive ground ice bodies in University Valley,McMurdo Dry Valleys of Antarctica: Using stable O–H isotope as tracers of sublimation in hyper-arid regions

To date, studies of the stability of subsurface ice in the McMurdo Dry Valleys of Antarctica have been mainly based on climate-based vapor diffusion models. In University Valley (1800 m), a small glacier is found at the base of the head of the valley, and adjacent to the glacier, a buried body of massive ice was uncovered beneath 20–40 cm of loose cryotic sediments and sandstone boulders. This study assesses the origin and stability of the buried body of massive ice by measuring the geochemistry and stable O–H isotope composition of the ice and applies a sublimation and molecular diffusion model that accounts for the observed trends. The results indicate that the buried massive ice body represents an extension of the adjacent glacier that was buried by a rock avalanche during a cold climate period. The contrasting δ¹⁸O profiles and regression slope values between the uppermost 6 cm of the buried massive ice (upward convex δ¹⁸O profile and S_D-18O = 5.1) and that below it (progressive increase in δ¹⁸O and S_D-18O = 6.4) suggest independent post-depositional processes affected the isotope composition of the ice. The upward convex δ¹⁸O profile in the uppermost 6 cm is consistent with the ice undergoing sublimation. Using a sublimation and molecular diffusion model, and assuming that diffusion occurred through solid ice, the sublimation rate needed to fit the measured δ¹⁸O profile is 0.2 ? 10^? 3 mm yr^? 1, a value that is more similar to net ice removal rates derived from ³He data from cobbles in Beacon Valley till (7.0 ? 10^? 3 mm yr^? 1) than sublimation rates computed based on current climate (0.1–0.2 mm yr^?1). We suggest that the climate-based sublimation rates are offset due to potential ice recharge mechanisms or to missing parameters, particularly the nature and thermo-physical properties of the overlying sediments (i.e., temperature, humidity, pore structure and ice content, grain size).     

Post-cyclic settlement and tilting potential of mat foundations

A simplified procedure for the estimation of cyclic shear stresses on foundation soil layers due to interactions of seismic excitation, foundation mat and overlying structural system, incorporating both kinematic and inertial aspects was described earlier in Cetin et al. [1], [2]. This simplified procedure is now used, along with laboratory-based cyclic shear and volumetric strain relationships for the assessment of cyclically induced settlement and tilting potential of mat foundations. The proposed methodology is calibrated by well-documented seismic foundation performance case histories of 3–6 story, relatively rigid residential structures with no basements. Immediately after 1999 Kocaeli and Duzce, Turkey earthquakes, foundation settlements of these case history buildings were carefully mapped relative to available elevation benchmarks, such as relatively rigid concrete pavements, drainage pipes and ditches, and entrance stairs, located in the immediate vicinity of the buildings. Tilt angles of these buildings were also mapped in orthogonal multi-directions. Relative settlement and tilt angle values vary in the range of 0–60 cm and 0–5.5 rad, respectively. As part of subsurface characterization studies, SPT and seismic CPTU were performed. For the retrieved disturbed and undisturbed samples, an intensive laboratory testing program including sieve analysis, Atterberg Limit, hydrometer, oedometer, static and cyclic triaxial and cyclic simple shear tests was executed. The foundation soil profiles generally consist of silty soils, sand–silt mixtures and silt–clay mixtures of SPT-N values varying in the range of 2–5 blows/30 cm in the upper 5 m and gradually increases up to a maximum value of 25 blows/30 cm beyond depths of 5–8 m. As part of the proposed framework, displacement potential indices, defined as the product of induced-cyclic strains and the thickness of soil sub-layers, were estimated, the sum of which produce settlement and tilting potentials. Then, these settlement and tilting potentials were calibrated against recorded settlement and tilt values through a statistically based calibration scheme. The proposed simplified procedure is shown to reliably and precisely capture both deviatoric and volumetric components of post-cyclic settlements, as well as tilting responses of mat foundations.     

Seismic response of deep Quaternary sediments in historical center of L’Aquila City (central Italy)

On 6 April 2009 a M_w=6.1 earthquake produced severe destruction and damage over the historic center of L’Aquila City (central Italy), in which the accelerometer stations AQK and AQU recorded a large amount of near-fault ground motion data. This paper analyzes the recorded ground motions and compares the observed peak accelerations and the horizontal to vertical response spectral ratios with those revealed from numerical simulations. The finite element method is considered herein to perform dynamic modeling on the soil profile underlying the seismic station AQU. The subsurface model, which is based on the reviewed surveys that were carried out in previous studies, consists of 200–400 m of Quaternary sediments overlying a Meso-Cenozoic carbonate bedrock. The Martin-Finn-Seed's pore-water pressure model is used in the simulations. The horizontal to vertical response spectral ratio that is observed during the weak seismic events shows three predominant frequencies at about 14 Hz, 3 Hz and 0.6 Hz, which may be related to the computed seismic motion amplification occurring at the shallow colluvium, at the top and base of the fluvial-lacustrine sequence, respectively. During the 2009 L’Aquila main shock the predominant frequency of 14 Hz shifts to lower values probably due to a peculiar wave-field incidence angle. The predominant frequency of 3 Hz shifts to lower values when the earthquake magnitude increases, which may be associated to the progressive softening of soil due to the excess pore-water pressure generation that reaches a maximum value of about 350 kPa in the top of fluvial-lacustrine sequence. The computed vertical peak acceleration underestimates the experimental value and the horizontal to vertical peak acceleration ratio that is observed at station AQU decreases when the earthquake magnitude increases, which reveals amplification of the vertical component of ground motion probably due to near-source effects.     

Transition zone structure under a stationary hot spot: Cape Verde

We report on a two-year seismic deployment in the Cape Verde Islands, one goal of which was to study the upper mantle to determine its structure under a hot spot that is stationary in the hot spot reference frame. We find from analysis of P-to-S receiver functions estimated from broadband seismic recordings that, within uncertainty, the time separation between the 410 and 660 km discontinuities is normal compared to radial earth models. Thus, to exist, even stationary hot spots do not require vertical thermal anomalies from deep melting sources anchored in the lower mantle or at the core–mantle boundary or their anomalies are narrower than ～ 250 km in the upper mantle.     

Monofrequency in situ damping measurements in Ottawa area soft soils

In Ottawa, Canada, unusually high amplification ratios have recently been measured in clayey silts (called ‘Leda Clays’) at low levels of earthquake-induced ground shaking. However, the contribution of seismic Q, or material damping (ξ=1/2Q), to the overall ground motion at soft soil sites across the city is not well understood. This research investigates attenuation measurements in soft soils (V_s<250 m/s) for ongoing seismic hazard evaluation in the Ottawa area. The work focuses on in situ measurements of damping in two deep boreholes drilled into Leda Clay. To investigate the possibility of frequency-dependent dynamic properties of these materials at low strains, a new approach to the spectral ratio technique has been developed for the measurement of Q_s in the field using a mono-frequency vibratory source (generating signals between 10 and 100 Hz), and two identical downhole 3-component geophones. Monofrequency signals also allowed for the measurement of dispersion (variation of velocity with frequency). Analysis of the data show that dynamic properties are, for the most part, independent of frequency in the homogenous silty soils, yielding negligible variation in shear wave velocity (<2 m/s) across the frequency test band, and small strain Q_s's ranging from 170 to 200 (damping of 0.25–0.30%) over soil thickness intervals ranging from 10 to 60 m. At intervals within 20 m of the ground surface, laminated silt and clay beds of elevated porosity are found to have slight influence on the frequency dependence of damping for frequencies greater than 70 Hz (damping increase to 0.6%).     

Geoelectrical,strain and tilt investigations of hydrological processes at the broadband Geodynamical Observatory Moxa,Germany

The Geodynamic Observatory Moxa, located in Thuringia/Germany, is dedicated to studies of temporal deformations of the earth's crust and of variations of the gravity field. One of the essential issues with respect to these investigations is the reduction of the hydrological impact on the data of the gravimeters, strainmeters and tiltmeters. In order to optimise the reductions, we investigated the changes in the hydrological conditions in the woody mountain slope above the observatory with time-lapse electrical resistivity tomography (ERT), and analysed the strain and tilt measurements for prominent signatures of pore pressure induced subsurface deformations.Here we present the results for two profiles – parallel and perpendicular to the slope – measured with ERT during 33 campaigns between June 2007 and April 2010. Resistivity changes and variations of apparent soil moisture, inferred from ERT sections, were found to primarily occur in the first two metres of the subsurface. These variations can be related to subsurface flow in the upper two metres induced by precipitation events and snowmelts. Trees close to the profiles only show a minimum impact on the resistivity and soil moisture changes.Furthermore, systematic hydrologically induced deformations can be observed in hodographs of strain and tilt measurements for large precipitation events (> 80 mm) and snowmelts. In the strain data a short-term (< 3 days) dilatational signal is found with an amplitude of 20 nstrain to 60 nstrain and a long-term (> 7 days) compressional signal between 40 nstrain and 180 nstrain. The preferential N–S direction of long-term deformational signals (> 1 week) is also observed in the tilt data. The direction of tilt changes (25 nrad–120 nrad) is nearly parallel to the drainage direction of the nearby Silberleite creek indicating variations of pore pressure gradients during hydrological events.The results of these hydrological studies at the Geodynamic Observatory Moxa can be used for removing the time dependent hydrological signal in strain and tilt data and, thus, better correction algorithms for hydrological impacts can be developed to enhance the value of the data for geodynamic studies.     

40Ar / 39Ar geochronology constraints on late miocene weathering rates in Minas Gerais,Brazil

⁴⁰Ar / ³⁹Ar incremental heating ages for twenty one grains of cryptomelane, collected at 0, 42, 45, and 60 m depths in the Cachoeira Mine weathering profile, Minas Gerais, permit calculating long-term (10 Ma time scale) weathering rate (saprolitization rate) in SE Brazil. Pure well-crystallized cryptomelane grains with high K contents (3–5 wt.%) yield reliable geochronological results. The ⁴⁰Ar / ³⁹Ar plateau ages obtained decrease from the top to the bottom of the profile (12.7 ± 0.1 to 7.6 ± 0.1 Ma at surface; 7.6  ± 0.2 to 6.1 ± 0.2 Ma at 42 m; and 7.1 ± 0.2 to 5.9 ± 0.1 Ma at 45 m; 6.6 ± 0.1 to 5.2 ± 0.1 Ma at 60 m), yielding a weathering front propagation rate of 8.9 ± 1.1 m/m.y. From the geochronological results and the mineral transformations implicit by the current mineralogy in the weathering profiles, it is possible to calculate the saprolitization rate for the Cachoeira Mine lithologies and for adjacent weathering profiles developed on granodiorites and schists. The measured weathering front propagation rate yields a saprolitization rate of 24.9 ± 3.1 t/km²/yr. This average long-term (> 10 Ma) saprolitization rate is consistent with mass balance calculations results for present saprolitization rates in weathering watersheds. These results are also consistent with long-term saprolitization rates estimated by combining cosmogenic isotope denudation rates with mass balance calculations.     

Luminescence dating of the top of a deep water core from the NESTOR site near the Hellenic Trench,east Mediterranean Sea

In the present study, the uppermost 3 cm of muddy, seabed sediment was collected from a deep-sea sediment core, drilled from a water depth of 4 km, near the deepest site of the Mediterranean Sea, outside Pylos, Greece. The core was divided into 7 layer samples, each 3–4 mm thick, in order to get an independent age assessment for each one using luminescence dating; from polymineral coarse grains in the range 30–60 μm. Between 11 and 22 aliquots were measured for each sample using the optically stimulated luminescence, single-aliquot regenerative-dose analysis. The estimated mean equivalent dose values had an uncertainty less than 3% and resulted in optical ages ranging from 3.5 to approximately 5 ka, with uncertainties lying between 5.5 and 7.4%. Multiple-aliquot, thermoluminescence-based, additive-dose, total bleaching approach provided equivalent doses with typical errors of 10–15% and ages in the range of 3.6 and 9.4 ka, with uncertainties up to 17%. Single-aliquot optical ages are shown to be relatively successful, due to their consistency with AMS ¹⁴C radiocarbon ages, obtained from Planktonic foraminifera from the same core. Luminescence dates for the topmost 1.5 cm indicate a substantial mix and post burial disturbance of the surface sediment. Below the topmost 1.5 cm, both luminescent approaches indicate ages which increase smoothly with depth. The concordant OSL and TL age estimates for the sample O5, in conjunction with specific luminescence properties and its major-element geochemical chemistry content, suggest that it was heated during the Santorini volcanic eruption. Deposition rates of 8.6–18.9 cm over 10 ka below the sample O5, provided by optical ages, exhibit an excellent level of agreement with the accumulation rate of 7–18 cm over 10 ka at the sea bottom, already reported for the site under study.     

A time probabilistic approach to seismic landslide hazard estimates in Iran

Understanding where seismically induced landslides are most likely to occur is crucial in land use planning and civil protection actions aimed at reducing property damage and loss of life in future earthquakes. For this purpose an approach proposed by Del Gaudio et al. [1] has been applied to the whole Iranian territory to provide the basis to assess location and temporal recurrence of conditions of seismic activation of slope failures, according to the Newmark's model [2]. Following this approach, occurrence probabilities for different levels of seismic shaking in a time interval of interest (50 years) were first obtained through a standard hazard estimate procedure. Then, empirical formulae in the form proposed by Jibson et al. [3] and calibrated for the main seismogenic Iranian regions were used to evaluate the slope critical acceleration (Ac)_x for which a prefixed probability exists that, under seismic shakings, Newmark's displacement D_N exceeds a threshold×corresponding to landslide triggering conditions. The obtained (Ac)_x values represent the minimum slope resistance required to limit the probability of landslide seismic triggering within the prefixed value. A map reporting the spatial distribution of these values gives comparative indications on regional different exposure of slopes to shaking capable of inducing failures and provides a reference for hazard estimate at local scale. The obtained results show that the exposure to landslide seismic induction is maximum in the Alborz Mountains region, where critical accelerations up to ∼0.1 g are required to limit the probability of seismic triggering of coherent type landslides within 10% in 50 years.     

Evaluation of site conditions for the Ankara Basin of Turkey based on seismic site characterization of near-surface geologic materials

Regional site conditions relevant for seismic hazard studies can be derived from various geologic, seismologic, and geotechnical sources. In this study, site conditions are derived for the Ankara Basin in Turkey by merging in-situ seismic measurements of dynamic properties, geologic information, and some geotechnical boring information. Field seismic refraction surveys were performed at 259 sites in the project area to classify and characterize Plio–Pleistocene fluvial deposits and Quaternary alluvial and terrace deposits. The shear-wave velocity profiles of the near-surface geologic units are used to characterize site classes according to the International Building Code [International Code Council, ICC. International Building Code. Structural and fire- and life-safety provisions (seismic, wind, accessibility, egress, occupancy and roof codes), 2006. Whittier, CA.] and the Turkish Seismic Code [Ministry of Public Works and Settlement, 1998. Turkish Seismic Code, Specification for Structures to be Built in Disaster Areas, Ankara, Turkey], and to develop a regional model for the average shear-wave velocity in the top 30 m. The resulting maps of site class indicate that the classification system for the Turkish Seismic Code results in a significant portion of the Ankara Basin being classified as Z4, the softest site class. The International Building Code site classification system results in most of the Ankara Basin being classified as D, stiff soil. These differences are caused by the Turkish Seismic Code incorporating information from only the surface layer, while the International Building Code incorporates information from the top 30 m.     

Velocity dispersion in fractured rocks in a wide frequency range

Experimental measurements of fracture-induced seismic waves velocity variations at frequencies ~ 1 kHz, ~ 40 kHz and ~ 1 MHz were performed directly in the field at the rocky outcrop and in the laboratory on specific rock samples collected from the outcrops. The peridotite–lherzolite outcrop appeared macroscopically uniform and contained three systems of visible parallel sub-vertical fractures. This rock has substantial bulk density and higher than average value of seismic wave velocity. The presence of fracture systems gives rise to its velocity anisotropy. The seismic waves passing through the rock fractures are subject to velocity dispersion and frequency dependent attenuation. Our data, obtained from field and laboratory measurements, were compared with theoretical model predictions. In this model we successfully used displacement discontinuity approach. For the velocity dispersion evaluation we used multi-frequency measurements. The a priori observation of orientations and densities of fracture sets allowed evaluation of their stiffness. Our approach revealed that the first arrivals of seismic waves can be used for evaluation of P-wave group velocities, the specific case, in which we expect anomalous velocity dispersion. Our observations contribute to the issue of up-scaling of well-log derived velocities in fractured rock to the scale of standard seismic exploration frequencies.     

Elasticity of antigorite,seismic detection of serpentinites,and anisotropy in subduction zones

Serpentinization of the mantle wedge is an important process that influences the seismic and mechanical properties in subduction zones. Seismic detection of serpentines relies on the knowledge of elastic properties of serpentinites, which thus far has not been possible in the absence of single-crystal elastic properties of antigorite. The elastic constants of antigorite, the dominant serpentine at high-pressure in subduction zones, were measured using Brillouin spectroscopy under ambient conditions. In addition, antigorite lattice preferred orientations (LPO) were determined using an electron back-scattering diffraction (EBSD) technique. Isotropic aggregate velocities are significantly lower than those of peridotites to allow seismic detection of serpentinites from tomography. The isotropic V_P/V_S ratio is 1.76 in the Voigt–Reuss–Hill average, not very different from that of 1.73 in peridotite, but may vary between 1.70 and 1.86 between the Voigt and Reuss bonds. Antigorite and deformed serpentinites have a very high seismic anisotropy and remarkably low velocities along particular directions. V_P varies between 8.9 km s^? 1 and 5.6 km s^? 1 (46% anisotropy), and 8.3 km s^? 1 and 5.8 km s^? 1 (37%), and V_S between 5.1 km s^? 1 and 2.5 km s^? 1 (66%), and 4.7 km s^? 1 and 2.9 km s^? 1 (50%) for the single-crystal and aggregate, respectively. The V_P/V_S ratio and shear wave splitting also vary with orientation between 1.2 and 3.4, and 1.3 and 2.8 for the single-crystal and aggregate, respectively. Thus deformed serpentinites can present seismic velocities similar to peridotites for wave propagation parallel to the foliation or lower than crustal rocks for wave propagation perpendicular to the foliation. These properties can be used to detect serpentinite, quantify the amount of serpentinization, and to discuss relationships between seismic anisotropy and deformation in the mantle wedge. Regions of high V_P/V_S ratios and extremely low velocities in the mantle wedge of subduction zones (down to about 6 and 3 km.s^?1 for V_P and V_S, respectively) are difficult to explain without strong preferred orientation of serpentine. Local variations of anisotropy may result from kilometer-scale folding of serpentinites. Shear wave splittings up to 1–1.5 s can be explained with moderately thick (10–20 km) serpentinite bodies.     

Optical signatures of dissolved organic matter in the watershed of a globally large river (Yangtze River,China)

Parallel factor analysis of fluorescence excitation emission matrices of surface water samples of a globally large river (Yangtze River, China) watershed identified three classes of fluorescent dissolved organic matter (FDOM) that had ex/em = 280/330 nm, 305/385 nm and 350/450 nm respectively, resembling “peak T”, “peak M” and “peak C” commonly identified in natural water, respectively. Peak T (a tyrosine/tryptophan-like FDOM) did not show correlations to peak M or C which were humic-like substances, while a positive correlation (r = 0.935, p < 0.001) was present between the natural log-transformed maximum fluorescence intensity (F_max) of peaks T and M indicating a tight link during their production and processing. F_max values (in Raman unit nm^?1) normalized to dissolved organic carbon (DOC) concentration were low, varying in ranges 15.93–85.95, 29.83–83.54 and 19.73–51.05 × 10^?5 nm^?1 (μmol/L)^?1 for peaks T, M and C, respectively, in line with the history of strong photobleaching of the water samples as indicated by fairly high absorption spectral slope ratios (0.75–1.53 with a mean 1.03). Intermediate fluorescence index (FI) (1.46–1.83 with a mean 1.61) and small specific absorption at 254 nm (0.64–1.93 with a mean 1.15 m^?1 mg^?1 L) of the water samples, indicated the presence of both aquatic microbial DOM (e.g. peak T) and soil DOM (e.g. peak C). Peak C could be substantially removed by UV-A (320–400 nm) irradiation, while peak M was slightly increased when a microbe-containing water was exposed to the same UV-A irradiation. Taken together, peak C was attributed to diffuse soil source while peak M was likely attributed to joint effects of microbial activities and solar irradiation on the chromophores in the sample.     

Temporal and spatial complexity in post-glacial sedimentation on the tectonically active,Poverty Bay continental margin of New Zealand

On the eastern Raukumara Ranges of the New Zealand East Coast, active tectonics, vigorous weather systems, and human colonisation have combined to cause widespread erosion of the mudstone- and sandstone-dominated hinterland. The Waipaoa River sedimentary dispersal system is an example that has responded to environmental change, and is now New Zealand's second largest river in terms of suspended sediment discharge. This paper presents new sediment accumulation rates for the continental shelf and slope that span century to post-glacial time scales. These data are derived from radiochemical tracer, palynological, tephrostratigraphic, and seismic methods. We hypothesise on the temporal and spatial complexity of post-glacial sedimentation across the margin and identify the broad extent of sediment dispersal from the Waipaoa system. The ∼15 km³ Poverty Bay mid-shelf basin lies adjacent to the mouth of the Waipaoa River, reaching a maximum thickness of ∼45 m. A post-glacial mud lobe of an additional ∼3 km³ extends through the Poverty Gap and out onto the uppermost slope, attaining 40 m thickness in a structurally controlled sub-basin. Here, an offset in the last-glacial erosion surface indicates that deposition was sympathetic with fault activity and the creation of accommodation space, implying that sedimentation was not supply limited. Contrary to classical shelf sedimentation models, the highest modern accumulation rate of 1 cm y⁻¹ occurs on the outer-shelf sediment lobe, approximately ∼2 times the rate recorded at the mid-shelf basin depocentre, and ∼10 times faster than the excess ²¹⁰Pb rates estimated from the slope. Pollen records from slope cores fingerprint Polynesian then European settlement, and broaden the spatial extent of post-settlement sedimentation initially documented from the Poverty Bay mid-shelf. Changes in sub-millennial sedimentation infer a 2–3-times increase in post-settlement accumulation on the shelf but a smaller 1–2 times increase on the slope. Over longer time scales, seismic evidence infers slower but steady sedimentation since the last transgression, and that significant cross-shelf sediment pathways pre-date the increase in sedimentation resulting from colonisation and deforestation. From a summation of coastal bedload, shelf and slope sediment mass accumulation, the total sediment budget for the Holocene is ∼1 Mt y⁻¹. Under modern conditions a larger proportion of the Waipaoa sediment dispersal system likely extends onto the slope and beyond.     

Prediction of ground motion due to the collapse of a large-scale cooling tower under strong earthquakes

The ground motion owing to the collapse of a large-scale cooling tower under strong earthquakes was appropriately predicted using a comprehensive approach. The predicted results can be used for the safety evaluation of nuclear-related facilities adjacent to the cooling tower as well as in the planning of nuclear power plant construction in China. In this study, a cooling tower–soil model was first developed based on a falling weight–soil model, which the authors verified by falling weight tests. Then the collapse process of a cooling tower was simulated, and the collapse-induced ground vibrations were assessed by using the proposed model. Finally, the ground motion, which was a combination of the earthquake-induced ground motion and the collapse-induced ground vibrations, was estimated based on the superposition principle of waves. It was found that the cooling tower may collapse under strong earthquakes with the peak ground accelerations (PGAs) in the range of 0.35–0.45 g in x (EW) and y (NS) directions, respectively. These PGAs are far beyond the PGA range of major earthquakes in the common seismic design in China. The types of the site geologies of towers can significantly affect the collapse-induced ground vibrations. For a typical hard soil consisting of strongly weathered sandy slate, moderate ground vibrations may occur in the considered region. The collapse-induced PGAs were in the range of 0.017–0.046 g for the observed points at distances of 350 m in radial direction. For a rock-like foundation, the collapse-induced radial PGAs may be as high as 0.08 g at distances of 350 m, indicating that the effect of the collapse-induced ground vibrations on the nuclear-related facilities should be seriously assessed in certain scenarios.