Seismic monitoring of EGS tests at the Coso Geothermal area, California, using accurate MEQ locations and full moment tensors

We studied high-resolution relative locations and full moment tensors of microearthquakes (MEQs) occurring before, during and following Enhanced Geothermal Systems (EGS) experiments in two wells at the Coso geothermal area, California. The objective was to map new fractures, determine the mode and sense of failure, and characterize the stress cycle associated with injection. New software developed for this work combines waveform crosscorrelation measurement of arrival times with relative relocation methods, and assesses confidence regions for moment tensors derived using linearprogramming methods. For moment tensor determination we also developed a convenient Graphical User Interface (GUI), to streamline the work. We used data from the U.S. Navy’s permanent network of three-component digital borehole seismometers and from 14 portable three-component digital instruments. The latter supplemented the permanent network during injection experiments in well 34A-9 in 2004 and well 34-9RD2 in 2005. In the experiment in well 34A-9, the co-injection earthquakes were more numerous, smaller, more explosive and had more horizontal motion, compared with the pre-injection earthquakes. In the experiment in well 34-9RD2 the relocated hypocenters reveal a well-defined planar structure, 700 m long and 600 m high in the depth range 0.8 to 1.4 km below sea level, striking N 20° E and dipping at 75° to the WNW. The moment tensors show that it corresponds to a mode I (opening) crack. For both wells, the perturbed stress state near the bottom of the well persisted for at least two months following the injection

Three-dimensional seismic ray tracing

Two methods for tracing seismic rays between 2 given end points through three dimensional, continuously varying velocity structures are available. This paper describes and compares them for problems of practical interest and for analytical ray paths through an idealized velocity structure. One method involves "shooting" the ray from one point with a given starting direction and then modifying this starting direction until the ray emerges at the desired target, while the other method involves "bending" an initial path between the end points until it satisfies the principle of stationary time. For most of the models investigated, "bending" is computationally faster than "shooting" by a factor of 10 or more. The "bending" method can be modified to deal with discontinuities in the velocity model, and can also be adapted for use in conjunction with a table of distances as a function of ray parameter when the three dimensional anomaly influences only a small fraction of the total ray path. The geometrical spreading effect on the amplitude of the ray may be retrieved easily from the "bending" solution.           ARK: https://n2t.net/ark:/88439/y041002 Permalink: https://geophysicsjournal.com/article/133 &nbsp

HiQuake: The Human-Induced Earthquake Database

HiQuake—The Human‐Induced Earthquake Database is the most complete database of anthropogenic projects proposed, on scientific grounds, to have induced earthquake sequences. It is freely available to download from the website given in Data and Resources. At the time this article was written, HiQuake contained ∼730∼730 anthropogenic projects proposed to have induced earthquakes, as well as associated project‐related and seismic data. The most commonly reported anthropogenic activities proposed to have induced earthquakes are mining and water reservoir impoundment. In recent years, the number of earthquake sequences proposed to have been induced by fluid‐injection activities has grown. The most commonly reported maximum observed magnitude in an induced earthquake sequence is 3≤MMAX<43≤MMAX<4 . The largest earthquake in HiQuake proposed to have been induced had a magnitude of MwMw 7.9 and occurred in China. Such large earthquakes release mostly stress of natural tectonic origin, but are conceivably triggered by small anthropogenic stress changes. The data in HiQuake are of variable quality because they are drawn from publications that span almost a century. We estimate underreporting to be ∼30%∼30% for M∼4M∼4 events, ∼60%∼60% for M∼3M∼3 events, and ∼90%∼90% for M∼2M∼2 events. The degree of certitude that the given earthquake sequences were anthropogenically induced is variable. HiQuake includes all earthquake sequences proposed on scientific grounds to have been human induced without regard to the strength of the case made. HiQuake is offered freely as a resource to interested parties, and judging the reliability of any particular case is the responsibility of the database user. HiQuake will be routinely updated to correct errors, update existing entries, and add new entries. It has the potential to help improve our understanding of induced earthquakes and to manage their impact on society

Time-dependent seismic tomography and its application to the Coso geothermal area, 1996-2006

