Geology topics

H. Zhang

Publications and source records attributed to H. Zhang.

13 recordsLinked to original sources

USGS2022

Using ensemble data assimilation to estimate transient hydrologic exchange flow under highly dynamic flow conditions

Quantifying dynamic hydrologic exchange flows (HEFs) within river corridors that experience high-frequency flow variations caused by dam regulations is important for understanding the biogeochemical processes at the river water and groundwater interfaces. Heat has been widely used as a tracer to infer steady-state flow velocities through analytical solutions of heat transport defined by the diurnal temperature signals. Under sub-daily dynamic flow conditions, however, such analytical solutions are not applicable due to the violation of their fundamental assumptions. In this study, we developed a data assimilation-based approach to estimate the sub-daily flux under highly dynamic flow conditions using multi-depth temperature observations at a 5-min resolution. If the hydraulic gradient is measured, Darcy's law was used to calculate the flux with permeability estimated from temperature responses below the riverbed. Otherwise, flux was estimated directly by assimilating multi-depth temperature data at 1- or 2-hr time intervals assuming one-dimensional flow and heat transport governing equation. By comparing estimated fluxes with model-generated synthetic truth, we demonstrated that both schemes have robust performance in estimating fluxes under highly dynamic flow conditions. This data assimilation-based flux estimation method was able to capture the vertical sub-daily fluxes using multi-depth high-resolution temperature data alone, even in the presence of multi-dimensional flow. This approach has been successfully applied to real field temperature data collected at the Hanford site, which experiences highly dynamic HEFs. Our study shows the promise of adopting distributed 1-D temperature monitoring to capture spatial and temporal exchange dynamics in river corridors at a watershed scale or beyond.

Washington↗

USGS2010

A California statewide three-dimensional seismic velocity model from both absolute and differential times

We obtain a seismic velocity model of the California crust and uppermost mantle using a regional-scale double-difference tomography algorithm. We begin by using absolute arrival-time picks to solve for a coarse three-dimensional (3D) P velocity (V P ) model with a uniform 30 km horizontal node spacing, which we then use as the starting model for a finer-scale inversion using double-difference tomography applied to absolute and differential pick times. For computational reasons, we split the state into 5 subregions with a grid spacing of 10 to 20 km and assemble our final statewide V P model by stitching together these local models. We also solve for a statewide S-wave model using S picks from both the Southern California Seismic Network and USArray, assuming a starting model based on the V P results and a V P =V S ratio of 1.732. Our new model has improved areal coverage compared with previous models, extending 570 km in the SW-NE directionand 1320 km in the NW-SE direction. It also extends to greater depth due to the inclusion of substantial data at large epicentral distances. Our V P model generally agrees with previous separate regional models for northern and southern California, but we also observe some new features, such as high-velocity anomalies at shallow depths in the Klamath Mountains and Mount Shasta area, somewhat slow velocities in the northern Coast Ranges, and slow anomalies beneath the Sierra Nevada at midcrustal and greater depths. This model can be applied to a variety of regional-scale studies in California, such as developing a unified statewide earthquake location catalog and performing regional waveform modeling.

Bulletin of the Seismological Society of America↗

USGS2009

Joint inversion for Vp, Vs, and Vp/Vs at SAFOD, Parkfield, California

We refined the three-dimensional (3-D) Vp, Vs and Vp/Vs models around the San Andreas Fault Observatory at Depth (SAFOD) site using a new double-difference (DD) seismic tomography code (tomoDDPS) that simultaneously solves for earthquake locations and all three velocity models using both absolute and differential P, S, and S-P times. This new method is able to provide a more robust Vp/Vs model than that from the original DD tomography code (tomoDD), obtained simply by dividing Vp by Vs. For the new inversion, waveform cross-correlation times for earthquakes from 2001 to 2002 were also used, in addition to arrival times from earthquakes and explosions in the region. The Vp values extracted from the model along the SAFOD trajectory match well with the borehole log data, providing in situ confirmation of our results. Similar to previous tomographic studies, the 3-D structure around Parkfield is dominated by the velocity contrast across the San Andreas Fault (SAF). In both the Vp and Vs models, there is a clear low-velocity zone as deep as 7 km along the SAF trace, compatible with the findings from fault zone guided waves. There is a high Vp/Vs anomaly zone on the southwest side of the SAF trace that is about 1-2 km wide and extends as deep as 4 km, which is interpreted to be due to fluids and fractures in the package of sedimentary rocks abutting the Salinian basement rock to the southwest. The relocated earthquakes align beneath the northeast edge of this high Vp/Vs zone. We carried out a 2-D correlation analysis for an existing resistivity model and the corresponding profiles through our model, yielding a classification that distinguishes several major lithologies. ?? 2009 by the American Geophysical Union.

