Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Topographic Map”

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
USGS2009

National assessment of historical shoreline change: a pilot study of historical coastal bluff retreat in the Great Lakes, Erie, Pennsylvania

Coastal bluff retreat is a chronic problem along many high-relief coastlines in the United States. As coastal populations continue to grow and community infrastructures are threatened by erosion, there is increased demand for accurate information regard-ing trends and rates of bluff retreat. There is also a need for a comprehensive analysis that is consistent from one coastal region to another. To address these national needs, the U.S. Geological Survey (USGS), as part of the National Assessment of Coastal Change Hazards Project, conducted a pilot study of bluff retreat along the Lake Erie, Pa., coastline to assess the feasibility of undertaking a larger, multi-state analysis in the Great Lakes region. This report provides an overview of the pilot-study location and bluff geomorphology, the data sources and methodology, results of the analysis, and a discussion of the feasibility of undertaking a similar analysis along eroding bluffs in other Great Lakes states. This pilot study is part of an ongoing effort by the USGS to provide a comprehensive analysis of historical shoreline change and cliff and bluff retreat along open-ocean coastlines of the conterminous United States and parts of Hawaii, Alaska, and the Great Lakes. One purpose of the work is to develop standard, repeatable methods for mapping and analyzing coastal change so that systematic and consistent periodic updates of coastal erosion can be made nationally. Bluff-retreat evaluations are conducted by comparing the location of a historical bluff edge digitized from aerial photographs with those of recent bluff edges interpreted from both aerial photographs and lidar topographic surveys. The historical bluff edge is from 1938, whereas the more recent bluff edges are from 1998 and 2006 lidar data. Long-term (68-year) rates of retreat are calculated using the available bluff-edge data. The rates of retreat presented in this report represent conditions from the 1930s to 1998/2006, and are not intended for predicting future bluff-edge positions or rates of retreat. The report presents bluff-retreat rates for 32 km of a 60-km stretch along the Lake Erie, Pa., coastline. Data are discontinuous due to gaps in source data and lack of continuous bluffs. The average rate of coastal bluff retreat for the Lake Erie, Pa., bluffs was -0.3 +- 0.1 m/yr (retreat rates are presented as negative numbers in this report), based on rates averaged from 1,595 individual transects. Retreat rates generally were lowest where bedrock outcrops are exposed as the basal unit in the bluff. The highest rates are associated with anthropogenic activities, including jetties that trap littoral sediment, depleting a source of material for the natural replenishment of protective beaches downcoast, and extensive irrigation of farmlands on the tops of the bluffs, which can destabilize bluffs by enhancing ground-water outflow.

Pennsylvania↗
USGS1989

Potentiometric surface of the Aquia Aquifer in southern Maryland during the fall of 1987

A map showing the potentiometric surface of the Aquia aquifer in the Paleocene Aquia Formation in southern Maryland during the fall of 1987 was prepared by using water level measurements from 80 observation wells. The potentiometric surface was above sea level near the northeastern boundary and outcrop area of the aquifer in topographically high areas of Anne Arundel and Prince Georges Counties. The potentiometric surface was below sea level in the remainder of the study area. The hydraulic gradient was generally to the southeast toward an extensive cone of depression, which is centered around well fields near Lexington Park. Much of the cone is more than 40 ft below sea level. (USGS)

Water-Resources Investigations Report↗
USGS2023

High-frequency time series comparison of Sentinel-1 and Sentinel-2 for open and vegetated water across the United States (2017-2021)

