Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Cells”

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,801 records · Page 100Linked to original sources

Simulation of ground water flow in the Glaciofluvial, Saginaw, Parma-Bayport, and Marshall Aquifers, Central Lower Peninsula of Michigan

A steady-state finite difference model was developed to simulate ground-water flow in four regional aquifers in Michigan’s Lower Peninsula. The Glaciofluvial, Saginaw, Parma-Bayport, and Marshall aquifers were simulated as layers 1 through 4, respectively, in the model. Separately calculated vertical conductances input to the model were used to simulate the intervening Till/“Red Beds”, Saginaw, and Michigan confining units, respectively. The model domain was laterally bound by a continuous specifiedhead boundary, formed from Lakes Michigan, Huron, St. Clair, and Erie, together with the St. Clair and Detroit River connecting channels. The model was developed to quantify regional ground-water flow in the aquifer systems using independently determined recharge estimates. The flow model showed that groundwater heads and flows in the Glaciofluvial aquifer are controlled by local stream stages and discharges, resulting in localized flow cells accounting for 95-percent of the overall model water budget. Simulation of recharge to an unspecified water table also enabled the estimation of ground-water discharge to three Great Lakes. A computer diskette contains all MODFLOW and MODFLOWP input files, as well as digital model surfaces and several Fortran processing routines used to construct the surfaces. The diskette also provides the data used for calibration and sensitivity analysis.

Michigan↗

Methods to achieve accurate projection of regional and global raster databases

Modeling regional and global activities of climatic and human-induced change requires accurate geographic data from which we can develop mathematical and statistical tabulations of attributes and properties of the environment. Many of these models depend on data formatted as raster cells or matrices of pixel values. Recently, it has been demonstrated that regional and global raster datasets are subject to significant error from mathematical projection and that these errors are of such magnitude that model results may be jeopardized (Steinwand, et al., 1995; Yang, et al., 1996; Usery and Seong, 2001; Seong and Usery, 2001). There is a need to develop methods of projection that maintain the accuracy of these datasets to support regional and global analyses and modeling

Open-File Report↗

Digital data used to relate nutrient input to water quality in the Chesapeake Bay watershed

Digital data sets compiled by the U.S. Geological Survey were used as input for a collection of Spatially Referenced Regressions On Watershed (SPARROW) attributes for the Chesapeake Bay region including parts of Delaware, Maryland, New York, Pennsylvania, Virginia, West Virginia, and the District of Columbia. These regressions use a nonlinear statistical approach to relate nutrient sources and land-surface characteristics to nutrient loads of streams throughout the Chesapeake Bay watershed. A digital segmented-watershed network serves as the primary framework for spatially referencing nutrient-source and land-surface characteristic data within a geographic information system. Flow direction and flow accumulation generated from a 30-meter cell-size Digital Elevation Model and attributes from 1:500,000-scale stream data were used to generate stream and watershed networks. Spatial data sets representing nutrient inputs of total nitrogen and total phosphorus from the early 1990's were created and compiled from numerous sources. Data include atmospheric deposition, septic systems, point-source locations, land use, land cover, and agricultural sources such as commercial fertilizer and manure. Some land-surface characteristic data sets representing factors that affect the transport of nutrients also were compiled. Data sets include land use, land cover, average-annual precipitation and temperature, slope, hydrogeomorphic regions, and soil permeability. Nutrient-input and land-surface characteristic data sets merged with the segmented-watershed network provide the spatial detail by watershed segment required by SPARROW. Stream-nutrient load estimates for 132 sampling sites representing the early 1990's (103 for total nitrogen and 121 for total phosphorus) serve as the dependent variables for the regressions. These estimates were used to calibrate models of total nitrogen and total phosphorus depicting 1992 land-surface conditions. Examples of model predictions consist of stream-nutrient load and source percentages contributed locally to each stream reach, as well as percentages of the load that reach Chesapeake Bay.

Open-File Report↗

Preliminary mineralogical characterization of weathered and less-weathered strata of the Meade Peak phosphatic shale member of the Permian Phosphoria Formation: Measured sections C and D, Dry Valley, Caribou County, Idaho

