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,153 records · Page 64Linked to original sources
USGS2008

Constraining fault constitutive behavior with slip and stress heterogeneity

We study how enforcing self-consistency in the statistical properties of the preshear and postshear stress on a fault can be used to constrain fault constitutive behavior beyond that required to produce a desired spatial and temporal evolution of slip in a single event. We explore features of rupture dynamics that (1) lead to slip heterogeneity in earthquake ruptures and (2) maintain these conditions following rupture, so that the stress field is compatible with the generation of aftershocks and facilitates heterogeneous slip in subsequent events. Our three-dimensional fmite element simulations of magnitude 7 events on a vertical, planar strike-slip fault show that the conditions that lead to slip heterogeneity remain in place after large events when the dynamic stress drop (initial shear stress) and breakdown work (fracture energy) are spatially heterogeneous. In these models the breakdown work is on the order of MJ/m2, which is comparable to the radiated energy. These conditions producing slip heterogeneity also tend to produce narrower slip pulses independent of a slip rate dependence in the fault constitutive model. An alternative mechanism for generating these confined slip pulses appears to be fault constitutive models that have a stronger rate dependence, which also makes them difficult to implement in numerical models. We hypothesize that self-consistent ruptures could also be produced by very narrow slip pulses propagating in a self-sustaining heterogeneous stress field with breakdown work comparable to fracture energy estimates of kJ/M2. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2010

Mineralogy and chemistry of cobbles at Meridiani Planum, Mars, investigated by the Mars Exploration Rover Opportunity

Numerous loose rocks with dimensions of a few centimeters to tens of centimeters and with no obvious physical relationship to outcrop rocks have been observed along the traverse of the Mars Exploration Rover Opportunity. To date, about a dozen of these rocks have been analyzed with Opportunity's contact instruments, providing information about elemental chemistry (Alpha Particle X‐ray Spectrometer), iron mineralogy and oxidation states (Mössbauer Spectrometer) and texture (Microscopic Imager). These “cobbles” appear to be impact related, and three distinct groups can be identified on the basis of chemistry and mineralogy. The first group comprises bright fragments of the sulfate‐rich bedrock that are compositionally and texturally indistinguishable from outcrop rocks. All other cobbles are dark and are divided into two groups, referred to as the “Barberton group” and the “Arkansas group,” after the first specimen of each that was encountered by Opportunity. Barberton group cobbles are interpreted as meteorites with an overall chemistry and mineralogy consistent with a mesosiderite silicate clast composition. Arkansas group cobbles appear to be related to Meridiani outcrop and contain an additional basaltic component. They have brecciated textures, pointing to an impact‐related origin during which local bedrock and basaltic material were mixed.

Journal of Geophysical Research E: Planets↗
USGS2017

A decade of induced slip on the causative fault of the 2015 Mw 4.0 Venus earthquake, northeast Johnson County, Texas

On 7 May 2015, a M w 4.0 earthquake occurred near Venus, northeast Johnson County, Texas, in an area of the Bend Arch-Fort Worth Basin that reports long-term, high-volume wastewater disposal and that has hosted felt earthquakes since 2009. In the weeks following the M w 4.0 earthquake, we deployed a local seismic network and purchased nearby active-source seismic reflection data to capture additional events, characterize the causative fault, and explore potential links between ongoing industry activity and seismicity. Hypocenter relocations of the resulting local earthquake catalog span ~4–6 km depth and indicate a fault striking ~230 ° , dipping to the west, consistent with a nodal plane of the M w 4.0 regional moment tensor. Fault plane solutions indicate normal faulting, with B axes striking parallel to maximum horizontal compressive stress. Seismic reflection data image the reactivated basement fault penetrating the Ordovician disposal layer and Mississippian production layer, but not displacing post-Lower Pennsylvanian units. Template matching at regional seismic stations indicates that low-magnitude earthquakes with similar waveforms began in April 2008, with increasing magnitude over time. Pressure data from five saltwater disposal wells within 5 km of the active fault indicate a disposal formation that is 0.9–4.8 MPa above hydrostatic. We suggest that the injection of 28,000,000 m 3 of wastewater between 2006 and 2015 at these wells led to an increase in subsurface pore fluid pressure that contributed to inducing this long-lived earthquake sequence. The 2015 M w 4.0 event represents the largest event in the continuing evolution of slip on the causative fault.

Texas↗
USGS2021

Setting and tracking suppression targets for sea lampreys in the Great Lakes

