Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Research and Exploration”

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.

At least 1,189 records · Page 66Linked to original sources
USGS2006

Active dust devils in Gusev crater, Mars: Observations from the Mars Exploration Rover Spirit

A full dust devil "season" was observed from Spirit from 10 March 2005 (sol 421, first active dust devil observed) to 12 December 2005 (sol 691, last dust devil seen); this corresponds to the period Ls 173.2?? to 339.5??, or the southern spring and summer on Mars. Thermal Emission Spectrometer data suggest a correlation between high surface temperatures and a positive thermal gradient with active dust devils in Gusev and that Spirit landed in the waning stages of a dust devil season as temperatures decreased. 533 active dust devils were observed, enabling new characterizations; they ranged in diameter from 2 to 276 m, with most in the range of 10-20 m in diameter, and occurred from about 0930 to 1630 hours local true solar time (with the maximum forming around 1300 hours) and a peak occurrence in southern late spring (Ls ??? 250??). Horizontal speeds of the dust devils ranged from <1 to 21 m/s, while vertical wind speeds within the dust devils ranged from 0.2 to 8.8 m/s. These data, when combined with estimates of the dust content within the dust devils, yield dust fluxes of 3.95 ?? 10-9 to 4.59-4 kg/m2/s. Analysis of the dust devil frequency distribution over the inferred dust devil zone within Gusev crater yields ???50 active dust devils/km2/sol, suggesting a dust loading into the atmosphere of ???19 kg/km2/sol. This value is less than one tenth the estimates by Cantor et al. (2001) for regional dust storms on Mars. Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research E: Planets↗
USGS2019

Investigating lake-area dynamics across a permafrost-thaw spectrum using airborne electromagnetic surveys and remote sensing time-series data in Yukon Flats, Alaska

Lakes in boreal lowlands cycle carbon and supply an important source of freshwater for wildlife and migratory waterfowl. The abundance and distribution of these lakes are supported, in part, by permafrost distribution, which is subject to change. Relationships between permafrost thaw and lake dynamics remain poorly known in most boreal regions. Here, new airborne electromagnetic (AEM) data collected during June 2010 and February 2016 were used to constrain deep permafrost distribution. AEM data were coupled with Landsat-derived lake surface-area data from 1979 through 2011 to inform temporal lake behavior changes in the 35 500- km 2 Yukon Flats ecoregion of Alaska. Together, over 1500 km of AEM data, and roughly 30 years of Landsat data were used to explore processes that drive lake dynamics across a variety of permafrost thaw states not possible in studies conducted with satellite imagery or field measurements alone. Clustered time-series data identified lakes with similar temporal dynamics. Clusters possessed similarities in lake permanence (i.e. ephemeral versus perennial), subsurface permafrost distribution, and proximity to rivers and streams. Of the clustered lakes, ~66% are inferred to have at least intermittent connectivity with other surface-water features, ~19% are inferred to have shallow subsurface connectivity to other surface water features that served as a low-pass filter for hydroclimatic fluctuations, and ~15% appear to be isolated by surrounding permafrost (i.e. no connectivity). Integrated analysis of AEM and Landsat data reveals a progression from relatively synchronous lake dynamics among disconnected lakes in the most spatially continuous, thick permafrost to quite high spatiotemporal heterogeneity in lake behavior among variably-connected lakes in regions with notably less continuous permafrost. Variability can be explained by the preferential development of thawed permeable gravel pathways for lateral water redistribution in this area. The general spatial progression in permafrost thaw state and lake area behavior may be extended to the temporal dimension. However, extensive permafrost thaw, beyond what is currently observed, is expected to promote ubiquitous subsurface connectivity, eventually evolving to a state of increased lake synchronicity.

Alaska↗
USGS2023

Uncertainties in intensity-based earthquake magnitude estimates

