Characterisation of subsurface spatial variability using a cone resistivity penetrometer

The subsurface spatial variation in clay soils, such as thin-layered sand seams, affects the mechanical strength and electrical resistivity. The objective of this study is the development and application of cone resistivity penetrometer (CRP), which measures the cone tip resistance, sleeve friction, and electrical resistivity to evaluate the subsurface spatial variability. The electrical resistivity is measured at the cone tip to increase its resolution. Two outer diameters of the cone resistivity penetrometers (CRPs) are developed: D=10 mm CRP with a projected area of 0.78 cm² and D=15 mm CRP with a projected area of 1.76 cm². The cone tip resistance is effectively separated using a friction sleeve. Strain gauges are used to measure the mechanical strength, and coaxial type electrodes monitor the electrical resistivity. The application tests in the laboratory are conducted using layered soils and saturated sands. In addition, the penetration tests in the field are carried out and compared with the standard piezocone test. The penetration tests show that the soil layers and the density changes are clearly detected by the electrical resistivity and mechanical strength. Field tests show that CRP clearly evaluates the subsurface profile. This study suggests that CRP may be a useful technique for the evaluation of subsurface spatial variability during penetration testing.     

Sea level rise in the north-western part of the Arabian Gulf

Relative sea level variations in the north-western part of the Arabian Gulf have been estimated in the past using no more than 10 to 15 years of observations. In this study, we have almost doubled the period to 28.7 years by examining all available tide gauge data in the area and constructing a mean gauge time-series from seven coastal tide gauges. We found for the period 1979–2007 a relative sea level rise of 2.2 ± 0.5 mm/year. Using the subsidence observed at 6 GPS stations within a radius of 100 km of the tide gauges as an indication of the vertical land motion, the corresponding absolute sea level rise is 1.5 ± 0.8 mm/year that is in agreement with the global estimate of 1.9 ± 0.1 mm/year (Church and White, 2011) for the same studied period. By taking into account the temporal correlations we conclude that previous published results underestimate the true sea level rate uncertainty in this area by a factor of 5–10.     

Sea level time series and ocean tide analysis from multipath signals at five GPS sites in different parts of the world

We present sea level observations derived from the analysis of signal-to-noise ratio (SNR) data recorded by five coastal GPS stations. These stations are located in different regions around the world, both in the northern and in the southern hemisphere, in different multipath environments, from rural coastal areas to busy harbors, and experience different tidal ranges.The recorded SNR data show periodic variations that originate from multipath, i.e. the interference of direct and reflected signals. The general assumption is that for satellite arcs facing the open sea, the rapid SNR variations are due to reflections off the sea surface. The SNR data recorded from these azimuth intervals were analyzed by spectral analysis with two methods: a standard analysis method assuming a static sea level during a satellite arc and an extended analysis method assuming a time dependent sea level during a satellite arc.The GPS-derived sea level results are compared to sea level records from co-located traditional tide gauges, both in the time and in the frequency domain. The sea level time series are highly correlated with correlation coefficients to the order of 0.89–0.99. The root-mean-square (RMS) difference is 6.2 cm for the station with the lowest tidal range of 165 cm and 43 cm for the station with the highest tidal range of 772 cm. The relative accuracy, defined as the ratio of RMS and tidal range, is between 2.4% and 10.0% for all stations.Comparing the standard analysis method and the extended analysis method, the results based on the extended analysis method agree better with the independent tide gauge records for the stations with a high tidal range. For the station with the highest tidal range (772 cm), the RMS is reduced by 47% when using the extended analysis method. Furthermore, the results also indicate that the standard analysis method, assuming a static sea level, can be used for stations with a tidal range of up to about 270 cm, without performing significantly worse than the extended analysis method.Tidal amplitudes and phases are derived by harmonic analysis of the sea level records. Again, a high level of agreement is observed between the tide gauge and the GPS-derived results. Comparing the GPS-derived results, the results based on the extended analysis method show a higher degree of agreement with the traditional tide gauge results for stations with larger tidal ranges. Spectral analysis of the residuals after the harmonic analysis reveals remaining signal power at multiples of the draconitic day. This indicates that the observed SNR data are to some level disturbed by additional multipath signals, in particular for GPS stations that are located in harbors.     

