Health and ecological risks associated with the use of mercury in gold mining are well known, with much recent attention focussed on contemporary small-scale artisanal mining. Legacy tailings from historical gold mining may also present ongoing risks, as the industry used large quantities of mercury with minimal environmental regulation to limit its discharge. This occurred in both alluvial (placer) mining and in processing auriferous ores. Analysis of historical data on mercury use in the mining industry in Victoria, Australia, indicates that at least 131 tonnes of elemental mercury were discharged into the environment as mine tailings between 1868–1888, with the total amount lost over the historic mining period likely to be much higher. The processing of pyritic ores also concentrated mercury losses in a small number of mining centres, including Bendigo, Ballarat, Castlemaine, Clunes, Maldon and Walhalla. This analysis provides a basis for further research needed to support improved management of legacy mine tailings.

Brazilian indigenous lands prevent the deforestation of the Amazon rainforest while protecting the land rights of indigenous peoples. However, they are at risk because they overlap with large areas of registered interest for mining. Indigenous lands have been in the spotlight of the pro-development wing of the parliament for decades, and the current president of Brazil, Jair Bolsonaro, promised that he would open up these territories for exploitation. Recently, bill PL191/2020 was released to downgrade the protection status of indigenous lands by regulating mining activities in these territories. Mining operations have an unavoidable socio-environmental impact on indigenous communities that is difficult to compensate. First, rapid demographic growth associated with the incoming migrant workforce often causes social disruption and threat indigenous societies. Moreover, sustained pollution related to mining procedures and accidental spills largely degrade the environment and imperil indigenous health. Finally, mining operations drive deforestation both within and beyond their operational boundaries. Mining is already an essential determinant of forest loss in the Amazon, where further deforestation may result in extended droughts with significant social and economic consequences. We conclude that, if mining operations were allowed in Brazilian indigenous lands, indigenous peoples would be imperiled along with regional and global climate and economies.

The deep-time dynamics of coupled socio-ecological systems at different spatial scales is viewed as a key framework to understand trends and mechanisms that have led to the Anthropocene. By integrating archeological and paleoenvironmental records, we test the hypothesis that Chilean societies progressively escalated their capacity to shape national biophysical systems as socio-cultural complexity and pressures on natural resources increased over the last three millennia. We demonstrate that Pre-Columbian societies intentionally transformed Chile’s northern and central regions by continuously adjusting socio-cultural practices and/or incorporating technologies that guaranteed resource access and social wealth. The fact that past human activities led to cumulative impacts on diverse biophysical processes, not only contradicts the notion of pristine pre-Industrial Revolution landscapes, but suggests that the Anthropocene derives from long-term processes that have operated uninterruptedly since Pre-Columbian times. Moreover, our synthesis suggests that most of present-day symptoms that describe the Anthropocene are rooted in pre-Columbian processes that scaled up in intensity over the last 3000 years, accelerating after the Spanish colonization and, more intensely, in recent decades. The most striking trend is the observed coevolution between the intensity of metallurgy and heavy-metal anthropogenic emissions. This entails that the Anthropocene cannot be viewed as a universal imprint of human actions that has arisen as an exclusive consequence of modern industrial societies. In the Chilean case, this phenomenon is intrinsically tied to historically and geographically diverse configurations in society-environment feedback relationships. Taken collectively with other case studies, the patterns revealed here could contribute to the discussion about how the Anthropocene is defined globally, in terms of chronology, stratigraphic markers and attributes. Furthermore, this deep-time narrative can potentially become a science-based instrument to shape better-informed discourses about the socio-environmental history in Chile. More importantly, however, this research provides crucial “baselines” to delineate safe operating spaces for future socio-ecological systems.