Frequent observations of surface water at fine spatial scales will provide critical data to support the management of aquatic habitat, flood risk and water quality. Sentinel-1 and Sentinel-2 satellites can provide such observations, but algorithms are still needed that perform well across diverse climate and vegetation conditions. We developed surface inundation algorithms for Sentinel-1 and Sentinel-2, respectively, at 12 sites across the conterminous United States (CONUS), covering a total of >536,000 km 2 and representing diverse hydrologic and vegetation landscapes. Each scene in the 5-year (2017–2021) time series was classified into open water, vegetated water, and non-water at 20 m resolution using variables from Sentinel-1 and Sentinel-2, as well as variables derived from topographic and weather datasets. The Sentinel-1 algorithm was developed distinct from the Sentinel-2 model to explore if and where the two time series could potentially be integrated into a single high-frequency time series. Within each model, open water and vegetated water (vegetated palustrine, lacustrine, and riverine wetlands) classes were mapped. The models were validated using imagery from WorldView and PlanetScope. Classification accuracy for open water was high across the 5-year period, with an omission and commission error of only 3.1% and 0.9% for the Sentinel-1 algorithm and 3.1% and 0.5% for the Sentinel-2 algorithm, respectively. Vegetated water accuracy was lower, as expected given that the class represents mixed pixels. The Sentinel-2 algorithm showed higher accuracy (10.7% omission and 7.9% commission error) relative to the Sentinel-1 algorithm (28.4% omission and 16.0% commission error). Patterns over time in the proportion of area mapped as open or vegetated water by the Sentinel-1 and Sentinel-2 algorithms were charted and correlated for a subset of all 12 sites. Our results showed that the Sentinel-1 and Sentinel-2 algorithm open water time series can be integrated at all 12 sites to improve the temporal resolution, but sensor-specific differences, such as sensitivity to vegetation structure versus pixel color, complicate the data integration for mixed-pixel, vegetated water. The methods developed here provide inundation at 5-day (Sentinel-2 algorithm) and 12-day (Sentinel-1 algorithm) time steps to improve our understanding of the short- and long-term response of surface water to climate and land use drivers in different ecoregions.

continental United States↗
USGS2006

Tectonic stressing in California modeled from GPS observations

What happens in the crust as a result of geodetically observed secular motions? In this paper we find out by distorting a finite element model of California using GPS-derived displacements. A complex model was constructed using spatially varying crustal thickness, geothermal gradient, topography, and creeping faults. GPS velocity observations were interpolated and extrapolated across the model and boundary condition areas, and the model was loaded according to 5-year displacements. Results map highest differential stressing rates in a 200-km-wide band along the Pacific-North American plate boundary, coinciding with regions of greatest seismic energy release. Away from the plate boundary, GPS-derived crustal strain reduces modeled differential stress in some places, suggesting that some crustal motions are related to topographic collapse. Calculated stressing rates can be resolved onto fault planes: useful for addressing fault interactions and necessary for calculating earthquake advances or delays. As an example, I examine seismic quiescence on the Garlock fault despite a calculated minimum 0.1-0.4 MPa static stress increase from the 1857 M???7.8 Fort Tejon earthquake. Results from finite element modeling show very low to negative secular Coulomb stress growth on the Garlock fault, suggesting that the stress state may have been too low for large earthquake triggering. Thus the Garlock fault may only be stressed by San Andreas fault slip, a loading pattern that could explain its erratic rupture history.

Journal of Geophysical Research B: Solid Earth↗
USGS1995

Potentiometric surface of the Aquia Aquifer in southern Maryland, September 1994

A map showing the potentiometric surface of the Aquia aquifer in the Aquia aquifer in the Aquia Formation of Paleocene age in southern Maryland during September 1994 was prepared from water levels measured in 74 wells. The potentiometric surface was above sea level near the northwestern boundary and outcrop area of the aquifer in topographically high areas of Anne Arundel and Prince Georges Counties, and was below sea level in the remainder of the study area. The hydraulic gradient was directed southeastward toward an extensive cone of depression around well fields at Lexington Park and Solomons Island. Ground-water levels were more than 70 feet below sea level in a 60-square-mile area surrounding the deepest part of the cone of depression and 131 feet below sea level in one well near the center of the cone.

Open-File Report↗
USGS1996

Potentiometric Surface of the Lower Patapsco Aquifer in Southern Maryland, September 1995

A map showing the potentiometric surface of the lower Patapsco aquifer in the Patapsco Formation of Cretaceous age in southern Maryland during September 1995 was prepared from water-level measurements in 69 wells. The potentiometric surface was 87 feet above sea level near the northwestern boundary and outcrop area of the aquifer in topographically high areas of Anne Arundel County. From this high area, the potentiometric surface declined to the south and southeast toward large well fields at Crofton, Severndale, Arnold, Annapolis, Waldorf, LaPlata, Indian Head, and Morgantown. The ground-water levels reached 18 feet below sea level at Crofton, 62 feet below sea level at Severndale, 26 feet below sea level at Arnold, more than 20 feet below sea level near Annapolis, 160 feet below sea level in the Waldorf area, nearly 100 feet below sea level at LaPlata, more than 110 feet below sea level at Indian Head, and 75 feet below sea level at the Morgantown powerplant.

Open-File Report↗
USGS2002

Basis and methods of NASA airborne topographic mapper lidar surveys for coastal studies