Geochemistry, Geophysics, Geosystems↗

USGS2009

Application of a rule-based model to estimate mercury exchange for three background biomes in the continental United States

Ecosystems that have low mercury (Hg) concentrations (i.e., not enriched or impactedbygeologic or anthropogenic processes) cover most of the terrestrial surface area of the earth yet their role as a net source or sink for atmospheric Hg is uncertain. Here we use empirical data to develop a rule-based model implemented within a geographic information system framework to estimate the spatial and temporal patterns of Hg flux for semiarid deserts, grasslands, and deciduous forests representing 45% of the continental United States. This exercise provides an indication of whether these ecosystems are a net source or sink for atmospheric Hg as well as a basis for recommendation of data to collect in future field sampling campaigns. Results indicated that soil alone was a small net source of atmospheric Hg and that emitted Hg could be accounted for based on Hg input by wet deposition. When foliar assimilation and wet deposition are added to the area estimate of soil Hg flux these biomes are a sink for atmospheric Hg. ?? 2009 American Chemical Society.

Environmental Science & Technology↗

USGS2009

Regional three-dimensional seismic velocity model of the crust and uppermost mantle of northern California

We present a three-dimensional (3D) tomographic model of the P wave velocity (Vp) structure of northern California. We employed a regional-scale double-difference tomography algorithm that incorporates a finite-difference travel time calculator and spatial smoothing constraints. Arrival times from earthquakes and travel times from controlled-source explosions, recorded at network and/or temporary stations, were inverted for Vp on a 3D grid with horizontal node spacing of 10 to 20 km and vertical node spacing of 3 to 8 km. Our model provides an unprecedented, comprehensive view of the regional-scale structure of northern California, putting many previously identified features into a broader regional context and improving the resolution of a number of them and revealing a number of new features, especially in the middle and lower crust, that have never before been reported. Examples of the former include the complex subducting Gorda slab, a steep, deeply penetrating fault beneath the Sacramento River Delta, crustal low-velocity zones beneath Geysers-Clear Lake and Long Valley, and the high-velocity ophiolite body underlying the Great Valley. Examples of the latter include mid-crustal low-velocity zones beneath Mount Shasta and north of Lake Tahoe. Copyright 2009 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

USGS2009

Web-based decision support and visualization tools for water quality management in the Chesapeake Bay watershed

Federal, State, and local water quality managers charged with restoring the Chesapeake Bay ecosystem require tools to maximize the impact of their limited resources. To address this need, the U.S. Geological Survey (USGS) and the Environmental Protection Agency's Chesapeake Bay Program (CBP) are developing a suite of Web-based tools called the Chesapeake Online Assessment Support Toolkit (COAST). The goal of COAST is to help CBP partners identify geographic areas where restoration activities would have the greatest effect, select the appropriate management strategies, and improve coordination and prioritization among partners. As part of the COAST suite of tools focused on environmental restoration, a water quality management visualization component called the Nutrient Yields Mapper (NYM) tool is being developed by USGS. The NYM tool is a web application that uses watershed yield estimates from USGS SPAtially Referenced Regressions On Watershed (SPARROW) attributes model (Schwarz et al., 2006) [6] to allow water quality managers to identify important sources of nitrogen and phosphorous within the Chesapeake Bay watershed. The NYM tool utilizes new open source technologies that have become popular in geospatial web development, including components such as OpenLayers and GeoServer. This paper presents examples of water quality data analysis based on nutrient type, source, yield, and area of interest using the NYM tool for the Chesapeake Bay watershed. In addition, we describe examples of map-based techniques for identifying high and low nutrient yield areas; web map engines; and data visualization and data management techniques.

Conference Paper↗

USGS2009

Using drift nets to capture early life stages and monitor spawning of the yangtze river chinese sturgeon (Acipenser sinensis)

