Wolfgang Rabbel

Deep low-frequency earthquakes in tremor localize to the plate interface in multiple subduction zones

Deep tremor under Shikoku, Japan, consists primarily, and perhaps entirely, of swarms of low-frequency earthquakes (LFEs) that occur as shear slip on the plate interface. Although tremor is observed at other plate boundaries, the lack of cataloged low-frequency earthquakes has precluded a similar conclusion about tremor in those locales. We use a network autocorrelation approach to detect and locate LFEs within tremor recorded at three subduction zones characterized by different thermal structures and levels of interplate seismicity: southwest Japan, northern Cascadia, and Costa Rica. In each case we find that LFEs are the primary constituent of tremor and that they locate on the deep continuation of the plate boundary. This suggests that tremor in these regions shares a common mechanism and that temperature is not the primary control on such activity. Citation: Brown, J.R., G. C. Beroza, S. Ide, K. Ohta, D. R. Shelly, S. Y. Schwartz, W. Rabbel, M. Thorwart, and H. Kao (2009), Deep low-frequency earthquakes in tremor localize to the plate interface in multiple subduction zones, Geophys. Res. Lett., 36, L19306, doi:10.1029/2009GL040027.

A tremor and slip event on the Cocos-Caribbean subduction zone as measured by a global positioning system (GPS) and seismic network on the Nicoya Peninsula, Costa Rica

In May 2007 a network of global positioning systems (GPS) and seismic stations on the Nicoya Peninsula, of northern Costa Rica, recorded a slow-slip event accompanied by seismic tremor. The close proximity of the Nicoya Peninsula to the seismogenic part of the Cocos-Caribbean subduction plate boundary makes it a good location to study such events. Several centimeters of southwest motion were recorded by the GPS stations over a period of several days to several weeks, and the seismic stations recorded three distinct episodes of tremor during the same time span. Inversion of the surface displacement data for the depth and pattern of slip on the plate interface shows peak slip at a depth of 25–30 km, downdip of the main seismogenic zone. Estimated temperatures here are ∼250°–300°C, lower than in other subduction zones where events of this nature have been previously identified. There may also be a shallower patch of slip at ∼6 km depth. These results are significant in that they are the first to suggest that slow slip can occur at the updip transition from stick slip to stable sliding, and that a critical temperature threshold is not required for slow slip. Tremor and low-frequency earthquake locations are more difficult to determine. Our results suggest they occur on or near the plate interface at the same depth range as the deep slow slip, but not spatially colocated.

Impacts of the use of the geological subsurface for energy storage: an investigation concept

New methods and technologies for energy storage are required to make a transitionto renewable energy sources; in Germany this transition is termed “Energiewende”. Subsurface georeservoirs, such as salt caverns for hydrogen, compressed air, and methane storage or porous formations for heat and gas storage, offer the possibility of hosting large amounts of energy. When employing these geological storage facilities, an adequate system and process understanding is essential in order to characterize and to predict the complex and interacting effects on other types of subsurface use and on protected entities. In order to make optimal use of georeservoirs, a comprehensive use planning of the subsurface is required that allocates specific uses to appropriate subsurface locations. This paper presents a generic methodology on how subsurface use planning can be conducted and how its scientific basis can be developed. Although synthetic, realistic scenarios for the use of the geological underground for energy storage are parameterized and numerically simulated, accounting for other kinds of subsurface use already in place. From these scenario analyses, the imposed coupled hydraulic, thermal, mechanical and chemical processes, as well as mutual effects and influences on protected entities are assessed and generalized. Based on these, a first methodology for large-scale planning of the geological subsurface considering different surface and subsurface usage scenarios may also be derived.

P and S velocity structure of the crust and the upper mantle beneath central Java from local tomography inversion

Here we present the results of local source tomographic inversion beneath central Java. The data set was collected by a temporary seismic network. More than 100 stations were operated for almost half a year. About 13,000 P and S arrival times from 292 events were used to obtain three-dimensional (3-D) Vp, Vs, and Vp/Vs models of the crust and the mantle wedge beneath central Java. Source location and determination of the 3-D velocity models were performed simultaneously based on a new iterative tomographic algorithm, LOTOS-06. Final event locations clearly image the shape of the subduction zone beneath central Java. The dipping angle of the slab increases gradually from almost horizontal to about 70°. A double seismic zone is observed in the slab between 80 and 150 km depth. The most striking feature of the resulting P and S models is a pronounced low-velocity anomaly in the crust, just north of the volcanic arc (Merapi-Lawu anomaly (MLA)). An algorithm for estimation of the amplitude value, which is presented in the paper, shows that the difference between the fore arc and MLA velocities at a depth of 10 km reaches 30% and 36% in P and S models, respectively. The value of the Vp/Vs ratio inside the MLA is more than 1.9. This shows a probable high content of fluids and partial melts within the crust. In the upper mantle we observe an inclined low-velocity anomaly which links the cluster of seismicity at 100 km depth with MLA. This anomaly might reflect ascending paths of fluids released from the slab. The reliability of all these patterns was tested thoroughly.

