Geology topics

C. Chen

Publications and source records attributed to C. Chen.

14 recordsLinked to original sources

USGS2011

Pseudospectral modeling and dispersion analysis of Rayleigh waves in viscoelastic media

Multichannel Analysis of Surface Waves (MASW) is one of the most widely used techniques in environmental and engineering geophysics to determine shear-wave velocities and dynamic properties, which is based on the elastic layered system theory. Wave propagation in the Earth, however, has been recognized as viscoelastic and the propagation of Rayleigh waves presents substantial differences in viscoelastic media as compared with elastic media. Therefore, it is necessary to carry out numerical simulation and dispersion analysis of Rayleigh waves in viscoelastic media to better understand Rayleigh-wave behaviors in the real world. We apply a pseudospectral method to the calculation of the spatial derivatives using a Chebyshev difference operator in the vertical direction and a Fourier difference operator in the horizontal direction based on the velocity-stress elastodynamic equations and relations of linear viscoelastic solids. This approach stretches the spatial discrete grid to have a minimum grid size near the free surface so that high accuracy and resolution are achieved at the free surface, which allows an effective incorporation of the free surface boundary conditions since the Chebyshev method is nonperiodic. We first use an elastic homogeneous half-space model to demonstrate the accuracy of the pseudospectral method comparing with the analytical solution, and verify the correctness of the numerical modeling results for a viscoelastic half-space comparing the phase velocities of Rayleigh wave between the theoretical values and the dispersive image generated by high-resolution linear Radon transform. We then simulate three types of two-layer models to analyze dispersive-energy characteristics for near-surface applications. Results demonstrate that the phase velocity of Rayleigh waves in viscoelastic media is relatively higher than in elastic media and the fundamental mode increases by 10-16% when the frequency is above 10. Hz due to the velocity dispersion of P and S waves. ?? 2011 Elsevier Ltd.

Soil Dynamics and Earthquake Engineering↗

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↗

USGS2005

Accurate elevation and normal moveout corrections of seismic reflection data on rugged topography

The application of the seismic reflection method is often limited in areas of complex terrain. The problem is the incorrect correction of time shifts caused by topography. To apply normal moveout (NMO) correction to reflection data correctly, static corrections are necessary to be applied in advance for the compensation of the time distortions of topography and the time delays from near-surface weathered layers. For environment and engineering investigation, weathered layers are our targets, so that the static correction mainly serves the adjustment of time shifts due to an undulating surface. In practice, seismic reflected raypaths are assumed to be almost vertical through the near-surface layers because they have much lower velocities than layers below. This assumption is acceptable in most cases since it results in little residual error for small elevation changes and small offsets in reflection events. Although static algorithms based on choosing a floating datum related to common midpoint gathers or residual surface-consistent functions are available and effective, errors caused by the assumption of vertical raypaths often generate pseudo-indications of structures. This paper presents the comparison of applying corrections based on the vertical raypaths and bias (non-vertical) raypaths. It also provides an approach of combining elevation and NMO corrections. The advantages of the approach are demonstrated by synthetic and real-world examples of multi-coverage seismic reflection surveys on rough topography. ?? The Royal Society of New Zealand 2005.

New Zealand Journal of Geology and Geophysics↗

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↗

USGS2004

Geophysical techniques in detection to river embankments - A case study: To locate sites of potential leaks using surface-wave and electrical methods

Geophysical technologies are very effective in environmental, engineering and groundwater applications. Parameters of delineating nature of near-surface materials such as compressional-wave velocity, shear-wave velocity can be obtained using shallow seismic methods. Electric methods are primary approaches for investigating groundwater and detecting leakage. Both of methods are applied to detect embankment in hope of obtaining evidences of the strength and moisture inside the body. A technological experiment has done for detecting and discovering the hidden troubles in the embankment of Yangtze River, Songzi, Hubei, China in 2003. Surface-wave and DC multi-channel array resistivity sounding techniques were used to detect hidden trouble inside and under dike like pipe-seeps. This paper discusses the exploration strategy and the effect of geological characteristics. A practical approach of combining seismic and electric resistivity measurements was applied to locate potential pipe-seeps in embankment in the experiment. The method presents a potential leak factor based on the shear-wave velocity and the resistivity of the medium to evaluate anomalies. An anomaly found in a segment of embankment detected was verified, where occurred a pipe-seep during the 98' flooding.

Conference Paper↗

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↗

USGS2001

Geological and geochemical characteristics of high arsenic coals from endemic arsenosis areas in southwestern Guizhou Province, China

Southwest Guizhou Province is one of the most important areas of disseminated, sediment-hosted-type Au deposits in China and is an important area of coal production. The chemistry of most of the coals in SW Guizhou is similar to those in other parts of China. Their As content is near the Chinese coal average, but some local, small coal mines contain high As coals. The highest As content is up to 3.5 wt.% in the coal. The use of high As coals has caused in excess of 3000 cases of As poisoning in several villages. The high As coals are in the Longtan formation, which is an alternating marine facies and terrestrial facies. The coals are distributed on both sides of faults that parallel the regional anticlinal axis. The As content of coal is higher closer to the fault plane. The As content of coal changes greatly in different coal beds and different locations of the same bed. Geological structures such as anticlines, faults and sedimentary strata control the distribution of high As coals. Small Au deposits as well as Sb, Hg, and Th mineralization, are found near the high As coals. Although some As-bearing minerals such as pyrite, arsenopyrite, realgar (?), As-bearing sulfate, As-bearing clays, and phosphate are found in the high As coals, their contents cannot account for the abundance of As in some coals. Analysis of the coal indicates that As mainly exists in the form of As5+ and As3+, perhaps, combined with organic compounds. The occurrence of such exceptionally high As contents in coal and the fact that the As is dominantly organically associated are unique observations. ?? 201 Elsevier Science Ltd.

Conference Paper↗

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↗