An Application of Hydraulic Tomography to a Large‐Scale Fractured Granite Site,Mizunami, Japan

While hydraulic tomography (HT) is a mature aquifer characterization technology, its applications to characterize hydrogeology of kilometer‐scale fault and fracture zones are rare. This paper sequentially analyzes datasets from two new pumping tests as well as those from two previous pumping tests analyzed by Illman et al. (2009) at a fractured granite site in Mizunami, Japan. Results of this analysis show that datasets from two previous pumping tests at one side of a fault zone as used in the previous study led to inaccurate mapping of fracture and fault zones. Inclusion of the datasets from the two new pumping tests (one of which was conducted on the other side of the fault) yields locations of the fault zone consistent with those based on geological mapping. The new datasets also produce a detailed image of the irregular fault zone, which is not available from geological investigation alone and the previous study. As a result, we conclude that if prior knowledge about geological structures at a field site is considered during the design of HT surveys, valuable non‐redundant datasets about the fracture and fault zones can be collected. Only with these non‐redundant data sets, can HT then be a viable and robust tool for delineating fracture and fault distributions over kilometer scales, even when only a limited number of boreholes are available. In essence, this paper proves that HT is a new tool for geologists, geophysicists, and engineers for mapping large‐scale fracture and fault zone distributions.     

Experimental verification of a wireless sensing and control system for structural control using MR dampers

The performance aspects of a wireless ‘active’ sensor, including the reliability of the wireless communication channel for real‐time data delivery and its application to feedback structural control, are explored in this study. First, the control of magnetorheological (MR) dampers using wireless sensors is examined. Second, the application of the MR‐damper to actively control a half‐scale three‐storey steel building excited at its base by shaking table is studied using a wireless control system assembled from wireless active sensors. With an MR damper installed on each floor (three dampers total), structural responses during seismic excitation are measured by the system's wireless active sensors and wirelessly communicated to each other; upon receipt of response data, the wireless sensor interfaced to each MR damper calculates a desired control action using an LQG controller implemented in the wireless sensor's computational core. In this system, the wireless active sensor is responsible for the reception of response data, determination of optimal control forces, and the issuing of command signals to the MR damper. Various control solutions are formulated in this study and embedded in the wireless control system including centralized and decentralized control algorithms. Copyright © 2007 John Wiley & Sons, Ltd.     

A GALERKIN FINITE-ELEMENT PROGRAM FOR SIMULATING UNSATURATED FLOW IN POROUS MEDIA

Abstract. A fully documented Galerkin finite-element FORTRAN program is presented for solving the one-dimensional, transient flow equation in unsaturated porous media. Material balance error summaries are presented to demonstrate accuracy of the numerical scheme. Comparison of our simulated results with other existing numerical solutions using the Galerkin scheme provided excellent agreement.     

GPS Height and Gravity Variations Due to Ocean Tidal Loading Around Taiwan

This study presents predicts ocean tidal loading (OTL) effects using a Green’s function approach and validates a novel tidal model for Taiwan. Numerical integration of OTL is performed using the Gauss quadrature method and a local tidal model for the inner zone and a global model for the outer zone. Observed time series of GPS-derived vertical displacements and gravity variations (3–7 days) at five co-located GPS-gravimeter stations along the South East China and Taiwan coasts were utilized to assess the accuracy of the proposed models and two other models. The OTL-induced gravity variations are 3–16 μgal and the vertical site displacements are 13–27 mm. Generally, an OTL model using a mixed global and local tidal model generates better agreement with the observations than an OTL model using a global tidal model only. However, containing a local model inside a global model does not always produce a good agreement with the observations. The relatively large discrepancies between modeled and observed OTL values at some stations indicate that there is a need for an improved local tidal model in the study area. An erratum to this article can be found at     

Effect of Near-fault Terrain upon Dislocation Modeling

—Coseismic surface deformation provides important information needed to determine source rupture geometry and slip distribution as well as to estimate seismic moment. In this study, numerical experiments were designed to analyze and classify how free-surface topography affects surface deformation. The investigation was performed by 3-D finite element modeling. Results of this study show that crustal deformation induced by near-fault terrain is significant and can be measured with present geodetic survey techniques. The characteristics of the terrain effects show that a hill structure produces more crustal deformation than a half-space model, and that the crustal deformation of a basin structure is less than that of the half-space model. The topographic correction is in the order of five percent of the fault dislocation. On the basis of the relationship between fault offset and earthquake magnitude, it is suggested that the terrain effects on the coseismic crustal deformation of shallow earthquakes with a magnitude greater than 5.6 should be considered as one of the major errors in coseismic deformation modeling which ignored the surface topography on the order of 300 meters.     

Tabulation of Mie properties for an effective microwave radiative model

Summary In microwave radiative transfer model simulations of atmospheric convective systems, the Mie calculations usually consume the majority of the computer time necessary for the calculations (70 to 90% for frequencies ranging from 6 to 300 GHz). For a large array of atmospheric profiles, the repeated Mie scattering calculations make the radiative transfer modeling not only expensive, but often impossible. A set of Mie tables, therefore, is developed to replace the actual Mie calculations in the microwave radiation model.The tables are divided into categories of liquid water (cloud water and rain), ice, snow, graupel and hail. The tables are free from restriction on variation of size distribution, particle density and mixing phases. Results from this study show that by using the Mie tables, the total CPU time is reduced by at least one order of magnitude, depending on the frequency. Compared to using full Mie calculation in the radiative model, the bias and root-mean-square (rms) deviation from the Mie tables are less than 1.2 K for any atmosphere with precipitation rates up to 125 mm hr^–1 at any frequency and any viewing angle.The Mie tables are particularly useful in simulating the microwave radiation field over a mesoscale meteorological system. Such a study would not be satisfactory using a Mie parameterization scheme, because one would have to trade off the accuracy and flexibility for modeling efficiency.With 7 FiguresThis research was done under contract at the Goddard Space Flight Center, Greenbelt, Maryland.     

