The impact of seismic noise produced by wind turbines on seismic borehole measurements

Seismic signals produced by wind turbines can have an adverse effect on seismological measurements up to distances of several kilometres. Based on numerical simulations of the emitted seismic wave field, we study the effectivity of seismic borehole installations as a way to reduce the incoming noise. We analyse the signal amplitude as a function of sensor depth and investigate effects of seismic velocities, damping parameters and geological layering in the subsurface. Our numerical approach is validated by real data from borehole installations affected by wind turbines. We demonstrate that a seismic borehole installation with an adequate depth can effectively reduce the impact of seismic noise from wind turbines in comparison to surface installations. Therefore, placing the seismometer at greater depth represents a potentially effective measure to improve or retain the quality of the recordings at a seismic station. However, the advantages of the borehole decrease significantly with increasing signal wavelength.</p

Anisotropy and XKS-splitting from geodynamic models of double subduction: Testing the limits of interpretation

In this study, we develop three-dimensional geodynamic models to predict XKS-splitting for double subduction scenarios characterized by two outward dipping slabs. These models are highly relevant in various realistic settings, such as the central Mediterranean. We focus on the analysis of XKS-splitting, a key geophysical observable used to infer seismic anisotropy and mantle flow patterns predicted from these geodynamic models. Our geodynamic models simulate the concurrent subduction of two identical oceanic plates which are separated by a continental plate. The variation of the separating plate strength, cause a transition from a retreating to a stationary trench. The models provide detailed insights into the temporal evolution of mantle flow patterns, especially the amount of trench parallel flow, induced by these double subduction scenarios. In a second step, we use the well-known D-Rex model (Kaminski et al., 2004) to efficiently estimate the CPO development in response to plastic deformation produced by mantle flow. Based on the results of the D-Rex model, which includes the full elastic tensor of a deformed multiphase polycrystalline mantle aggregate within the three-dimensional model, we obtain synthetic apparent splitting parameters at receivers placed at the surface by applying multiple-layer anisotropic waveform modeling. Employing analytical techniques, we show the ambiguous nature of apparent splitting parameters, as already suggested by previous studies based on numerical modeling. In the light of the results, we postulate that a meaningful inversion, based on the commonly applied 2-layer anisotropic model, requires additional constraints on fast-axis orientation or strength of anisotropy (delay time). Finally, we show that constraints from our texture simulations (i.e., the integrated delay time) can be used to achieve unique 2-layer models that perfectly fit the synthetic observables. Such models could serve as reference for the interpretation of the observations. Our study highlights the necessity of combining geodynamic modeling and XKS-splitting analysis to shed light on complex upper mantle flow patterns such as those that might occur around subduction zones

Shear wave splitting across the Iceland hot spot: Results from the ICEMELT experiment

We report on observations of upper mantle anisotropy from the splitting of teleseismic shear waves (SKS, SKKS, and PKS) recorded by the ICEMELT broadband seismometer network in Iceland. In a ridge-centered hot spot locale, mantle anisotropy may be generated by flow-induced lattice-preferred orientation of olivine grains or the anisotropic distribution of magma. Splitting measurements of teleseismic shear waves may thus provide diagnostic information on upper mantle flow and/or the distribution of retained melt associated with the Iceland mantle plume. In eastern Iceland, fast polarization directions lie between N10°W and N45°W and average N24°W; delay times between the fast and slow shear waves are generally 0.7–1.35 s. In western Iceland, in contrast, the fast polarization directions, while less well constrained, yield an average value of N23°E and delay times are smaller (0.2–0.95 s). We propose that splitting in eastern Iceland is caused by a 100- to 200-km-thick anisotropic layer in the upper mantle. The observed fast directions in eastern Iceland, however, do not correspond either to the plate spreading direction or to a pattern of radial mantle flow from the center of the Iceland hot spot. We suggest that the relatively uniform direction and magnitude of splitting in eastern Iceland, situated on the Eurasian plate, may therefore reflect the large-scale flow field of the North Atlantic upper mantle. We hypothesize that the different pattern of anisotropy beneath western Iceland, part of the North American plate, is due to the different absolute motions of the two plates. By this view, splitting in eastern and western Iceland is the consequence of shear by North American and Eurasian plate motion relative to the background mantle flow. From absolute plate motion models, in which the Eurasian plate is approximately stationary and the North American plate is moving approximately westward, the splitting observations in both eastern and western Iceland can be satisfied by a background upper mantle flow in the direction N34°W and a velocity of 3 cm/yr in a hot spot reference frame. This inference can be used to test mantle flow models. In particular, it is inconsistent with kinematic flow models, which predict southward flow, or models where flow is dominated by subduction-related sources of mantle buoyancy, which predict westward flow. Our observations are more compatible with the flow field predicted from global seismic tomography models, which in particular include the influence of the large-scale lower mantle upwelling beneath southern Africa. While the hypothesized association between our observations and this upwelling is presently speculative, it makes a very specific and testable prediction about the flow field and hence anisotropy beneath the rest of the Atlantic basin.This work was supported by the National Science Foundation under grants EAR-9316137, OCE-9402991, and EAR-9707193.Peer Reviewe

Crustal and uppermost mantle shear wave velocity structure beneath the Middle East from surface wave tomography

