Precursors to earthquakes: Seismoelectromagnetic signals

Field measurements in the past several years have documented electromagnetic signals which are attributed to precursory stress and strain changes which were followed ultimately by earthquakes. Precursory electric field changes observed in Greece on multiple dipoles have been used to issue earthquake predictions. While the source of these signals is still unknown, a sufficient number of predictions has been issued to allow some, but not all, statistical analyses to show this method is better than randomly sampling the earthquake catalog. Ongoing efforts to identify the sources of both these signals and the magnetic field variations prior to the Loma Prieta earthquake are focusing on electrokinetic coupling of fluid flow and transient electric fields. A mechanism related to local fluid flow appears to be best suited at this time of explaining the variety of purported precursors. However, much more work is needed to improve the observations and refine the models of precursor generation. Efforts to monitor magnetotelluric transfer functions at longer periods (T>10s) have been hampered by variability of the functions. The use of modern noise reduction techniques such as remote referencing should reduce this variability, but may not reduce errors to a level needed for monitoring. Monitoring of high frequency (81 kHz) seismoelectric emissions may be promising, but lack of simultaneous observations on multiple instruments hinders the utility of this technique.     

Temporal variability in summertime bottom hypoxia in shallow areas of Mobile Bay,Alabama

This paper addresses temporal variability in bottom hypoxia in broad shallow areas of Mobile Bay, Alabama. Time-series data collected in the summer of 2004 from one station (mean depth of 4 m) exhibit bottom dissolved oxygen (DO) variations associated with various time scales of hours to days. Despite a large velocity shear, stratification was strong enough to suppress vertical mixing most of the time. Bottom DO was closely related to the vertical salinity gradient (ΔS). Hypoxia seldom occurred when ΔS (over 2.5 m) was <2 psu and occurred almost all the time when ΔS was >8 psu in the absence of extreme events like hurricanes. Oxygen balance between vertical mixing and total oxygen demand was considered for bottom water from which oxygen demand and diffusive oxygen flux were estimated. The estimated decay rates at 20°C ranging between 0.175–0.322 d⁻¹ and the corresponding oxygen consumption as large as 7.4 g O₂ m⁻² d⁻¹ fall at the upper limit of previously reported ranges. The diffusive oxygen flux and the corresponding vertical diffusivity estimated for well mixed conditions range between 8.6–9.5 g O₂ m⁻² d⁻¹ and 2.6–2.9 m² d⁻¹, respectively. Mobile Bay hypoxia is likely to be associated with a large oxygen demand, supported by both water column and sediment oxygen demands, so that oxygen supply from surface water during destratification events would be quickly exhausted to return to hypoxic conditions within a few hours to days after destratification events are terminated.     

Discrimination of small earthquakes and artificial explosions in the Korean Peninsula using Pg/Lg ratios

The purpose of this study is to develop a technique to discriminate artificial explosions from local small earthquakes ( M ≤ 4.0) in the time–frequency domain. In order to obtain spectral features of artificial explosions and earthquakes, 3-D spectrograms (frequency, time and amplitude) have been used. They represent a useful tool for studying the frequency content of entire seismic waveforms observed at local and regional distances (Kim, Simpson & Richards 1994). P and S(L g ) waves from quarry blasts show that the frequency content associated with the dominant amplitude appears above 10  Hz and Rg phases are observed at close distances. P and S(L g ) waves from the Tongosan earthquake have strong amplitudes below 10  Hz. For the Munkyong earthquake, however, a broader frequency content up to 20  Hz is found.
  For discrimination between small earthquakes and explosions, Pg/L g spectral ratios are used below 10  Hz, and through spectrogram analysis we can see different frequency contents of explosions and earthquakes. Unfortunately, because explosion data recorded at KSRS array are digitized at 20  sps, we cannot avoid analysing below 10  Hz because of the Nyquist frequency. In order to select time windows, the group velocity was computed using multiple-filter analysis (MFA), and free-surface effects have been removed from all three-component data in order to improve data quality. Using FFT, a log-average spectral amplitude is calculated over seven frequency bands: 0.5 to 3, 2 to 4, 3 to 5, 4 to 6, 5 to 7, 6 to 8 and 8 to 10  Hz. The best separation between explosions and earthquakes is observed from 6 to 8  Hz. In this frequency band we can separate explosions with log ( Pg/L g ) above −0.5, except EXP1 recorded at SIHY1-1, and earthquakes below −0.5, except the Munkyong earthquake record at station KMH.     

Landslide susceptibility mapping using frequency ratio, analytic hierarchy process, logistic regression, and artificial neural network methods at the Inje area, Korea

Every year, the Republic of Korea experiences numerous landslides, resulting in property damage and casualties. This study compared the abilities of frequency ratio (FR), analytic hierarchy process (AHP), logistic regression (LR), and artificial neural network (ANN) models to produce landslide susceptibility index (LSI) maps for use in predicting possible landslide occurrence and limiting damage. The areas under the relative operating characteristic (ROC) curves for the FR, AHP, LR, and ANN LSI maps were 0.794, 0.789, 0.794, and 0.806, respectively. Thus, the LSI maps developed by all the models had similar accuracy. A cross-tabulation analysis of landslide occurrence against non-occurrence areas showed generally similar overall accuracies of 65.27, 64.35, 65.51, and 68.47 % for the FR, AHP, LR, and ANN models, respectively. A correlation analysis between the models demonstrated that the LR and ANN models had the highest correlation (0.829), whereas the FR and AHP models had the lowest correlation (0.619).     