We use remote sensing to enhance the interpretation of the first baseline dataset of hydrologic, isotopic and hydrochemical variables spanning 620 km of the upper Marañón River, in Andean Peru, from the steep alpine canyons to the lower lying jungle. Remote, data-scarce river systems are under increased hydropower development pressure to meet rising energy demands. The upstream-downstream river continuum, which serves as a conduit for resource exchange across ecosystems, is at risk, potentially endangering the people, environments, and economies that rely on river resources. The Marañón River, one of the final free-flowing headwater connections between the Andes and the Amazon, is the subject of myriad large-scale hydropower proposals. Due to challenging access, environmental data are scarce in the upper Marañón, limiting our ability to do system-wide river basin planning. We capture key processes and transitions in the context of hydropower development. Two hydrologic regimes control the Marañón dry-season flow: in the higher-elevation upper reaches, a substantial baseflow is fed by groundwater recharged from wet season rains, in contrast to the lower reaches where the mainstem discharge is controlled by rain-fed tributaries that receive rain from lowland Amazon moisture systems. Sustainability of the upper corridor’s dry-season baseflow appears to be more highly connected to the massive natural storage capacity of extensive wetlands in the puna (alpine grasslands) than with cryospheric water inputs. The extent and conservation of puna ecosystems and glacier reservoirs may be interdependent, bringing to bear important conservation questions in the context of changing climate and land use in the region. More generally, this case study demonstrates an efficient combined remote sensing and field observation approach to address data scarcity across regional scales in mountain basins facing imminent rapid change.

Roots and associated microbes generate acid-forming CO 2 and organic acids and accelerate mineral weathering deep within Earth’s critical zone (CZ). At the Calhoun CZ Observatory in the USA’s Southern Piedmont, we tested the hypothesis that deforestation-induced deep root losses reduce root- and microbially-mediated weathering agents well below maximum root density (to 5 m), and impart land-use legacies even after ~70 y of forest regeneration. In forested plots, root density declined with depth to 200 cm; in cultivated plots, roots approached zero at depths >70 cm. Below 70 cm, root densities in old-growth forests averaged 2.1 times those in regenerating forests. Modeled root distributions suggest declines in density with depth were steepest in agricultural plots, and least severe in old-growth forests. Root densities influenced biogeochemical environments in multiple ways. Microbial community composition varied with land use from surface horizons to 500 cm; relative abundance of root-associated bacteria was greater in old-growth soils than in regenerating forests, particularly at 100–150 cm. At 500 cm in old-growth forests, salt-extractable organic C (EOC), an organic acid proxy, was 8.8 and 12.5 times that in regenerating forest and agricultural soils, respectively. The proportion of soil organic carbon comprised of EOC was greater in old-growth forests (20.0 ± 2.6%) compared to regenerating forests (2.1 ± 1.1) and agricultural soils (1.9 ± 0.9%). Between 20 and 500 cm, [EOC] increased more with root density in old-growth relative to regenerating forests. At 300 cm, in situ growing season [CO 2 ] was significantly greater in old-growth forests relative to regenerating forests and cultivated plots; at 300 and 500 cm, cultivated soil [CO 2 ] was significantly lower than in forests. Microbially-respired δ 13 C-CO 2 suggests that microbes may rely partially on crop residue even after ~70 y of forest regeneration. We assert that forest conversion to frequently disturbed ecosystems limits deep roots and reduces biotic generation of downward-propagating weathering agents.

Unconventional oil and gas exploration in the US has become a significant new source of atmospheric hydrocarbons. Field measurements and monitoring have been initiated to determine integral effects from this geographically dispersed source in and downwind of shale areas, driven mostly by concerns related to photochemical ozone production. The Texas Commission on Environmental Quality (TCEQ) deployed its first air quality monitor near the Eagle Ford shale in south Texas in summer 2013, followed by a more centrally located monitor in winter 2014/15. Here, we report on the latter monitor’s 2015 data, showing at times extraordinarily high levels of saturated hydrocarbons, similar to earlier findings in this area. Using hydrocarbon ratios, we establish that the dominant sources at this site appear to be oil and gas exploration. A non-negative matrix factorization (NMF) analysis revealed six consistent source factors, of which two were associated with pre-existing local sources from car traffic and industry, three with regional oil and gas exploration, and one with diesel emissions. The dominant source factors were associated with evaporative and fugitive emissions, and with flaring and (diesel-powered) compressor engine emissions. The former is a major source of saturated hydrocarbons while the latter is a major source of NO x and unsaturated hydrocarbons, confirming earlier findings. Due to the rural nature of the site, road traffic is a minor NO x source in this area, and the NMF results support inventory estimates showing oil and gas exploration to be the dominant regional source of NO x emissions. The NMF based source apportionment results also suggests that benzene levels in this rural area in 2015, while comparable to levels in Houston now, were probably three to five times lower before the shale boom.

