Multi-aquifer pumping test to determine cutoff wall length for groundwater flow control during site excavation in Tokyo, Japan

Multi-aquifer pumping tests, using a multi-screen pumping well and multi-level piezometers, were carried out for groundwater flow control in a large-scale excavation site in Tokyo, Japan. The site was underlain by multi-layered confined aquifers. In the tests, pumping was carried out using a multi-aquifer pumping well in which a screen depth was chosen arbitrarily. Changes in groundwater pressure heads in each aquifer were measured at each screen position of the multi-aquifer pumping well. Hydraulic conductivity (K) and specific storage (S _s) of not only aquifers, but also for low permeability layers between the aquifers, were estimated using the Cooper-Jacob method, and calibrated by a finite element method (FEM) groundwater model. Four different cutoff wall lengths were assumed for final excavation depth, and correlations among wall length, pumping discharge and drawdown at the back of the cutoff wall were obtained from simulations using the K and S _s parameters in the FEM model. Then, the most suitable wall length was selected based on the simulated correlations considering environmental condition, construction period and cost of the cutoff wall.     

Maintaining real-time precise point positioning during outages of orbit and clock corrections

The precise point positioning (PPP) is a popular positioning technique that is dependent on the use of precise orbits and clock corrections. One serious problem for real-time PPP applications such as natural hazard early warning systems and hydrographic surveying is when a sudden communication break takes place resulting in a discontinuity in receiving these orbit and clock corrections for a period that may extend from a few minutes to hours. A method is presented to maintain real-time PPP with 3D accuracy less than a decimeter when such a break takes place. We focus on the open-access International GNSS Service (IGS) real-time service (RTS) products and propose predicting the precise orbit and clock corrections as time series. For a short corrections outage of a few minutes, we predict the IGS-RTS orbits using a high-order polynomial, and for longer outages up to 3 h, the most recent IGS ultra-rapid orbits are used. The IGS-RTS clock corrections are predicted using a second-order polynomial and sinusoidal terms. The model parameters are estimated sequentially using a sliding time window such that they are available when needed. The prediction model of the clock correction is built based on the analysis of their properties, including their temporal behavior and stability. Evaluation of the proposed method in static and kinematic testing shows that positioning precision of less than 10 cm can be maintained for up to 2 h after the break. When PPP re-initialization is needed during the break, the solution convergence time increases; however, positioning precision remains less than a decimeter after convergence.     

The degree of Cl-loss for the particulate matter (PM) and fog water sampled at the same air mass at the Hachimantai mountain range in northern Japan

It is well known that sea-salt aerosols in particulate matter (PM) react with acids such as H₂SO₄ and HNO₃ during transportation and thereby lose chloride ions (Cl-loss). The PM and fog were sampled concurrently at different altitudes in the Hachimantai mountain range, northern Japan. The PM and fog sampled at different altitudes had nearly identical properties for the ion components. However, the PM was in a Cl-depleted state (more than 80% of all samples), but the fog water was not in so Cl-depleted state (less than 29%). As a result, it could be explained that this phenomenon caused because the fog droplets took up the gaseous state HCl other than sea-salt PM. After all Cl^- in the fog water recovered and was rather rich compared with the sea-salt or the PM by the uptake of the gaseous state HCl. Moreover, it was found that for PM_coarse(2.5 < D < 10), 86% of the acid (H₂SO₄ and HNO₃) was consumed for Cl-loss reactions and/or for dissolution of Ca and Mg in soil particles.     

A heating process of Kuchi-erabu-jima volcano,Japan, as inferred from geomagnetic field variations and electrical structure