This study used powder X-ray diffraction, with Rietveld quantification software to quantify and characterize the mineralogy of the 83 samples collected from two stratigraphic sections measured by the U.S. Geological Survey at the Dry Valley mine in the Meade Peak Member of the Phosphoria Formation. Analyses show extensive variability of carbonate substitution in the fluorapatite structure, determined by measuring the apatite a-cell dimension, as well as patterns of correlation between mineralogy and the stratigraphy.

Idaho↗

User's guide to Model Viewer, a program for three-dimensional visualization of ground-water model results

Model Viewer is a computer program that displays the results of three-dimensional groundwater models. Scalar data (such as hydraulic head or solute concentration) may be displayed as a solid or a set of isosurfaces, using a red-to-blue color spectrum to represent a range of scalar values. Vector data (such as velocity or specific discharge) are represented by lines oriented to the vector direction and scaled to the vector magnitude. Model Viewer can also display pathlines, cells or nodes that represent model features such as streams and wells, and auxiliary graphic objects such as grid lines and coordinate axes. Users may crop the model grid in different orientations to examine the interior structure of the data. For transient simulations, Model Viewer can animate the time evolution of the simulated quantities. The current version (1.0) of Model Viewer runs on Microsoft Windows 95, 98, NT and 2000 operating systems, and supports the following models: MODFLOW-2000, MODFLOW-2000 with the Ground-Water Transport Process, MODFLOW-96, MOC3D (Version 3.5), MODPATH, MT3DMS, and SUTRA (Version 2D3D.1). Model Viewer is designed to directly read input and output files from these models, thus minimizing the need for additional postprocessing. This report provides an overview of Model Viewer. Complete instructions on how to use the software are provided in the on-line help pages.

Open-File Report↗

Fracture-correlated lineaments at Great Bay, southeastern New Hampshire

Analysis by remote-sensing techniques and observations of exposed bedrock structure were preliminary steps taken in a study to locate potential bedrock-fracture zones that may store and transmit ground water near Great Bay, N.H. To help correlate lineaments on the surface with fractures, structural measurements were made at exposed bedrock, largely along the shoreline of the bay, and analyzed to identify fracture trends and fracture characteristics. With these fracture data, lineament-filtering techniques, such as (1) buffer analysis around individual lineaments, (2) discrete-measurement analysis by domain, and (3) spacing-normalized analysis by domain, identified 'fracture-correlated lineaments.' Of the 927 lineaments identified in the study area (180 square kilometers), 406 (44 percent) were evaluated because they either were located within 305 meters of an outcrop with fracture data or intersected one of five 3,300-meter-square grid domain cells that encompassed the fracture data. Of the 406 lineaments, 190 (47 percent) are fracture correlated, although only 15 percent were correlated by more than one filtering technique. The large number of lineaments found in areas of thin glacial overburden and high densities of fractured outcrops suggests that filtering techniques are useful in these areas to selectively identify fracture-correlated lineaments. Fractures parallel to bedding in the Kittery Formation are open locally and often associated with vugs, with up to 1-centimeter aperture, and may provide appreciable secondary porosity in this rock unit. Discrete-measurement analysis by domain identified fracture-correlated lineaments with orientations parallel to these open and vug-filled fractures. Fracture-correlated lineaments related to closely spaced fractures were identified by the spacing-normalized analysis by domain. Analysis results may be used to indicate the potential bedrock pathways for ground-water-discharge points along the shoreline of Great Bay.

Open-File Report↗

Southwest Washington coastal erosion workshop report 2000

This report is a compilation of abstracts that correspond to oral presentations and posters presented at the fifth principal investigators workshop of the Southwest Washington Coastal Erosion Study. The workshop was held November 15 - 17, 2000 at the Department of Ecology headquarters building in Olympia, WA. For the fourth consecutive year in November, the workshop convened the entire multi-disciplinary group of scientists and engineers working on the Study or on related projects within the Columbia River littoral cell (CRLC) (Figures 1 and 2). The workshop participants are listed in the List of Attendees section towards the end of this report. The purpose of this workshop was to bring all Study investigators and associated engineers and scientists together to discuss recent work, ongoing tasks, and future research plans in the CRLC. Investigators were asked to present recent data, preliminary interpretations, and research results to invoke discussion and correlation with parallel scientific efforts. The abstracts compiled in this report represent a wealth of information on the CRLC, but because much of the work is in progress, the reader is advised that the information provided herein is preliminary and subject to change.

