Examining Landscape Dynamics at a Watershed Scale Using Landsat TM Imagery for Detection of Wintering Hooded Crane Decline in Yashiro,Japan

It is usually inappropriate to define rectangular land areas or administrative units as the extent for quantifying landscapes that possess hierarchical structure. As a functional unit established by geophysical relationships, the watershed is one of many natural scales in the hierarchical landscape. We examined the dynamics of the Yashiro watershed of Japan at the landscape level using pattern metrics based on Landsat thematic mapper (TM) imagery from 1985 to 1998. This watershed provides important habitats for the hooded crane (Grus monachus), a vulnerable species. While its world population has remained stable, the number wintering at Yashiro has declined in recent years. Changes in landscape metrics reveal that the spatial pattern within the watershed underwent homogenization due to depopulation of local people and shifts in local energy requirements and forest management policy at Yashiro. Specific changes include: a decrease in bare land area from 6.2% to 1.0% of the landscape, increased forest cover from 69.2% to 76.1%, reduction in patch number from 1194 to 616 and enlarged mean patch size, and a decrease in total edge from 223,740 m to 158,040 m. The rate of change in landscape metrics indicates a rapid change towards homogeneity in the landscape since 1990. The temporal changes in hooded crane populations corresponded to the changes in landscape. An alternative explanation has been proposed that decline of the species is influenced by landscape dynamics affecting both habitat selection and food resources. Conservation at the watershed scale is suggested to be complementary to the current conservation measures of the species.     

Landscape Changes in the Mullica River Basin of the Pinelands National Reserve,New Jersey,USA

In 1979, the Pinelands Commission was established as a regional land-use planning and regulatory agency charged with the implementation of a Comprehensive Management Plan (CMP) for the Pinelands National Reserve (New Jersey, USA). The CMP was created to balance land preservation and development interests in the Reserve. Because water-quality degradation from developed and agricultural landscapes is associated with changes in the composition of biological communities, monitoring landscape changes provides one of the most direct measures of the impact of land-use policies on the Pinelands ecosystem. In this study, we prepared detailed, land-cover maps within randomly selected aerial-photograph plots for a major watershed in the Reserve. We used these land-cover maps to quantify changes in landscape composition and structure (i.e., patch size, patch area, and number of patches) and characterize land-cover transitions in the basin between 1979 and 1991. Because the study period represented the first 12 years of the regional-planning effort in the Reserve, we evaluated the relationship between land-cover transitions and Commission management-area designations which permit different land-use intensities. Although the landscape composition was similar in 1979 and 1991, we found an increase in the total area and number of developed-land, managed-grassland, and barren-land patches. An increase in the number of patches and a decrease in the total area and median and third-quartile patch sizes for forest land and for all patches regardless of cover type indicated that fragmentation of forest land and the landscape as a whole occurred during the study period. The major land-cover transitions that occurred during the period were the loss of forest land to development and associated cover types and the conversion of one agricultural type to another. Overall, land-cover transitions during the period were found to be consistent with the Commission management-area designations, which indicated that the regional-planning effort has been successful in directing human activities to appropriate areas of the basin.     

The Control of Land-Use Patterns for Stormwater Management at Multiple Spatial Scales

While most research about the relationship between land use and watershed hydrological output has focused primarily on land-use types and their impact on hydrological processes, the relationship between characteristics of land-use patterns (such as pattern fragmentation, connectivity, and coherence) and hydrological processes has not been well examined. Using historical stormwater data, this study evaluates the hydrological effects of different land-use scenarios in the Qing-shui watershed in Beijing, China, at a variety of spatial scales. This study demonstrates that planning and managing land-use patterns can significantly reduce runoff under different scales, particularly for small storm events. In contrast to other aspects of land-use structure characteristics, such as the shape complexity of land-use patches, fragmented level of the patches of land-use types appear as dominant drivers of runoff. The results of the study suggest that land-use pattern management should be an important component of Best Management Practices to reduce the impacts of urbanization on natural hydrological processes.     

Multi-scale analysis of oxygen demand trends in an urbanizing Oregon watershed, USA