Measurements of temporal changes in Earth structure are commonly determined using localearthquake tomography computer programs that invert multiple seismic-wave arrival time data sets separately and assume that any differences in the structural results arise from real temporal variations. This assumption is dangerous because the results of repeated tomography experiments would differ even if the structure did not change, simply because of variation in the seismic ray distribution caused by the natural variation in earthquake locations. Even if the source locations did not change (if only explosion data were used, for example), derived structures would inevitably differ because of observational errors. A better approach is to invert multiple data sets simultaneously, which makes it possible to determine what changes are truly required by the data. This problem is similar to that of seeking models consistent with initial assumptions, and techniques similar to the “damped least squares” method can solve it. We have developed a computer program, dtomo, that inverts multiple epochs of arrival-time measurements to determine hypocentral parameters and structural changes between epochs. We shall apply this program to data from the seismically active Coso geothermal area, California, in the near future. The permanent network operated there by the US Navy, supplemented by temporary stations, has provided excellent earthquake arrival-time data covering a span of more than a decade. Furthermore, structural change is expected in the area as a result of geothermal exploitation of the resource. We have studied the period 1996 through 2006. Our results to date using the traditional method show, for a 2-km horizontal grid spacing, an irregular strengthening with time of a negative VP/VS anomaly in the upper ~ 2 km of the reservoir. This progressive reduction in VP/VS results predominately from an increase of VS with respect to VP. Such a change is expected to result from effects of geothermal operations such as decreasing fluid pressure and the drying of argillaceous minerals such as illite

Microearthquake moment tensors from the Coso Geothermal area

The Coso geothermal area, California, has produced hot water and steam for electricity generation for more than 20 years, during which time there has been a substantial amount of microearthquake activity in the area. Seismicity is monitored by a high-quality permanent network of 16 three-component digital borehole seismometers operated by the US Navy and supplemented by a ~ 14-station portable array of surface three-component digital instruments. The portable stations improve focal sphere coverage, providing seismic-wave polarity and amplitude data sets sufficient for determining full moment-tensor microearthquake mechanisms by the linearprogramming inversion method. We have developed a GUI-based interface to this inversion software that greatly increases its ease of use and makes feasible analyzing larger numbers of earthquakes than previously was practical. We show examples from an injection experiment conducted in well 34-9RD2, on the East Flank of the Coso geothermal area. This tight well was re-drilled February – March 2005 with the intention of hydrofracturing it, but instead, pervasive porosity and fractures were encountered at about 2660 m depth. Total drilling mud losses occurred, obviating the need to stimulate the well. These mud losses induced a 50-minute swarm of 44 microearthquakes, with magnitudes in the range -0.3 to 2.6. Most of the largest microearthquakes occurred in the first 2 minutes. Accurate relative relocations and moment tensors for the best-recorded subset reveal fine details of the fracture that was stimulated. This comprised a fault striking at N 20° E and dipping at 75° to the WNW, which propagated to the NNE and upward. Co-injection focal mechanisms reveal combined crack-opening and shear motion. Stress release and mode of failure differed between the pre-, co- and post-swarm periods. Some post-swarm events involved cavity collapse, suggesting that some of the cavity opening caused by the fluid injection was quickly reversed. Stress & mode of failure had not returned to pre-swarm conditions within 1 month following the injection, posing the question of how long stress perturbations persist following a stimulation experiment. This question may be answered by processing data spanning a longer postinjection period, work that is currently underway and will be reported in this presentation

Time-dependent seismic tomography of the Coso geothermal area, 1996-2004

Local-earthquake tomographic images were calculated for each of the years 1996 - 2004 using arrival times from the U.S. Navy’s permanent seismometer network at the Coso geothermal area, California. The results show irregular strengthening with time of the wave-speed ratio VP/VS at shallow depths. These changes result predominately from progressive relative increase in VS with respect to VP, and could result from processes associated with geothermal operations such as decrease in fluid pressure and the drying of argillaceous minerals such as illite

Crustal structure beneath western and eastern Iceland from surface waves and receiver functions