Estimating the magnitude of historical earthquakes is crucial for assessing seismic hazard. Magnitudes of early‐instrumental earthquakes can be inferred using a combination of instrumental records, field observations, and the observed distribution of shaking intensity determined from macroseismic observations. For earthquakes before 1900, shaking intensity distributions often provide the only information to constrain earthquake magnitude. Considerable effort has been made to develop methods to estimate the magnitude of moderate‐to‐large historical earthquakes using shaking intensities derived from macroseismic data. In this study, we consider earthquakes in California with known instrumental magnitudes to explore uncertainties in estimating the magnitude of historical earthquakes from intensity information alone. We use three California‐specific intensity prediction equations (IPEs) and an IPE based on a global ground‐motion model (GMM) to determine optimum intensity‐based magnitudes for 33 moderate‐to‐large California earthquakes between 1979 and 2021. Intensity‐based magnitudes are close to instrumental magnitudes on average. However, intensity‐based magnitudes for individual events differ by as much as 2.2 magnitude units from instrumental magnitudes. This result reflects the weak dependence of ground motions and shaking intensities on moment magnitude and their strong dependence on stress drop. Considering the intensity distributions of the 1906 San Francisco and 1989 Loma Prieta earthquakes, we show that information that could constrain rupture length is discarded when considering only the 2D decay of intensity with distance. We also show that ground‐motion intensity conversion equations used in a GMM‐based approach may cause a systematic overestimation of large historical earthquake magnitudes. This study underscores both the reducible and potentially irreducible uncertainties associated with using intensity data to estimate magnitudes of historical earthquakes using IPEs and highlights the value of using additional information to constrain rupture dimensions. Using intensity observations alone, moment magnitude uncertainties are typically on the order of a full unit.

California↗
USGS2020

Exploring the exceptional performance of a deep learning stream temperature model and the value of streamflow data

Stream water temperature ( T s ) is a variable of critical importance for aquatic ecosystem health. T s is strongly affected by groundwater-surface water interactions which can be learned from streamflow records, but previously such information was challenging to effectively absorb with process-based models due to parameter equifinality. Based on the long short-term memory (LSTM) deep learning architecture, we developed a basin-centric lumped daily mean T s model, which was trained over 118 data-rich basins with no major dams in the conterminous United States, and showed strong results. At a national scale, we obtained a median root-mean-square error of 0.69°C, Nash–Sutcliffe model efficiency coefficient of 0.985, and correlation of 0.994, which are marked improvements over previous values reported in literature. The addition of streamflow observations as a model input strongly elevated the performance of this model. In the absence of measured streamflow, we showed that a two-stage model could be used, where simulated streamflow from a pre-trained LSTM model ( Q sim ) still benefited the T s model even though no new information was brought directly into the inputs of the T s model. The model indirectly used information learned from streamflow observations provided during the training of Q sim , potentially to improve internal representation of physically meaningful variables. Our results indicate that strong relationships exist between basin-averaged forcing variables, catchment attributes, and T s that can be simulated by a single model trained by data on the continental scale.

Environmental Research Letters↗
USGS2014

Preslip and cascade processes initiating laboratory stick slip

Recent modeling studies have explored whether earthquakes begin with a large aseismic nucleation process or initiate dynamically from the rapid growth of a smaller instability in a “cascade-up” process. To explore such a case in the laboratory, we study the initiation of dynamic rupture (stick slip) of a smooth saw-cut fault in a 76mm diameter cylindrical granite laboratory sample at 40–120MPa confining pressure. We use a high dynamic range recording system to directly compare the seismic waves radiated during the stick-slip event to those radiated from tiny (M _6) discrete seismic events, commonly known as acoustic emissions (AEs), that occur in the seconds prior to each large stick slip. The seismic moments, focal mechanisms, locations, and timing of the AEs all contribute to our understanding of their mechanics and provide us with information about the stick-slip nucleation process. In a sequence of 10 stick slips, the first few microseconds of the signals recorded from stick-slip instabilities are nearly indistinguishable from those of premonitory AEs. In this sense, it appears that each stick slip begins as an AE event that rapidly (~20 μs) grows about 2 orders of magnitude in linear dimension and ruptures the entire 150mm length of the simulated fault. We also measure accelerating fault slip in the final seconds before stick slip. We estimate that this slip is at least 98% aseismic and that it both weakens the fault and produces AEs that will eventually cascade-up to initiate the larger dynamic rupture.