Human alteration of the landscape has an extensive influence on the biogeochemical processes that drive oxygen cycling in streams. We estimated trends from the mid-1990s to 2003, using the seasonal Mann-Kendall's test, for percent saturation dissolved oxygen (DO), chemical oxygen demand (COD), total Kjeldahl nitrogen (TKN), and ammonia-nitrogen (NH(3)-N) for 12 sites in the Rock Creek watershed, northwest Oregon, USA. In order to understand the influence of landscape change, scale, and stormwater runoff management on dissolved oxygen trends, we calculated land cover change through aerial photo interpretation at full-basin, local (near sample point) basin, and 100m stream buffer scales, for the years 1994 and 2000. Significant (p < or = 0.05) trends occurred in DO (increasing at five sites), COD (decreasing at seven sites), TKN (decreasing at five sites, increasing at one site), and NH(3)-N (decreasing at one site, increasing at one site). Significant land cover change occurred in agricultural land cover (-8% for the entire basin area) and residential land cover (+10% for the entire basin area) (p < or = 0.05). Correlation results indicated that: (1) forest cover negatively influenced COD at the full basin scale and positively influences NH(3)-N at local scales, (2) residential land cover influenced oxygen demand variables at local scales, (3) agricultural land cover did not influence oxygen demand, (4) local topography negatively influenced TKN and NH(3)-N, and (5) stormwater runoff management infrastructure correlated positively with COD at the local scale. This study indicates that landscape factors influencing DO conditions for the study streams act at multiple scales, suggesting that better knowledge of scale-process interactions can guide watershed managers' decision making in order to maintain improving water quality conditions.     

Modeling the spatial dynamics of regional land use: the CLUE-S model

Land-use change models are important tools for integrated environmental management. Through scenario analysis they can help to identify near-future critical locations in the face of environmental change. A dynamic, spatially explicit, land-use change model is presented for the regional scale: CLUE-S. The model is specifically developed for the analysis of land use in small regions (e.g., a watershed or province) at a fine spatial resolution. The model structure is based on systems theory to allow the integrated analysis of land-use change in relation to socio-economic and biophysical driving factors. The model explicitly addresses the hierarchical organization of land use systems, spatial connectivity between locations and stability. Stability is incorporated by a set of variables that define the relative elasticity of the actual land-use type to conversion. The user can specify these settings based on expert knowledge or survey data. Two applications of the model in the Philippines and Malaysia are used to illustrate the functioning of the model and its validation.     

INTEGRATING LANDSCAPE ASSESSMENT AND HYDROLOGIC MODELING FOR LAND COVER CHANGE ANALYSIS1

ABSTRACT: Significant land cover changes have occurred in the watersheds that contribute runoff to the upper San Pedro River in Sonora, Mexico, and southeast Arizona. These changes, observed using a series of remotely sensed images taken in the 1970s, 1980s, and 1990s, have been implicated in the alteration of the basin hydrologic response. The Cannonsville subwatershed, located in the Catskill/Delaware watershed complex that delivers water to New York City, provides a contrast in land cover change. In this region, the Cannonsville watershed condition has improved over a comparable time period. A landscape assessment tool using a geographic information system (GIS) has been developed that automates the parameterization of the Soil and Water Assessment Tool (SWAT) and KINEmatic Runoff and EROSion (KINEROS) hydrologic models. The Automated Geospatial Watershed Assessment (AGWA) tool was used to prepare parameter input files for the Upper San Pedro Basin, a subwatershed within the San Pedro undergoing significant changes, and the Cannonsville watershed using historical land cover data. Runoff and sediment yield were simulated using these models. In the Cannonsville watershed, land cover change had a beneficial impact on modeled watershed response due to the transition from agriculture to forest land cover. Simulation results for the San Pedro indicate that increasing urban and agricultural areas and the simultaneous invasion of woody plants and decline of grasslands resulted in increased annual and event runoff volumes, flashier flood response, and decreased water quality due to sediment loading. These results demonstrate the usefulness of integrating remote sensing and distributed hydrologic models through the use of GIS for assessing watershed condition and the relative impacts of land cover transitions on hydrologic response.     

EFFECTS OF SCALE ON LAND USE AND WATER QUALITY RELATIONSHIPS: A LONGITUDINAL BASIN‐WIDE PERSPECTIVE1