In response to invasive species, the course of action taken by management agencies often evolves over a range of options from a do-nothing approach to suppression to complete eradication. As a case study of suppression targets, we explore the history of approaches used by the Great Lakes Fishery Commission in response to the invasion of the Laurentian Great Lakes by sea lampreys ( Petromyzon marinus ). With the early realization that eradication was not possible using available techniques, focus shifted to suppression of sea lampreys to support fish community objectives for the lakes. As a surrogate for damage to the fishery, a suppression target was defined for the maximum acceptable marking rate, indicated by the average number of sea lamprey wounds observed on every 100 large lake trout ( Salvelinus namaycush ) surveyed. Historic marking rates were used to delineate time periods of acceptable sea lamprey predation levels as an approximate measure of sea lamprey-induced mortality. A second target, independent of lake trout population status, was based on the average sea lamprey abundance estimated during the specified time periods. These intuitive targets have served as suppression benchmarks for the Sea Lamprey Control Program in the Great Lakes, allowing decision makers to gauge progress towards targets, refine control strategies and prioritize geographic areas for increased control effort. Here we document the development and subsequent changes to targets, summarize the methods used to implement these changes, and provide considerations for the future.

Great Lakes↗
USGS2021

Geodetic constraints on a 25-year magmatic inflation episode near Three Sisters, central Oregon

Crustal inflation near the Three Sisters volcanic center documented since the mid-1990s has persisted for more than two decades. We update past analyses of the event through 2020 by simultaneously inverting InSAR interferograms, GPS time series, and leveling data for time-dependent volcanic deformation source parameters. We explore several source models to estimate how the deformation rate varied through time and to identify parameters that can reproduce measured deformation. Our preferred model is a Mogi source 4.1 km below sea level (5.9 km below the surface) about 5 km west of the summit of South Sister. Inflation started in late 1995 or 1996; the rate increased rapidly during 1998–1999, and peaked in late 1999, resulting in maximum surface uplift of about 30 cm by mid-2020. Since 2000, the inflation rate generally declined exponentially with a time constant of about 6 years. Two source inflation scenarios fit the data equally well. In the first, the crust surrounding the source is elastic and the net source-volume increase, which we attribute to persistent magma input, has been about 49 × 10 6 m 3 . The second scenario adds a viscoelastic shell surrounding the Mogi source. In that case, an injection of about 21 × 10 6 m 3 of magma prior to 2000, followed by continuing relaxation of the viscoelastic shell, can account for most of the observed surface deformation. In both scenarios, modeling reveals quasiperiodic increases in the inflation rate (pulses) with a recurrence interval of 3–4 years, both before and after 2000.

Oregon↗
USGS2023

Diversification of forest management can mitigate wind damage risk and maintain biodiversity

Mitigating future forest risks, safeguarding timber revenues and improving biodiversity are key considerations for current boreal forest management. Alternatives to rotation forestry likely have an important role, but how they will perform under a changing climate remains unclear. We used a boreal forest growth simulator to explore how variations on traditional clear-cutting, in rotation length, thinning intensity, and increasing number of remaining trees after final harvest (green tree retention), and on extent of continuous cover forestry will affect stand-level probability of wind damage, timber production, deadwood volume, and habitats for forest species. We used business-as-usual rotation forestry as a baseline and compared alternative management adaptations under the reference and two climate change scenarios. Climate change increased overall timber production and had lower impacts on biodiversity compared to management adaptations. Shortening the rotation length reduced the probability of wind damage compared to business-as-usual, but also decreased both deadwood volume and suitable habitats for our focal species. Continuous cover forestry, and management with refraining from thinnings, and extension of rotation length represent complementary approaches benefiting biodiversity, with respective effects of improving timber revenues, reducing wind damage risk, and benefiting old-growth forest structures. However, extensive application of rotation length shortening to mitigate wind damage risk may be detrimental for forest biodiversity. To safeguard forest biodiversity over the landscape, shortening of the rotation length could be complemented with widespread application of regimes promoting old-growth forest structures.

European Journal of Forest Research↗
USGS2015

Shale-gas assessment: Comparison of gas-in-place versus performance-based approaches