Journal of Geophysical Research B: Solid Earth↗
USGS2017

Effects of food web changes on Mysis diluviana diet in Lake Ontario

Mysids are important benthic-pelagic omnivores in many deep-lake food webs, yet quantitative data on their diet are limited. We explored the trophic role of Mysis diluviana in offshore Lake Ontario using samples collected in May, July, and September 2013 with a focus on seasonal and ontogenetic patterns in herbivory and zooplanktivory using two approaches. We hypothesized that Mysis diet in 2013 differs from the last investigation in 1995 in response to changes in pelagic prey over 1995 to 2013. Gut fluorescence indicated high grazing by adult and juvenile Mysis in May 2013. In July, smaller mysids were more herbivorous than larger individuals, a pattern that was less pronounced in September. Microscopic gut analysis showed copepods, including Limnocalanus , were common in diets of both size groups in May. In July, mainly cladocerans were consumed, including Cercopagis pengoi which represents a change from a past investigation that preceded Cercopagis invasion in the lake. Our results are consistent with earlier observations of a larger proportion of algae in mysid diets in spring, transitioning to relatively more zooplanktivory and use of cladocerans in the summer and fall. Higher chlorophyll content in small mysids in July than in September may be associated with the presence of a deep chlorophyll layer in July that had largely dissipated by September. Overall, Mysis in Lake Ontario continues to be a generalist omnivore, incorporating new prey items and exhibiting higher herbivory in spring.

Lake Ontario↗
USGS2019

Exploring trends in wet-season precipitation and drought indices in wet, humid and dry regions

This study examines wet season droughts using eight products from the FROGS database. The study begins by evaluating wet season precipitation totals and wet day counts at seasonal and decadal time scales. While we find a high level of agreement among the products at a seasonal timescale, evaluations of 10-year variability indicate substantial non-stationary inter-product differences that make the assessment of low-frequency changes difficult, especially in data-sparse regions. Some products, however, appear more reliable than others on decadal time scales. Global time series of dry, middle, and wet region standardized precipitation index (SPI) time series indicate little coherent change. There is substantial coherence in year-to-year variations in these time series, for the better performing products, likely indicative of skill for monitoring variations at large spatial scales. During the wet season, the data do not appear to indicate wide spread global increases in precipitation, RefET or Standardized Precipitation Evapotranspiration Index (SPEI) values. Neither the precipitation, RefET or SPEI indicate a wide-spread regional shift towards more arid conditions. Focusing on SPEI values for dry regions during droughts, however, indicate substantial increases in dry region aridity when wet season precipitation is below normal. Dry region SPEI values during droughts have decreased by -0.2 since the 1990s. More detailed analysis in further studies will be needed to confirm this result. For wet regions, however, the majority of products appear to indicate increases in wet season precipitation, although many products perform poorly in these regions due to limited observation networks, and estimated increases vary substantially.

Environmental Research Letters↗
USGS2019

Hydrodynamic and morphologic response of a back-barrier estuary to an extratropical storm

We investigated the hydrodynamic and morphologic response of Barnegat Bay-Little Egg Harbor, New Jersey, USA to Hurricane Sandy. We implemented a three-dimensional, coupled ocean-wave-sediment transport model of the estuary and explored the role of offshore water levels, offshore waves, local winds and waves by systematically removing forcings from a series of simulations. Offshore water levels had the largest impact on water levels in the bay while waves and local wind forcing created substantial spatial variation along the longitudinal axis of the bay. The shape of the bay and its orientation relative to the storm track influenced the response to winds and restricted the maximum water levels in the northern bay and reduced the maximum volume of surge. Basin-average hydrodynamic residence time was reduced by 40%, though its typical spatial distribution remained during the storm. Wave and current-induced bed shear stress resuspended fine sediment resulting in net erosion from the shoals with ensuing net deposition over fringing low-lying land. The net sediment exchange between the bay and the ocean was several times smaller than the exchange at the peak of the storm resulting in negligible net change in the bay volume. Overall, our results suggest that water level responses are highly sensitive to the specific orientation of storm winds relative to the estuary, thereby limiting the utility of simple inundation models. The sediment transport patterns indicate that storms are an important mechanism for redistributing sediment from shoals to fringing wetlands, while net change to sediment budget can be negligible.

