Geology topics

C. Chen

Publications and source records attributed to C. Chen.

10 recordsLinked to original sources

USGS2009

High-frequency Rayleigh-wave method

High-frequency (≥2 Hz) Rayleigh-wave data acquired with a multichannel recording system have been utilized to determine shear (S)-wave velocities in near-surface geophysics since the early 1980s. This overview article discusses the main research results of high-frequency surface-wave techniques achieved by research groups at the Kansas Geological Survey and China University of Geosciences in the last 15 years. The multichannel analysis of surface wave (MASW) method is a non-invasive acoustic approach to estimate near-surface S-wave velocity. The differences between MASW results and direct borehole measurements are approximately 15% or less and random. Studies show that simultaneous inversion with higher modes and the fundamental mode can increase model resolution and an investigation depth. The other important seismic property, quality factor ( Q ), can also be estimated with the MASW method by inverting attenuation coefficients of Rayleigh waves. An inverted model (S-wave velocity or Q ) obtained using a damped least-squares method can be assessed by an optimal damping vector in a vicinity of the inverted model determined by an objective function, which is the trace of a weighted sum of model-resolution and model-covariance matrices. Current developments include modeling high-frequency Rayleigh-waves in near-surface media, which builds a foundation for shallow seismic or Rayleigh-wave inversion in the time-offset domain; imaging dispersive energy with high resolution in the frequency-velocity domain and possibly with data in an arbitrary acquisition geometry, which opens a door for 3D surface-wave techniques; and successfully separating surface-wave modes, which provides a valuable tool to perform S-wave velocity profiling with high-horizontal resolution.

Earth Science↗

USGS2009

Discovery of lake-effect clouds on Titan

Images from instruments on Cassini as well as from telescopes on the ground reveal the presence of sporadic small-scale cloud activity in the cold late-winter north polar region of Saturn's large moon Titan. These clouds lie underneath the previously discovered uniform polar cloud attributed to a quiescent ethane cloud at ???40 km and appear confined to the same latitudes as those of the largest known hydrocarbon lakes at the north pole of Titan. The physical properties of these clouds suggest that they are due to methane convection and condensation. Such convection could be caused by a process in some ways analogous to terrestrial lake-effect clouds. The lakes on Titan could be a key connection between the surface and the meteorological cycle. ?? 2009 by the American Geophysical Union.

Geophysical Research Letters↗

USGS2009

Investigation of river eutrophication as part of a low dissolved oxygen total maximum daily load implementation

In the United States, environmentally impaired rivers are subject to regulation under total maximum daily load (TMDL) regulations that specify watershed wide water quality standards. In California, the setting of TMDL standards is accompanied by the development of scientific and management plans directed at achieving specific water quality objectives. The San Joaquin River (SJR) in the Central Valley of California now has a TMDL for dissolved oxygen (DO). Low DO conditions in the SJR are caused in part by excessive phytoplankton growth (eutrophication) in the shallow, upstream portion of the river that create oxygen demand in the deeper estuary. This paper reports on scientific studies that were conducted to develop a mass balance on nutrients and phytoplankton in the SJR. A mass balance model was developed using WARMF, a model specifically designed for use in TMDL management applications. It was demonstrated that phytoplankton biomass accumulates rapidly in a 88 km reach where plankton from small, slow moving tributaries are diluted and combined with fresh nutrient inputs in faster moving water. The SJR-WARMF model was demonstrated to accurately predict phytoplankton growth in the SJR. Model results suggest that modest reductions in nutrients alone will not limit algal biomass accumulation, but that combined strategies of nutrient reduction and algal control in tributaries may have benefit. The SJR-WARMF model provides stakeholders a practical, scientific tool for setting remediation priorities on a watershed scale.

Water Science and Technology↗

USGS2007

Wave field features of shallow vertical discontinuity and their application in non-destructive detection

The geotechnical integrity of critical infrastructure can be seriously compromised by the presence of fractures or crevices. Non-destructive techniques to accurately detect fractures in critical infrastructure such as dams and highways could be of significant benefit to the geotechnical industry. This paper investigates the application of shallow seismic and georadar methods to the detection of a vertical discontinuity using numerical simulations. The objective is to address the kinematical analysis of a vertical discontinuity, determine the resulting wave field characteristics, and provide the basis for determining the existence of vertical discontinuities based on the recorded signals. Simulation results demonstrate that: (1) A reflection from a vertical discontinuity produces a hyperbolic feature on a seismic or georadar profile; (2) In order for a reflection from a vertical discontinuity to be produced, a reflecting horizon below the discontinuity must exist, the offset between source and receiver (x0) must be non-zero, on the same side of the vertical discontinuity; (3) The range of distances from the vertical discontinuity where a reflection event is observed is proportional to its length and to x0; (4) Should the vertical crevice (or fracture) pass through a reflecting horizon, dual hyperbolic features can be observed on the records, and this can be used as a determining factor that the vertical crevice passes through the interface; and (5) diffractions from the edges of the discontinuity can be recorded with relatively smaller amplitude than reflections and their ranges are not constrained by the length of discontinuity. If the length of discontinuity is short enough, diffractions are the dominant feature. Real-world examples show that the shallow seismic reflection method and the georadar method are capable of recording the hyperbolic feature, which can be interpreted as vertical discontinuity. Thus, these methods show some promise as effective non-destructive detection methods for locating vertical discontinuities (e.g., fractures or crevices) in infrastructure such as dams and highway pavement. ?? 2007 Elsevier B.V. All rights reserved.