Changes in mixed layer depth under climate change projections in two CGCMs

Two coupled general circulation models, i.e., the Meteorological Research Institute (MRI) and Geophysical Fluid Dynamics Laboratory (GFDL) models, were chosen to examine changes in mixed layer depth (MLD) in the equatorial tropical Pacific and its relationship with ENSO under climate change projections. The control experiment used pre-industrial greenhouse gas concentrations whereas the 2 × CO₂ experiment used doubled CO₂ levels. In the control experiment, the MLD simulated in the MRI model was shallower than that in the GFDL model. This resulted in the tropical Pacific’s mean sea surface temperature (SST) increasing at different rates under global warming in the two models. The deeper the mean MLD simulated in the control simulation, the lesser the warming rate of the mean SST simulated in the 2 × CO₂ experiment. This demonstrates that the MLD is a key parameter for regulating the response of tropical mean SST to global warming. In particular, in the MRI model, increased stratification associated with global warming amplified wind-driven advection within the mixed layer, leading to greater ENSO variability. On the other hand, in the GFDL model, wind-driven currents were weak, which resulted in mixed-layer dynamics being less sensitive to global warming. The relationship between MLD and ENSO was also examined. Results indicated that the non-linearity between the MLD and ENSO is enhanced from the control run to the 2 × CO₂ run in the MRI model, in contrast, the linear relationship between the MLD index and ENSO is unchanged despite an increase in CO₂ concentrations in the GFDL model.     

Linear spatial interpolation: Analysis with an application to San Joaquin Valley

The properties of linear spatial interpolators of single realizations and trend components of regionalized variables are examined in this work. In the case of the single realization estimator explicit and exact expressions for the weighting vector and the variances of estimator and estimation error were obtained from a closed-form expression for the inverse of the Lagrangian matrix. The properties of the trend estimator followed directly from the Gauss-Markoff theorem. It was shown that the single realization estimator can be decomposed into two mutually orthogonal random functions of the data, one of which is the trend estimator. The implementation of liear spatial estimation was illustrated with three different methods, i.e., full information maximum likelihood (FIML), restricted maximum likelihood (RML), and Rao's minimum norm invariant quadratic unbiased estimation (MINQUE) for the single realization case and via generalized least squares (GLS) for the trend. The case study involved large correlation length-scale in the covariance of specific yield producing a nested covariance structure that was nearly positive semidefinite. The sensitivity of model parameters, i.e., drift and variance components (local and structured) to the correlation length-scale, choice of covariance model (i.e., exponential and spherical), and estimation method was examined. the same type of sensitivity analysis was conducted for the spatial interpolators. It is interesting that for this case study, characterized by a large correlation length-scale of about 50 mi (80 km), both parameter estimates and linear spatial interpolators were rather insensitive to the choice of covariance model and estimation method within the range of credible values obtained for the correlation length-scale, i.e., 40–60 mi (64–96 km), with alternative estimates falling within ±5% of each other.     

Combined effect of El Niño-Southern Oscillation and Pacific Decadal Oscillation on the East Asian winter monsoon

Using long-term observational data and numerical model experiments, the combined effect of the El Niño-Southern Oscillation (ENSO) and Pacific Decadal Oscillation (PDO) on the variability of the East Asian winter monsoon is examined. In the observations, it is found that when the ENSO and PDO are in-phase combinations (i.e., El Niño/positive PDO phase and La Niña/negative PDO phase), a negative relationship between ENSO and East Asian winter monsoon is significantly intensified. In other words, when El Niño (La Niña) occurs with positive (negative) PDO phase, anomalous warm (cold) temperatures are dominant over the East Asian winter continent. On the other hand, there are no significant temperature anomalies when the ENSO and PDO are out-of-phase combinations (i.e., El Niño/negative PDO phase and La Niña/positive PDO phase). Further analyses indicate that the anticyclone over the western North Pacific including the East Asian marginal seas plays an essential role in modulating the intensity of the East Asian winter monsoon under the changes of ENSO–PDO phase relationship. Long-lasting high pressure and warm sea surface temperature anomalies during the late fall/winter and following spring over the western North Pacific, which appear as the El Niño occurs with positive PDO phase, can lead to a weakened East Asian winter monsoon by transporting warm and wet conditions into the East Asian continent through the southerly wind anomalies along the western flank of the anomalous high pressure, and vice versa as the La Niña occurs with negative PDO phase. In contrast, the anomalous high pressure over the western North Pacific does not show a prominent change under the out-of-phase combinations of ENSO and PDO. Numerical model experiments confirm the observational results, accompanying dominant warm temperature anomalies over East Asia via strong anticyclonic circulation anomalies near the Philippine Sea as the El Niño occurs with positive PDO phase, whereas such warming is weakened as the El Niño occurs with negative PDO phase. This result supports the argument that the changes in the East Asian winter monsoon intensity with ENSO are largely affected by the strength of the anticyclone over the western North Pacific, which significantly changes according to the ENSO–PDO phase relationship.