Wave reflection from partially perforated-wall caisson breakwater

In 1995, Suh and Park developed a numerical model that computes the reflection of regular waves from a fully perforated-wall caisson breakwater. This paper describes how to apply this model to a partially perforated-wall caisson and irregular waves. To examine the performance of the model, existing experimental data are used for regular waves, while a laboratory experiment is conducted in this study for irregular waves. The numerical model based on a linear wave theory tends to over-predict the reflection coefficient of regular waves as the wave nonlinearity increases, but such an over-prediction is not observed in the case of irregular waves. For both regular and irregular waves, the numerical model slightly over- and under-predicts the reflection coefficients at larger and smaller values, respectively, because the model neglects the evanescent waves near the breakwater.     

Microbial activity and chemical weathering in the Middendorf aquifer,South Carolina

We use reactive transport modeling to better understand the kinetics of chemical weathering in the Cretaceous Middendorf aquifer of South Carolina, USA, and the relationship of this process to subsurface microbial activity. We constructed a model accounting for the kinetics of mineral dissolution and precipitation, ion exchange, and the CO₂ and bicarbonate produced by iron reducing and sulfate reducing bacteria in the aquifer. We then fit the model to observed trends in the chemical composition of groundwater along the aquifer by adjusting the rate constants for the kinetic reactions considered. The modeling portrays weathering in the Middendorf as a slow process by which groundwater gradually reacts toward equilibrium with minerals in the aquifer. The rate constants predicted are 6 to 7 orders of magnitude smaller than measured in laboratory experiments and 3 to 4 orders of magnitude less than those inferred from weathering rates in soils. The rate constants are smaller even than expected by projecting observed trends with the duration of weathering to the geologic age of the Middendorf. Weathering is driven largely by biological activity: about half the acid consumed is CO₂ derived from the recharge area, and about half is supplied by iron reducing bacteria in the aquifer; only about 1% of the acid is of atmospheric origin, from CO₂ dissolved in rainwater.     

Tephrostratigraphy and paleo-environmental implications of Late Quaternary sediments and interstitial water in the western Ulleung Basin, East/Japan Sea

Two piston cores, collected from the western Ulleung Basin of the East/Japan Sea, were used to investigate the Late Quaternary tephrostratigraphy, lithology, and mineralogy of the sediments, as well as the elemental composition of both sediments and interstitial water, and their paleo-environmental implications. The cores show two lapilli tephra layers and one rhyolitic ash layer at the boundaries between sedimentary lithofacies units I, II, and III. These layers can be correlated with the well-known Ulleung-Oki (ca. 9.3 ka; boundary of units I/II), Ulleung-Yamato (ca. 25–33 ka), and Aira-Tanzawa (ca. 25.1 ka; boundary of units II/III) layers, respectively. These data suggest that the cores cover the period from the middle stage of marine isotope stage (MIS) 3 to the Holocene. In each core, a so-called dark laminated mud (DLM) layer in unit II commonly has high contents of Si and Al, suggesting that the DLM layer generally contains significant amounts of fine-grained silicates and/or aluminosilicates. In addition, these DLM layers are enriched in Fe relative to Mn, iron being predominantly bound in framboidal pyrite. The size distribution of pyrite in the DLM layers suggests that these have formed under anoxic (euxinic) conditions at times of reduced water circulation in the basin during sea-level lowstand (euxinic environment). The C/N ratios [5–12] suggest that sedimentary organic matter in the cores is predominantly of marine origin. Total organic carbon (TOC) contents increased shortly before Termination I. This could have been caused by an increased flux of marine organic matter in association with sea-level rise. The C and S values of units II and III (Late Pleistocene sediments; C: <2.0%; S: 0.5–2.1%) suggest a more euxinic environment than that of unit I (Holocene sediments; C: 1.0–3.5%; S: <1.0%). Concentrations of SO₄^2– in the interstitial water decrease with increasing burial depth, whereas CH₄ concentrations show the reverse trend. Therefore, it seems that sulfate reduction, probably related to microbial activity, predominates in the upper core sections (<5 m), shifting to methanogenesis in the lower core sections.     

Direct shortwave radiative forcing of the Asian dust aerosol on dust emission

The physical processes of the feedback mechanism of direct shortwave radiative forcing of the Asian dust aerosol on dust emission has been examined using simulated results with the coupled (with dust shortwave radiative forcing) and the non-coupled model (without dust shortwave radiative forcing) based on the MM5 model and the Asian Dust Aerosol Model on 19 March 2002. The results indicate that a significant dust emission reduction occurs in the high dust concentration (HDC) region of the dust source region whereas an enhanced dust emission appears in the downstream of the dust source region. It is found that Asian dust aerosols raised during the daytime by the strong surface wind cause negative shortwave radiative flux near the surface, which in turn reduces the sensible heat flux causing the cooling of the air, thereby enhancing stable stratification. The dynamic adjustment of the negative radiative flux of the dust induces a positive pressure anomaly over the HDC region and a negative pressure anomaly toward the synoptic low pressure center, resulting in a dipole shape of pressure anomaly field near the surface. The associated secondary circulation of this pressure anomaly together with the reduction of turbulent intensity due to the reduced sensible heat flux reduces the low-level wind speed thereby reducing dust emission in the upstream of the HDC region of the dust source region (Region I), while enhancing the low-level wind speed in the downstream region (Region II), which in turn enhances dust emission. This enhanced dust emission is smaller than the emission reduction in the upstream, resulting in overall dust emission reduction during the daytime.