New Jersey↗
USGS2009

Using nitrate dual isotopic composition (δ15N and δ18O) as a tool for exploring sources and cycling of nitrate in an estuarine system: Elkhorn Slough, California

Nitrate (NO 3 − ) concentrations and dual isotopic composition ( δ 15 N and δ 18 O) were measured during various seasons and tidal conditions in Elkhorn Slough to evaluate mixing of sources of NO 3 − within this California estuary. We found the isotopic composition of NO 3 − was influenced most heavily by mixing of two primary sources with unique isotopic signatures, a marine (Monterey Bay) and terrestrial agricultural runoff source (Old Salinas River). However, our attempt to use a simple two end‐member mixing model to calculate the relative contribution of these two NO 3 − sources to the Slough was complicated by periods of nonconservative behavior and/or the presence of additional sources, particularly during the dry season when NO 3 − concentrations were low. Although multiple linear regression generally yielded good fits to the observed data, deviations from conservative mixing were still evident. After consideration of potential alternative sources, we concluded that deviations from two end‐member mixing were most likely derived from interactions with marsh sediments in regions of the Slough where high rates of NO 3 − uptake and nitrification result in NO 3 − with low δ 15 N and high δ 18 O values. A simple steady state dual isotope model is used to illustrate the impact of cycling processes in an estuarine setting which may play a primary role in controlling NO 3 − isotopic composition when and where cycling rates and water residence times are high. This work expands our understanding of nitrogen and oxygen isotopes as biogeochemical tools for investigating NO 3 − sources and cycling in estuaries, emphasizing the role that cycling processes may play in altering isotopic composition.

California↗
USGS2009

Stimulating a Great Lakes coastal wetland seed bank using portable cofferdams: implications for habitat rehabilitation

Coastal wetland seed banks exposed by low lake levels or through management actions fuel the reestablishment of emergent plant assemblages (i.e., wetland habitat) critical to Great Lakes aquatic biota. This project explored the effectiveness of using portable, water-filled cofferdams as a management tool to promote the natural growth of emergent vegetation from the seed bank in a Lake Erie coastal wetland. A series of dams stretching approximately 450??m was installed temporarily to isolate hydrologically a 10-ha corner of the Crane Creek wetland complex from Lake Erie. The test area was dewatered in 2004 to mimic a low-water year, and vegetation sampling characterized the wetland seed bank response at low, middle, and high elevations in areas open to and protected from bird and mammal herbivory. The nearly two-month drawdown stimulated a rapid seed-bank-driven response by 45 plant taxa. Herbivory had little effect on plant species richness, regardless of the location along an elevation gradient. Inundation contributed to the replacement of immature emergent plant species with submersed aquatic species after the dams failed and were removed prematurely. This study revealed a number of important issues that must be considered for effective long-term implementation of portable cofferdam technology to stimulate wetland seed banks, including duration of dewatering, product size, source of clean water, replacement of damaged dams, and regular maintenance. This technology is a potentially important tool in the arsenal used by resource managers seeking to rehabilitate the functions and values of Great Lakes coastal wetland habitats.

Journal of Great Lakes Research↗
USGS2017

Considerations in comparing the U.S. Geological Survey one‐year induced‐seismicity hazard models with “Did You Feel It?” and instrumental data