The recent interest in exploration for shale gas increases the demand for a reliable, compatible resource assessment. Many different assessment methods are used, commonly depending on types and quantity of data available, which may lead to significantly divergent results for the same shale-gas play. This study compares results obtained using performance-based and gas-in-place methodologies to assess a well-developed and active shale-gas play (Woodford Shale, Arkoma Basin, USA) and two untested, hypothetical shale-gas plays (Shublik and Brookian, Alaska North Slope, USA). Results show that the two assessment methods produce comparable results when assessment units are identically defined and similar geological constraints are used as input parameters. Inherent uncertainties are associated with both assessment methods, and these are related to aspects of shale-gas production that are not well understood. The performance-based method relies on decline trend analysis to generate distributions of estimated ultimate recovery (EUR), and uncertainty increases in cases of short production history. The gas-in-place method requires the application of a recovery factor to estimate technically recoverable resources, and both absolute values of recovery factors and their spatial variability are poorly documented, and therefore a source of uncertainty.

Natural Resources Research↗
USGS2011

Sea otter mortality in fish and shellfish traps: Estimating potential impacts and exploring possible solutions

Sea otters Enhydra lutris can be bycaught and drowned in fishing pots and traps, which may pose a threat to the welfare of otter populations. We explored this potential problem and its solutions using a wide variety of analyses. We exposed live California (USA) sea otters to finfish traps, lobster traps, and mock Dungeness crab traps in captive trials and found that the animals attempted to enter the circular and rectangular fyke openings, with some becoming entrapped. Using both live and dead sea otters, we found that a 3 × 9 inch (7.6 × 22.9 cm) fyke opening (1 inch narrower than the 4 × 9 inch [10.2 × 22.9 cm] openings currently used in California’s commercial Dungeness crab fishery) would exclude most free-living (i.e. weaned from their mothers) otters while permitting the undiminished capture of crabs. Observer programs do not currently exist in California for these fisheries, so we calculated the effort required by an observer program to document sea otter bycatch over a range of hypothetical levels and evaluated the impact of those mortality rates on population growth. These analyses demonstrate that significant mortality from bycatch might easily go undetected, even with seemingly high levels of observer effort. As sea otters reoccupy portions of their former habitat in California, co-occurrence with finfish and shellfish traps with relatively large fyke openings will increase.

Endangered Species Research↗
USGS2020

Assessing the influence of natural copper-nickel-bearing bedrocks of the Duluth Complex on water quality in Minnesota, 2013–15

The U.S. Geological Survey, in cooperation with the University of Minnesota-Duluth Natural Resources Research Institute, completed an assessment of regional water quality in areas of potential base-metal mining in Minnesota. Bedrock, soil, streambed sediment, and surface-water samples were collected in three watersheds that cross the basal part of the Duluth Complex with different mineral-deposit settings: (1) copper-nickel-platinum group element mineralization (Filson Creek), (2) iron-titanium-oxide mineralization (headwaters of the St. Louis River), and (3) no identified mineralization (Keeley Creek). At least 10 bedrock, 30 soil (2 each from 15 sites), and as many as 13 streambed sediment samples were collected in each watershed and analyzed for 44 major and trace elements, total and inorganic carbon, and 10 loosely bound metals (when possible). Surface-water samples were collected at four to nine locations in each watershed three to four times per year for 2 years (total of 141 environmental samples). Surface-water samples were analyzed for 10 trace metals (total and dissolved concentrations), 8 trace elements, 8 major ions (dissolved concentrations), alkalinity, and total and dissolved organic carbon. Metal and element concentrations in solid media varied by watershed, representing local geology. Copper-nickel sulfide mineralization in the Filson Creek watershed was evidenced in bedrock, soil, and streambed sediments. In the Keeley Creek watershed, silicate mineralogy of underlying bedrock contributed metals to streambed sediments. Thick glacial cover masked potential bedrock contributions to solid media in the St. Louis River watershed. Water-quality data indicate that waters in all three watersheds are dilute. Water quality is more similar between the Filson and Keeley Creek watersheds, compared to the St. Louis River watershed, because of the difference in glacial cover. Metal concentrations (copper and nickel, in particular) in surface-water samples follow similar patterns of concentrations in solid media, indicating the influence of bedrock on water quality in Filson and Keeley Creeks. Data from this study provide a baseline of metal concentrations and general water quality within an area of active mineral exploration.

Minnesota↗
USGS2023

Stream temperature prediction in a shifting environment: The influence of deep learning architecture