Oregon;Washington↗

Coastal vulnerability assessment of Cape Cod National Seashore to sea-level rise

A coastal vulnerability index (CVI) was used to map the relative vulnerability of the coast to future sea-level rise within the Cape Cod National Seashore (CACO). The CVI ranks the following in terms of their physical contribution to sea-level rise-related coastal change: geomorphology, regional coastal slope, rate of relative sea-level rise, shoreline change rates, mean tidal range and mean wave height. The rankings for each variable were combined and an index value calculated for 1-minute grid cells covering the park. The CVI highlights those regions where the physical effects of sea-level rise might be the greatest. This approach combines the coastal system's susceptibility to change with its natural ability to adapt to changing environmental conditions, yielding a quantitative, although relative, measure of the park's natural vulnerability to the effects of sea-level rise. CACO consists of high glacial cliffs, beaches, sand spits, and salt marsh wetlands. The areas most vulnerable to sea-level rise are those with the lowest regional coastal slopes, geomorphologic types that are susceptible to inundation, and the highest rates of shoreline change. Most of CACO's infrastructure lies on high elevation uplands away from the shore; most high use areas are accessible by foot only. The CVI provides an objective technique for evaluation and long-term planning by scientists and park managers.

Massachusetts↗

ERF1_2 -- Enhanced River Reach File 2.0

The digital segmented network based on watershed boundaries, ERF1_2, includes enhancements to the U.S. Environmental Protection Agency's River Reach File 1 (RF1) (USEPA, 1996; DeWald and others, 1985) to support national and regional-scale surface water-quality modeling. Alexander and others (1999) developed ERF1, which assessed the hydrologic integrity of the digital reach traces and calculated the mean water time-of-travel in river reaches and reservoirs. ERF1_2 serves as the foundation for SPARROW (Spatially Referenced Regressions (of nutrient transport) on Watershed) modeling. Within the context of a Geographic Information System, SPARROW estimates the proportion of watersheds in the conterminous U.S. with outflow concentrations of several nutrients, including total nitrogen and total phosphorus, (Smith, R.A., Schwarz, G.E., and Alexander, R.B., 1997). This version of the network expands on ERF1 (Version 1.2; Alexander, et al., 1999) and includes the incremental and total drainage area derived from 1-kilometer (km) elevation data for North America. Previous estimates of the water time-of-travel were recomputed for reaches with water-quality monitoring sites that included two reaches. The mean flow and velocity estimates for these split reaches are based on previous estimation methods (Alexander et al., 1999) and are unchanged in ERF1_2. Drainage area calculations provide data used to estimate the contribution of a given nutrient to the outflow. Data estimates depend on the accuracy of node connectivity. Reaches split at water-quality or pesticide-monitoring sites indicate the source point for estimating the contribution and transport of nutrients and their loads throughout the watersheds. The ERF1_2 coverage extends the earlier drainage area founded on the 1-kilometer data for North America (Verdin, 1996; Verdin and Jenson, 1996). A 1-kilometer raster grid of ERF1_2 projected to Lambert Azimuthal Equal Area, NAD 27 Datum (Snyder, 1987), was merged with the HYDRO1K flow direction data set (Verdin and Jenson, 1996) to generate a DEM-based watershed grid, ERF1_2WS_LG. The watershed boundaries are maintained in a raster (grid cell) format as well as a vector (polygon) format for subsequent model analysis. Both the coverage, ERF1_2, and the grid, ERF1_2WS_LG, are available at: URL:http://water.usgs.gov/lookup/getspatial?erf1_2

Open-File Report↗

Modeling GPR data to interpret porosity and DNAPL saturations for calibration of a 3-D multiphase flow simulation