This paper provides an overview of the basic principles of airborne laser altimetry for surveys of coastal topography, and describes the methods used in the acquisition and processing of NASA Airborne Topographic Mapper (ATM) surveys that cover much of the conterminous US coastline. This form of remote sensing, also known as "topographic lidar", has undergone extremely rapid development during the last two decades, and has the potential to contribute within a wide range of coastal scientific investigations. Various airborne laser surveying (ALS) applications that are relevant to coastal studies are being pursued by researchers in a range of Earth science disciplines. Examples include the mapping of "bald earth" land surfaces below even moderately dense vegetation in studies of geologic framework and hydrology, and determination of the vegetation canopy structure, a key variable in mapping wildlife habitats. ALS has also proven to be an excellent method for the regional mapping of geomorphic change along barrier island beaches and other sandy coasts due to storms or long-term sedimentary processes. Coastal scientists are adopting ALS as a basic method in the study of an array of additional coastal topics. ALS can provide useful information in the analysis of shoreline change, the prediction and assessment of landslides along seacliffs and headlands, examination of subsidence causing coastal land loss, and in predicting storm surge and tsunami inundation.

Journal of Coastal Research↗
USGS2005

Climatic and topographic controls on the style and timing of Late Quaternary glaciation throughout Tibet and the Himalaya defined by 10Be cosmogenic radionuclide surface exposure dating

Temporal and spatial changes in glacier cover throughout the Late Quaternary in Tibet and the bordering mountains are poorly defined because of the inaccessibility and vastness of the region, and the lack of numerical dating. To help reconstruct the timing and extent of glaciation throughout Tibet and the bordering mountains, we use geomorphic mapping and 10Be cosmogenic radionuclide (CRN) surface dating in study areas in southeastern (Gonga Shan), southern (Karola Pass) and central (Western Nyainqentanggulha Shan and Tanggula Shan) Tibet, and we compare these with recently determined numerical chronologies in other parts of the plateau and its borderlands. Each of the study regions receives its precipitation mainly during the south Asian summer monsoon when it falls as snow at high altitudes. Gonga Shan receives the most precipitation (>2000 mm a-1) while, near the margins of monsoon influence, the Karola Pass receives moderate amounts of precipitation (500-600 mm a-1) and, in the interior of the plateau, little precipitation falls on the western Nyainqentanggulha Shan (???300 mm a -1) and the Tanggula Shan (400-700 mm a-1). The higher precipitation values for the Tanggula Shan are due to strong orographic effects. In each region, at least three sets of moraines and associated landforms are preserved, providing evidence for multiple glaciations. The 10Be CRN surface exposure dating shows that the formation of moraines in Gonga Shan occurred during the early-mid Holocene, Neoglacial and Little Ice Age, on the Karola Pass during the Lateglacial, Early Holocene and Neoglacial, in the Nyainqentanggulha Shan date during the early part of the last glacial cycle, global Last Glacial Maximum and Lateglacial, and on the Tanggula Shan during the penultimate glacial cycle and the early part of the last glacial cycle. The oldest moraine succession in each of these regions varies from the early Holocene (Gonga Shan), Lateglacial (Karola Pass), early Last Glacial (western Nyainqentanggulha Shan), and penultimate glacial cycle (Tanggula Shan). We believe that the regional patterns and timing of glaciation reflect temporal and spatial variability in the south Asian monsoon and, in particular, in regional precipitation gradients. In zones of greater aridity, the extent of glaciation has become increasingly restricted throughout the Late Quaternary leading to the preservation of old (???100 ka) glacial landforms. In contrast, in regions that are very strongly influenced by the monsoon (???1600 mm a-1), the preservation potential of pre-Lateglacial moraine successions is generally extremely poor. This is possibly because Lateglacial and Holocene glacial advances may have been more extensive than early glaciations and hence may have destroyed any landform or sedimentary evidence of earlier glaciations. Furthermore, the intense denudation, mainly by fluvial and mass movement processes, which characterize these wetter environments, results in rapid erosion and re-sedimentation of glacial and associated landforms, which also contributes to their poor preservation potential. ?? 2004 Elsevier Ltd. All rights reserved.

Quaternary Science Reviews↗
USGS1996

Potentiometric Surface of the Upper Patapsco Aquifer in Southern Maryland, September 1995