Stream temperature is a fundamental control on ecosystem health. Recent efforts incorporating process guidance into deep learning models for predicting stream temperature have been shown to outperform existing statistical and physical models. This performance is in part because deep learning architectures can actively learn spatiotemporal relationships that govern how water and energy propagate through a river network. However, exploration of how spatiotemporal awareness and process guidance influence a model's generalizability under shifting environmental conditions such as climate change is limited. Here, we use Explainable Artificial Intelligence (XAI) to interrogate how differing deep learning architectures affect a model's learned spatial and temporal dependencies, and how those learned dependencies affect a model's ability to maintain high accuracy when applied to unseen environmental conditions. Using the Delaware River Basin in the northeastern United States as a test case, we compare two spatiotemporally aware process-guided deep learning models for predicting stream temperature (a recurrent graph convolution network—RGCN, and a temporal convolution graph model—Graph WaveNet). Both models achieve equally high predictive performance when testing data are well represented in the training data (test root mean squared errors of 1.64°C and 1.65°C); however, Graph WaveNet significantly outperforms RGCN in 4 out of 5 experiments where test partitions represent different types of unseen environmental conditions. XAI results show that the architecture of Graph WaveNet leads to learned spatial relationships with greater fidelity to physical processes, and that this fidelity improves the generalizability of the model when applied to shifting and/or unseen environmental conditions.

Delaware River Basin↗
USGS2017

Delayed seismicity rate changes controlled by static stress transfer

On 15 June 2010, a M w 5.7 earthquake occurred near Ocotillo, California, in the Yuha Desert. This event was the largest aftershock of the 4 April 2010 M w 7.2 El Mayor-Cucapah (EMC) earthquake in this region. The EMC mainshock and subsequent Ocotillo aftershock provide an opportunity to test the Coulomb failure hypothesis (CFS). We explore the spatiotemporal correlation between seismicity rate changes and regions of positive and negative CFS change imparted by the Ocotillo event. Based on simple CFS calculations we divide the Yuha Desert into three subregions, one triggering zone and two stress shadow zones. We find the nominal triggering zone displays immediate triggering, one stress shadowed region experiences immediate quiescence, and the other nominal stress shadow undergoes an immediate rate increase followed by a delayed shutdown. We quantitatively model the spatiotemporal variation of earthquake rates by combining calculations of CFS change with the rate-state earthquake rate formulation of Dieterich (1994), assuming that each subregion contains a mixture of nucleation sources that experienced a CFS change of differing signs. Our modeling reproduces the observations, including the observed delay in the stress shadow effect in the third region following the Ocotillo aftershock. The delayed shadow effect occurs because of intrinsic differences in the amplitude of the rate response to positive and negative stress changes and the time constants for return to background rates for the two populations. We find that rate-state models of time-dependent earthquake rates are in good agreement with the observed rates and thus explain the complex spatiotemporal patterns of seismicity.

California↗
USGS2020

Submarine landslide kinematics derived from high-resolution imaging in Port Valdez, Alaska

Submarine landslides caused by strong ground shaking during the M9.2 1964 Great Alaska earthquake generated a tsunami that destroyed much of the old town of Valdez, Alaska, and was responsible for 32 deaths at that location. We explore structural details of the 1964 landslide deposit, as well as landslide deposits from earlier events, in order to characterize kinematics of the landslide process. We present a new high‐resolution seismic reflection data set that images the 1964 landslide deposit and six pre‐1964 deposits with great detail. These deposits are represented by thick packages (~7–23 m) of debris within >500 m of fjord sedimentation above basement. Internal slide structures are associated with distinctive landslide failure mechanisms, including detailed erosional and depositional features and structures resolved within both landslide blocks and distal debris flow layers. Based on comparisons of deposit volume from subbottom structure and differenced bathymetry, we refine prior interpretations of the source of failed material. New data show evidence for basal erosion and reworking of fjord‐floor sedimentation. Additionally, material comprising the 1964 landslide appears to have been translated and deformed by lateral thrusting, rather than having been sourced entirely from upslope evacuation zones. Taking into account these complexities in depositional patterns, we show variations in slide size through Holocene time and relate the history of landslides to the paleoseismic record. Collectively, these new observations demonstrate that Port Valdez has a repeated history of large submarine landslides, which are likely associated with large megathrust earthquakes.

Alaska↗
USGS2021

Warmer winters increase the biomass of phytoplankton in a large floodplain river