The recent steep increase in seismicity rates in Oklahoma, southern Kansas, and other parts of the central United States led the U.S. Geological Survey (USGS) to develop, for the first time, a probabilistic seismic hazard forecast for one year (2016) that incorporates induced seismicity. In this study, we explore a process to ground‐truth the hazard model by comparing it with two databases of observations: modified Mercalli intensity (MMI) data from the “Did You Feel It?” (DYFI) system and peak ground acceleration (PGA) values from instrumental data. Because the 2016 hazard model was heavily based on earthquake catalogs from 2014 to 2015, this initial comparison utilized observations from these years. Annualized exceedance rates were calculated with the DYFI and instrumental data for direct comparison with the model. These comparisons required assessment of the options for converting hazard model results and instrumental data from PGA to MMI for comparison with the DYFI data. In addition, to account for known differences that affect the comparisons, the instrumental PGA and DYFI data were declustered, and the hazard model was adjusted for local site conditions. With these adjustments, examples at sites with the most data show reasonable agreement in the exceedance rates. However, the comparisons were complicated by the spatial and temporal completeness of the instrumental and DYFI observations. Furthermore, most of the DYFI responses are in the MMI II–IV range, whereas the hazard model is oriented toward forecasts at higher ground‐motion intensities, usually above about MMI IV. Nevertheless, the study demonstrates some of the issues that arise in making these comparisons, thereby informing future efforts to ground‐truth and improve hazard modeling for induced‐seismicity applications.

Kansas, Oklahoma↗
USGS2022

Orbital and in-situ investigation of periodic bedrock ridges in Glen Torridon, Gale Crater, Mars

Wind has been the dominant agent of landscape modification on Mars for the past ~3 billion years. Among the diversity of features formed by aeolian abrasion on the surface of Mars are periodic bedrock ridges (PBRs), a relatively recently recognized class of erosional bedforms on Mars for which Earth analogues are rare. Gale crater, the field site for NASA’s Mars Science Laboratory Curiosity rover since it landed there in 2012, contains a diverse and extensive record of aeolian deposition and erosion. This study focuses on a series of periodic, linear bedrock ridges that occur within the Fe/Mg-smectite clay-bearing Glen Torridon region of Aeolis Mons (informally Mount Sharp). During Curiosity’s exploration of the Glen Torridon region between sols ~2300-3080, the rover drove through this field of ridges, providing the first opportunity for the in situ observation of these enigmatic erosional features. This study characterizes the Glen Torridon ridges using orbiter and rover data to determine their morphology, spatial distribution, compositional and material properties, and association with other aeolian features in the area. Based on these observations, the Glen Torridon ridges are interpreted to be consistent with an origin as wind-eroded periodic bedrock ridges carved during the most recent exhumation of Mount Sharp into the present-day mound. Although there is evidence for multidirectional winds in the Glen Torridon region based on the orientation of modern ripples, megaripples, TARs and other bedrock indicators, the consistent orientation of the Glen Torridon ridges, coupled with morphologic asymmetries within the ridges, support formation and elongation of the Glen Torridon PBRs forms parallel to a net regional northerly wind direction in and around Gale crater.

Journal of Geophysical Research E: Planets↗
USGS2007

Evidence for montmorillonite or its compositional equivalent in Columbia Hills, Mars

During its exploration of the Columbia Hills, the Mars Exploration Rover "Spirit" encountered several similar samples that are distinctly different from Martian meteorites and known Gusev crater soils, rocks, and sediments. Occurring in a variety of contexts and locations, these "Independence class" samples are rough-textured, iron-poor (equivalent FeO ??? 4 wt%), have high Al/Si ratios, and often contain unexpectedly high concentrations of one or more minor or trace elements (including Cr, Ni, Cu, Sr, and Y). Apart from accessory minerals, the major component common to these samples has a compositional profile of major and minor elements which is similar to the smectite montmorillonite, implicating this mineral, or its compositional equivalent. Infrared thermal emission spectra do not indicate the presence of crystalline smectite. One of these samples was found spatially associated with a ferric sulfate-enriched soil horizon, possibly indicating a genetic relationship between these disparate types of materials. Compared to the nearby Wishstone and Watchtower class rocks, major aqueous alteration involving mineral dissolution and mobilization with consequent depletions of certain elements is implied for this setting and may be undetectable by remote sensing from orbit because of the small scale of the occurrences and obscuration by mantling with soil and dust.

Journal of Geophysical Research E: Planets↗
USGS2020

Formation criteria for hyporheic anoxic microzones: Assessing interactions of hydraulics, nutrients and biofilms