A map showing the potentiometric surface of the upper Patapsco aquifer in the Patapsco Formation of Cretaceous age in southern Maryland during September 1995 was prepared from water-level measurements in 42 wells. The potentiometric surface was nearly 120 feet above sea level near the northwestern boundary and outcrop area of the aquifer in topographically high areas of Anne Arundel County, and 55 feet above sea level in a similar setting in Prince Georges County. From these high areas, the potentiometic surface declined to the south and southeast toward large well fields in the Annapolis and Waldorf areas and at the Chalk Point powerplant. Ground-water levels reached nearly 30 feet below sea level in the Annapolis area, 113 feet below sea level southwest of Waldorf, and more than 30 feet below sea level at the Chalk Point powerplant.

Open-File Report↗
USGS1973

Map showing depth to top of the principal aquifer, Sugar House quadrangle, Salt Lake County, Utah, February 1972

The depth to the top of the principal aquifer in the Sugar House quadrangle ranges from about 50 feet to more than 400 feet below land surface. The principal aquifer supplied about 4 percent, or 9,000 acre-feet, of the municipal and industrial water used annually in Salt Lake County curing 1964-68. The least depths occur in topographically low areas of the Jordan Valley, such as near Murray in the southwest corner of the Sugar House quadrangle. The greatest depths occur near the mountain front in the vicinity of Mill creek and Neffs Canyons where thick alluvial-fan deposits overlie the principal aquifer.

Utah↗
USGS2012

Mapping plant species ranges in the Hawaiian Islands: developing a methodology and associated GIS layers

This report documents a methodology for projecting the geographic ranges of plant species in the Hawaiian Islands. The methodology consists primarily of the creation of several geographic information system (GIS) data layers depicting attributes related to the geographic ranges of plant species. The most important spatial-data layer generated here is an objectively defined classification of climate as it pertains to the distribution of plant species. By examining previous zonal-vegetation classifications in light of spatially detailed climate data, broad zones of climate relevant to contemporary concepts of vegetation in the Hawaiian Islands can be explicitly defined. Other spatial-data layers presented here include the following: substrate age, as large areas of the island of Hawai'i, in particular, are covered by very young lava flows inimical to the growth of many plant species; biogeographic regions of the larger islands that are composites of multiple volcanoes, as many of their species are restricted to a given topographically isolated mountain or a specified group of them; and human impact, which can reduce the range of many species relative to where they formerly were found. Other factors influencing the geographic ranges of species that are discussed here but not developed further, owing to limitations in rendering them spatially, include topography, soils, and disturbance. A method is described for analyzing these layers in a GIS, in conjunction with a database of species distributions, to project the ranges of plant species, which include both the potential range prior to human disturbance and the projected present range. Examples of range maps for several species are given as case studies that demonstrate different spatial characteristics of range. Several potential applications of species-range maps are discussed, including facilitating field surveys, informing restoration efforts, studying range size and rarity, studying biodiversity, managing invasive species, and planning of conservation efforts.

Hawai'i↗
USGS2026

Compounding of 100-year coastal floods by rainfall in an urban environment

Coastal and pluvial flooding are both becoming more prevalent and severe due to climate change and urbanization in floodplains. The co-occurrence of these flood drivers is generally assumed to exacerbate the resulting flood impacts, a result referred to as compound flooding. However, few observational or modeling studies have investigated the circumstances under which this occurs. Here, we study the impacts of these combined flood drivers and evaluate the implicit hypothesis of official flood maps, which is that rainfall has a negligible impact on the flood depth and flooded area due to a 100 year coastal flood. A coastal system model, configured to capture coastal and pluvial flood drivers, is used. We evaluate the flooding for different urban landform types, including coastal landfill (human-made land), convergent areas (topographic depressions) and other urban terrain, within a model domain covering the Jamaica Bay watershed of New York City. A scenario-based strategy is adopted with a 100 year coastal flood as a control simulation, to which we add a set of realistic scenarios of rainfall data from historical tropical cyclones. We also apply a joint probability analysis framework with historical data to evaluate the probability of these compound coastal-pluvial scenarios. Results reveal cases where the pluvial driver compounds the coastal flood through expansion of the flood zone, with a 17% chance of rainfall increasing the flood area by 6%–38%, and a 5% chance of an increase of 61%–73%. It is rare that floods are significantly deepened but when deepening occurs, it is more common for the convergent zone than for the coastal landfill. These findings quantitatively assess the potential of the pluvial driver to exacerbate flooding, which may influence emergency management strategies such as evacuation plans, shelter arrangements, and related preparedness measures.

New York↗
USGS2016

Creating a Coastal National Elevation Database (CoNED) for science and conservation applications