Modeling, parameterization and evaluation of monitoring methods for CO2 storage in deep saline formations: the CO2-MoPa project

Capture and geological sequestration of CO2 from large industrial sources is considered a measure for reducing anthropogenic emissions of CO2 and thus mitigating climate change. One of the main storage options proposed are deep saline formations, as they provide the largest potential storage capacities among the geologic options. A thorough assessment of this type of storage site therefore is required. The CO2-MoPa project aims at contributing to the dimensioning of CO2 storage projects and to evaluating monitoring methods for CO2 injection by an integrated approach. For this, virtual, but realistic test sites are designed geometrically and fully parameterized. Numerical process models are developed and then used to simulate the effects of a CO2 injection into the virtual test sites. Because the parameterization of the virtual sites is known completely, investigation as well as monitoring methods can be closely examined and evaluated by comparing the virtual monitoring result with the simulation. To this end, the monitoring or investigation method is also simulated, and the (virtual) measurements are recorded and evaluated like real data. Application to a synthetic site typical for the north German basin showed that pressure response has to be evaluated taking into account the layered structure of the storage system. Microgravimetric measurements are found to be promising for detecting the CO2 phase distribution. A combination of seismic and geoelectric measurements can be used to constrain the CO2 phase distribution for the anticline system used in the synthetic site.

Shear wave anisotropy of laminated lower crust beneath Urach (SW Germany): a comparison with xenoliths and with exposed lower crustal sections

In order to investigate the laminated lower crust beneath Urach (SW Germany), we analysed (1) the anisotropy and reflectivity of seismic shear waves, and (2) the composition of xenoliths samples. We conclude that the directional dependence of shear waves of the lower crust is effectively transversely isotropic with a maximum velocity difference of 3 to 6%. Both seismic signature and petrological information are compared with reference profiles deduced from the exposed lower crustal sections of Ivrea and Calabria (Italy). The comparison, based on synthetic seismograms, shows that the observed combination of shear wave anisotropy and reflectivity pattern corresponds to alternating layers containing a high amount of metapelites, such as found in the Calabria profile and in the Val Strona subsection of Ivrea. The main portion of the whole rock anisotropy is caused by the lattice-preferred orientation of anisotropic minerals (so-called intrinsic anisotropy). The effect of alternating rock layering (so-called layering anisotropy) contributes only a minor portion to the whole rock anisotropy. Layered mafic intrusions, such as found in the Val Sesia subsection of Ivrea, are nearly isotropic. Both seismic anisotropy and xenoliths indicate that mafic layers do not play a major role for the laminated lower crust at Urach. They may be restricted to a thin layer at the crust-mantle boundary.

Superdeep vertical seismic profiling at the KTB deep drill hole (Germany): Seismic close-up view of a major thrust zone down to 8.5 km depth

The lowermost section of the continental superdeep drill hole German Continental Deep Drilling Program (KTB) ( south Germany) has been investigated for the first time by vertical seismic profiling (VSP). The new VSP samples the still accessible range of 6 - 8.5 km depth. Between 7 and 8.5 km depth, the drill hole intersects a major cataclastic fault zone which can be traced back to the Earths surface where it forms a lineament of regional importance, the Franconian line. To determine the seismic properties of the crust in situ, in particular within and around this deep fault zone, was one of the major goals of the VSP. For the measurements a newly developed high-pressure/high-temperature borehole geophone was used that was capable of withstanding temperatures and pressures up to 260 degreesC and 140 MPa, respectively. The velocity-depth profiles and reflection images resulting from the VSP are of high spatial resolution due to a small geophone spacing of 12.5 m and a broad seismic signal spectrum. Compared to the upper part of the borehole, we found more than 10% decrease of the P wave velocity in the deep, fractured metamorphic rock formations. P wave velocity is similar to 5.5 km/s at 8.5 km depth compared to 6.0 - 6.5 km/s at more shallow levels above 7 km. In addition, seismic anisotropy was observed to increase significantly within the deep fracture zone showing more than 10% shear wave splitting and azimuthal variation of S wave polarization. In order to quantify the effect of fractures on the seismic velocity in situ we compared lithologically identical rock units at shallow and large depths: Combining seismic velocity and structural logs, we could determine the elastic tensors for three gneiss sections. The analysis of these tensors showed that we need fracture porosity in the percent range in order to explain seismic velocity and anisotropy observed within the fault zone. The opening of significant pore space around 8 km depth can only be maintained by differential tectonic stress combined with intense macroscopic fracturing. VSP reflection imaging based on PP and PS converted reflected waves showed that the major fault system at the KTB site is wider and more complex than previously known. The so-called SE1 reflection previously found in two- and three-dimensional surface seismic surveys corresponds to the top of an similar to1 km wide fault system. Its lower portion was not illuminated by surface seismic acquisition geometry. VSP imaging shows that the fault zone comprises two major and a number of smaller SE dipping fault planes and several conjugate fracture planes. The previously recognized upper fault plane is not associated with a strong velocity anomaly but indicates the depth below which the dramatic velocity decrease starts. Regarding the complexly faulted crustal section of the KTB site as a whole, we found that fluctuation spectra of rock composition and seismic velocity show similar patterns. We could verify that a significant amount of P wave energy is continuously converted into shear energy by forward scattering and that multipathing plays an important role in signal formation. The media behaves effectively smoothly only at wavelength larger than 150 m. It was shown by moving source profiling that the media is orthorhombic on are regional scale. The tilt of the symmetry axes of anisotropy varies with depth following the dip of the geological structure.