On the electric field at the tip of dart leaders in lightning flashes

The results obtained in this study show that as the dart leader tip passes a given point on the defunct return stroke channel the electric field increases within a fraction of a microsecond to values larger than the critical electric field necessary for the initiation of cold electron runaway in low-density air comprising the channel. These results are in support of the hypothesis that cold runaway electron breakdown may play a role in the emission of X-ray bursts by dart leaders. The calculations also show that the peak power dissipated by a typical dart leader is about 300–500 MW/m and the energy dissipated within the first 10 μs or so is about 500–600 J/m. Furthermore, the minimum resistance and the maximum radius of the core of a typical dart leader are estimated to be about 3 Ω/m and 0.003 m, respectively.     

In situ cosmogenic 10Be production-rate calibration from the Southern Alps,New Zealand

We present a ¹⁰Be production-rate calibration derived from an early Holocene debris-flow deposit at about 1000 m above sea level in the central Southern Alps, New Zealand, in the mid-latitude Southern Hemisphere. Ten radiocarbon ages on macrofossils from a soil horizon buried by the deposit date the deposit to 9690 ± 50 calendar years before AD2008. Surface ¹⁰Be concentrations of seven large boulders partially embedded in the stable surface of the deposit are tightly distributed, yielding a standard deviation of ～2%. Conversion of the ¹⁰Be measurements to sea level/high-latitude values using each of five standard scaling methods indicates ¹⁰Be production rates of 3.84 ± 0.08, 3.87 ± 0.08, 3.83 ± 0.08, 4.15 ± 0.09, and 3.74 ± 0.08 atoms g^?1 a^?1, relative to the ‘07KNSTD’ ¹⁰Be AMS standard, and including only the local time-integrated production-rate uncertainties. When including a sea level high-latitude scaling uncertainty the overall error is ～2.5% (1σ) for each rate. To test the regional applicability of this production-rate calibration, we measured ¹⁰Be concentrations in a set of nearby moraines deposited before 18 060 ± 200 years before AD2008. The ¹⁰Be ages are only consistent with minimum-limiting ¹⁴C age data when calculated using the new production rates. This also suggests that terrestrial in situ cosmogenic-nuclide production did not change significantly from Last Glacial Maximum to Holocene time in New Zealand. Our production rates agree well with those of a recent calibration study from northeastern North America, but are 12–14% lower than other commonly adopted values. The production-rate values presented here can be used elsewhere in New Zealand for rock surfaces exposed during or since the last glacial period.     

First results on determination of cosmogenic 36Cl in limestone from the Yenicekale Complex in the Hittite capital of Hattusha (Turkey)