Deltas are highly sensitive to erosion, flooding, and salinization with consequential agricultural productivity losses and out-migration, which is a preferred adaptive measure for the inhabitants of deltaic islands. This study investigates the associations between agricultural productivity decrease, household poverty and the probability of out-migration in the Indian Sundarban Delta (ISD). Using newly collected survey data from randomly selected households within the ISD, we analysed these relationships by means of descriptive statistics and regression modeling. Results suggest the significant positive association between a decrease in agricultural productivity and out-migration. The results further show that ceteris paribus, out-migration is negatively associated with household poverty, which is likely to be explained by the effect of remittances. The results yield important policy implications at the local level and can contribute to the progress towards sustainable livelihoods in these deltaic islands.

Salinity intrusion in coastal Bangladesh has serious population health implications, which are yet to be clearly understood. The study was undertaken through the ‘Assessing Health, Livelihoods, Ecosystem Services and Poverty Alleviation in Populous Deltas’ project in coastal Bangladesh. Drinking water salinity and blood pressure measurements were carried out during the household survey campaign. The study explored association among Socio-Ecological Systems (SESs), drinking water salinity and blood pressure. High blood pressure (prehypertension and hypertension) was found significantly associated with drinking water salinity. People exposed to slightly saline (1000–2000 mg/l) and moderately saline (≥2000 mg/l) concentration drinking water had respectively 17% (p < 0.1) and 42% (p < 0.05) higher chance of being hypertensive than those who consumed fresh water (<1000 mg/l). Women had 31% higher chance of being hypertensive than men. Also, respondents of 35 years and above were about 2.4 times more likely to be hypertensive compared to below 35 years age group. For the 35 years and above age group, both prehypertension and hypertension were found higher than national rural statistics (50.1%) for saline water categories (53.8% for slightly and 62.5% for moderate saline). For moderate salinity exposure, hypertension prevalence was found respectively 21%, 60% and 48% higher than national statistics (23.6%) in consecutive survey rounds among the respondents. Though there was small seasonal variation in drinking water salinity, however blood pressure showed an increasing trend and maximum during the dry season. Mean salinity and associated hypertension prevalence were found higher for deep aquifer (21.6%) compared to shallow aquifer (20.8%). Localized increase in soil and groundwater salinity was predicted over the study area. Shallow aquifer salinity increase was projected based on modelled output of soil salinity. Rather than uniform increase, there were localized extreme values. Deep aquifer salinity was also predicted to exhibit increasing trend over the period. Study findings and recommendations are suggested for immediate and planned intervention.

This study examines the morphodynamic response of a deltaic system to extreme weather events. The Wax Lake Delta (WLD) in Louisiana, USA, is used to illustrate the impact of extreme events (hurricanes) on a river-dominated deltaic system. Simulations using the open source Delft3D model reveal that Hurricane Rita, which made landfall 120 km to the west of WLD as a Category 3 storm in 2005, caused erosion on the right side and deposition on the left side of the hurricane eye track on the continental shelf line (water depth 10 m to 50 m). Erosion over a wide area occurred both on the continental shelf line and in coastal areas when the hurricane moved onshore, while deposition occurred along the Gulf coastline (water depth < 5 m) when storm surge water moved back offshore. The numerical model estimated that Hurricane Rita’s storm surge reached 2.5 m, with maximum currents of 2.0 m s –1 , and wave heights of 1.4 m on the WLD. The northwestern-directed flow and waves induced shear stresses, caused erosion on the eastern banks of the deltaic islands and deposition in channels located west of these islands. In total, Hurricane Rita eroded more than 500,000 m 3 of sediments on the WLD area. Including waves in the analysis resulted in doubling the amount of erosion in the study area, comparing to the wave-excluding scenario. The exclusion of fluvial input caused minor changes in deltaic morphology during the event. Vegetation cover was represented as rigid rods in the model which add extra source terms for drag and turbulence to influence the momentum and turbulence equations. Vegetation slowed down the floodwater propagation and decreased flow velocity on the islands, leading to a 47% reduction in the total amount of erosion. Morphodynamic impact of the hurricane track relative to the delta was explored. Simulations indicate that the original track of Hurricane Rita (landfall 120 km west of the WLD) produced twice as much erosion and deposition at the delta compared to a hurricane of a similar intensity that made landfall directly on the delta. This demonstrates that the wetlands located on the right side of a hurricane track experience more significant morphological changes than areas located directly on the hurricane track.