Dense nonaqueous phase liquids (DNAPLs) are a pervasive and persistent category of groundwater contamination. In an effort to better understand their unique subsurface behavior, a controlled and carefully monitored injection of PCE (perchloroethylene), a typical DNAPL, was performed in conjunction with the University of Waterloo at Canadian Forces Base Borden in 1991. Of the various geophysical methods used to monitor the migration of injected PCE, the U.S. Geological Survey collected 500-MHz ground penetrating radar (GPR) data. These data are used in determining calibration parameters for a multiphase flow simulation. GPR data were acquired over time on a fixed two-dimensional surficial grid as the DNAPL was injected into the subsurface. Emphasis is on the method of determining DNAPL saturation values from this time-lapse GPR data set. Interactive full-waveform GPR modeling of regularized field traces resolves relative dielectric permittivity versus depth profiles for pre-injection and later-time data. Modeled values are end members in recursive calculations of the Bruggeman-Hanai-Sen (BHS) mixing formula, yielding interpreted pre-injection porosity and post-injection DNAPL saturation values. The resulting interpreted physical properties of porosity and DNAPL saturation of the Borden test cell, defined on a grid spacing of 50 cm with 1-cm depth resolution, are used as observations for calibration of a 3-D multiphase flow simulation. Calculated values of DNAPL saturation in the subsurface at 14 and 22 hours after the start of injection, from both the GPR and the multiphase flow modeling, are interpolated volumetrically and presented for visual comparison.

Open-File Report↗

Coastal vulnerability assessment of Gulf Islands National Seashore (GUIS) to sea-level rise

A coastal vulnerability index (CVI) was used to map the relative vulnerability of the coast to future sea-level rise within Gulf Islands National Seashore (GUIS) in Mississippi and Florida. The CVI ranks the following in terms of their physical contribution to sea-level rise-related coastal change: geomorphology, regional coastal slope, rate of relative sea-level rise, shoreline change rates, mean tidal range and mean wave height. The rankings for each variable were combined and an index value calculated for 1-minute grid cells covering the park. The CVI highlights those regions where the physical effects of sea-level rise might be the greatest. This approach combines the coastal system's susceptibility to change with its natural ability to adapt to changing environmental conditions, yielding a quantitative, although relative, measure of the park's natural vulnerability to the effects of sea-level rise. The Gulf Islands in Mississippi and Florida consist of stable and washover dominated portions of barrier beach backed by wetland and marsh. The areas likely to be most vulnerable to sea-level rise are those with the highest occurrence of overwash, the highest rates of shoreline change, the gentlest regional coastal slope, and the highest rates of relative sea-level rise. The CVI provides an objective technique for evaluation and long-term planning by scientists and park managers.

Florida, Mississippi↗

Coastal vulnerability assessment of Fire Island National Seashore to sea-level rise

A coastal vulnerability index (CVI) was used to map the relative vulnerability of the coast to future sea-level rise within Fire Island National Seashore (FIIS), New York. The CVI ranks the following in terms of their physical contribution to sea-level rise-related coastal change: geomorphology, regional coastal slope, rate of relative sea-level rise, shoreline change rates, mean tidal range and mean wave height. The rankings for each variable were combined and an index value calculated for 1-minute grid cells covering the park. The CVI highlights those regions where the physical effects of sea-level rise might be the greatest. This approach combines the coastal system's susceptibility to change with its natural ability to adapt to changing environmental conditions, yielding a quantitative, although relative, measure of the park's natural vulnerability to the effects of sea-level rise. Fire Island consists of stable and washover dominated portions of barrier beach backed by lagoons, tidal wetlands and marsh. The areas most vulnerable to sea-level rise are those with the highest historic occurrence of overwash and the highest rates of shoreline change. Implementation of large-scale beach nourishment and other coastal engineering alternatives being considered for Fire Island could alter the CVI computed here. The CVI provides an objective technique for evaluation and long-term planning by scientists and park managers.

New York↗

Mineral Commodity Profiles -- Rubidium

Overview -- Rubidium is a soft, ductile, silvery-white metal that melts at 39.3 ?C. One of the alkali metals, it is positioned in group 1 (or IA) of the periodic table between potassium and cesium. Naturally occurring rubidium is slightly radioactive. Rubidium is an extremely reactive metal--it ignites spontaneously in the presence of air and decomposes water explosively, igniting the liberated hydrogen. Because of its reactivity, the metal and several of its compounds are hazardous materials, and must be stored and transported in isolation from possible reactants. Although rubidium is more abundant in the earth?s crust than copper, lead, or zinc, it forms no minerals of its own, and is, or has been, produced in small quantities as a byproduct of the processing of cesium and lithium ores taken from a few small deposits in Canada, Namibia, and Zambia. In the United States, the metal and its compounds are produced from imported raw materials by at least one company, the Cabot Corporation (Cabot, 2003). Rubidium is used interchangeably or together with cesium in many uses. Its principal application is in specialty glasses used in fiber optic telecommunication systems. Rubidium?s photoemissive properties have led to its use in night-vision devices, photoelectric cells, and photomultiplier tubes. It has several uses in medical science, such as in positron emission tomographic (PET) imaging, the treatment of epilepsy, and the ultracentrifugal separation of nucleic acids and viruses. A dozen or more other uses are known, which include use as a cocatalyst for several organic reactions and in frequency reference oscillators for telecommunications network synchronization. The market for rubidium is extremely small, amounting to 1 to 2 metric tons per year (t/yr) in the United States. World resources are vast compared with demand.