We have measured ³⁶Cl in three rock surfaces of the Yenicekale building complex in Hattusha (Bo?azköy, Turkey). Hattusha was the capital of Hittite Empire which lasted from about 1650/1600 to 1200 BC. At Yenicekale, Hittite masons flattened the summit of an outcropping limestone knoll to form an artificial platform as the foundation for a building. Next they built a circuit wall along the lateral precipices of the flattened bedrock platform. We took one sample from the limestone bedrock platform and two samples from limestone building blocks of the circuit wall for cosmogenic ³⁶Cl analysis. Calculated exposure ages are 20 ± 1 ka for the sample from the bedrock platform and 24 ± 1 ka and 52 ± 2 ka for the circuit wall blocks. These exposure ages are significantly older than the age expected based on the estimated time of construction between 3.2 ka and 3.7 ka. We conclude that the sampled surfaces contain significant inherited cosmogenic ³⁶Cl. We cannot directly determine exposure ages for the building complex based on these three samples. On the other hand we may use the measured concentrations to determine how much of the rock was removed from the platform during flattening. To this end we modeled the variation of ³⁶Cl production with depth at Yenicekale using the results from the bedrock sample. We conclude that the Hittite masons removed only around 3 m from top of the limestone block. This means that the volume of rock removed from the bedrock platform is significantly less than the volume in the circuit wall atop the platform. They did not gain enough rock from this flattening to make the building. In agreement with this, the first results of our detailed microfacies analysis indicate that many of the building blocks are not of the same facies as the underlying limestone and must have been quarried elsewhere. Although we were not able to exposure date the Yenicekale complex due to the presence of inherited ³⁶Cl, our data suggest that Hittite masons excavated (most of) the building stones not at Yenicekale, but in quarries outside of Hattusha and then transported them to the construction site. These quarries have not yet been identified.     

Optically stimulated luminescence dating of fluvioglacial (sandur) sediments from north-eastern Germany

This study presents first results of a project which aims to develop an absolute chronology of the Weichselian ice marginal positions in north-eastern Germany using Optically Stimulated Luminescence (OSL) dating. Owing to the special geomorphological and stratigraphical setting, many of the sampling sites were situated close to the former ice margin. At some sections the sediments showed an incomplete resetting of the OSL signal. For this study the representative sandur sediment section from the Beelitz outwash cone was chosen to elaborate a strategy to reliably detect insufficient bleaching and to choose an appropriate age model for the calculation of the true burial dose of each sample. The age determined for the formation of the Beelitz outwash cone is 34.1 ± 3.0 ka. Although this age might still be a maximum age, this study shows that by applying the proposed methodological strategy, plausible ages can be derived from insufficiently bleached, ice proximal sandur sediments.     

Coastal and global sea level change

Both coastal and global mean sea level rise by about 3.0 ± 0.5 mm/year from January 1993 to December 2004. Over shorter intervals the coastal sea level rises faster and over longer intervals slowly than the global mean, which trend is almost constant for each interval and is equal to 2.9 ± 0.5 mm/year in 1993–2008. The different trends are due to the higher interannual variability of coastal sea level, caused by the sea level regional variability, that is further averaged out when computing the global mean.Coastal sea level rise is well represented by a selected set of 267 stations of the Permanent Service for Mean Sea Level and by the corresponding co-located altimeter points. Its departure from coastal sea level computed from satellite altimetry in a 150 km distance from coast, dominated by a large rise in the Eastern Pacific, is due to the regional interannual variability.Regionally the trends of the coastal and open-ocean sea level variability are in good agreement and the main world basins have a positive averaged trend. The interannual variability is highly correlated with the El Nino Southern Oscillation (SO) and the North Atlantic Oscillation (NAO) climatic indices over both the altimeter period and the interval 1950–2001. Being the signal of large scale a small number of stations with good spatial coverage is needed. The reconstruction of the interannual variability using the spatial pattern from altimetry and the temporal patterns from tide gauges correlated to NAO and SOI restitutes about 50% of the observed interannual variability over 1993–2001.     

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.     

Infrastructure and evolution of ocean-ridge discontinuities in Iceland