The population of the Ganges Brahmaputra Meghna (GBM) delta is highly vulnerable to food insecurity and malnutrition due to the specific environmental, climatic and human development factors affecting agricultural production and fisheries. To better understand the impacts of climate and environmental change on food security and nutrition in this delta, this study combines spatially explicit data from the 2007 and 2011 Bangladesh Demographic and Health Survey (BDHS) with a standard satellite remotely sensed vegetation greenness index (Normalised Difference Vegetation Index, NDVI), used as a proxy for rice production. The strength of association between NDVI and child nutrition in this tropical mega-delta were tested, showing correlations between two widely used indicators of child malnutrition; stunting and wasting, and deviations from a 10 year mean NDVI (anomalies) for rice crop growing seasons – regarded as critical to individual children’s early lives. For children surveyed in 2007 we found that the likelihood of being stunted decreased with increased NDVI as a measure of food production. Similarly, for children surveyed in 2011, the likelihood of being wasted reduced with increased NDVI. However, regression results for stunting in 2011 and wasting in 2007 were not statistically significant. Our findings suggest that NDVI can be regarded as indicative of climatic variability and periods of low food production but is only partly successful as an indicator of climate related impacts on child nutrition in the GBM delta. Furthermore, our study highlights some of the uncertainties and challenges with linking environmental indicators such as the NDVI with household survey data across spatial and temporal scales.

The Mediterranean basin (including the Black Sea) is characterized by a plethora of deltas that have developed in a wave-influenced setting. Many of these deltas are sourced in sediments by river catchments that have been variably dammed. The vulnerability status of a selection of ten deltas subject to different levels of reduction in fluvial sediment supply following damming was analysed by quantifying changes in delta protrusion area and protrusion angle over the last 30 years. The rationale for choosing these two metrics, which do not require tricky calculations of longshore bedload transport volumes and river ‘influence’, is that as sediment supply wanes, increasing relative efficiency of waves leads to longshore redistribution of reworked sediments and progressive ‘flattening’ of the delta protrusion. The results show that eight of the ten deltas (Nile, Rhône, Ebro, Ceyhan, Arno, Ombrone, Moulouya, Medjerda) are in erosion, whereas two (Danube, Po) show stability, but the statistical relationship between change in delta protrusion area and sediment flux reduction is poor, thus suggesting that the role of dams in causing delta shoreline erosion may have been over-estimated. But this poor relationship could also be due to a long temporal lag between dam construction and bedload removal and transport to the coast downstream of dams, and, where the delta protrusion is being eroded, to bedload trapping by shoreline engineering structures and by elongating delta-flank spits. Other potential influential factors in shoreline change include subsidence, sea-level rise, storminess, exceptional river floods, and managed sediment releases downstream of dams. A longer observation period and high-resolution sediment-budget studies will be necessary to determine more definitively to which extent continued trapping of sediment behind dams will impact overall delta stability in the Mediterranean and Black Seas. Mitigation of delta erosion is likely to become costlier under continued sediment starvation and sea-level rise.

Deltaic environments are often densely populated with high socio-economic values, and thus are hotspots of climatic, environmental and anthropogenic change. Large scale engineered structures, such as dike systems, have played an important role in shaping both environmental and socio-economic conditions in deltas, with such interventions more likely where there is a high population and a wealthy economy. Engineered interventions interact with the morphological evolution of the delta, reducing or removing sedimentation and accelerating subsidence, increasing the consequences of flooding and necessitating further adaptation. They also encourage further development, reinforcing this feedback. Thus, in these cases, the deltaic landscape and associated livelihoods can be considered to be the result of a coevolution process between natural delta processes and human engineered interventions. This paper explores this hypothesis. It analyses the history of large scale engineering interventions and their implications in five representative, large, populated deltas across the globe (Ganges-Brahmaputra-Meghna, Yangtze, Rhine-Meuse-Scheldt, Mekong and Nile). The results demonstrate coevolution has occurred and indicate that the response type and the management approach of these engineered structures have significant implications for future delta development. To understand and manage unintended consequences and the development of lock-in trajectories in deltas, a systematic understanding of delta development, including these coevolution processes is essential.