The U.S. Geological Survey is creating the Coastal National Elevation Database, an expanding set of topobathymetric elevation models that extend seamlessly across coastal regions of high societal or ecological significance in the United States that are undergoing rapid change or are threatened by inundation hazards. Topobathymetric elevation models are raster datasets useful for inundation prediction and other earth science applications, such as the development of sediment-transport and storm surge models. These topobathymetric elevation models are being constructed by the broad regional assimilation of numerous topographic and bathymetric datasets, and are intended to fulfill the pressing needs of decision makers establishing policies for hazard mitigation and emergency preparedness, coastal managers tasked with coastal planning compatible with predictions of inundation due to sea-level rise, and scientists investigating processes of coastal geomorphic change. A key priority of this coastal elevation mapping effort is to foster collaborative lidar acquisitions that meet the standards of the USGS National Geospatial Program's 3D Elevation Program, a nationwide initiative to systematically collect high-quality elevation data. The focus regions are located in highly dynamic environments, for example in areas subject to shoreline change, rapid wetland loss, hurricane impacts such as overwash and wave scouring, and/or human-induced changes to coastal topography.

Journal of Coastal Research↗
USGS1995

Potentiometric surface of the upper Patapsco Aquifer in southern Maryland, September 1994

A map showing the potentiometric surface of the upper Patapsco aquifer in the Patapsco Formation of Cretaceous age in southern Maryland during September 1994 was prepared from water levels measured in 43 wells. The potentiometric surface was at least 70 feet above sea level near the northwestern boundary and outcrop area of the aquifer in a topographically high area of Anne Arundel County, and nearly 60 feet above sea level in a similar setting in Prince Georges County. From these high areas, the potentiometric surface declined to the south and southeast toward large well fields in the Annapolis and Waldorf areas and at the Chalk Point powerplant. Ground-water levels reached nearly 30 feet below sea level in the Annapolis area, nearly 110 feet below sea level southwest of Waldorf, and more than 25 feet below sea level at the Chalk Point powerplant.

Open-File Report↗
USGS1995

Potentiometric surface of the lower Patapsco Aquifer in southern Maryland, September 1994

A map showing the potentiometric surface of the lower Patapsco aquifer in the Patapsco Formation of Cretaceous age in southern Maryland during September 1994 was prepared from water levels measured in 56 wells. The potentiometric surface was at least 104 feet above sea level near the northwestern boundary and outcrop area of the aquifer in a topographically high area of Prince Georges County, and at least 89 feet above sea level in a similar setting in Anne Arundel County. From these high areas, the potentiometric surface declined to the south and southeast toward large well fields at Crofton, Severndle, Arnold, Annapolis, Waldorf, La Plata, Indian Head, and Morgantown. The ground-water levels reached 20 feet below sea level at Crofton, 80 feet below sea level at Severndale, 25 feet below sea level at Arnold, more than 10 feet below sea level near Annapolis, more than 130 feet below sea level in the Waldorf area, more than 100 feet below sea level at La Plata an Indian Head, an nearly 80 feet below sea level at the Morgantown powerplant.

Open-File Report↗
USGS1988

Methods for delineating flood-prone areas in the Great Basin of Nevada and adjacent states

The Great Basin is a region of about 210,000 square miles having no surface drainage to the ocean; it includes most of Nevada and parts of Utah, California, Oregon, Idaho, and Wyoming. The area is characterized by many parallel mountain ranges and valleys trending north-south. Stream channels usually are well defined and steep within the mountains, but on reaching the alluvial fan at the canyon mouth, they may diverge into numerous distributary channels, be discontinuous near the apex of the fan, or be deeply entrenched in the alluvial deposits. Larger rivers normally have well-defined channels to or across the valley floors, but all terminate at lakes or playas. Major floods occur in most parts of the Great Basin and result from snowmelt, frontal-storm rainfall, and localized convective rainfall. Snowmelt floods typically occur during April-June. Floods resulting from frontal rain and frontal rain on snow generally occur during November-March. Floods resulting from convective-type rainfall during localized thunderstorms occur most commonly during the summer months. Methods for delineating flood-prone areas are grouped into five general categories: Detailed, historical, analytical, physiographic, and reconnaissance. The detailed and historical methods are comprehensive methods; the analytical and physiographic are intermediate; and the reconnaissance method is only approximate. Other than the reconnaissance method, each method requires determination of a T-year discharge (the peak rate of flow during a flood with long-term average recurrence interval of T years) and T-year profile and the development of a flood-boundary map. The procedure is different, however, for each method. Appraisal of the applicability of each method included consideration of its technical soundness, limitations and uncertainties, ease of use, and costs in time and money. Of the five methods, the detailed method is probably the most accurate, though most expensive. It is applicable to hydraulic and topographic conditions found in many parts of the Great Basin. The historical method is also applicable over a wide range of conditions and is less expensive than the detailed method. However, it requires more historical flood data than are usually available, and experience and judgement are needed to obtain meaningful results. The analytical method is also less expensive than the detailed method and can be used over a wide range of conditions in which the T-year discharge can be determined directly. Experience, good judgement, and thorough knowledge of hydraulic principles are required to obtain adequate results, and the method has limited application in other than rigid-channel situations. The physiographic method is applicable to rigid-boundary channels and is less accurate than the detailed method. The reconnaissance method is relatively imprecise, but it may be the most rational method to use on alluvial fans or valley floors with discontinuous channels. In general, a comprehensive method is most suitable for use with rigid-bank streams in urban areas; only an approximate method seems justified in undeveloped areas.