Journal of Applied Geophysics↗

USGS2006

Estimation of elastic moduli in a compressible Gibson half-space by inverting Rayleigh-wave phase velocity

A Gibson half-space model (a non-layered Earth model) has the shear modulus varying linearly with depth in an inhomogeneous elastic half-space. In a half-space of sedimentary granular soil under a geostatic state of initial stress, the density and the Poisson's ratio do not vary considerably with depth. In such an Earth body, the dynamic shear modulus is the parameter that mainly affects the dispersion of propagating waves. We have estimated shear-wave velocities in the compressible Gibson half-space by inverting Rayleigh-wave phase velocities. An analytical dispersion law of Rayleigh-type waves in a compressible Gibson half-space is given in an algebraic form, which makes our inversion process extremely simple and fast. The convergence of the weighted damping solution is guaranteed through selection of the damping factor using the Levenberg-Marquardt method. Calculation efficiency is achieved by reconstructing a weighted damping solution using singular value decomposition techniques. The main advantage of this algorithm is that only three parameters define the compressible Gibson half-space model. Theoretically, to determine the model by the inversion, only three Rayleigh-wave phase velocities at different frequencies are required. This is useful in practice where Rayleigh-wave energy is only developed in a limited frequency range or at certain frequencies as data acquired at manmade structures such as dams and levees. Two real examples are presented and verified by borehole S-wave velocity measurements. The results of these real examples are also compared with the results of the layered-Earth model. ?? Springer 2006.

Surveys in Geophysics↗

USGS2006

Nonlinear inversion of potential-field data using a hybrid-encoding genetic algorithm

Using a genetic algorithm to solve an inverse problem of complex nonlinear geophysical equations is advantageous because it does not require computer gradients of models or "good" initial models. The multi-point search of a genetic algorithm makes it easier to find the globally optimal solution while avoiding falling into a local extremum. As is the case in other optimization approaches, the search efficiency for a genetic algorithm is vital in finding desired solutions successfully in a multi-dimensional model space. A binary-encoding genetic algorithm is hardly ever used to resolve an optimization problem such as a simple geophysical inversion with only three unknowns. The encoding mechanism, genetic operators, and population size of the genetic algorithm greatly affect search processes in the evolution. It is clear that improved operators and proper population size promote the convergence. Nevertheless, not all genetic operations perform perfectly while searching under either a uniform binary or a decimal encoding system. With the binary encoding mechanism, the crossover scheme may produce more new individuals than with the decimal encoding. On the other hand, the mutation scheme in a decimal encoding system will create new genes larger in scope than those in the binary encoding. This paper discusses approaches of exploiting the search potential of genetic operations in the two encoding systems and presents an approach with a hybrid-encoding mechanism, multi-point crossover, and dynamic population size for geophysical inversion. We present a method that is based on the routine in which the mutation operation is conducted in the decimal code and multi-point crossover operation in the binary code. The mix-encoding algorithm is called the hybrid-encoding genetic algorithm (HEGA). HEGA provides better genes with a higher probability by a mutation operator and improves genetic algorithms in resolving complicated geophysical inverse problems. Another significant result is that final solution is determined by the average model derived from multiple trials instead of one computation due to the randomness in a genetic algorithm procedure. These advantages were demonstrated by synthetic and real-world examples of inversion of potential-field data. ?? 2005 Elsevier Ltd. All rights reserved.

Computers & Geosciences↗

USGS2006

Simple equations guide high-frequency surface-wave investigation techniques