There are two main ocean-ridge discontinuities in Iceland: the Tjörnes Fracture Zone (TFZ) and the South Iceland Seismic Zone (SISZ). The TFZ is a 120-km-long and as much as 70-km-wide WNW-trending zone of high seismicity. It has three main seismic lineaments: the Husavik-Flatey Fault (HFF), the Dalvik lineament, and the Grimsey lineament. The HFF, a dextral strike-slip fault and active as a transform fault for about 9 Ma, has a cumulative transform-parallel displacement of some 60 km. Offshore, the HFF is marked by a transform (fracture-zone) valley, 5–10 km wide and 3–4 km deep. Onshore the Flateyjarskagi Peninsula the HFF is marked by a 3–5-km-wide zone of intense crustal deformation with numerous strike-slip and normal faults, transform-parallel dykes, dense sets of mineral veins, and subzones of completely crushed rocks, that is, fault cores. Where the HFF comes on land on Tjörnes there is a similar, but much thinner, zone of crushed rocks. The seismic lineaments are located a few tens of kilometres south (Dalvik) and north (Grimsey) of, and run subparallel with, the HFF. Both lineaments are composed of sets of NNW-trending sinistral faults arranged en echelon.The SISZ is a 70-km-long and 10–20-km wide zone of almost continuous seismicity located between the overlapping West and East Volcanic Zones. It produces the largest earthquakes in Iceland, some of which exceed M7, during which the N–S width of the zone may be as great as 50–60 km. The SISZ is partly covered with Holocene lava flows where the seismogenic faults occur as dextral NNE-trending and sinistral ENE-trending conjugate arrays with push-ups between their nearby ends. The same fault-segment trends occur in the Pleistocene pile north of the Holocene lava flows.The HFF is neither perpendicular to the nearby ridge segments nor parallel with the spreading vector. As a consequence, the North Volcanic Zone has propagated to the north and the Kolbeinsey Ridge to the south during the past 1 Ma, resulting in the development of the Grimsey and Dalvik lineaments. Similarly, the tip of the East Volcanic Zone has been propagating rapidly to the southwest during the past 3 Ma. The tip has been at its present location for no more than several hundred thousand years, thus making the SISZ less stable than the HFF. If the propagation of the tip of the East Volcanic Zone continues, it will eventually reach the Reykjanes Ridge, whereby either the West or the East Volcanic Zone becomes extinct. Then the SISZ dies out as a major seismic zone.     

Liquefaction assessment of cohesionless soils in the vicinity of large embankments

In a typical seismic dam safety evaluation, standard penetration, cone penetration, Becker penetration, or shear wave velocity (V_s) tests are often first conducted near the toe of an earth dam to infer if any liquefiable soil exists in the foundation of the dam footprint. In current practice, a level-ground condition is commonly assumed when normalizing penetration resistance and V_s, and may be assumed (particularly in preliminary assessments) in applying the cyclic stress method (with or without the K_α correction) to evaluate liquefaction. However, the presence of an earth dam, or any other large embankment or structure, significantly alters the normal and shear stresses in the foundation. This paper identifies and quantifies potential errors in ignoring altered stresses near heavy structures, and presents a methodology to incorporate these effects within the framework of the simplified procedure. Specifically, the effects of these altered stresses (in comparison to the level-ground assumption with and without K_α correction) on the: (1) normalization of field measurements such as penetration resistance and V_s; (2) cyclic stress ratio (CSR); (3) cyclic resistance ratio (CRR); and (4) factor of safety against liquefaction triggering (FS_liq), are evaluated by considering static and dynamic analyses of a generic earthen embankment (60 m high) resting on a saturated, cohesionless foundation (30 m deep). Our analyses indicated that ignoring the presence of induced static shear stresses can result in potentially unconservative errors in overburden correction factors of 30% to 60% at shallow depth (although this error is greatly muted at depths exceeding about 15 m), while errors in CSR potentially can range from about 20% too conservative to 40% unconservative. Potential errors in CRR can approach 50% unconservative at shallow depths, but again, this error is muted at depths exceeding about 15 m. Combining these factors, potentially unconservative errors in computing FS_liq could exceed 100% at shallow depths (less than 15 m to 20 m) while at greater depth (exceeding 20 m) errors approach 20% on the conservative side.     

Post-seismic relaxation: An example of earthquake triggering in the Apennine belt (1915–1920)