ABSTRACT: Human land use is a major source of change in catchments in developing areas. To better anticipate the long‐term effects of growth, land use planning requires estimates of how changes in land use will affect ecosystem processes and patterns across multiple scales of space and time. The complexity of biogeochemical and hydrologic interactions within a basin makes it difficult to scale up from process‐based studies of individual reaches to watershed scales over multiple decades. Empirical models relating land use/land cover (LULC) to water quality can be useful in long‐term planning, but require an understanding of the effects of scale on apparent land use‐water quality relationships. We empirically determined how apparent relationships between water quality and LULC data change at different scales, using LIJLC data from the Willapa Bay watershed (Washington) and water quality data collected along the Willapa and North Rivers. Spatial scales examined ranged from the local riparian scale to total upstream catchment. The strength of the correlations between LTJLC data and longitudinal water quality trends varied with scale. Different water quality parameters also varied in their response to changes in scale. Intermediate scales of land use data generally were better predictors than local riparian or total catchment scales. Additional data from the stream network did not increase the strength of relationships significantly. Because of the likelihood of scale‐induced artifacts, studies quantifying land use‐water quality relationships performed at single scales should be viewed with great caution.     

Landscape Character and Fish Assemblage Structure and Function in Western Lake Superior Streams: General Relationships and Identification of Thresholds

As part of a comparative watershed project investigating land-cover/land-use disturbance gradients for streams in the western Lake Superior Basin, we examined general relationships between landscape character and fish assemblage structure and function. We also examined the shape of those relationships to identify discontinuity thresholds where small changes in landscape character were associated with marked shifts in the fish assemblages. After completing a geographic analysis of second- and third-order watersheds in the western Lake Superior drainage, we randomly selected 48 streams along mature forest and watershed storage gradients in 2 hydrogeomorphic regions as our study sites. During the summers of 1997 and 1998, we used electrofishing to sample fish assemblages from each stream. Each of the landscape factors was significantly associated with fish assemblage structure and function based on analysis of covariance. Watershed storage was related to the greatest number of fish assemblage characteristics, but hydrogeopmorphic region and mature forest cover were strongly associated as well. The hydrogeomorphic region also mediated relationships between watershed character and fish assemblages. Discontinuity thresholds for our fish assemblages averaged 11% for watershed storage and 50% for watershed mature forest cover based on piecewise regression analysis. Although many of the landscape–fish relationships might have been manifest through effects on in-stream habitat, our results highlight the importance of management and land-use planning decisions at the watershed and landscape scales.Published online     

Exploring changes in Ecuadorian land use for food production and their effects on natural resources

A model is used for the dynamic and spatially explicit exploration of near future agricultural land-use changes. In a case study for Ecuador, different plausible scenarios are formulated, taking into account possible developments in national food demand until the year 2010. The protection of nature parks and restrictions due to land degradation are evaluated with respect to their possible spatial impacts on the land-use change dynamics within the country. Under the assumptions of the demand scenarios, agricultural land-use expands significantly, resulting in more use of land in existing agricultural areas and frontier-type expansion into relatively undisturbed natural areas. The patterns of change depend on the increase in demand, competition between land-use types, changes in driving factors of land use, and the area and characteristics of land that is excluded from agricultural use. The modelled land-use dynamics are related to their possible impacts on the natural resource base, specifically soil fertility. The results indicate potential negative effects of land-use changes on the soil nutrient balance and biodiversity. It is argued that spatial and temporal quantification of land-use dynamics at the landscape level can support research and policies aimed at understanding the driving factors of land-use change and the behaviour of complex agro-ecosystems under changing conditions at different scales. In this way, issues dealing with sustainable food production and the management of the natural resource base can be addressed in a more integrated and quantitative manner.     

FORECASTING FUTURE LAND USE FOR WATERSHED ASSESSMENT1

ABSTRACT: Protecting surface water quality in watersheds undergoing demographic change requires both the management of existing threats and planning to address potential future stresses arising from changing land use. Many reservoirs and threatened waterbodies are located in areas undergoing rapid population growth, and increases in density of residential and commercial land use, accompanied by increased amount of impervious surface area, can result in increased pollutant loading and degradation of water quality. Effective planning to address potential threats, including zoning and growth management, requires analytical tools to predict and compare the impacts of different management options. The focus of this paper is not on developing demographic projections, but rather the translation of such projections into changes in land use which form the basis for assessment of future watershed loads. Land use change can be forecast at a variety of spatial and temporal scales. A semi-lumped, GIS-based, transition matrix approach is recommended as consistent with the level of complexity achievable in most watershed models. Practical aspects of forecasting future land use for watershed assessment are discussed. Several recent reservoir water supply projection studies are used to demonstrate a general framework for simulating changes in land use and resulting impacts on water quality. In addition to providing a technical basis for selecting optimal management alternatives, such a tool is invaluable for demonstrating to different stakeholder groups the trade-offs among management alternatives, both in terms of water quality and future land use patterns within the watershed.     

