SEARCH · Search USGS
Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Explore the source record for details and available documents.
This paper describes the importance of wave-current interaction in an inlet-estuary system. The three-dimensional, fully coupled, Coupled Ocean-Atmosphere-Wave-Sediment Transport (COAWST) modeling system was applied in Willapa Bay (Washington State) from 22 to 29 October 1998 that included a large storm event. To represent the interaction between waves and currents, the vortex-force method was used. Model results were compared with water elevations, currents, and wave measurements obtained by the U.S. Army Corp of Engineers. In general, a good agreement between field data and computed results was achieved, although some discrepancies were also observed in regard to wave peak directions in the most upstream station. Several numerical experiments that considered different forcing terms were run in order to identify the effects of each wind, tide, and wave-current interaction process. Comparison of the horizontal momentum balances results identified that wave-breaking-induced acceleration is one of the leading terms in the inlet area. The enhancement of the apparent bed roughness caused by waves also affected the values and distribution of the bottom shear stress. The pressure gradient showed significant changes with respect to the pure tidal case. During storm conditions the momentum balance in the inlet shares the characteristics of tidal-dominated and wave-dominated surf zone environments. The changes in the momentum balance caused by waves were manifested both in water level and current variations. The most relevant effect on hydrodynamics was a wave-induced setup in the inner part of the estuary.
We use three-dimensional elastic models to help guide the kinematic interpretation of crustal deformation associated with strike-slip faults. Deformation of the brittle upper crust in the vicinity of strike-slip fault systems is modeled with the assumption that upper crustal deformation is driven by the relative plate motion in the upper mantle. The driving motion is represented by displacement that is specified on the bottom of a 15-km-thick elastic upper crust everywhere except in a zone of finite width in the vicinity of the faults, which we term the “shear zone.” Stress-free basal boundary conditions are specified within the shear zone. The basal driving displacement is either pure strike slip or strike slip with a small oblique component, and the geometry of the fault system includes a single fault, several parallel faults, and overlapping en echelon faults. We examine the variations in deformation due to changes in the width of the shear zone and due to changes in the shear strength of the faults. In models with weak faults the width of the shear zone has a considerable effect on the surficial extent and amplitude of the vertical and horizontal deformation and on the amount of rotation around horizontal and vertical axes. Strong fault models have more localized deformation at the tip of the faults, and the deformation is partly distributed outside the fault zone. The dimensions of large basins along strike-slip faults, such as the Rukwa and Dead Sea basins, and the absence of uplift around pull-apart basins fit models with weak faults better than models with strong faults. Our models also suggest that the length-to-width ratio of pull-apart basins depends on the width of the shear zone and the shear strength of the faults and is not constant as previously suggested. We show that pure strike-slip motion can produce tectonic features, such as elongate half grabens along a single fault, rotated blocks at the ends of parallel faults, or extension perpendicular to overlapping en echelon faults, which can be misinterpreted to indicate a regional component of extension. Zones of subsidence or uplift can become wider than expected for transform plate boundaries when a minor component of oblique motion is added to a system of parallel strike-slip faults.
The pattern of historical ground deformation, seismicity, and crustal structure near Medicine Lake volcano illustrates a close relation between magmatism and tectonism near the margin of the Cascade volcanic chain and the Basin and Range tectonic province. Between leveling surveys in 1954 and 1989 the summit of Medicine Lake volcano subsided 389±43 mm with respect to a reference bench mark 40 km to the southwest (average rate = 11.1±1.2 mm/yr). A smaller survey across the summit caldera in 1988 suggests that the subsidence rate was 15–28 mm/yr during 1988–1989. Swarms of shallow earthquakes ( M ≤ 4.6) occurred in the region during August 1978, January–February 1981, and September 1988. Except for the 1988 swarm, which occurred beneath Medicine Lake caldera, most historical earthquakes were located at least 25 km from the summit. The spatial relation between subsidence and seismicity indicates (1) radially symmetric downwarping of the volcano's summit and flanks centered near the caldera and (2) downfaulting of the entire edifice along regional faults located 25–30 km from the summit. We propose that contemporary subsidence, seismicity, and faulting are caused by (1) loading of the crust by more than 600 km 3 of erupted products plus a large volume of mafic intrusives; (2) east‐west extension in the western Basin and Range province; and, to a lesser extent, (3) crystallization or withdrawal of magma beneath the volcano. Thermal weakening of the subvolcanic crust by mafic intrusions facilitates subsidence and influences the distribution of earthquakes. Subsidence occurs mainly by aseismic creep within 25 km of the summit, where the crust has been heated and weakened by intrusions, and by normal faulting during episodic earthquake swarms in surrounding, cooler terrain.
In this paper we examine the variations of the boreal summer season sea breeze circulation along the Florida panhandle coast from relatively high resolution (10 km) regional climate model integrations. The 23 year climatology (1979–2001) of the multidecadal dynamically downscaled simulations forced by the National Centers for Environmental Prediction–Department of Energy (NCEP‐DOE) Reanalysis II at the lateral boundaries verify quite well with the observed climatology. The variations at diurnal and interannual time scales are also well simulated with respect to the observations. We show from composite analyses made from these downscaled simulations that sea breezes in northwestern Florida are associated with changes in the size of the Atlantic Warm Pool (AWP) on interannual time scales. In large AWP years when the North Atlantic Subtropical High becomes weaker and moves further eastward relative to the small AWP years, a large part of the southeast U.S. including Florida comes under the influence of relatively strong anomalous low‐level northerly flow and large‐scale subsidence consistent with the theory of the Sverdrup balance. This tends to suppress the diurnal convection over the Florida panhandle coast in large AWP years. This study is also an illustration of the benefit of dynamic downscaling in understanding the low‐frequency variations of the sea breeze.
Analysis of two absorption features near 3 μ m in the lunar reflectance spectrum, observed by the orbiting M 3 spectrometer and interpreted as being due to OH and H 2 O, is presented, and the results are used to discuss the processes producing these molecules. This analysis focuses on the dependence of the absorptions on lunar physical properties, including composition, illumination, latitude, and temperature. Solar wind proton-induced hydroxylation is proposed as the creation process, and its products could be a source for other reported types of hydrogen-rich material and water. The irregular and damaged fine-grained lunar soil seems especially adapted for trapping solar wind protons and forming OH owing to abundant dangling oxygen bonds. The M 3 data reveal that the strengths of the two absorptions are correlated and widespread, and both are correlated with lunar composition but in different ways. Feldspathic material seems richer in OH. These results seem to rule out water from the lunar interior and cometary infall as major sources. There appear to be correlations of apparent band strengths with time of day and lighting conditions. However, thermal emission from the Moon reduces the apparent strengths of the M 3 absorptions, and its removal is not yet completely successful. Further, many of the lunar physical properties are themselves intercorrelated, and so separating these dependencies on the absorptions is difficult, due to the incomplete M 3 data set. This process should also operate on other airless silicate surfaces, such as Mercury and Vesta, which will be visited by the Dawn spacecraft in mid-2011.
Rapid changes in ground tilt and GPS positions on Kīlauea Volcano, Hawai'i, are interpreted as resulting from a shallow, two‐segment dike intrusion into the east rift zone that began at 1217 UTC (0217 HST) on 17 June 2007 and lasted almost 3 days. As a result of the intrusion, a very small volume of basalt (about 1500 m 3 ) erupted on 19 June. Northward tilt at a coastal tiltmeter, subsidence of south flank GPS sites, southeastward displacements at southwestern flank GPS sites, and a swarm of flank earthquakes suggest that a slow slip event occurred on the décollement beneath Kīlauea's south flank concurrent with the rift intrusion. We use 4 min GPS positions that include estimates of time‐dependent tropospheric gradients and ground tilt data to study the spatial and temporal relationships between the two inferred shallow, steeply dipping dike segments extending from the surface to about 2 km depth and décollement slip at 8 km depth. We invert for the temporal evolution of distributed dike opening and décollement slip in independent inversions at each time step using a nonnegative least squares algorithm. On the basis of these inversions, the intrusion occurred in two stages that correspond spatially and temporally with concentrated rift zone seismicity. The dike opening began on the western of the two segments before jumping to the eastern segment, where the majority of opening accumulated. Dike opening preceded the start of décollement slip at an 84% confidence level; the latter is indicated by the onset of northward tilt of a coastal tiltmeter. Displacements at southwest flank GPS sites began about 18 h later and are interpreted as resulting from slow slip on the southwestern flank. Additional constraints on the evolution of the intrusion and décollement slip come from inversion of an Envisat interferogram that spans the intrusion until 0822 UTC on 18 June 2007, combined with GPS and tilt data. This inversion shows that up to 0822 UTC on 18 June, décollement slip is only required in a limited region offshore of Ka'ena Point. A similar inversion of the complete event, which includes GPS and tilt data up to 21 June and a second Envisat interferogram spanning the complete intrusion until 21 June, shows décollement slip spread westward across the south flank. This may suggest westward migration of the décollement slip as the event progressed.
New seismic and geodetic data from Costa Rica provide insight into seismogenic zone processes in Central America, where the Cocos and Caribbean plates converge. Seismic data are from combined land and ocean bottom deployments in the Nicoya peninsula in northern Costa Rica and near the Osa peninsula in southern Costa Rica. In Nicoya, inversion of GPS data suggests two locked patches centered at 14 ± 2 and 39 ± 6 km depth. Interplate microseismicity is concentrated in the more freely slipping intermediate zone, suggesting that small interseismic earthquakes may not accurately outline the updip limit of the seismogenic zone, the rupture zone for future large earthquakes, at least over the short (∼1 year) observation period. We also estimate northwest motion of a coastal “sliver block” at 8 ± 3 mm/yr, probably related to oblique convergence. In the Osa region to the south, convergence is orthogonal to the trench. Cocos-Caribbean relative motion is partitioned here, with ∼8 cm/yr on the Cocos-Panama block boundary (including a component of permanent shortening across the Fila Costeña fold and thrust belt) and ∼1 cm/yr on the Panama block–Caribbean boundary. The GPS data suggest that the Cocos plate–Panama block boundary is completely locked from ∼10–50 km depth. This large locked zone, as well as associated forearc and back-arc deformation, may be related to subduction of the shallow Cocos Ridge and/or younger lithosphere compared to Nicoya, with consequent higher coupling and compressive stress in the direction of plate convergence.
Prominent low‐albedo wind streaks issue from embayments at the north end of Victoria crater in Meridiani Planum, the site of surface investigations by the Mars Exploration Rover Opportunity. Opportunity was sent to examine the darkest of these wind streaks and compare it to the adjacent soil in the brighter corridor between the streaks. Two nearby sites inside and outside of the streak were selected for visible and infrared remote sensing and for elemental abundance measurements and images of the microscopic soil morphology. The rover was next sent to study the interaction of the wind streak with an obstacle, a small rock near the source of the streak. Opportunity's observations suggest that the Victoria wind streaks are deposits of basaltic sand blown out of the crater from the dark dunes nestled below the crater rim, particularly at the base of the alcove leading up to the darkest streak. No local sources of sand have been identified within the Victoria crater wall rock, suggesting that the sand is supplied from outside the crater and is presently escaping its temporary topographic trap. This process presents a possible explanation for the serrated margin of Victoria crater through abrasion of the soft rock as trapped sands are blown out of the crater and carve alcoves under various seasonal winds. Some indication of the rate of activity of these wind streaks will be provided by continued monitoring of the erasure of the rover's tracks.
The emission rate of carbon dioxide escaping from the summit of Kīlauea Volcano, Hawaiʻi, proved highly variable, averaging 4900 ± 2000 metric tons per day (t/d) in June–July 2003 during a period of summit inflation. These results were obtained by combining over 90 measurements of COSPEC-derived SO 2 emission rates with synchronous CO 2 /SO 2 ratios of the volcanic gas plume along the summit COSPEC traverse. The results are lower than the CO 2 emission rate of 8500 ± 300 t/d measured by the same method in 1995–1999 during a period of long-term summit deflation [Gerlach, T.M., McGee, K.A., Elias, T., Sutton, A.J. and Doukas, M.P., 2002. Carbon dioxide emission rate of Kīlauea Volcano: Implications for primary magma and the summit reservoir. Journal of Geophysical Research-Solid Earth, 107(B9): art. no.-2189.]. Analysis of the data indicates that the emission rates of the present study likely reflect changes in the magma supply rate and residence time in the summit reservoir. It is also likely that emission rates during the inflation period were heavily influenced by SO 2 pulses emitted adjacent to the COSPEC traverse, which biased CO 2 /SO 2 ratios towards low values that may be unrepresentative of the global summit gas plume. We conclude that the SO 2 pulses are consequences of summit re-inflation under way since 2003 and that CO 2 emission rates remain comparable to, but more variable than, those measured prior to re-inflation.
Geologists are in a quandary over the correct interpretation of paleomagnetic data for Jurassic rocks of the North American plate. Conflicting reference paleopoles and alternate methods of constructing apparent polar wander (APW) paths have led to a controversy regarding the configuration of Jurassic APW for cratonic North America. These differences have been recently disputed at meetings (see Eos , Spring Meeting Supplement, April 7, 1992, p. 94) and in an exchange of letters in the Journal of Geophysical Research-Solid Earth [ Butler et al ., 1992; Van Fossen and Kent , 1992a]. At stake is important information concerning the Jurassic paleogeography of North America and the whole of Pangea, as well as the nature and driving mechanism of plate movements. In addition, cratonic paleopoles provide a reference frame for measuring relative displacements of tectonostratigraphic terranes.
A multitemporal 1 km advanced very high resolution radiometer (AVHRR) land cover analysis approach was used as the basis for regional land cover mapping, fire disturbance-regeneration, and multiresolution land cover scaling studies in the boreal forest ecosystem of central Canada. The land cover classification was developed by using regional field observations from ground and low-level aircraft transits to analyze spectral-temporal clusters that were derived from an unsupervised cluster analysis of monthly normalized difference vegetation index (NDVI) image composites (April-September 1992). Quantitative areal proportions of the major boreal forest components were determined for a 821 km ?? 619 km region, ranging from the southern grasslands-boreal forest ecotone to the northern boreal transitional forest. The boreal wetlands (mostly lowland black spruce, tamarack, mosses, fens, and bogs) occupied approximately 33% of the region, while lakes accounted for another 13%. Upland mixed coniferous-deciduous forests represented 23% of the ecosystem. A SW-NE productivity gradient across the region is manifested by three levels of tree stand density for both the boreal wetland conifer and the mixed forest classes, which are generally aligned with isopleths of regional growing degree days. Approximately 30% of the region was directly affected by fire disturbance within the preceding 30-35 years, especially in the Canadian Shield Zone where large fire-regeneration patterns contribute to the heterogeneous boreal landscape. Intercomparisons with land cover classifications derived from 30-m Landsat Thematic Mapper (TM) data provided important insights into the relative accuracy of the 1 km AVHRR land cover classification. Primarily due to the multitemporal NDVI image compositing process, the 1 km AVHRR land cover classes have an effective spatial resolution in the 3-4 km range; therefore fens, bogs, small water bodies, and small patches of dry jack pine cannot be resolved within the wet conifer mosaic. Major differences in the 1-km AVHRR and 30-m Landsat TM-derived land cover classes are most likely due to differences in the spatial resolution of the data sets. In general, the 1 km AVHRR land cover classes are vegetation mosaics consisting of mixed combinations of the Landsat classes. Detailed mapping of the global boreal forest with this approach will benefit from algorithms for cloud screening and to atmospherically correct reflectance data for both aerosol and water vapor effects. We believe that this 1 km AVHRR land cover analysis provides new and useful information for regional water, energy, carbon, and trace gases studies in BOREAS, especially given the significant spatial variability in land cover type and associated biophysical land cover parameters (e.g., albedo, leaf area index, FPAR, and surface roughness). Multiresolution land cover comparisons (30 m, 1 km, and 100 km grid cells) also illustrated how heterogeneous landscape patterns are represented in land cover maps with differing spatial scales and provided insights on the requirements and challenges for parameterizing landscape heterogeneity as part of land surface process research.
The foundational ecosystem processes of gross primary production (GPP) and ecosystem respiration (ER) cannot be measured directly but can be modeled in aquatic ecosystems from subdaily patterns of oxygen (O 2 ) concentrations. Because rivers and streams constantly exchange O 2 with the atmosphere, models must either use empirical estimates of the gas exchange rate coefficient ( K 600 ) or solve for all three parameters (GPP, ER, and K 600 ) simultaneously. Empirical measurements of K 600 require substantial field work and can still be inaccurate. Three-parameter models have suffered from equifinality, where good fits to O 2 data are achieved by many different parameter values, some unrealistic. We developed a new three-parameter, multiday model that ensures similar values for K 600 among days with similar physical conditions (e.g., discharge). Our new model overcomes the equifinality problem by (1) flexibly relating K 600 to discharge while permitting moderate daily deviations and (2) avoiding the oft-violated assumption that residuals in O 2 predictions are uncorrelated. We implemented this hierarchical state-space model and several competitor models in an open-source R package, streamMetabolizer . We then tested the models against both simulated and field data. Our new model reduces error by as much as 70% in daily estimates of K 600 , GPP, and ER. Further, accuracy benefits of multiday data sets require as few as 3 days of data. This approach facilitates more accurate metabolism estimates for more streams and days, enabling researchers to better quantify carbon fluxes, compare streams by their metabolic regimes, and investigate controls on aquatic activity.
Very low frequency (VLF, 0.001–0.005 Hz) waves are important drivers of flooding of low-lying coral reef-islands. In particular, VLF wave resonance is known to drive large wave runup and subsequent overwash. Using a 5 month data set of water levels and waves collected along a cross-reef transect on Roi-Namur Island in the Republic of the Marshall Islands, the observed VLF motions were categorized into four different classes: (1) resonant, (2) (nonresonant) standing, (3) progressive-growing, and (4) progressive-dissipative waves. Each VLF class is set by the reef flat water depth and, in the case of resonance, the incident-band offshore wave period. Using an improved method to identify VLF wave resonance, we find that VLF wave resonance caused prolonged (∼0.5–6.0 h), large-amplitude water surface oscillations at the inner reef flat ranging in wave height from 0.14 to 0.83 m. It was induced by relatively long-period, grouped, incident-band waves, and occurred under both storm and nonstorm conditions. Moreover, observed resonant VLF waves had nonlinear, bore-like wave shapes, which likely have a larger impact on the shoreline than regular, sinusoidal waveforms. As an alternative technique to the commonly used Fast Fourier Transformation, we propose the Hilbert-Huang Transformation that is more computationally expensive but can capture the wave shape more accurately. This research demonstrates that understanding VLF waves on reef flats is important for evaluating coastal flooding hazards.
The New Madrid seismic zone (NMSZ), which experienced severe liquefaction during the great New Madrid, Missouri, earthquakes of 1811 and 1812 as well as during several prehistoric earthquakes, is a superb laboratory for the study of world-class, earthquake-induced liquefaction features and their use in paleoseismology. In seismically active regions like the NMSZ, frequent large earthquakes can produce a complex record of liquefaction events that is difficult to interpret. Lessons learned studying liquefaction features in the NMSZ may help to unravel the paleoseismic record in other seismically active regions. Soil characteristics of liquefaction features, as well as their structural and stratigraphic relations to Native American occupation horizons and other cultural features, can help to distinguish prehistoric liquefaction features from historic features. In addition, analyses of artifact assemblages and botanical content of cultural horizons can help to narrow the age ranges of liquefaction features. Future research should focus on methods for defining source areas and estimating magnitudes of prehistoric earthquakes from liquefaction features. Also, new methods for dating liquefaction features are needed.
Measurements of aeolian sediment transport support our understanding of mineral dust impacts on Earth and human systems and assessments of aeolian process sensitivities to global environmental change. However, sample design principles are often overlooked in aeolian research. Here, we use high‐density field measurements of sediment mass flux across land use and land cover types to examine sample size and power effects on detecting change in aeolian transport. Temporal variances were 1.6 to 10.1 times the magnitude of spatial variances in aeolian transport for six study sites. Differences in transport were detectable for >67% of comparisons among sites using ~27 samples. Failure to detect change with smaller sample sizes suggests that aeolian transport measurements and monitoring are much more uncertain than recognized. We show how small and selective sampling, common in aeolian research, gives the false impression that differences in aeolian transport can be detected, potentially undermining inferences about process and impacting reproducibility of aeolian research.
This synthesis addresses the vulnerability of the North American high-latitude soil organic carbon (SOC) pool to climate change. Disturbances caused by climate warming in arctic, subarctic, and boreal environments can result in significant redistribution of C among major reservoirs with potential global impacts. We divide the current northern high-latitude SOC pools into (1) near-surface soils where SOC is affected by seasonal freeze-thaw processes and changes in moisture status, and (2) deeper permafrost and peatland strata down to several tens of meters depth where SOC is usually not affected by short-term changes. We address key factors (permafrost, vegetation, hydrology, paleoenvironmental history) and processes (C input, storage, decomposition, and output) responsible for the formation of the large high-latitude SOC pool in North America and highlight how climate-related disturbances could alter this pool's character and size. Press disturbances of relatively slow but persistent nature such as top-down thawing of permafrost, and changes in hydrology, microbiological communities, pedological processes, and vegetation types, as well as pulse disturbances of relatively rapid and local nature such as wildfires and thermokarst, could substantially impact SOC stocks. Ongoing climate warming in the North American high-latitude region could result in crossing environmental thresholds, thereby accelerating press disturbances and increasingly triggering pulse disturbances and eventually affecting the C source/sink net character of northern high-latitude soils. Finally, we assess postdisturbance feedbacks, models, and predictions for the northern high-latitude SOC pool, and discuss data and research gaps to be addressed by future research.
We investigate the dynamic interactions between the Southern San Andreas Fault (SSAF) and a proximal normal fault (NF) beneath the Salton Sea in southern California. The NF, positioned near the SSAF terminus at Bombay Beach, exhibits 11–15 displacement events across 14 stratigraphic sequences, with a range of 0.2–1.4 m of vertical offset since ∼2–3 ka. Notably, four of these events may align temporally with SSAF earthquakes, raising questions about the possible interplay between the two faults. Utilizing dynamic rupture models, we analyze the coseismic interactions between the SSAF and NF, addressing under what conditions the SSAF induces slip on the NF. Our findings reveal that a suite of SSAF ruptures, particularly those propagating from north to south, can trigger slip on the normal fault and replicate observed vertical offsets. If the SSAF extends beneath the Salton Sea, earthquakes originating south of the NF intersection are less likely to trigger normal fault slip, although we cannot exclude this possibility. Some SSAF ruptures do not trigger discernible slip on the NF, rendering such events undetectable in the stratigraphic record. Our research contributes toward discussions regarding the seismic hazard in southern California, shedding light on the interplay between the SSAF and NF.