The strong earthquake (M = 7) that occurred in the Fucino basin (central Italy) on January 13, 1915 was followed by six earthquakes of M > 5.5 and several other shocks of M > 5 in the major seismic zones of the northern Apennines from 1916 to 1920. This seismicity pattern is consistent with the implications of the present tectonic setting in the study area, which suggests that strong decoupling earthquakes in the central Apennines cause a significant increase of tectonic load, and possibly of seismicity, in the northern Apennines. A numerical simulation, carried out by an elastic-viscous model, of the stress diffusion induced by the Fucino and successive largest earthquakes, shows that each of the above shocks occurred when the respective zone was reached by the highest values of the strain and strain rate perturbation triggered by the previous events. Furthermore, the computed strain regime at each earthquake site is consistent with the known faulting pattern. The results provide important insights into the physical mechanism that controls the interaction of seismic sources in the central and northern Apennines.     

How well do we understand production of 36Cl in limestone and dolomite?

We have evaluated all parameters for the calculation of cosmogenic ³⁶Cl production rates and thus surface exposure ages in dolomite and limestone. We found that we can use either of both published negative muon stopping rates until more information is available. The largest uncertainty of the age estimation in the upper meter of rock comes from the ³⁶Cl production rate from Ca spallation and, in the case of 50–100 ppm Cl content, from the production rate of epithermal neutrons, which we estimate at 760 ± 150 neutrons/g_air/yr (1σ). For a sample with representative amounts of Ca and Cl (20 wt% Ca and 50 ppm Cl, or 40 wt% Ca and 100 ppm Cl), the age can be calculated with a precision of 7–10% in the top 1.5 m of the depth profile. Further improvement of ³⁶Cl calculations depends on new calibration of ³⁶Cl production from Ca spallation, re-evaluation of ³⁶Cl production by low-energy neutron capture on ³⁵Cl, as well as of the muon flux and muon capture based on the most recent measurement data.     

Dike emplacement and flank instability at Mount Etna: Constraints from a poro-elastic-model of flank collapse

Many volcanic edifices are subject to flank failure, usually produced by a combination of events, rather than any single process. From a dynamic point of view, the cause of collapse can be divided into factors that contribute to an increase in shear stress, and factors that contribute to the reduction in the friction coefficient μ of a potential basal failure plane. We study the potential for flank failure at Mount Etna considering a schematic section of the eastern flank, approximated by a wedge-like block. For such geometry, we perform a (steady state) limit equilibrium analysis: the resolution of the forces parallel to the possible basal failure plane allows us to determine the total force acting on the potentially unstable wedge. An estimate of the relative strength of these forces suggests that, in first approximation, the stability is controlled primarily by the balance between block weight, lithostatic load and magmatic forces. Any other force (sea load, hydrostatic uplift, and the uplift due to mechanical and thermal pore-fluid pressure) may be considered of second order. To study the model sensitivity, we let the inferred slope α of the basal surface failure vary between ?10° and 10°, and consider three possible scenarios: no magma loading, magmastatic load, and magmastatic load with magma overpressure. We use error propagation to include in our analysis the uncertainties in the estimates of the mechanics and geometrical parameters controlling the block equilibrium. When there is no magma loading, the ratio between destabilizing and stabilizing forces is usually smaller than the coefficient of friction of the basal failure plane. In the absence of an initiating mechanism, and with the nominal values of the coefficient of friction μ = 0.7 ± 0.1 proposed, the representative wedge will remain stable or continue to move at constant speed. In presence of magmastatic forces, the influence of the lateral restraint decreases. If we consider the magmastatic load only, the block will remain stable (or continue to move at constant speed), unless the transient mechanical and thermal pressurization significantly decrease the friction coefficient, increasing the instability of the flank wedge for α > 5° (seaward dipping decollement). When the magma overpressure contribution is included in the equilibrium analysis, the ratio between destabilizing and stabilizing forces is of the same order or larger than the coefficient of friction of the basal failure plane, and the block will become unstable (or accelerate), especially in the case of the reduction in friction coefficient. Finally, our work suggests that the major challenge in studying flank instability at Mount Etna is not the lack of an appropriate physical model, but the limited knowledge of the mechanical and geometrical parameters describing the block equilibrium.     