Winters are changing rapidly across the globe but the implications for aquatic productivity and food webs are not well understood. In addition, the degree to which winter dynamics in aquatic systems respond to large-scale climate versus ecosystem-level factors is unclear but important for understanding and managing potential changes. We used a unique winter data set from the Upper Mississippi River System to explore spatial and temporal patterns in phytoplankton biomass (chlorophyll a , CHL) and associated environmental covariates across 25 years and ∼1,500 river km. To assess the role of regional climate versus site-specific drivers of winter CHL, we evaluated whether there were coherent long-term CHL dynamics from north to south and across lotic-lentic areas. We then estimated the degree to which these patterns were associated with climate variability (i.e., the Multivariate El Nino-Southern Oscillation Index), winter severity (freezing degree days), river discharge, or site-specific environmental variables (ice depth, snow depth, and nutrient concentrations). We found that winter CHL was typically highest in ice-free reaches and backwater lakes, occasionally exceeding summer values. We did not find highly synchronous CHL dynamics across the basin, but instead show that temporal trends were independent among river reaches and lotic-lentic areas of the river. Moreover, after accounting for these spatial dynamics, we found that CHL was most responsive to winter air temperature, being consistently higher in years with warmer winters across the basin. These results indicate that although productivity dynamics are highly dynamic within large river ecosystems, changes in the duration and severity of winter may uniformly increase wintertime productivity.

Minnesota, Missouri, Illinois, Iowa, Wisconsin↗
USGS2017

Exploring potential effects of cormorant predation on the fish community in Saginaw Bay, Lake Huron

Stakeholders and fishery managers expressed concern that double-crested cormorant Phalacrocorax auritus predation may be a factor in the recent poor survival of yellow perch Perca flavescens in Saginaw Bay. We quantified cormorant diets from two nesting colonies in Saginaw Bay during April–September in 2013 and 2014, with special emphasis on impacts to yellow perch. Cormorants ( n = 691) were collected when returning to colonies after foraging. Stomachs were removed and preserved in the field. Diet items were identified, enumerated, and measured ( n = 23.373). Cormorant diets from Saginaw Bay indicate a heavy reliance on round goby and Notropis species as prey during the breeding season, consistent with other areas of the Great Lakes where round goby and cormorants coincide. Respectively, the three most common prey species observed by number (%) and biomass (%) pooled across years and sites were round goby Neogobius melanostomus (56.6%, 42.1%), emerald shiner Notropis antherinoides (25.2%, 12.5%), and yellow perch (8.0%, 14.1%). Diet composition was more variable at Spoils Island than at Little Charity Island. Overall cormorant consumption (estimated using cormorant consumption demand rates) of yellow perch was compared to walleye consumption. Cormorant consumption of age-1 yellow perch was 13–17% as much as mean walleye consumption of yellow perch in 2013 and 8–11% in 2014. The cumulative effects of walleye and spring cormorant predation likely represent a recruitment bottleneck for yellow perch in Saginaw Bay. Future studies determining age-specific abundance of yellow perch would facilitate better determination of cormorant predation significance.

Lake Huron, Saginaw Bay↗
USGS2016

Inhibition of Akt enhances the chemopreventive effects of topical rapamycin in mouse skin

The PI3Kinase/Akt/mTOR pathway has important roles in cancer development for multiple tumor types, including UV-induced non-melanoma skin cancer. Immunosuppressed populations are at increased risk of aggressive cutaneous squamous cell carcinoma (SCC). Individuals who are treated with rapamycin, (sirolimus, a classical mTOR inhibitor) have significantly decreased rates of developing new cutaneous SCCs compared to those that receive traditional immunosuppression. However, systemic rapamycin use can lead to significant adverse events. Here we explored the use of topical rapamycin as a chemopreventive agent in the context of solar simulated light (SSL)-induced skin carcinogenesis. In SKH-1 mice, topical rapamycin treatment decreased tumor yields when applied after completion of 15 weeks of SSL exposure compared to controls. However, applying rapamycin during SSL exposure for 15 weeks, and continuing for 10 weeks after UV treatment, increased tumor yields. We also examined whether a combinatorial approach might result in more significant tumor suppression by rapamycin. We validated that rapamycin causes increased Akt (S473) phosphorylation in the epidermis after SSL, and show for the first time that this dysregulation can be inhibited in vivo by a selective PDK1/Akt inhibitor, PHT-427. Combining rapamycin with PHT-427 on tumor prone skin additively caused a significant reduction of tumor multiplicity compared to vehicle controls. Our findings indicate that patients taking rapamycin should avoid sun exposure, and that combining topical mTOR inhibitors and Akt inhibitors may be a viable chemoprevention option for individuals at high risk for cutaneous SCC.