Open-File Report↗

Sea floor topography and backscatter intensity of the Historic Area Remediation Site (HARS), offshore of New York, based on multibeam surveys conducted in 1996, 1998, and 2000

An area offshore of Sandy Hook, New Jersey, has been used extensively for disposal of dredged and other materials, derived from the New York/New Jersey Harbor and surrounding areas, since the late 1800's (Figure 1). Between 1976 and 1995, the New York Bight Dredged Material Disposal Site, also known as the Mud Dump Site (Figure 2), received on average about 6 million cubic yards of material each year from federal and private maintenance dredging and from harbor deepening activities (Massa and others, 1996). In September 1997 the Mud Dump Site (MDS) was closed as an official ocean disposal site by the U.S. Environmental Protection Agency (http://www.epa.gov/), and the MDS and surrounding areas were designated as the Historic Area Remediation Site (HARS). The HARS is subdivided into a Primary Remediation Area (PRA, subdivided into 9 cells), a Buffer Zone, and a No-Discharge Zone (Figure 2). The sea floor of the HARS, approximately 9 square nautical miles in area, is being remediated by placing at least a one-meter cap of Category I (clean) dredged material on top of the existing surface sediments that exhibit varying degrees of degradation. (See http://www.nan.usace.army.mil/business/prjlinks/dmmp/benefic/hars.htm)(Category I sediments have no potential short or long-term impacts and are acceptable for unrestricted ocean disposal (EPA, 1996)). About 1.1 million cubic yards of dredged material for remediation was placed in the HARS in 1999, and 2.5 million cubic yards in 2000. Three multibeam echosounder surveys were carried out to map the topography and surficial geology of the HARS. The surveys were conducted November 23 - December 3, 1996, October 26 - November 11, 1998, and April 6 - 30, 2000. The surveys were carried out as part of a larger survey of the Hudson Shelf Valley and adjacent shelf (Butman and others, 1998, (http://pubs.usgs.gov/openfile/of98-616/). This report presents maps showing topography, shaded relief, and backscatter intensity (a measure of sea floor texture and roughness) at a scale of 1:25,000. Comparison of the topography and backscatter intensity from the three surveys show changes in topography and surficial sediment properties resulting from placement of dredged material in 1996 and 1997 prior to closure of the Mud Dump Site, as well as placement of capping material for remediation of the HARS. The surficial geology and sediments of the HARS and the surrounding region are described in Butman and others (1998), (http://pubs.usgs.gov/openfile/of98-616/, Schwab and others (1997, 2000, http://pubs.usgs.gov/openfile/of00-295/). A history of waste disposal in the New York Bight region is presented in Massa and others (1996).

Open-File Report↗

User's Guide for the Agricultural Non-Point Source (AGNPS) Pollution Model Data Generator

BACKGROUND Throughout this user guide, we refer to datasets that we used in conjunction with developing of this software for supporting cartographic research and producing the datasets to conduct research. However, this software can be used with these datasets or with more 'generic' versions of data of the appropriate type. For example, throughout the guide, we refer to national land cover data (NLCD) and digital elevation model (DEM) data from the U.S. Geological Survey (USGS) at a 30-m resolution, but any digital terrain model or land cover data at any appropriate resolution will produce results. Another key point to keep in mind is to use a consistent data resolution for all the datasets per model run. The U.S. Department of Agriculture (USDA) developed the Agricultural Nonpoint Source (AGNPS) pollution model of watershed hydrology in response to the complex problem of managing nonpoint sources of pollution. AGNPS simulates the behavior of runoff, sediment, and nutrient transport from watersheds that have agriculture as their prime use. The model operates on a cell basis and is a distributed parameter, event-based model. The model requires 22 input parameters. Output parameters are grouped primarily by hydrology, sediment, and chemical output (Young and others, 1995.) Elevation, land cover, and soil are the base data from which to extract the 22 input parameters required by the AGNPS. For automatic parameter extraction, follow the general process described in this guide of extraction from the geospatial data through the AGNPS Data Generator to generate input parameters required by the pollution model (Finn and others, 2002.)