Recent experimental studies have detected the presence of anoxic microzones in hyporheic sediments. These microzones are small‐scale anoxic pores, embedded within oxygen‐rich porous media and can act as anaerobic reaction sites producing reduction compounds such as nitrous oxide, a potent greenhouse gas. Microbes are a key control on nutrient transformation in hyporheic sediment, but their associated biomass growth is also capable of altering hydraulic flux, leading to potential bioclogging. Here, we developed one of the first computational modeling approaches that combined hydraulics and microbial conditions to explore the continuous evolution of microzones in stream sediments. The model assessed stream and sediment conditions with different hydraulic flux (0.1–1.0 m/day Darcy flux), nutrient concentrations (O 2 = 8 mg/L, OrgC = 20 mg/L, NO − 3 = 1.5–3 mg/L, and NH 3 = 0.5–1 mg/L), and biomass scenarios (with and without). The model domain is a pore network model with random sized pore‐throat radii creating heterogeneous and anisotropic flow that is representative of a natural streambed composed of medium sand with a hydraulic conductivity of 0.8 m/day. Results from 30 day simulations indicate that hyporheic microzone formation will occur and microzone distributions are not simply controlled by residence time alone, rather by the complex interactions of hydraulic flux, nutrient concentrations, and biomass, with bioclogging having strong feedbacks on both hydraulics and nutrients. Under all conditions with biomass growth, anoxic microzones were unstable, perishing a few days after formation, because bioclogging primarily occurs near the influent (downwelling) area of the hyporheic zone. In turn, this bioclogging shifts transport conditions from advection‐dominated to diffusion‐dominated transport, removing all oxic regions in the hyporheic zone. Overall, results from the modeling show that anoxic microzones are likely to form under many hyporheic zone conditions, and be dynamic through space and time as they are dependent on both hydraulic flux and nutrient transport.

Water Resources Research↗
USGS2024

Brittle regime slip partitioned damage and deformation mechanisms along the eastern Denali fault zone in southwestern, Yukon

Rare bedrock exposures of the eastern Denali fault zone in southwestern Yukon allow for the measurement, sampling, and analyses of brittle regime fault slip data and deformation mechanisms to explore relations to far field, oblique plate motions. Host rock lithologies and associated slip surfaces show episodic damage zone‐related deformation and calcite ± hematite ± chlorite related hydrothermal fluid flow. This regional scale network of asymmetric fault damage is spatially and kinematically linked to a discrete and narrow fault core. Fault network observations, orientations, slip data, and strain inversions document a slip partitioned strike‐slip fault system with locally and mutually overprinting strike‐, oblique‐, and dip‐slip components. Microstructural analyses reveal crystal plastic and co‐seismic brittle deformation mechanisms active in a narrow range of upper crustal temperature, pressure, fluid, and chemical conditions. The net damage related slip is not exclusively formed by a single kinematic system, but rather a fully partitioned, time integrated system likely operative for much of the fault's brittle regime evolution temporally constrained by previously published thermochronometric data. Although the fault slip data was collected from outcrop‐scale exposures at sites tens of kilometers apart, results show remarkable correlation between fault kinematics and plate motions along the ∼580 km long eastern Denali fault segment. End member, subhorizontal, northeast directed reverse and north directed dextral strike slip fault strain axes closely reflect relative plate motion interactions over at least the last 30 m.y. and act as a proxy for far‐field stresses compatible with the kinematics of the damage zone network.

Alaska↗
USGS2020

Mapping the 3-D extent of the Stillwater Complex, Montana—Implications for potential platinum group element exploration and development