Recurrent, tidally driven coastal flooding is one of the most visible signs of sea level rise. Recent studies have shown that such flooding will become more frequent and extensive as sea level continues to rise, potentially altering the landscape and livability of coastal communities decades before sea level rise causes coastal land to be permanently inundated. In this study, we identify US communities that will face effective inundation—defined as having 10% or more of livable land area flooded at least 26 times per year—with three localized sea level rise scenarios based on projections for the 3 rd US National Climate Assessment. We present these results in a new, online interactive tool that allows users to explore when and how effective inundation will impact their communities. In addition, we identify communities facing effective inundation within the next 30 years that contain areas of high socioeconomic vulnerability today using a previously published vulnerability index. With the Intermediate-High and Highest sea level rise scenarios, 489 and 668 communities, respectively, would face effective inundation by the year 2100. With these two scenarios, more than half of communities facing effective inundation by 2045 contain areas of current high socioeconomic vulnerability. These results highlight the timeframes that US coastal communities have to respond to disruptive future inundation. The results also underscore the importance of limiting future warming and sea level rise: under the Intermediate-Low scenario, used as a proxy for sea level rise under the Paris Climate Agreement, 199 fewer communities would be effectively inundated by 2100.

Wetlands occupy a position in the landscape that makes them vulnerable to the effects of current land use and the legacies of past land use. Many wetlands in agricultural regions like the North American Midwest are strongly affected by elevated nutrient inputs as well as high rates of invasion by the hybrid cattail Typha × glauca . These two stressors also exacerbate each other: increased nutrients increase invasion success, and invasions increase nutrient retention and nutrient loads in the wetland. This interaction could create a positive feedback that would inhibit efforts to manage and control invasions, but little is known about the effects of past or present nutrient inputs on wetland invasive plant management. We augmented a previously-published community-ecosystem model (M ONDRIAN ) to simulate the most common invasive plant management tools: burning, mowing, and herbicide application. We then simulated different management strategies and 3 different durations in low and high nutrient input conditions, and found that the most effective management strategy and duration depends strongly on the amount of nutrients entering the wetland. In high-nutrient wetlands where invasions were most successful, a combination of herbicide and fire was most effective at reducing invasion. However, in low-nutrient wetlands this approach did little to reduce invasion. A longer treatment duration (6 years) was generally better than a 1-year treatment in high-nutrient wetlands, but was generally worse than the 1-year treatment in low-nutrient wetlands. At the ecosystem level, we found that management effects were relatively modest: there was little effect of management on ecosystem C storage, and while some management strategies decreased wetland nitrogen retention, this effect was transient and disappeared shortly after management ceased. Our results suggest that considering nutrient inputs in invaded wetlands can inform and improve management, and reducing nutrient inputs is an important component of an effective management strategy.

As climate-science graduate students at the University of Washington, we had the opportunity to engage in a political process focused on implementing legislation to reduce greenhouse gas emissions in Washington State. Our insights gained from this rare, first-hand, experience may be particularly relevant to other climate scientists. We argue that inflexible research goals within the United States climate-science community limit the relevance of the knowledge our community creates. The mismatch between climate-science research and the information needs of policy makers, while widely acknowledged in certain domains, has yet to be fully appreciated within many earth science disciplines. Broadening the climate-science training of graduate students to include education on the uses of climate information outside of academic settings would both inform and motivate new research directions, and engender validation of non-traditional research within disciplinary cultures.

The ubiquitous marine bacterium Vibrio parahaemolyticus is a leading cause of illness associated with seafood consumption. The emergence of two genetically distinct ecotypes (ST3 and ST36) has led to an alarming increase in the size and frequency of disease outbreaks. We conducted a genomic comparison of 30 V. parahaemolyticus genomes that represent a diverse collection of 15 genetically distinct ecotypes, including newly sequenced representatives of ST3 and ST36, isolated from both clinical and environmental sources. A multistep evolutionary analysis showed that genes associated with sensing and responding to environmental stimuli have evolved under positive selection, identifying examples of convergent evolution between ST3 and ST36. A comparison of predicted proteomes indicated that ST3 and ST36 ecotypes laterally acquired tens of novel genes associated with a variety of functions including dormancy, homeostasis and membrane transport. Genes identified in this study play an apparent role in environmental fitness and may confer cross protection against stressors encountered in the human host. Together, these results show the evolution of stress response is an important genetic mechanism correlated with the recent emergence of the ST3 and ST36 ecotypes.

