Quasi-2-Day Signals Observed in the Warm Pool Region during TOGA/COARE IOP

During Tropical Ocean and Global Atmosphere (TOGA)/Coupled Ocean and Atmosphere Research Experiment (COARE) Intensive Observing Period (IOP), upward-looking acoustic Doppler current profilers (ADCP) and current meters were moored at two equatorial sites (147°E and 154°E) and two off-equatorial sites (2°N and 2°S, 156°E) in the warm pool region of the western equatorial Pacific. Using current data obtained by these moorings, we have shown that there is a dominant signal with a period of about 2 days from the end of November to the middle of December in 1992, except at the equatorial site on 147°E (Ueki et al., 1998). The energy of this quasi-2-day signal for the meridional current is larger than that for the zonal one and the signal has a high coherence between two off-equatorial sites. In this paper, using band-passed time series of the meridional curerent, we investigated characters of the quasi-2-day signal and attempted to interpret this signal as an equatorially trapped wave. Complex empirical orthogonal function (CEOF) analysis reveals two different phase propagating features between the equatorial and off-equatorial site. One is an upward propagating signal, which is dominant near the surface at two off-equatorial sites, and the other is a downward propagating signal, which is dominant near 200 m at the equatorial site. If one interprets the quasi-2-day signal as an equatorially trapped wave, it is suggested that it cannot be explained as a single wave but can be described as the superimposition of several wave signals. The main part of these signals consists of two signals, one caused by a meteorological forcing and another by another factor in the ocean field.     

Mechanism of the M<Subscript><Emphasis Type="Italic">L</Emphasis></Subscript>4.0 25 April 2016 earthquake in southwest of France in the vicinity of the Lacq gas field

The source mechanism of the M_L 4.0 25 April 2016 Lacq earthquake (Aquitaine Basin, South-West France) is analyzed from the available public data and discussed with respect to the geometry of the nearby Lacq gas field. It is one of the biggest earthquakes in the area in the past few decades of gas extraction and the biggest after the end of gas exploitation in 2013. The routinely obtained location shows its hypocenter position inside the gas reservoir. We first analyze its focal mechanism through regional broad-band seismograms recorded in a radius of about 50 km epicentral distances and obtain EW running normal faulting above the reservoir. While the solution is stable using regional data only, we observe a large discrepancy between the recorded data on nearby station URDF and the forward modeling up to 1 Hz. We then look for the best epicenter position through performing wave propagation simulations and constraining the potential source area by the peak ground velocity (PGV). The resulting epicentral position is a few to several km away to the north or south direction with respect to station URDF such that the simulated particle motions are consistent with the observation. The initial motion of the seismograms shows that the epicenter position in the north from URDF is preferable, indicating the north-east of the Lacq reservoir. This study is an application of full waveform simulations and characterization of near-field ground motion in terms of an engineering factor such as PGV. The finally obtained solution gives a moment magnitude of M_w 3.9 and the best focal depth of 4 km, which corresponds to the crust above the reservoir rather than its interior. This position is consistent with the tendency of Coulomb stress change due to a compaction at 5 km depth in the crust. Therefore, this earthquake can be interpreted as a relaxation of the shallow crust due to a deeper gas reservoir compaction so that the occurrence of similar events cannot be excluded in the near future. It would be necessary to continue monitoring such local induced seismicity in order to better understand the reservoir/overburden behavior and better assess the local seismic hazard even after the end of gas exploitation.     

Using Estimated Risk to Develop Stimulation Strategies for Enhanced Geothermal Systems

Enhanced geothermal systems (EGS) are an attractive source of low-carbon electricity and heating. Consequently, a number of tests of this technology have been made during the past couple of decades, and various projects are being planned or under development. EGS work by the injection of fluid into deep boreholes to increase permeability and hence allow the circulation and heating of fluid through a geothermal reservoir. Permeability is irreversibly increased by the generation of microseismicity through the shearing of pre-existing fractures or fault segments. One aspect of this technology that can cause public concern and consequently could limit the widespread adoption of EGS within populated areas is the risk of generating earthquakes that are sufficiently large to be felt (or even to cause building damage). Therefore, there is a need to balance stimulation and exploitation of the geothermal reservoir through fluid injection against the pressing requirement to keep the earthquake risk below an acceptable level. Current strategies to balance these potentially conflicting requirements rely on a traffic light system based on the observed magnitudes of the triggered earthquakes and the measured peak ground velocities from these events. In this article we propose an alternative system that uses the actual risk of generating felt (or damaging) earthquake ground motions at a site of interest (e.g. a nearby town) to control the injection rate. This risk is computed by combining characteristics of the observed seismicity of the previous 6 h with a (potentially site-specific) ground motion prediction equation to obtain a real-time seismic hazard curve; this is then convolved with the derivative of a (potentially site-specific) fragility curve. Based on the relation between computed risk and pre-defined acceptable risk thresholds, the injection is increased if the risk is below the amber level, decreased if the risk is between the amber and red levels, or stopped completely if the risk is above the red level. Based on simulations using a recently developed model of induced seismicity in geothermal systems, which is checked here using observations from the Basel EGS, in this article it is shown that the proposed procedure could lead to both acceptable levels of risk and increased permeability.     

Effects of load irregularity on the cyclic behaviour of sand