We determine the crustal structures beneath 14 broad-band seismic stations, deployed in western, eastern, central and southern Iceland, using surface wave dispersion curves and receiver functions. We implement a method to invert receiver functions using constraints obtained from genetic algorithm inversion of surface waves. Our final models satisfy both data sets. The thickness of the upper crust, as defined by the velocity horizon Vs= 3.7 km s−1, is fairly uniform at ∼6.5–9 km beneath the Tertiary intraplate areas of western and eastern Iceland, and unusually thick at 11 km beneath station HOT22 in the far south of Iceland. The depth to the base of the lower crust, as defined by the velocity horizon Vs= 4.1 km s−1 is ∼20–26 km in western Iceland and ∼27–33 km in eastern Iceland. These results agree with those of explosion profiles that detect a thinner crust beneath western Iceland than beneath eastern Iceland. An earlier report of a substantial low-velocity zone beneath the Middle Volcanic Zone in the lower crust is confirmed by a similar observation beneath an additional station there. As was found in previous receiver function studies, the most reliable feature of the results is the clear division into an upper sequence that is a few kilometres thick where velocity gradients are high, and a lower, thicker sequence where velocity gradients are low. The transition to typical mantle velocities is variable, and may range from being very gradational to being relatively sharp and clear. A clear Moho, by any definition, is rarely seen, and there is thus uncertainty in estimates of the thickness of the crust in many areas. Although a great deal of seismic data are now available constraining the structures of the crust and upper mantle beneath Iceland, their geological nature is not well understood

Seismic evidence for a tilted mantle plume and north–south mantle flow beneath Iceland

Post-print / lokagerð höfundaShear waves converted from compressional waves at mantle discontinuities near 410- and 660-km depth recorded by two broadband seismic experiments in Iceland reveal that the center of an area of anomalously thin mantle transition zone lies at least 100 km south of the upper-mantle low-velocity anomaly imaged tomographically beneath the hotspot. This offset is evidence for a tilted plume conduit in the upper mantle, the result of either northward flow of the Icelandic asthenosphere or southward flow of the upper part of the lower mantle in a no-net-rotation reference frame.Peer Reviewe

A History of Discrete Event Simulation Programming Languages

The history of simulation programming languages is organized as a progression in periods of similar developments. The five periods, spanning 1955-1986, are labeled: The Period of Search (1955-1960); The Advent (1961-1965); The Formative Period (1966-1970); The Expansional Period (1971-1978); and The Period of Consolidation and Regeneration (1979-1986). The focus is on recognizing the people and places that have made important contributions in addition to the nature of the contribution. A balance between comprehensive and in-depth treatment has been reached by providing more detailed description of those languages which have or have had major use. Over 30 languages are mentioned, and numerous variations are described in the major contributors. A concluding summary notes the concepts and techniques either originating with simulation programming languages or given significant visibility by them

APOE ɛ4 exacerbates age-dependent deficits in cortical microstructure

The apolipoprotein E ɛ4 allele is the primary genetic risk factor for the sporadic type of Alzheimer’s disease. However, the mechanisms by which apolipoprotein E ɛ4 are associated with neurodegeneration are still poorly understood. We applied the Neurite Orientation Dispersion Model to characterize the effects of apolipoprotein ɛ4 and its interactions with age and education on cortical microstructure in cognitively normal individuals. Data from 1954 participants were included from the PREVENT-Dementia and ALFA (ALzheimer and FAmilies) studies (mean age = 57, 1197 non-carriers and 757 apolipoprotein E ɛ4 carriers). Structural MRI datasets were processed with FreeSurfer v7.2. The Microstructure Diffusion Toolbox was used to derive Orientation Dispersion Index maps from diffusion MRI datasets. Primary analyses were focused on (i) the main effects of apolipoprotein E ɛ4, and (ii) the interactions of apolipoprotein E ɛ4 with age and education on lobar and vertex-wise Orientation Dispersion Index and implemented using Permutation Analysis of Linear Models. There were apolipoprotein E ɛ4 × age interactions in the temporo-parietal and frontal lobes, indicating steeper age-dependent Orientation Dispersion Index changes in apolipoprotein E ɛ4 carriers. Steeper age-related Orientation Dispersion Index declines were observed among apolipoprotein E ɛ4 carriers with lower years of education. We demonstrated that apolipoprotein E ɛ4 worsened age-related Orientation Dispersion Index decreases in brain regions typically associated with atrophy patterns of Alzheimer’s disease. This finding also suggests that apolipoprotein E ɛ4 may hasten the onset age of dementia by accelerating age-dependent reductions in cortical Orientation Dispersion Index