Open-File Report↗

Coastal vulnerability assessment of Assateague Island National Seashore (ASIS) to sea-level rise

A coastal vulnerability index (CVI, http://pubs.usgs.gov/of/2004/1020/html/cvi.htm) was used to map relative vulnerability of the coast to future sea-level rise within Assateague Island National Seashore (ASIS) in Maryland and Virginia. The CVI ranks the following in terms of their physical contribution to sea-level rise-related coastal change: geomorphology, regional coastal slope, rate of relative sea-level rise, shoreline change rates, mean tidal range and mean wave height. Rankings for each variable were combined and an index value calculated for 1-minute grid cells covering the park. The CVI highlights those regions where the physical effects of sea-level rise might be the greatest. This approach combines the coastal system's susceptibility to change with its natural ability to adapt to changing environmental conditions, yielding a quantitative, although relative, measure of the park's natural vulnerability to the effects of sea-level rise. The CVI provides an objective technique for evaluation and long-term planning by scientists and park managers. Assateague Island consists of stable and washover dominated portions of barrier beach backed by wetland and marsh. The areas within Assateague that are likely to be most vulnerable to sea-level rise are those with the highest occurrence of overwash and the highest rates of shoreline change.

Maryland, Virginia↗

Coastal vulnerability assessment of Olympic National Park to sea-level rise

A coastal vulnerability index (CVI) was used to map the relative vulnerability of the coast to future sea-level rise within Olympic National Park (OLYM), Washington. The CVI scores the following in terms of their physical contribution to sea-level rise-related coastal change: geomorphology, regional coastal slope, rate of relative sea-level rise, shoreline change rates, mean tidal range and mean wave height. The rankings for each variable were combined and an index value calculated for 1-minute grid cells covering the park. The CVI highlights those regions where the physical effects of sea-level rise might be the greatest. This approach combines the coastal system's susceptibility to change with its natural ability to adapt to changing environmental conditions, yielding a quantitative, although relative, measure of the park's natural vulnerability to the effects of sea-level rise. The CVI provides an objective technique for evaluation and long-term planning by scientists and park managers. The Olympic National Park coast consists of rocky headlands, pocket beaches, glacial-fluvial features, and sand and gravel beaches. The Olympic coastline that is most vulnerable to sea-level rise are beaches in gently sloping areas.

Washington↗

Quaternary stratigraphy and tectonics, and late prehistoric agriculture of the Safford Basin (Gila and San Simon river valleys), Graham County, Arizona

This guidebook accompanied the 46th annual meeting of the Rocky Mountain Cell of the Friends of the Pleistocene (FOP) and the 2002 Fall Field Trip of the Arizona Geological Society. The meeting and field trip were held in the Safford Basin, southeastern Arizona. The Friends of the Pleistocene is an informal gathering of Quaternary geologists, geomorphologists, and pedologists who meet annually for a field conference. The first part of the guidebook consists of road logs with descriptions of stops covering the three days of the field trip. An overview of the geology of the Safford Basin is given in Stop 1-1. The second part of the guidebook consists of four short papers that discuss adjacent areas or that expand upon the road log descriptions of the field trip stops. The first paper by Reid and Buffler is a summary of upper Cenozoic depositional facies in the Duncan Basin, the first basin to the east of the Safford Basin. The next three papers expand upon (1) the soil study of the gridded field agricultural complex (Stop 2-3, Homburg and Sandor), (2) the vertebrate fossils of the San Simon Valley in the southeastern part of the Safford Basin (Stop 3-1, Thrasher), and (3) paleoIndian irrigation systems and settlements in Lefthand Canyon at the foot of the Pinaleno Mountains (Stop 3-2, Neely and Homburg).

Open-File Report↗

Refine your search to explore more results.