The K–Ar Cassignol–Gillot technique applied to western Martinique lavas: A record of Lesser Antilles arc activity from 2 Ma to Mount Pelée volcanism

We present 23 new ages from three volcanic complexes of the Lesser Antilles arc in Martinique Island (French West Indies). These ages obtained with the K–Ar Cassignol–Gillot technique are distributed within the whole Quaternary. They allowed us to reconstruct a detailed history of successive volcanic growth and flank collapse stages. Trois Ilets Volcanism has been active during at least 2 Ma, between 2.35 ± 0.03 Ma and 346 ± 27 ka, with monogenetic volcanoes of basaltic-andesite to andesitic compositions. We here propose that magma mixing, which characterizes this volcanism, could have been initiated between 617 and 346 ka by the activation of arc-parallel and arc-transverse fault systems. Meanwhile, the Carbet complex was active 25 km to the north from 998 ± 14 to 322 ± 6 ka, and was partially destroyed by a flank collapse after 602 ± 10 ka. Together with geochemical data, our ages show that Mount Conil and Mount Pelée volcanoes are parts of the same edifice sharing a single magmatic reservoir. Mount Conil started to emerge before 543 ± 8 ka, and andesites erupted until 127 ± 2 ka, when a flank collapse destroyed the western flank of the edifice, probably triggering the emplacement of Piton Marcel, the last eruption of this first stage. We note that this collapse occurred during the transition from oxygen stages 6 to 5, i.e. during glacial to interglacial change, when eustatic level rapidly increased. After that, and until present, Mount Pelée volcano was built with periods of cone growth intercalated by flank collapse events. We here show that a peak of activity occurred between 550 and 330 ka in western Martinique within the three complexes, which are spaced of 15–25 km. Since 330 ka volcanic activity is limited to the northernmost Mount Conil–Mount Pelée complex. Our data are in agreement with the regional scale observations that the whole recent Lesser Antilles arc was subject to a high volcanic activity since 600 ka, probably linked to an increase in magma production. This permanent establishment of rising magma in regularly spaced batches and tectonically controlled, could explain the individual chemical evolution of each edifice and the different eruptive dynamisms occurring at the same time along the recent arc.     

Crustal deformation before the 2008 Wenchuan MS8.0 earthquake studied using GPS data

We used GPS velocities from approximately 700 stations in western China to study the crustal deformation before the Wenchuan M_S8.0 earthquake. The processing methods included analyses of the strain rate field, inversion of fault locking and the GPS velocity profiles. The GPS strain rate in the E-W direction in the Qinghai-Tibet block shows that extensional deformation was dominant in the western region of the block (west of 92.5° E), while compressive deformation predominated in the eastern region of the block (from 92.5° E to 100° E). On a regional scale, the hypocentral region of the Wenchuan earthquake was located at the edge of an intense compression deformation zone of about 1.9 × 10⁻⁸/a in an east-west direction. The characteristic deformation in the seismogenic fault was compressive with a dextral component. The compression deformation rate was greater in the fault's western region than in its eastern region, and the strain accumulation was very slow on the fault scale. The results of a fault locking inversion show that the locking fraction and slip deficit was greater in the middle-northern section of the seismogenic fault than in the southern section. The GPS velocity profile before the Wenchuan earthquake shows that the compression deformation was smaller than the dextral deformation, which is asymmetrical with respect to the distribution of co-seismic displacement. These deformation characteristics should provide some clues to the Wenchuan earthquake which occurred in the later period of the earthquake cycle.     

Settlements of saturated clean sand deposits in shaking table tests