Seismic anisotropy at the continental deep drilling site (Germany)

Abstract The existence of seismic anisotropy at the Continental Deep Drilling Site (KTB, Oberpfalz, Germany) has been established by the DEKORP ISO 89 borehole experiments. The anisotropy, quantified in-situ for a sequence of gneissic rocks, comprises two parts: a background component reflecting mineral composition and its state of ordering (foliation), and a superimposed component originating from fractures. The anisotropy values obtained in the direct vicinity of the KTB pilot hole, however, cannot be extrapolated into its surroundings without modification. Whereas metamorphic rocks in the immediate vicinity of the KTB exhibit a relatively strong anisotropy (> 10%) with an approximately NW-SE-striking axis of high velocity, the anisotropy of an adjacent granite area is considerably lower and of different strike direction (roughly W-E). The directions of high seismic velocities can be correlated with strike directions of geologically mapped fracture zones and rock units as well as to estimates of the direction of maximum horizontal tectonic stress.

Rayleigh-to-shear wave conversion at the tunnel face — From 3D-FD modeling to ahead-of-drill exploration

For safe tunnel excavation, it is important to predict lithologic and structural heterogeneities ahead of construction. Conventional tunnel seismic prediction systems utilize body waves (P- and S-waves) that are directly generated at the tunnel walls or near the cutter head of the tunnel boring machine (TBM). We propose a new prediction strategy that has been discovered by 3D elastic finite-difference (FD) modeling: Rayleigh waves arriving at the front face of the tunnel are converted into high-amplitude S-waves propagating further ahead. Reflected or backscattered S-waves are converted back into Rayleigh waves which can be recorded along the sidewalls. We name these waves RSSR waves. In our approach, the front face acts as an S-wave transceiver. One technical advantage is that both the sources and the receivers may be placed behind the cutter head of the TBM. The modeling reveals that the RSSR waves exhibit significantly higher amplitudes than the directly reflected body waves. The excavation damage zone causes dispersion of the RSSR wave leading to multimodal reflection response. For the detection of geologic interfaces ahead, RSSR waves recorded along the sidewalls are corrected for dispersion and stacked. From the arrival times, the distance to the S-S reflection point can be estimated. A recurrent application, while the tunnel approaches the interface, allows one to quantify the orientation of the reflecting interfaces as well. Our approach has been verified successfully in a field experiment at the Piora adit of the Gotthard base tunnel. The distance to the Piora fault zone estimated from stacked RSSR events agrees well with the information obtained by geologic surveying and exploratory drilling.

The steeply subducting edge of the Cocos Ridge: Evidence from receiver functions beneath the northern Talamanca Range, south-central Costa Rica

The deep structure of the south-central Costa Rican subduction zone has not been studied in great detail so far because large parts of the area are virtually inaccessible. We present a receiver function study along a transect of broadband seismometers through the northern flank of the Cordillera de Talamanca (south Costa Rica). Below Moho depths, the receiver functions image a dipping positive conversion signal. This is interpreted as the subducting Cocos Plate slab, compatible with the conversions in the individual receiver functions. In finite difference modeling, a dipping signal such as the one imaged can only be reproduced by a steeply (80°) dipping structure present at least until a depth of about 70–100 km; below this depth, the length of the slab cannot be determined because of possible scattering effects. The proposed position of the slab agrees with previous results from local seismicity, local earthquake tomography, and active seismic studies, while extending the slab location to greater depths and steeper dip angle. Along the trench, no marked change is observed in the receiver functions, suggesting that the steeply dipping slab continues until the northern flank of the Cordillera de Talamanca, in the transition region between the incoming seamount segment and Cocos Ridge. Considering the time predicted for the establishment of shallow angle underthrusting after the onset of ridge collision, the southern Costa Rican subduction zone may at present be undergoing a reconfiguration of subduction style, where the transition to shallow underthrusting may be underway but still incomplete.