We discuss five useful equations related to high-frequency surface-wave techniques and their implications in practice. These equations are theoretical results from published literature regarding source selection, data-acquisition parameters, resolution of a dispersion curve image in the frequency-velocity domain, and the cut-off frequency of high modes. The first equation suggests Rayleigh waves appear in the shortest offset when a source is located on the ground surface, which supports our observations that surface impact sources are the best source for surface-wave techniques. The second and third equations, based on the layered earth model, reveal a relationship between the optimal nearest offset in Rayleigh-wave data acquisition and seismic setting - the observed maximum and minimum phase velocities, and the maximum wavelength. Comparison among data acquired with different offsets at one test site confirms the better data were acquired with the suggested optimal nearest offset. The fourth equation illustrates that resolution of a dispersion curve image at a given frequency is directly proportional to the product of a length of a geophone array and the frequency. We used real-world data to verify the fourth equation. The last equation shows that the cut-off frequency of high modes of Love waves for a two-layer model is determined by shear-wave velocities and the thickness of the top layer. We applied this equation to Rayleigh waves and multi-layer models with the average velocity and obtained encouraging results. This equation not only endows with a criterion to distinguish high modes from numerical artifacts but also provides a straightforward means to resolve the depth to the half space of a layered earth model. ?? 2005 Elsevier Ltd. All rights reserved.

Soil Dynamics and Earthquake Engineering↗

USGS2004

Delineation of a collapse feature in a noisy environment using a multichannel surface wave technique

A collapse developed at Calvert Cliffs Nuclear Power Plant, Maryland, in early 2001. The location of the collapse was over a groundwater drainage system pipe buried at an elevation of +0??9 m (reference is to Chesapeake Bay level). The cause of the collapse was a subsurface drain pipe that collapsed because of saltwater corrosion of the corrugated metal pipe. The inflow/outflow of sea water and groundwater flow caused soil to be removed from the area where the pipe collapsed. To prevent damage to nearby structures, the collapse was quickly filled with uncompacted sand and gravel (???36000 kg). However, the plant had an immediate need to determine whether more underground voids existed. A high-frequency multichannel surface-wave survey technique was conducted to define the zone affected by the collapse. Although the surface-wave survey at Calvert Cliffs Nuclear Power Plant was conducted at a noise level 50-100 times higher than the normal environment for a shallow seismic survey, the shear (S)-wave velocity field calculated from surface-wave data delineated a possible zone affected by the collapse. The S-wave velocity field showed chimney-shaped low-velocity anomalies that were directly related to the collapse. Based on S-wave velocity field maps, a potential zone affected by the collapse was tentatively defined.

Geotechnique↗

USGS2003

Concentrations and co-occurrence correlations of 88 volatile organic compounds (VOCs) in the ambient air of 13 semi-rural to urban locations in the United States

The ambient air concentrations of 88 volatile organic compounds were determined in samples taken at 13 semi-rural to urban locations in Maine, Massachusetts, New Jersey, Pennsylvania, Ohio, Illinois, Louisiana, and California. The sampling periods ranged from 7 to 29 months, yielding a large data set with a total of 23,191 individual air concentration values, some of which were designated "ND" (not detected). For each compound at each sampling site, the air concentrations (ca, ppbV) are reported in terms of means, medians, and means of the detected values. The analytical method utilized adsorption/thermal desorption with air-sampling cartridges. The analytes included numerous halogenated alkanes, halogenated alkenes, ethers, alcohols, nitriles, esters, ketones, aromatics, a disulfide, and a furan. At some sites, the air concentrations of the gasoline-related aromatic compounds and the gasoline additive methyl tert-butyl ether were seasonally dependent, with concentrations that maximized in the winter. For each site studied here, the concentrations of some compounds were highly correlated one with another (e.g., the BTEX group (benzene, toluene, ethylbenzene, and the xylenes). Other aromatic compounds were also all generally correlated with one another, while the concentrations of other compound pairs were not correlated (e.g., benzene was not correlated with CFC-12). The concentrations found for the BTEX group were generally lower than the values that have been previously reported for urbanized and industrialized areas of other nations. ?? 2003 Elsevier Ltd. All rights reserved.

Atmospheric Environment↗

USGS1999

Preliminary study on the mode of occurrence of arsenic in high arsenic coals from southwest Guizhou Province

Coal samples from high arsenic coal areas have been analyzed by electron micropbe analyzer (EMPA), scanning electron microscopy with an energy dispersive X-ray analyzer (SEM-EDX), X-ray diffinction analysis (XRD), low temperature ashing (LTA), transmission electron microscopy (TEM), X-ray absorption fine structure (XAFS), instnunent neutron activation analysis (INAA) and wet chemical analysis. Although some As-bearing minerals such as pyrite, arsenopyrite, realgar (?), As-bearing sulfate, and As-bearing clays are found in the high arsenic coals, their contents do not account for the abundance of arsenic in the some coals. Analysis of the coal indicates that arsenic exists mainly in the form of A s 5+ and A s 3+ , combined with compounds in the organic matrix. The occurrence of such exceptionally high arsenic contents in coal and the fact that the arsenic is dominantly organically associated are unique observations. The modes of occurrence of arsenic in high As-coals are discussed.

Science in China, Series D: Earth Sciences↗