Cancer Prevention Research↗
USGS2020

Resolving small-scale forest snow patterns using an energy-balance snow model with a 1-layer canopy

Modelling spatiotemporal dynamics of snow in forests is challenging, as involved processes are strongly dependent on small-scale canopy properties. In this study, we explore how local canopy structure information can be integrated in a medium-complexity energy-balance snow model to replicate observed snow patterns at very high spatial resolutions. Snow depth distributions simulated with the Flexible Snow Model (FSM2) were tested against extensive experimental data acquired in discontinuous subalpine forest stands in Eastern Switzerland over three winters. While the default canopy implementation in FSM2 fails to capture the observed snow depth variability, performance is considerably improved when local canopy cover fraction and hemispherical sky view fraction are additionally accounted for (30% reduction in RMSE). However, realistic snow depth distribution patterns throughout the season are only achieved if effective temperatures of near and distant canopy elements are discerned, and if a mechanism to mimic preferential deposition of snow in canopy gaps is included. We demonstrate that by diversifying the canopy structure input in order to reflect respective portions of the canopy relevant to different processes, even a simple model based on widely used process parametrizations and canopy metrics can be applied for high-resolution simulations of the sub-canopy snow cover with just a few modifications. The presented approaches could be implemented in commonly used land surface models, allowing upscaling experiments and development of sub-grid parametrizations without necessitating complex high-resolution models.

Water Resources Research↗
USGS1998

Two-way coupling between Vesuvius eruptions and southern Apennine earthquakes, Italy, by elastic stress transfer

During the past 1000 years, eruptions of Vesuvius have often been accompanied by large earthquakes in the Apennines 50–60 km to the northeast. Statistical investigations had shown that earthquakes often preceded eruptions, typically by less than a decade, but did not provide a physical explanation for the correlation. Here, we explore elastic stress interaction between earthquakes and eruptions under the hypothesis that small stress changes can promote events when the Apennine normal faults and the Vesuvius magma body are close to failure. We show that earthquakes can promote eruptions by compressing the magma body at depth and opening suitably oriented near-surface conduits. Voiding the magma body in turns brings these same normal faults closer to Coulomb failure, promoting earthquakes. Such a coupling is strongest if the magma reservoir is a dike oriented normal to the regional extension axis, parallel to the Apennines, and the near-surface conduits and fissures are oriented normal to the Apennines. This preferred orientation suggests that the eruptions issuing from such fissures should be most closely linked in time to Apennine earthquakes. Large Apennine earthquakes since 1400 are calculated to have transferred more stress to Vesuvius than all but the largest eruptions have transferred to Apennine faults, which may explain why earthquakes more commonly lead than follow eruptions. A two-way coupling may thus link earthquakes and Vesuvius eruptions along a 100-km-long set of faults. We test the statistical significance of the earthquake-eruption correlation in the two-way coupling zone, and find a correlation significant at the 95% confidence level.

Journal of Geophysical Research B: Solid Earth↗
USGS1976

Exploration for porphyry copper deposits in Pakistan using digital processing of Landsat-1 data

Rock-type classification by digital-computer processing of Landsat-1 (formerly ERTS-1) multispectral scanner data has been used to select 23 prospecting targets in the Chagai District, Pakistan, 5 of which have proved to be large areas of hydrothermally altered porphyry containing pyrite. The known porphyry copper deposit at Saindak was used as a control area. Empirical maximum and minimum apparent reflectance limits were selected for each multispectral scanner band in each rock type classified, and a relatively unrefined classification table was prepared. Where the values for all four bands fitted within the limits designated for a particular class, a symbol for the presumed rock type was printed by the computer at the appropriate location; these symbols formed a classification map. Drainage channels, areas of mineralized quartz diorite, areas of pyrite-rich rock, and the approximate limit of propylitic alteration were very well delineated on the computer-generated map of the test area. The classification method was then used to evaluate 2,100 km 2 in the Mashki Chah region, east of the test area, and a partial check of the results was made in the field. The results of the experiment show that outcrops of hydrothermally altered and mineralized rock can be identified from Landsat-1 data under favorable conditions. The empirical method of digital computer classification of the multispectral scanner data was relatively unrefined and rapid. The five mineralized prospecting sites identified are in locations that I would not have selected as favorable on the basis of geologic knowledge or photogeologic interpretation at the time of the investigation.

Journal of Research of the U.S. Geological Survey↗