Since August 2000, we have recorded the total intensity of the geomagnetic field at the summit area of Kuchi-erabu-jima volcano, where phreatic eruptions have repeatedly occurred. A time series analysis has shown that the variations in the geomagnetic field since 2001 have a strong relationship to an increase in volcanic activity. These variations indicate thermal demagnetization of the subsurface around the presently active crater. The demagnetization source for the early variations, until summer 2002, was estimated at about 200 m below sea level. For the variations since 2003, the source was modeled on the basis of the expansion of a uniformly magnetized ellipsoid. The modeling result showed that the source is located at 300 m above sea level beneath the crater. We carried out an audio-frequency magnetotelluric survey with the aim of obtaining a relation between the demagnetization source and the shallow structure of the volcano. A two-dimensional inversion applied to the data detected two good conductors, a shallow thin one which is restricted to a region around the summit area, while the other extends over the edifice at depths between 200 and 800 m. These conductors are regarded as clay-rich layers with low permeability, which were assumed to be generated through hydrothermal alteration. The demagnetization source for the early variations was possibly located at the lower part of the deep conductor and the source after 2003 lies between the two conductors, where groundwater is considered to be abundant. Based on these results, as well as on seismological, geodetic, and geochemical information, we propose a heating process of the Kuchi-erabu-jima volcano. In the initial stage, high-temperature volcanic gases supplied from the deep-seated magma remained temporarily at the level around the lower part of the less permeable deep conductor since the ascent path had not yet been established. Then, when the pathway developed as a result of repeated earthquakes, it became possible for a massive flux of volcanic gases to ascend through the conductor. The high temperature gases reached the aquifer located above the conductor and the heat was efficiently transported to the surrounding rocks through the groundwater. As a consequence, an abrupt increase of the gas flux and diffusion of the heat through the aquifer occurred and the high-temperature zone expanded. Since the high-temperature zone is located beneath another conductor, which acts as caprock, we assume that the energy of the phreatic explosion is accumulated there.     

Fractal size and spatial distributions of fault zones: An investigation into the seismic Chelungpu Fault, Taiwan

Abstract   The fractality of fault zones of thickness T for the 248-m-long core from a borehole penetrating the northern segment of the seismic Chelungpu Fault, Taiwan, was analyzed. The frequency curve of T shows that the fractal dimension is normal when T is smaller than a characteristic thickness T _c, and it becomes abnormally large when T exceeds T _c. The fractal dimensions of size and spatial distributions of T increase as the mean distribution density of T increases, which is inconsistent with the evolution laws for ordinary brittle faults. This discrepancy implies that the thickening rate of T when T is more than T _c is not constant, but a decreasing function of fault displacement. The slow thickening rate is related to the elastohydrodynamic lubrication which was effective on the fault when T exceeds T _c. This slip instability mechanism can explain the large, fast and smooth slip on the northern segment of the Chelungpu Fault during the 1999 Chi-Chi earthquake.     

Clayey injection veins and pseudotachylyte from two boreholes penetrating the Chelungpu Fault, Taiwan: Their implications for the contrastive seismic slip behaviors during the 1999 Chi-Chi earthquake

Abstract   The seismic slip that occurred during the 1999 Chi-Chi earthquake in Taiwan showed contrastive behaviors in different regions along the Chelungpu Fault: A large and smooth slip occurred in the north, while a relatively small slip associated with high-frequency seismic wave radiation occurred in the south. The core samples from shallow boreholes at northern (Fengyuan) and southern (Nantou) sites penetrating the seismic Chelungpu Fault were analyzed. The fault zones at the northern site are characterized by soft clayey material associated with clayey injection veins. This suggests that the fault zones were pressurized during ancient seismic slip events, and hence the elastohydrodynamic lubrication occurred effectively. In contrast, the fault rock from the southern site is old pseudotachylyte that has been shattered by repeated ancient seismic slip events. Statistical analysis of many pseudotachylyte fragments reveals that the degree of frictional melting tended to be low. In this case, the seismic slip is restrained by the mechanical barrier of a highly viscous melt layer. These contrastive fault rocks were produced by repeated ancient seismic slip events, but the two corresponding mechanisms of friction are likely to have also occurred during the 1999 Chi-Chi earthquake, thus causing the contrastive slip behaviors in the north and south.     

Analysis of a guided wave along a conducting structure in a borehole

A guided wave along a borehole is often observed in borehole radar measurements. These guided waves deform the antenna pattern and can cause artefacts in radar measurements. A water-filled borehole or a conducting logging cable can function as a waveguide for electromagnetic waves under some conditions. We describe the theoretical characteristics of such a guided wave in a borehole and compare them with our experiments. The measured signal discussed was obtained with a directional borehole radar. This radar uses a cylindrical conformal array antenna as receiver and is a model of a conducting structure in a borehole. The induced field around the borehole was compared with the theory. The most fundamental symmetrical and asymmetrical modes were TM₀₁ and HE₁₁, and they were identified in the measured signals using time–frequency distribution analysis and by observation of the azimuthal field distribution of the magnetic field.     