Water Supply Paper↗
USGS2023

Spatially explicit models of seed availability improve predictions of conifer regeneration following the 2018 Carr Fire in northern California

For many conifer species in dry conifer forests of North America, seeds must be present for postfire regeneration to occur, suggesting that seed dispersal from surviving trees plays a critical role in postfire forest recovery. However, the application of tree fecundity and spatial arrangement to postfire conifer recovery predictions have only recently become more common, and is often included at relatively coarse scales (i.e., 30 meters). In this study, we mapped surviving trees using lidar and created a spatially explicit estimate of seed density (seed shadows) with 10 m, 50 m, and 100 m median dispersal distances. We estimated the number of seeds produced by each tree using allometric relationships between tree size and fecundity. Along with the seed shadows, we used a suite of topographic variables as inputs to negative binomial hurdle models to predict conifer seedling abundance in 131 plots following the 2018 Carr Fire in northern California, USA. We compared models using each of the seed shadows to each other as well as to a model using the distance to the nearest surviving tree, which served as a baseline. All model formulations indicated that estimated seed availability was positively associated with conifer regeneration. Despite the importance of seed availability plays in regeneration and the substantial differences in seed availability represented by the different seed shadows in our analysis, we found surprisingly little difference in model performance regardless of which seed shadow was used. However, the models employing seed shadows outperformed the models with distance to the nearest live tree. Although we have demonstrated a modest improvement in predicting postfire conifer regeneration, the uncertainty in our results highlights the importance of tree detection and classification in future studies of this kind. Future studies may find it useful to consider other factors such as predation, site suitability, and seed mortality as potential drivers of discrepancies between total and realized dispersal kernels.

California↗
USGS2025

Assessing the topographic distribution of legacy soil phosphorus in agricultural fields of the Delmarva Peninsula, Mid-Atlantic Coastal Plain, USA

Phosphorus (P) management remains a challenge in agricultural watersheds. The Choptank River Conservation Effects Assessment Project watershed, located in Maryland and Delaware and draining to the Chesapeake Bay, contains legacy soil P from historical dairy and poultry manure applications. These practices elevated soil P beyond crop needs, contributing to persistent P export to aquatic ecosystems. We assessed spatial P distribution and analyzed GIS (Geographic Information Systems)-derived landscape features driving legacy P movement on a farm (47 ha). We hypothesized that P accumulates in drained lowlands and depressional areas due to gravity-driven processes that accelerate P-enriched water to receiving waters via overland flow. In collaboration with the US Department of Agriculture Legacy P project, we collected 105 soil samples (0- to 5-cm and 5- to 15-cm depths) and 14 ditch sediment samples across five topographic openness classes from a farm with >100 years of dairy manure application. Average Mehlich-III P concentrations were 218 and 179 mg kg −1 at 0- to 5-cm and 5- to 15-cm depths, respectively, with legacy areas defined by P content > 100 mg kg −1 . Soil P and clay particle size were positively correlated ( r = 0.42, p < 0.05), increased as landscape openness decreased, and were negatively correlated with topographic openness (ranging from −0.2 to −0.4, p < 0.05), indicating accumulation of P and clay in low-lying areas. These patterns suggest that historical field-level managements have primarily shaped P distribution, while hydrologic and landscape properties further influence its redistribution via transport pathways and drainage. These findings support the development of landscape models to map critical source areas in low-relief watersheds and guide targeted mitigation in high-risk P export zones.

Maryland↗