Shaking table tests were conducted on saturated clean Vietnam sand in the large biaxial laminar shear box (1880 mm×1880 mm×1520 mm) at the National Center for Research on Earthquake Engineering (NCREE), Taiwan. The settlement of sand specimens was measured and evaluated during and after each shaking test. Without liquefaction, the settlement of sand caused by shaking is very small. Significant volume changes occur only when there is liquefaction of sand. The volumetric strain of liquefied sand was calculated according to the measured settlement and the thickness of liquefied sand in the specimen. Relations between volumetric strain after liquefaction and the relative density of saturated clean sand were developed for various shaking durations and earthquake magnitudes. They are not affected by the shaking amplitude, frequency, and direction (one- or multidirectional shaking).     

Elastic anisotropy of ferromagnesian post-perovskite in Earth's D” layer

We report experimental observation of a sizable elastic anisotropy in a polycrystalline sample of ferromagnesian silicate in post-perovskite (ppv) structure. Using a novel composite X-ray transparent gasket to contain and synthesize ppv in a panoramic diamond-anvil cell along with oblique X-ray diffraction geometry, we observed the anisotropic lattice strain and {1 0 0} or {1 1 0} slip-plane texture in the sample at 140 GPa. We deduced the elasticity tensor (c_ij), orientation-dependent compressional wave velocities, polarization-dependent shear-wave velocities, and the velocity anisotropy of the silicate ppv. Our results are consistent with calculations and indicate that with sufficient preferred orientation, the elastic anisotropy of this phase can produce large shear-wave splitting.     

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).     

Astrochronology of the late Jurassic Kimmeridge Clay (Dorset,England) and implications for Earth system processes

The Late Jurassic Kimmeridge Clay Formation (KCF) is an economically important, organic-rich source rock of Kimmeridgian–Early Tithonian age. The main rock types of the KCF in Dorset, UK, include grey to black laminated shale, marl, coccolithic limestone, and dolostone, which occur with an obvious cyclicity at astronomical timescales. In this study, we examine two high-resolution borehole records (Swanworth Quarry 1 and Metherhills 1) obtained as part of a Rapid Global Geological Events (RGGE) sediment drilling project. Datasets examined were total organic carbon (TOC), and borehole wall microconductivity by Formation Microscanner (FMS). Our intent is to assess the rhythmicity of the KCF with respect to the astronomical timescale, and to discuss the results with respect to other key Late Jurassic geological processes. Power spectra of the untuned data reveal a hierarchy of cycles throughout the KCF with ～ 167 m, ～ 40 m, 9.1 m, 3.8 m and 1.6 m wavelengths. Tuning the ～ 40 m cycles to the 405-kyr eccentricity cycle shows the presence of all the astronomical parameters: eccentricity, obliquity, and precession index. In particular, ～ 100-kyr and 405-kyr eccentricity cycles are strongly expressed in both records. The 405-kyr eccentricity cycle corresponds to relative sea-level changes inferred from sequence stratigraphy. Intervals with elevated TOC are associated with strong obliquity forcing. The 405-kyr-tuned duration of the lower KCF (Kimmeridgian Stage) is 3.47 Myr, and the upper KCF (early part of the Tithonian Stage, elegans to fittoni ammonite zones) is 3.32 Myr. Two other chronologies test the consistency of this age model by tuning ～ 8–10 m cycles to 100-kyr (short eccentricity), and ～ 3–5 m cycles to 36-kyr (Jurassic obliquity). The ‘obliquity-tuned’ chronology resolves an accumulation history for the KCF with a variation that strongly resembles that of Earth's orbital eccentricity predicted for 147.2 Ma to 153.8 Ma. There is evidence for significant non-deposition (up to 1 million years) in the lowermost KCF (baylei–mutabilis zones), which would indicate a Kimmeridgian/Oxfordian boundary age of 154.8 Ma. This absolute calibration allows assignment of precise numerical ages to zonal boundaries, sequence surfaces, and polarity chrons of the lower M-sequence.