Relationships between landscape pattern, wetland characteristics, and water quality in agricultural catchments

Water quality in streams is dependent on landscape metrics at catchment and wetland scales. A study was undertaken to evaluate the correlation between landscape metrics, namely patch density and area, shape, heterogeneity, aggregation, connectivity, land-use ratio, and water quality variables (salinity, nutrients, sediments, alkalinity, other potential pollutants and pH) in the agricultural areas of a semiarid Mediterranean region dominated by irrigated farmlands (NE Spain). The study also aims to develop wetland construction criteria in agricultural catchments. The percentage of arable land and landscape homogeneity (low value of Simpson index) are significantly correlated with salinity (r(2) = 0.72) and NO(3)-N variables (r(2) = 0.49) at catchment scale. The number of stock farms was correlated (Spearman's corr. = 0.60; p < 0.01) with TP concentration in stream water. The relative abundance of wetlands and the aggregation of its patches influence salinity variables at wetland scale (r(2) = 0.59 for Na(+) and K(+) concentrations). The number and aggregation of wetland patches are closely correlated to the landscape complexity of catchments, measured as patch density (r(2) = 0.69), patch size (r(2) = 0.53), and landscape heterogeneity (r(2) = 0.62). These results suggest that more effective results in water quality improvement would be achieved if we acted at both catchment and wetland scales, especially reducing landscape homogeneity and creating numerous wetlands scattered throughout the catchment. A set of guidelines for planners and decision makers is provided for future agricultural developments or to improve existing ones.     

基于RS与GIS的福建省土地利用动态变化研究

利用1996年和2001年2期1：10万TM假彩色影像资料，采用人机交互目视解译方法，分析了福建省5年的土地利用动态特征。结果表明，土地利用的变化主要表现在林地、耕地的减少和草地、城镇建设用地的增加；5年间，城镇建设用地的年变化率最大，因此。加强林地保护，提高耕地质量，合理控制建设用地规模是切实可行的土地利用对策。     

INFLUENCES OF WATERSHED,RIPARIAN‐CORRIDOR,AND REACH‐SCALE CHARACTERISTICS ON AQUATIC BIOTA IN AGRICULTURAL WATERSHEDS1

ABSTRACT: Multivariate analyses and correlations revealed strong relations between watershed and riparian‐corridor land cover, and reach‐scale habitat versus fish and macroinvertebrate assemblages in 38 warmwater streams in eastern Wisconsin. Watersheds were dominated by agricultural use, and ranged in size from 9 to 71 km² Watershed land cover was summarized from satellite‐derived data for the area outside a 30‐m buffer. Riparian land cover was interpreted from digital orthophotos within 10‐, 10‐to 20‐, and 20‐to 30‐m buffers. Reach‐scale habitat, fish, and macroinvertebrates were collected in 1998 and biotic indices calculated. Correlations between land cover, habitat, and stream‐quality indicators revealed significant relations at the watershed, riparian‐corridor, and reach scales. At the watershed scale, fish diversity, intolerant fish and EPT species increased, and Hilsenhoff biotic index (HBI) decreased as percent forest increased. At the riparian‐corridor scale, EPT species decreased and HBI increased as riparian vegetation became more fragmented. For the reach, EPT species decreased with embeddedness. Multivariate analyses further indicated that riparian (percent agriculture, grassland, urban and forest, and fragmentation of vegetation), watershed (percent forest) and reach‐scale characteristics (embeddedness) were the most important variables influencing fish (IBI, density, diversity, number, and percent tolerant and insectivorous species) and macroinvertebrate (HBI and EPT) communities.     

URBAN IMPACTS ON PHYSICAL STREAM CONDITION: EFFECTS OF SPATIAL SCALE,CONNECTIVITY, AND LONGITUDINAL TRENDS1

ABSTRACT: An assessment of physical conditions in urban streams of the Puget Sound region, coupled with spatially explicit watershed characterizations, demonstrates the importance of spatial scale, drainage network connectivity, and longitudinal downstream trends when considering the effects of urbanization on streams. A rapid stream assessment technique and a multimetric index were used to describe the physical conditions of multiple reaches in four watersheds. Watersheds were characterized using geographic information system (GIS) derived landscape metrics that represent the magnitude of urbanization at three spatial scales and the connectivity of urban land. Physical conditions, as measured by the physical stream conditions index (PSCI), were best explained for the watersheds by two landscape metrics: quantity of intense and grassy urban land in the subwatershed and quantity of intense and grassy urban land within 500 m of the site (R²= 0.52, p > 0.0005). A multiple regression of PSCI with these metrics and an additional connectivity metric (proximity of a road crossing) provided the best model for the three urban watersheds (R²= 0.41, p > 0.0005). Analyses of longitudinal trends in PSCI within the three urban watersheds showed that conditions improved when a stream flowed through an intact riparian buffer with forest or wetland vegetation and without road crossings. Results demonstrate that information on spatial scale and patterns of urbanization is essential to understanding and successfully managing urban streams.     

Topographic Effects on Soil Organic Carbon in Louisiana Watersheds

Terrestrial carbon storage is influenced by a number of environmental factors, among which topographic and geomorphological features are of special significance. This study was designed to examine the relationships of soil organic carbon (SOC) density to various terrain parameters and watershed characteristics across Louisiana, USA. A polygon data set of 484 watersheds and 12 river drainage basins for Louisiana was used to form the landscape units. SOC densities were calculated for each soil map unit using the State Soil Geographic (STATSGO) database. Average drainage densities and average slopes at watershed and basin scales were quantified with the 1:24 K Digital Elevation Models (DEM) data, and the Louisiana hydrographic water features. Correlation and regression analyses were performed to determine relationships among drainage density, slope, elevation, and SOC. The study found an average watershed drainage density of 1.6 km/km² and an average watershed slope of 2.9 degrees in Louisiana. The results revealed that SOC density at both watershed and basin scales was closely related to drainage density, slope, and elevation. SOC density was positively correlated with watershed drainage density, but negatively correlated with watershed slope gradient and elevation. Regression models were developed for predicting SOC density at watershed and basin scales, obtaining regression coefficients (r ²) ranging from 0.43 to 0.83. The study showed that estimation of SOC at watershed and drainage basin scales combining DEM data can be a feasible approach to improve the understanding of the relationships among SOC, topographic, and geomorphological features.     

基于县域尺度的土地利用景观格局特征及其影响因素研究——以江苏省丰县为例

从类型水平、景观水平两个维度选取8个景观指数定量分析了江苏省丰县土地利用景观格局特征,探讨其影响因素.研究表明,耕地作为丰县景观基底,其景观结合度水平较高;耕地和园地景观的聚集度高、破碎度低、连接度高,而其他农用地和林地则相反;距离交通线路越近,斑块密度、景观多样性、均匀度等景观指数越大;经济较发达镇的斑块密度较大,景观分散度和分离度指数较高,聚集度指数较低;距县城越远,景观类型分布越集中且较单一.     

Modeling the relationships between land use and land cover on private lands in the Upper Midwest, USA

This paper presents an approach to modeling land-cover change as a function of land-use change. We argue that, in order to model the link between socio-economic change and changes in forest cover in a region that is experiencing residential and recreational development and agricultural abandonment, land-use and land-cover change need to be represented as separate processes. Forest-cover change is represented here using two transition probabilities that were calculated from Landsat imagery and that, taken together, describe a Markov transition matrix between forest and non-forest over a 10-year period. Using a three-date land-use data set, compiled and interpreted from digitized parcel boundaries, and scanned aerial photography for 136 sites (c. 2500 ha) sampled from the Upper Midwest, USA, we test functional relationships between forest-cover transition probabilities, standardized to represent changes over a decade, and land-use conditions and changes within sample sites. Regression models indicated that about 60% of the variation in the average forest-cover transition probabilities (i.e. from forest to non-forest and vice versa) can be predicted using three variables: amount of agricultural land use in a site; amount of developed land use; and the amount of area increasing in development. In further analysis, time lags were evaluated, showing that agricultural abandonment had a relatively strong time-lag effect but development did not. We demonstrate an approach to using forest-cover transition probabilities to develop spatially-constrained simulations of forest-cover change. Because the simulations are based on transition probabilities that are indexed to a particular time and place, the simulations are improved over previous applications of Markov transition models. This modeling approach can be used to predict forest-cover changes as a result of socio-economic change, by linking to models that predict land-use change on the basis of exogenous human-induced drivers.     

A Comprehensive Land-Use/Hydrological Modeling System for Scenario Simulations in the Elbow River Watershed, Alberta, Canada

The Elbow River watershed in Alberta covers an area of 1,238 km² and represents an important source of water for irrigation and municipal use. In addition to being located within the driest area of southern Canada, it is also subjected to considerable pressure for land development due to the rapid population growth in the City of Calgary. In this study, a comprehensive modeling system was developed to investigate the impact of past and future land-use changes on hydrological processes considering the complex surface–groundwater interactions existing in the watershed. Specifically, a spatially explicit land-use change model was coupled with MIKE SHE/MIKE 11, a distributed physically based catchment and channel flow model. Following a rigorous sensitivity analysis along with the calibration and validation of these models, four land-use change scenarios were simulated from 2010 to 2031: business as usual (BAU), new development concentrated within the Rocky View County (RV-LUC) and in Bragg Creek (BC-LUC), respectively, and development based on projected population growth (P-LUC). The simulation results reveal that the rapid urbanization and deforestation create an increase in overland flow, and a decrease in evapotranspiration (ET), baseflow, and infiltration mainly in the east sub-catchment of the watershed. The land-use scenarios affect the hydrology of the watershed differently. This study is the most comprehensive investigation of its nature done so far in the Elbow River watershed. The results obtained are in accordance with similar studies conducted in Canadian contexts. The proposed modeling system represents a unique and flexible framework for investigating a variety of water related sustainability issues.     

APPLICATION OF THE SNOWMELT RUNOFF MODEL USING MULTIPLE-PARAMETER LANDSCAPE ZONES ON THE TOWANDA CREEK BASIN,PENNSYLVANIA1

ABSTRACT: The Snowmelt Runoff Model (SRM) is designed to compute daily stream discharge using satellite snow cover data for a basin divided into elevation zones. For the Towanda Creek basin, a Pennsylvania watershed with relatively little relief, analysis of snow cover images revealed that both elevation and land use affected snow accumulation and melt on the landscape. The distribution of slope and aspect on the watershed was also considered; however, these landscape features were not well correlated with the available snow cover data. SRM streamflow predictions for 1990, 1993 and 1994 snowmelt seasons for the Towanda Creek basin using a combination of elevation and land use zones yielded more precise streamflow estimates than the use of standard elevation zones alone. The use of multiple-parameter zones worked best in non-rain-on-snow conditions such as in 1990 and 1994 seasons where melt was primarily driven by differences in solar radiation. For seasons with major rain-on-snow events such as 1993, only modest improvements were shown since melt was dominated by rainfall energy inputs, condensation and sensible heat convection. Availability of GIS coverages containing satellite snow cover data and other landscape attributes should permit similar reformulation of multiple-parameter watershed zones and improved SRM streamflow predictions on other basins.     

Boosted Regression Tree Models to Explain Watershed Nutrient Concentrations and Biological Condition

Boosted regression tree (BRT) models were developed to quantify the nonlinear relationships between landscape variables and nutrient concentrations in a mesoscale mixed land cover watershed during base‐flow conditions. Factors that affect instream biological components, based on the Index of Biotic Integrity (IBI), were also analyzed. Seasonal BRT models at two spatial scales (watershed and riparian buffered area [RBA]) for nitrite‐nitrate (NO₂‐NO₃), total Kjeldahl nitrogen, and total phosphorus (TP) and annual models for the IBI score were developed. Two primary factors — location within the watershed (i.e., geographic position, stream order, and distance to a downstream confluence) and percentage of urban land cover (both scales) — emerged as important predictor variables. Latitude and longitude interacted with other factors to explain the variability in summer NO₂‐NO₃ concentrations and IBI scores. BRT results also suggested that location might be associated with indicators of sources (e.g., land cover), runoff potential (e.g., soil and topographic factors), and processes not easily represented by spatial data indicators. Runoff indicators (e.g., Hydrological Soil Group D and Topographic Wetness Indices) explained a substantial portion of the variability in nutrient concentrations as did point sources for TP in the summer months. The results from our BRT approach can help prioritize areas for nutrient management in mixed‐use and heavily impacted watersheds.