Evolution of the Komiji Syncline in the North Fossa Magna, central Japan: Paleomagnetic and K–Ar age insights

Abstract We carried out paleomagnetic measurements and K–Ar dating on Neogene andesitic lavas and sills of the Shigarami Formation in North Fossa Magna, central Japan. The Shigarami Formation is distributed in the axial part of the Komiji Syncline in the folding zone of the southwestern North Fossa Magna. Results of the present study indicate that the Komiji Syncline was formed shortly after 4.42 ± 0.12 Ma during the Pliocene. The sedimentary rocks of the Shigarami Formation consist of shallow marine and fluvial deposits. Intrusions of andesitic sills are found in the shallow marine deposits and two andesitic lava flows are present in the fluvial deposits. Oriented samples were taken from the sills at four sites and from the lavas at three sites. The samples produced stable remanent magnetization through stepwise alternating field and thermal demagnetizations. Results of a positive fold test indicate that the stable remanent magnetizations concentrate around a mean reversed polarity of declination = 169.0°, inclination = ?58.5° and 95% confidence limit = 9.0° after corrections have been made according to the direction of the bedding of the sedimentary rocks. Four fresh samples were selected for K–Ar dating from the samples used for paleomagnetic measurements. The groundmass of three samples taken from the sills yield ages of 4.42 ± 0.12, 4.49 ± 0.22 and 4.69 ± 0.13 Ma, whereas the one taken from the lower lava has an age of 5.91 ± 0.26 Ma. We believe that the Komiji Syncline was formed after the emplacement of lavas and sills in the area, because the descending Miocene strata were folded concordantly with the Shigarami Formation. The Pliocene and Pleistocene strata rest unconformably on the folded strata. The deformation might have progressed during the Pliocene, then slowed down in the Early Pleistocene. Our results suggest that the northwestward motion of the Philippine Sea Plate and the collision of the Tanzawa Block affected not only the South Fossa Magna, but also the North Fossa Magna.     

Far-Field Tsunami Potential and a Real-Time Forecast System for the Pacific Using the Inversion Method

Estimating tsunami potential is anessential part of mitigating tsunami disasters. Weproposed a new method to estimate the far-fieldtsunami potential by assuming faultmodels on the Pacific Rim. We find thata tsunami that generates in the areas wherethere is no tsunami in the history can damagethe Japanese coast. This shows that it isimportant to estimate tsunami potential byassuming fault models other than the pastearthquake data.Another important activity to mitigate tsunamidisasters is to provide appropriatewarnings to coastal communities when dangerfrom a tsunami is imminent. We applied anew inversion method using wavelet transformto a part of the real-time tsunami forecastsystem for the Pacific. Because this inversionmethod does not require fault location, it ispossible to analyze a tsunami in real timewithout all seismic information. In order tocheck the usability of the system, anumerical simulation was executed assuming anearthquake at sea off Taiwan. The correlationcoefficient for the estimated initialwaveform to the assumed one was calculatedto be 0.78. It takes 90 min to capturetime-series waveform data from tsunamigauges and 5 sec to estimate the 2-D initialwaveform using the inversion method. After that,it takes 2 minutes to forecast thetsunami heights at the Japanese coast. Since thesum of these times is less than the 105minutes transit time of the tsunami fromTaiwan to Japan, it is possible to give a warningto the residents before the tsunami attacksthe Japanese coast. Comparing the tsunamiheights forecasted by this system with thosecalculated by the fault model, the averageerror was 0.39 m. The average error ofthe arrival time was 0.007 min.     

Detailed Fracture System of the Soultz-sous-Forêts HDR Field Evaluated Using Microseismic Multiplet Analysis

—?The reservoir structure of the Soultz HDR field has been investigated by examining induced microearthquake multiplets. Microseismic events with similar waveforms have been selected from microseismic data obtained during a 1993 hydraulic fracturing experiment. Precise relative arrival times and source locations have been determined by cross-spectrum analysis. The cross-spectrum analysis decreased the residual from 0.75?ms to 0.1?ms. The estimated orientations of the multiplet planes are consistent with fracture orientations detected in Soultz boreholes. A comparison between the stress field and the orientation of structural planes suggests that the structural planes were under a critical condition of frictional slip.