The cycling of mercury between ocean and atmosphere is an important part of the global Hg cycle. Here we study the regional contribution of the air-sea exchange in the North- and Baltic Sea region. We use a newly developed coupled regional chemistry transport modeling (CTM) system to determine the flux between atmosphere and ocean based on the meteorological model COSMO-CLM, the ocean-ecosystem model ECOSMO, the atmospheric CTM CMAQ and a newly developed module for mercury partitioning and speciation in the ocean (MECOSMO). The model was evaluated using atmospheric observations of gaseous elemental mercury (GEM), surface concentrations of dissolved gaseous mercury (DGM), and air-sea flux (ASF) calculations based on observations made on seven cruises in the western and central Baltic Sea and three cruises in the North Sea performed between 1991 and 2006. It was shown that the model is in good agreement with observations: DGM (Normalized Mean Bias NMB=-0.27 N=413), ASF (NMB=-0.32, N=413), GEM (NMB=0.07, N=2359). Generally, the model was able to reproduce the seasonal DGM cycle with the best agreement during winter and autumn (NMB Winter =-0.26, NMB Spring =-0.41, NMB Summer =-0.29, NMB Autumn =-0.03). The modelled mercury evasion from the Baltic Sea ranged from 3400 to 4000 kg/a for the simulation period 1994–2007 which is on the lower end of previous estimates. Modelled atmospheric deposition, river inflow and air-sea exchange lead to an annual net Hg accumulation in the Baltic Sea of 500 to 1000 kg/a. For the North Sea the model calculates an annual mercury flux into the atmosphere between 5700 and 6000 kg/a. The mercury flux from the ocean influenced coastal atmospheric mercury concentrations. Running CMAQ coupled with the ocean model lead to better agreement with GEM observations. Directly at the coast GEM concentrations could be increased by up to 10% on annual average and observed peaks could be reproduced much better. At stations 100km downwind the impact was still observable but reduced to 1–3%.

Atmosphere–surface exchange of mercury, although a critical component of its global cycle, is currently poorly constrained. Here we use the GEOS-Chem chemical transport model to interpret atmospheric Hg 0 (gaseous elemental mercury) data collected during the 2013 summer Nitrogen, Oxidants, Mercury and Aerosol Distributions, Sources and Sinks (NOMADSS) aircraft campaign as well as ground- and ship-based observations in terms of their constraints on the atmosphere–surface exchange of Hg 0 over eastern North America. Model–observation comparison suggests that the Northwest Atlantic may be a net source of Hg 0 , with high evasion fluxes in summer (our best sensitivity simulation shows an average oceanic Hg 0 flux of 3.3 ng m -2 h -1 over the Northwest Atlantic), while the terrestrial ecosystem in the summer of the eastern United States is likely a net sink of Hg 0 (our best sensitivity simulation shows an average terrestrial Hg 0 flux of -0.6 ng m -2 h -1 over the eastern United States). The inferred high Hg 0 fluxes from the Northwest Atlantic may result from high wet deposition fluxes of oxidized Hg, which are in turn related to high precipitation rates in this region. We also find that increasing simulated terrestrial fluxes of Hg 0 in spring compared to other seasons can better reproduce observed seasonal variability of Hg 0 concentration at ground-based sites in eastern North America.

I reassess a recent analysis of uncertainty in estimates of nitrogen export from stormwater control measures, using structured expert judgment, which concluded that nitrogen export from a watershed in the Piedmont physiographic province of the Chesapeake Bay basin was an order of magnitude greater than from a watershed in the adjacent the Coastal Plain province. Re-analysis of expert responses suggests that hydrographic measurement error is a likely large source of uncertainty in N export from one of the watersheds. Mass-balance estimates of impervious runoff into stormwater drainage systems suggest that nitrogen export from the Coastal Plain watershed is an order of magnitude larger than estimated. This analysis highlights the importance of stormwater drainage infrastructure in driving the hydrology of streams in urban catchments by quarantining impervious runoff from watershed soils.

We reply to a comment on our recent structured expert judgment analysis of stormwater nitrogen retention in suburban watersheds. Low relief, permeable soils, a dynamic stream channel, and subsurface flows characterize many lowland Coastal Plain watersheds. These features result in unique catchment hydrology, limit the precision of streamflow measurements, and challenge the assumptions for calculating runoff from rainfall and catchment area. We reiterate that the paucity of high-resolution nitrogen loading data for Chesapeake Bay watersheds warrants greater investment in long-term empirical studies of suburban watershed nutrient budgets for this region.