SUMMARY We have constructed a 3-D shear wave velocity (Vs) model for the crust and uppermost mantle beneath the Middle East using Rayleigh wave records obtained from ambient-noise cross-correlations and regional earthquakes. We combined one decade of data collected from 852 permanent and temporary broad-band stations in the region to calculate group-velocity dispersion curves. A compilation of &amp;gt;54 000 ray paths provides reliable group-velocity measurements for periods between 2 and 150 s. Path-averaged group velocities calculated at different periods were inverted for 2-D group-velocity maps. To overcome the problem of heterogeneous ray coverage, we used an adaptive grid parametrization for the group-velocity tomographic inversion. We then sample the period-dependent group-velocity field at each cell of a predefined grid to generate 1-D group-velocity dispersion curves, which are subsequently inverted for 1-D Vs models beneath each cell and combined to approximate the 3-D Vs structure of the area. The Vs model shows low velocities at shallow depths (5–10 km) beneath the Mesopotamian foredeep, South Caspian Basin, eastern Mediterranean and the Black Sea, in coincidence with deep sedimentary basins. Shallow high-velocity anomalies are observed in regions such as the Arabian Shield, Anatolian Plateau and Central Iran, which are dominated by widespread magmatic exposures. In the 10–20 km depth range, we find evidence for a band of high velocities (&amp;gt;4.0 km s–1) along the southern Red Sea and Arabian Shield, indicating the presence of upper mantle rocks. Our 3-D velocity model exhibits high velocities in the depth range of 30–50 km beneath western Arabia, eastern Mediterranean, Central Iranian Block, South Caspian Basin and the Black Sea, possibly indicating a relatively thin crust. In contrast, the Zagros mountain range, the Sanandaj-Sirjan metamorphic zone in western central Iran, the easternmost Anatolian plateau and Lesser Caucasus are characterized by low velocities at these depths. Some of these anomalies may be related to thick crustal roots that support the high topography of these regions. In the upper mantle depth range, high-velocity anomalies are obtained beneath the Arabian Platform, southern Zagros, Persian Gulf and the eastern Mediterranean, in contrast to low velocities beneath the Red Sea, Arabian Shield, Afar depression, eastern Turkey and Lut Block in eastern Iran. Our Vs model may be used as a new reference crustal model for the Middle East in a broad range of future studies.</jats:p

Multiple mantle upwellings in the transition zone beneath the northern East-African Rift system from relative P-wave travel-time tomography

Mantle plumes and consequent plate extension have been invoked as the likely cause of East African Rift volcanism. However, the nature of mantle upwelling is debated, with proposed configurations ranging from a single broad plume connected to the large low-shear-velocity province beneath Southern Africa, the so-called African Superplume, to multiple lower-mantle sources along the rift. We present a new P-wave travel-time tomography model below the northern East-African, Red Sea, and Gulf of Aden rifts and surrounding areas. Data are from stations that span an area from Madagascar to Saudi Arabia. The aperture of the integrated data set allows us to image structures of 100 km length-scale down to depths of 700– 800 km beneath the study region. Our images provide evidence of two clusters of low-velocity structures consisting of features with diameter of 100–200 km that extend through the transition zone, the first beneath Afar and a second just west of the Main Ethiopian Rift, a region with off-rift volcanism. Considering seismic sensitivity to temperature, we interpret these features as upwellings with excess temperatures of 100 6 50 K. The scale of the upwellings is smaller than expected for lower mantle plume sources. This, together with the change in pattern of the low-velocity anomalies across the base of the transition zone, suggests that ponding or flow of deep-plume material below the transition zone may be spawning these upper mantle upwellings

Seismic radiation from wind turbines: observations and analytical modeling of frequency-dependent amplitude decays

In this study, we determine spectral characteristics and amplitude decays of wind turbine induced seismic signals in the far field of a wind farm (WF) close to Uettingen, Germany. Average power spectral densities (PSDs) are calculated from 10 min time segments extracted from (up to) 6 months of continuous recordings at 19 seismic stations, positioned along an 8 km profile starting from the WF. We identify seven distinct PSD peaks in the frequency range between 1 and 8 Hz that can be observed to at least 4 km distance; lower-frequency peaks are detectable up to the end of the profile. At distances between 300 m and 4 km the PSD amplitude decay can be described by a power law with exponent b. The measured b values exhibit a linear frequency dependence and range from b=0.39 at 1.14 Hz to b=3.93 at 7.6 Hz. In a second step, the seismic radiation and amplitude decays are modeled using an analytical approach that approximates the surface wave field. Since we observe temporally varying phase differences between seismograms recorded directly at the base of the individual wind turbines (WTs), source signal phase information is included in the modeling approach. We show that phase differences between source signals have significant effects on the seismic radiation pattern and amplitude decays. Therefore, we develop a phase shift elimination method to handle the challenge of choosing representative source characteristics as an input for the modeling. To optimize the fitting of modeled and observed amplitude decay curves, we perform a grid search to constrain the two model parameters, i.e., the seismic shear wave velocity and quality factor. The comparison of modeled and observed amplitude decays for the seven prominent frequencies shows very good agreement and allows the constraint of shear velocities and quality factors for a two-layer model of the subsurface. The approach is generalized to predict amplitude decays and radiation patterns for WFs of arbitrary geometry