In an effort to investigate the effects of irregular nature of load application on the behaviour of soils during earthquakes, a series of uni-directional irregular loading tests were performed on saturated Fuji river snad under undrained condition using a simple shear test apparatus. The samples prepared to three different densities in this test device were subjected to twelve irregular time histories of shear stress alteration. For comparison purposes, the conventional type of constant amplitude cyclic simple shear tests were also conducted on samples prepared under identical conditions. On the basis of the test results, values of coefficients allowing for the effect of load irregularity were elucidated for the samples with three different densities. It was then suggested that the values of coefficient thus determined be used to facilitate the evaluation of the soil strength under actual seismic loading conditions, on the basis of the cyclic shear strength obtained from the conventional type of constant-amplitude triaxial tests.     

Deep Flows on the Slope Inshore of the Kuril-Kamchatka Trench Southeast off Cape Erimo, Hokkaido

Deep flows on the slope inshore of the Kuril-Kamchatka Trench southeast off Cape Erimo, Hokkaido were observed for about five years from June 1989 to March 1995, using a mooring system with two current meters. In 1991 and 1993 directionally stable southwestward flows were observed at the upper layer (1000 m). These appear to be typical of the Oyashio because the characteristics of the flows were high mean kinetic energy, low eddy energy and high stability. However, the magnitudes of other mean flows at the upper layer, except for 1991 and 1993, were less than their standard deviations. This suggests that the Oyashio was observed for only a limited period of time. On the other hand, at the lower layer (3000 m) the magnitudes of the mean flows for 10–11 months were 1–3 cm s^-1 and ellipses of their eddy kinetic energy were extremely flattened in the direction of the local isobath. The directions of the mean flows in 1990, 1991 and 1993 were southwestward along the local isobath. The relationships between the upper and the lower flows are discussed in terms of monthly change of kinetic energy, since the low-frequency fluctuations longer than 30-day are predominant from the eddy kinetic energy spectra. The results show that there are cases when the kinetic energy of the monthly mean flows at the lower layers are larger than those at the upper layers. This suggests the possibility that the lower flows are in part a southward deep western boundary current.     

Wakefieldite‐(Nd), a New Neodymium Vanadate Mineral in the Arase Stratiform Ferromanganese Deposit,Kochi Prefecture,Japan

Wakefieldite‐(Nd), NdVO₄, is a new mineral found from the Arase stratiform ferromanganese deposit in Kochi Prefecture, Shikoku Island, Japan. It is the Nd‐dominant analogue of wakefieldite‐(Y) and wakefieldite‐(Ce). The ferromanganese ore specimen mainly consists of hematite and caryopilite, and wakefieldite‐(Nd) is typically enclosed in caryopilite. Wakefieldite‐(Nd) is tetragonal, I4₁/amd, a = 7.338(16) Å, c = 6.509(19) Å, V = 350.5(18) ų, Z = 4. The four strongest lines in the X‐ray diffraction pattern [d(Å), I/I₀, hkl] using a Gandolfi camera are (3.67, 100, 200); (2.74, 51, 112); (4.84, 27, 101) and (1.89, 25, 312). Chemical composition of wakefieldite‐(Nd) are V₂O₃ 35.25, As₂O₃ 0.93, SiO₂ 0.14, MnO 1.45, Fe₂O₃ 0.41, Y₂O₃ 2.87, La₂O₃ 7.61, Ce₂O₃ 7.37, Pr₂O₃ 6.04, Nd₂O₃ 26.79, Sm₂O₃ 4.41, Eu₂O₃ 1.36, Gd₂O₃ 3.41, Tb₂O₃ 0.22, Dy₂O₃ 1.41, Er₂O₃ 0.10, total 99.77 wt.%. The empirical formula is (Nd_0.403La_0.118Ce_0.114Pr_0.093Y_0.064Sm_0.064Mn_0.052Gd_0.048Eu_0.020Dy_0.019Fe_0.013Tb_0.003Er_0.001)_1.012(V_0.981As_0.020Si_0.006)_1.007O₄ on the basis of O = 4. The calculated density is 4.782 g/cm³. Microtexture and co‐existing relationship between wakefieldite‐(Nd) and caryopilite suggest that recrystallization and dehydration of Fe‐ and Mn‐oxyhydroxide led to the generation of hematite, caryopilite, rhodochrosite and wakefieldite‐(Nd) by the metamorphism during the accretion of the host unit of the Arase deposit. Chondrite‐normalized REE pattern of the host ferromanganese ore, which is regarded as oceanic metalliferous sediment in origin, shows negative Ce anomaly. Chemical composition of wakefieldite‐(Nd) reflects Ce‐depleted bulk composition of REE‐enriched ferromanganese ore.     

Finite Difference Simulations of Seismic Wave Propagation for the 2007 Mw 6.6 Niigata-ken Chuetsu-Oki Earthquake: Validity of Models and Reliable Input Ground Motion in the Near-Field

Finite difference simulations of seismic wave propagation are performed in the Niigata area, Japan, for the 2007 Mw 6.6 Niigata-ken Chuetsu-Oki earthquake at low frequencies. We test three 3D structural models built independently in various studies. First aftershock simulations are carried out. The model based on 3D tomography yields correct body waves in the near field, but later phases are imperfectly reproduced due to the lack of shallow sediment layers; other models based on various 1D/2D profiles and geological interpretation provide good site responses but generate seismic phases that may be shifted from those actually observed. Next, for the mainshock simulations, we adopt two different finite source models that differ in the near-field ground motion, especially above the fault plane (but under the sea) and then along the coastline. Each model is found to be calibrated differently for the given stations. For engineering purposes, the variations observed in simulated ground motion are significant, but for seismological purposes, additional parameter calibrations would be possible for such a complex 3D case.