Geophysical models characterize the exposed and interpreted buried extent of the Stillwater Complex, critical for understanding the origin of the layered mafic intrusion and its associated high-grade platinum group element resources. The 3D models, constrained by gravity, magnetic, xenolith, seismic, borehole, and rock property data indicate that the likely maximum extent of the Stillwater Complex beneath Phanerozoic cover is ~10 times greater than its outcrop, ~2240 km 2 . The thickness values are poorly constrained but vary from ~7000 to 12,000 m, depending on crustal and mantle density variations and depths to the top of the lower crust and mantle. This thickness may include dense metasedimentary units of the basin into which the Stillwater Complex intruded. Using the modeled thickness results in a volume estimate of ~24,700 km 3 , albeit poorly constrained. New analyses of xenoliths from the Cretaceous Sliderock and Suzie Peak intrusions produce ages of 2706–2716 Ma, corresponding to the age of the Stillwater Complex, and 2813 Ma, corresponding to the age of Archean gneissic basement. Seismic reflectors in inferred Archean crystalline basement, possibly including the Stillwater Complex, dip ~25–30° north, with segments dipping as much as 70° north. Layered reflectors beneath the Phanerozoic sedimentary section and above the inferred Archean crystalline basement may represent metasedimentary units, perhaps a southern extension of the Mesoproterozoic Belt Basin. The potential field models and the seismic reflection data suggest that the Stillwater Complex was dipping northward prior to deposition of Cambrian strata, perhaps uplifted in the late Archean or Proterozoic as previously proposed, and that during Laramide times, faulting and intrusions highly disrupted the complex. Temperature measurements from boreholes help constrain the depths of feasible mining of the complex.

Montana↗
USGS2006

Erosion rates at the Mars Exploration Rover landing sites and long-term climate change on Mars

Erosion rates derived from the Gusev cratered plains and the erosion of weak sulfates by saltating sand at Meridiani Planum are so slow that they argue that the present dry and desiccating environment has persisted since the Early Hesperian. In contrast, sedimentary rocks at Meridiani formed in the presence of groundwater and occasional surface water, and many Columbia Hills rocks at Gusev underwent aqueous alteration during the Late Noachian, approximately coeval with a wide variety of geomorphic indicators that indicate a wetter and likely warmer environment. Two-toned rocks, elevated ventifacts, and perched and undercut rocks indicate localized deflation of the Gusev plains and deposition of an equivalent amount of sediment into craters to form hollows, suggesting average erosion rates of ???0.03 nm/yr. Erosion of Hesperian craters, modification of Late Amazonian craters, and the concentration of hematite concretions in the soils of Meridiani yield slightly higher average erosion rates of 1-10 nm/yr in the Amazonian. These erosion rates are 2-5 orders of magnitude lower than the slowest continental denudation rates on Earth, indicating that liquid water was not an active erosional agent. Erosion rates for Meridiani just before deposition of the sulfate-rich sediments and other eroded Noachian areas are comparable with slow denudation rates on Earth that are dominated by liquid water. Available data suggest the climate change at the landing sites from wet and likely warm to dry and desiccating occurred sometime between the Late Noachian and the beginning of the Late Hesperian (3.7-3.5 Ga). Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research E: Planets↗
USGS1978

Porphyry copper exploration model for northern Sonora, Mexico

The regional tectonic pattern of the porphyry copper province of northern Sonora, Mexico, is similar to, but more complex than, the Colorado mineral belt. Four northeast-trending shear zones, spaced at 30- to 50-kilometer intervals from Hermosillo, Mexico, north to Nogales, Mexico, are interpreted from analysis of lineament data from Landsat-1 images. From comparison with other areas, these zones are believed to have Precambrian ancestry. North-northwest-trending lineaments, which commonly occur along range fronts, are interpreted to be the principal basin and range faults. North-south- and west-northwest-trending lineaments are believed to be related to lesser faults and fractures. Areas of hydrothermal alteration have been mapped on a regional scale by use of digital image-processing techniques. These altered areas occur on the flanks of the northeast-trending shear zones, generally near intersections with north-northwest-trending lineaments. This tectonic framework with associated alteration is similar to patterns in the Colorado Plateau, the porphyry copper province of Arizona, and central Mexico and is supported by relationships found in limited local mapping. The pattern of northeast-trending shear zones and north-northwest-trending faults in northern Sonora closely resembles the tectonic framework of the Colorado mineral belt. The younger sedimentary and volcanic cover, however, has not been eroded in Sonora as it has been in the Colorado mineral belt. This comparison with the Colorado mineral belt suggests that mineral deposits are most likely to be found in altered areas where a northeast-trending shear zone is intersected by north-northwest-trending faults. Seven areas with these characteristics have been defined in Sonora and are being evaluated by means of regional geophysical and geochemical data.

Sonora↗