A sampling system for capturing sturgeon eggs using a D-shaped bottom anchored drift net was used to capture early life stages (ELS) of Chinese sturgeon, Acipenser sinensis, and monitor annual spawning success at Yichang on the Yangtze River, 1996-2004, before and just after the Three Gorges Dam began operation. Captured were 96 875 ELS (early life stages: eggs, yolk-sac larvae = eleuthero embryos, and larvae); most were eggs and only 2477 were yolk-sac larvae. Most ELS were captured in the main river channel and inside the bend at the Yichang spawning reach. Yolk-sac larvae were captured for a maximum of 3 days after hatching began, indicating quick dispersal downstream. The back-calculated day of egg fertilization over the eight years indicated a maximum spawning window of 23 days (20 October-10 November). Spawning in all years was restricted temporally, occurred mostly at night and during one or two spawning periods, each lasting several days. The brief temporal spawning window may reduce egg predation by opportunistic predators by flooding the river bottom with millions of eggs. During 1996-2002, the percentage of fertilized eggs in an annual 20-egg sample was between 63.5 to 94.1%; however, in 2003 the percentage fertilized was only 23.8%. This sudden decline may be related to the altered environmental conditions at Yichang caused by operation of the Three Gorges Dam. Further studies are needed to monitor spawning and changes in egg fertilization in this threatened population. ?? 2009 Blackwell Verlag GmbH.

Conference Paper↗

USGS2007

Three-dimensional P wave velocity model for the San Francisco Bay region, California

[1] A new three-dimensional P wave velocity model for the greater San Francisco Bay region has been derived using the double-difference seismic tomography method, using data from about 5,500 chemical explosions or air gun blasts and approximately 6,000 earthquakes. The model region covers 140 km NE-SW by 240 km NW-SE, extending from 20 km south of Monterey to Santa Rosa and reaching from the Pacific coast to the edge of the Great Valley. Our model provides the first regional view of a number of basement highs that are imaged in the uppermost few kilometers of the model, and images a number of velocity anomaly lows associated with known Mesozoic and Cenozoic basins in the study area. High velocity ( V p > 6.5 km/s) features at ∼15-km depth beneath part of the edge of the Great Valley and along the San Francisco peninsula are interpreted as ophiolite bodies. The relocated earthquakes provide a clear picture of the geometry of the major faults in the region, illuminating fault dips that are generally consistent with previous studies. Ninety-five percent of the earthquakes have depths between 2.3 and 15.2 km, and the corresponding seismic velocities at the hypocenters range from 4.8 km/s (presumably corresponding to Franciscan basement or Mesozoic sedimentary rocks of the Great Valley Sequence) to 6.8 km/s. The top of the seismogenic zone is thus largely controlled by basement depth, but the base of the seismogenic zone is not restricted to seismic velocities of ≤6.3 km/s in this region, as had been previously proposed.

California↗

USGS2007

Factors that influence the extraction of polycyclic aromatic hydrocarbons from coal

Coal samples and carbonaceous mudstone were collected from the Huaibei coalfield, China, and experiments investigating the factors influencing the extraction of the sixteen US EPA (Environmental Protection Agency) priority polycyclic aromatic hydrocarbons (PAHs) were carried out. Different extraction times, solvents, and methods were used. Major interest was focused on finding optimum conditions for extracting the PAHs from coal. We conclude that (1) coal composition, including the H/C and O/C ratios, is an important factor for the distribution of PAHs in coals; (2) the total amount of EPA priority PAHs increases with increasing extraction time, 30 min being suitable for ultrasonic-assisted extraction and 24 h for Soxhlet extraction; (3) CS2 is effective in extracting low molecular weight PAHs, while CH2Cl2 is better for extracting high molecular weight PAHs (both are excellent extraction solvents vs hexane); (4) both Soxhlet and ultrasonic extraction showed a similar PAH concentration profile, but the ultrasonic method is less efficient. ?? 2007 American Chemical Society.

Energy and Fuels↗

USGS2006

Three-dimensional compressional wavespeed model, earthquake relocations, and focal mechanisms for the Parkfield, California, region

We present a new three-dimensional (3D) compressional vvavespeed (V p) model for the Parkfield region, taking advantage of the recent seismicity associated with the 2003 San Simeon and 2004 Parkfield earthquake sequences to provide increased model resolution compared to the work of Eberhart-Phillips and Michael (1993) (EPM93). Taking the EPM93 3D model as our starting model, we invert the arrival-time data from about 2100 earthquakes and 250 shots recorded on both permanent network and temporary stations in a region 130 km northeast-southwest by 120 km northwest-southeast. We include catalog picks and cross-correlation and catalog differential times in the inversion, using the double-difference tomography method of Zhang and Thurber (2003). The principal Vp features reported by EPM93 and Michelini and McEvilly (1991) are recovered, but with locally improved resolution along the San Andreas Fault (SAF) and near the active-source profiles. We image the previously identified strong wavespeed contrast (faster on the southwest side) across most of the length of the SAF, and we also improve the image of a high Vp body on the northeast side of the fault reported by EPM93. This narrow body is at about 5- to 12-km depth and extends approximately from the locked section of the SAP to the town of Parkfield. The footwall of the thrust fault responsible for the 1983 Coalinga earthquake is imaged as a northeast-dipping high wavespeed body. In between, relatively low wavespeeds (<5 km/sec) extend to as much as 10-km depth. We use this model to derive absolute locations for about 16,000 earthquakes from 1966 to 2005 and high-precision double-difference locations for 9,000 earthquakes from 1984 to 2005, and also to determine focal mechanisms for 446 earthquakes. These earthquake locations and mechanisms show that the seismogenic fault is a simple planar structure. The aftershock sequence of the 2004 mainshock concentrates into the same structures defined by the pre-2004 seismicity, confirming earlier observations (Waldhauser et al., 2004) that the seismicity pattern at Parkfield is long lived and persists through multiple cycles of mainshocks.

Bulletin of the Seismological Society of America↗

USGS2004

Fine-scale structure of the San Andreas fault zone and location of the SAFOD target earthquakes

We present results from the tomographic analysis of seismic data from the Parkfield area using three different inversion codes. The models provide a consistent view of the complex velocity structure in the vicinity of the San Andreas, including a sharp velocity contrast across the fault. We use the inversion results to assess our confidence in the absolute location accuracy of a potential target earthquake. We derive two types of accuracy estimates, one based on a consideration of the location differences from the three inversion methods, and the other based on the absolute location accuracy of “virtual earthquakes.” Location differences are on the order of 100–200 m horizontally and up to 500 m vertically. Bounds on the absolute location errors based on the “virtual earthquake” relocations are ∼50 m horizontally and vertically. The average of our locations places the target event epicenter within about 100 m of the SAF surface trace.

California↗

USGS2003

Atmospheric mercury emissions from mine wastes and surrounding geologically enriched terrains

Waste rock and ore associated with Hg, precious and base metal mining, and their surrounding host rocks are typically enriched in mercury relative to natural background concentrations (<0.1 ??g Hg g-1). Mercury fluxes to the atmosphere from mineralized areas can range from background rates (0-15 ng m-2 h-1) to tens of thousands of ng m-2 h-1. Mercury enriched substrate constitutes a long-term source of mercury to the global atmospheric mercury pool. Mercury emissions from substrate are influenced by light, temperature, precipitation, and substrate mercury concentration, and occur during the day and night. Light-enhanced emissions are driven by two processes: desorption of elemental mercury accumulated at the soil:air interface, and photo reduction of mercury containing phases. To determine the need for and effectiveness of regulatory controls on short-lived anthropogenic point sources the contribution of mercury from geologic non-point sources to the atmospheric mercury pool needs to be quantified. The atmospheric mercury contribution from small areas of mining disturbance with relatively high mercury concentrations are, in general, less than that from surrounding large areas of low levels of mercury enrichment. In the arid to semi-arid west-ern United States volatilization is the primary means by which mercury is released from enriched sites.

Conference Paper↗

USGS1999

Nevada STORMS project: Measurement of mercury emissions from naturally enriched surfaces

Diffuse anthropogenic and naturally mercury-enriched areas represent long-lived sources of elemental mercury to the atmosphere. The Nevada Study and Tests of the Release of Mercury From Soils (STORMS) project focused on the measurement of mercury emissions from a naturally enriched area. During the project, concurrent measurements of mercury fluxes from naturally mercury-enriched substrate were made September 1–4, 1997, using four micrometeorological methods and seven field flux chambers. Ambient air mercury concentrations ranged from 2 to nearly 200 ng m −3 indicating that the field site is a source of atmospheric mercury. The mean daytime mercury fluxes, during conditions of no precipitation, measured with field chambers were 50 to 360 ng m −2 h −1 , and with the micrometeorological methods were 230 to 600 ng m −2 h −1 . This wide range in mercury emission rates reflects differences in method experimental designs and local source strengths. Mercury fluxes measured by many field chambers were significantly different ( p <0.05) but linearly correlated. This indicates that field chambers responded similarly to environmental conditions, but differences in experimental design and site heterogeneity had a significant influence on the magnitude of mercury fluxes. Data developed during the field study demonstrated that field flux chambers are ideal for assessment of the physicochemical processes driving mercury flux and development of an understanding of the magnitude of the influence of individual factors on flux. In general, mean mercury fluxes measured with micrometeorological methods during daytime periods were nearly 3 times higher than mean fluxes measured with field flux chambers. Micrometeorological methods allow for derivation of a representative mercury flux occurring from an unconstrained system and provide an assessment of the actual magnitude and variability of fluxes occurring from an area.

Journal of Geophysical Research D: Atmospheres↗