Introduction

Aquatic scientists have long noted major differences in nutrient control of net primary productivity across different types of ecosystems, with phosphorus (P) limitation in many lakes (Schindler 1974; Vollenweider 1975), nitrogen (N) limitation in many estuaries and coastal marine ecosystems (Ryther 1954), and co-limitation by N and P in many open ocean systems such as the subtropical gyres (Redfield 1958). Exceptions to these generalities exist, and some coastal marine ecosystems are limited by P rather than N while some lakes are limited by N (Howarth 1988; Bergström and Jansson 2006; Elser et al. 2009; Hayes et al. 2018). Time scale matters, and ecosystems can switch between N and P limitation seasonally, including both lakes and estuaries (Fisher et al. 1992; Dolman et al. 2016).

The topic of nutrient limitation in aquatic ecosystems is controversial, with long-standing debates about how one best determines which nutrient is most limiting. Our goal in this paper is not to further debate the evidence by which one determines which nutrient is more limiting, which we have previously discussed in some detail (Howarth 1988; Howarth and Marino 2006); see also NRC (2000). We agree with the suggestion of Moss et al. (2013) that limitation by a single nutrient is more likely in human-impacted, more nutrient-enriched systems, and that co-limitation by N and P is more likely in oligotrophic ecosystems (see also Begon et al. 2014). In our view, P limitation is more prevalent than N limitation in mesotrophic and eutrophic lakes, and N limitation is more prevalent in temperate-zone estuaries. Our personal research experience informs this belief, including our long-term study of a coastal marine lagoon that is clearly N limited (Howarth et al. 2014) and research on a lake that is clearly P limited (Peterson et al. 1974; Roberts and Howarth 2006). The coastal system, West Falmouth Harbor on Cape Cod, MA, is characterized by low dissolved inorganic N concentrations (generally < 1 µM) and moderate concentrations of soluble reactive P (0.25–1 µM; Hayn et al. 2014; Howarth et al. 2014). The lake, Cayuga Lake in the Finger Lakes region of central NY, is characterized by similar concentrations of soluble reactive P (SRP) but concentrations of dissolved inorganic N that are more than two magnitudes of order greater (mean value of 380 µM; Halfman 2017). Although the molar N:P ratio of dissolved inorganic nutrient concentrations is an imperfect indicator of nutrient limitation, particularly when near the Redfield ratio of 16:1 or in ecosystems where inorganic N and P are extremely low, the ratio near 1:1 in West Falmouth Harbor and the average ratio of more than 1000:1 in Cayuga Lake clearly indicate N and P limitation on primary production, respectively (NRC 2000).

The focus of this paper is to examine the role of external nutrient inputs to aquatic ecosystems as one of the factors that may tend to make primary production more likely to be limited by N, P, or both. We and others have written extensively about the role of biogeochemical processes such as N fixation, denitrification, and phosphorus sorption and desorption as factors that can lead more toward limitation of N or P in an aquatic ecosystem (see for example Howarth and Marino 2006 and references therein). However, the ratio of N:P for the external nutrient inputs to the ecosystem sets the stage against which these internal biogeochemical processes can act. For example, N fixation by planktonic cyanobacteria can help maintain P limitation in lakes when the N deficit is relatively small (Schindler 1977; Flett et al. 1980; Schindler et al. 2008). But in lakes where the N:P ratio of external inputs is sufficiently low, N fixation may not be able to make up the N deficit compared to P availability, and N limitation can occur (Moss et al. 2013; Grantz et al. 2014; Hayes et al. 2018). Significant planktonic N fixation is largely absent from saline estuaries, which is one of the reasons we believe coastal marine ecosystems are more likely than lakes to be N limited (Howarth et al. 1988; Marino et al. 2002; Marino and Howarth 2016).

For the surface waters of much of the world’s oceans away from shore and continental shelves, including the subtropical gyres that cover half of the area of the oceans and 35% of the total surface of the Earth, the largest source of nutrient inputs by far is mixing of deeper ocean waters across the thermocline (Michaels et al. 1996). The mixing rate is low, but the nutrient content of this water is quite high, so the input to the surface oceans is substantial. As we present below, the ratio of N to P in these deep waters is quite close to the Redfield ratio (Karl et al. 1997; Michaels et al. 1996). This leads to co-limitation of net primary productivity by both nutrients in the subtropical gyres, consistent with the thought that co-limitation by N and P may be the norm in oligotrophic ecosystems (Begon et al. 2014). This balance between N and P cycling in open ocean waters is at the center of the Redfield (1958) concept.

Most lakes receive the vast majority of their nutrient inputs from their watersheds, although atmospheric inputs of both N and P can be important, particularly in lakes with relatively pristine landscapes and consequently small rates of watershed loading. Estuaries and coastal marine ecosystems also receive substantial inputs of nutrients from watersheds and atmospheric deposition, but they additionally receive inputs of nutrients from ocean sources (Nixon et al. 1996; Howarth et al. 2011). This major difference in the external nutrient sources between lakes and coastal marine systems often escapes attention and is one of the major foci of this paper. In this paper, we first explore the average inputs of nutrients to aquatic ecosystems from their watersheds in the context of the United States, using a mass balance modeling approach at the regional scale, and then we examine dissolved N and P data from across the global oceans as sources of nutrients to the coastal zone.

Since this paper is an invited contribution for the special issue of Biogeochemistry celebrating more than 35 years of the journal, we note that the watershed nutrient inputs we discuss next are based on our 1996 paper published in Biogeochemistry as part of another earlier special issue (Howarth et al. 1996). This paper won the John Martin ASLO award in 2018 in recognition of being a highly influential paper. The 1996 special issue focused on papers developing an N budget for the North Atlantic Ocean as part of the International SCOPE Nitrogen Project. Other papers from that issue (Nixon et al. 1996; Michaels et al. 1996; Seitzinger and Giblin 1996) are also central to the development of the logic behind this current paper.

Nutrient inputs from the landscape

The rate of nutrient loading from the landscape to lakes and estuaries has been estimated for thousands of systems globally. However, there is no uniformly accepted methodology for estimating N and P loading (NRC 2000), complicating comparisons across systems. Here, we build on our previous work and use a regional nutrient accounting approach to estimate the average N and P loads from watersheds across the United States by region and over time, presenting 5-year averages from 1987 to 2012. The approach is based on the Net Anthropogenic Nitrogen Inputs (NANI) and Net Anthropogenic Phosphorus Inputs (NAPI) to the landscape (Howarth et al. 1996, 2006, 2012; Hong et al. 2011, 2012, 2013; Han et al. 2011).

NANI is the sum of N inputs to a region from synthetic N fertilizer, N fixation by agricultural crops, atmospheric deposition of N, and the N in food and feed for humans and livestock. If there is a net export of N in food and feed, this is treated as a negative input and is subtracted from the other inputs (Howarth et al. 1996; Hong et al. 2011). The approach for NAPI is very similar, but the net import of P in detergents and other non-food uses is included, and atmospheric deposition of P is assumed to be negligible (Han et al. 2011; Hong et al. 2012). We present data for several large regions of the United States based on county-scale data (Fig. 1a). As of the writing of this manuscript, county-scale data for fertilizer sales after 2012 are still not available from the USDA agricultural census (Slater, AAPFCO, pers. comm.).

Fig. 1
figure 1figure 1

a Regions of the United States for which we report data on NANI and NAPI. b NANI for each major region of the United States over time from 1987 to 2017. c NAPI for each major region of the United States over time from 1987 to 2017. d Molar ratio of nitrogen to phosphorus exported from each major region of the United States over time from 1987 to 2017. The Redfield ratio of 16:1 is shown as a dashed red line for context

The input terms for NAPI are as described in Swaney and Howarth (2019). For NANI, the input terms are as described in Swaney et al. (2018) except for N deposition. In our previously published work, atmospheric N deposition estimates for NANI were based on the U.S. Environmental Protection Agency CMAQ model. Here, we use updated estimates of atmospheric N deposition from a hybrid approach now favored by the US EPA (Schwede and Lear 2014), combining CMAQ output with National Atmospheric Deposition Program observational data to make deposition estimates (NADP 2018). The refinement to the approach shows lower values of NOy deposition for 2012 than our previous estimates (Howarth et al. 2012; Hong et al. 2013), consistent with observed declines in emissions (Butler et al. 2011).

Net Anthropogenic Nitrogen Inputs per watershed area are the highest in the Mid-Atlantic region but are also high in California, the Great Lakes area, the Mississippi River basin, and the South Atlantic/Gulf Coast (Fig. 1b). Values are considerably lower in New England, the Pacific Northwest, and the Texas Gulf coast area. In all of the regions, changes in NANI over time have been relatively small, sometimes increasing and sometimes decreasing. Net Anthropogenic Phosphorus Inputs per area of watershed are highest in the Mid-Atlantic and South Atlantic/Gulf Coast regions, and are far lower in most other regions, particularly New England and the Pacific Northwest (Fig. 1c). As with NANI, time trends for NAPI over time have been relatively small, with the exception of the Great Lakes region which saw a steady decline from 1987 to 2012.

To examine the ratio of total N to total P likely exported from each region to aquatic ecosystems, we first calculate the portions of NANI and NAPI that move downstream. For NANI, an average of 24% is exported from large watersheds as total N, based on a data set for 154 watersheds in the United States, Europe, and the United Kingdom (Howarth et al. 2012; Hong et al. 2013). Note that this is virtually identical to the 26% export of NANI that we estimated from the analysis of large regions in our original NANI paper (Howarth et al. 1996). While there is some variation in the percentage of NANI that is exported across these watersheds, it is relatively small. Surprisingly, the variation does not seem tied to the specific N inputs, the intensity of agriculture, topography, soil types, or other such factors (Howarth et al. 1996, 2012; Hong et al. 2013). Most of the differences in export across regions are related to climate, with greater percentage exports in regions where precipitation and river discharge per area are higher. For simplicity in this paper, we use the average export of 24%. Most of the N that is not exported downstream is denitrified or stored in the watershed (van Breemen et al. 2002).

For NAPI, we assume that 5.9% is exported downstream as total P. This is the average observed for 57 moderately large watersheds in the United States (Zhang et al. 2015). Somewhat lower export percentages are observed in China (2.9%) and Europe (3.9%), and Zhang et al. (2015) report an overall average export of 3.4% for 158 watersheds in the United States, China, and Europe. For this paper, we are using the 5.9% value since we are considering only regions in the United States. We recognize that this may bias our P export numbers towards the high end.

Using NANI and NAPI export percentages of 24% and 5.9% respectively, we can calculate the molar N:P ratio of nutrients leaving each region of the United States over time (Fig. 1d). In all regions and at all times since 1987 this N:P ratio is high, well above the Redfield ratio of 16:1. As we state in the paragraph above, our P export values may be biased high, and if that is the case and actual P fluxes are less than we estimate, the N:P ratios shown in Fig. 1d would be even higher. The N:P ratios of exports from the landscape are particularly high in the Great Lakes, the Mississippi River basin, New England, and the Pacific Northwest, ranging from somewhat > 50:1 to well over 100:1. The N:P ratio is lowest in the Texas Gulf region. Changes in the N:P ratio of exported nutrients over time have been small in most regions, but the Great Lakes region saw a large increase in this N:P ratio from 1987 to 2012 (Fig. 1d), largely driven by the decrease in NAPI in this region over time (Fig. 1c).

There is some debate over when and if nutrients stored in the landscape (“legacy nutrients”) are ultimately released, thereby enriching receiving waters, and how significant this is as a proportion of total nutrient loading. Some modeling efforts suggest current nutrient loads are due in part to release of legacy nutrients (Van Meter et al. 2017; McCrackin et al. 2018). However, these legacy nutrient sources are ultimately the result of excess NANI and NAPI storage in the landscape, and are not an additional input to the system. Van Meter et al. (2017) suggest that current nutrient exports have a larger proportion of legacy N than legacy P, indicating that using the long-term average export proportions as in the calculation above would produce a conservative estimate of N relative to P, and the N:P ratio could be somewhat higher.

It is important to reiterate that NANI and NAPI have several input terms and generally are not driven by any particular term such as fertilizer use (Howarth et al. 1996, 2012; Hong et al. 2011, 2012). This is illustrated by contrasting two different regions, New England and the Mississippi River basin. For NANI in New England, the biggest input terms are atmospheric deposition and net import of food and animal feeds, with fertilizer use and N fixation by agricultural crops being considerably smaller (Fig. 2a). On the other hand, the largest N inputs in the Mississippi River basin are from fertilizer use and agricultural N fixation. Further, the food and animal feeds term is negative in the Mississippi River basin, indicating a net export rather than an import, as in New England. In fact, more N leaves the Mississippi River basin via export of food and feed than flows down the Mississippi River (Fig. 2a). In both New England and the Mississippi River basin, it is NANI that determines riverine N flows, not any single input term. For NAPI, fertilizer is the largest input for both New England and the Mississippi River basin, followed by import in food and feeds (Fig. 2b). While the import of P in food and feeds in the Mississippi River basin is small, it is a positive value on average, in contrast to the net export of N from the basin. Mean values of the input terms for NANI and NAPI for the 1987–2012 period for each of the regions are shown in tables in on-line only supplemental material.

Fig. 2
figure 2

a Comparison on input terms for NANI in New England (bottom) and the Mississippi River basin (top) in kg N per square kilometer per year. Also shown are exports of N in water flowing to coastal ecosystems based on 24% of NANI. See text. b Comparison on input terms for NAPI in New England (bottom) and the Mississippi River basin (top) in kg P per square kilometer per year. Also shown are exports of P in water flowing to coastal ecosystems based on 5.9% of NAPI. See text

The regional-scale watershed fluxes we estimate are for total N and total P. Not all of this N and P will be biologically available to primary producer organisms in recipient ecosystems. Of the organic N and P, some will be mineralized to readily available forms on ecologically meaningful time scales, but some will simply be buried and not become available, particularly for total N (Howarth et al. 1996; Nixon et al. 1996). It is important to note that as human activity increases overall nutrient loads, inorganic N and P increasingly come to dominate total N and P fluxes (Meybeck 1982; Howarth et al. 1996; Seitzinger et al. 2005). All of the inorganic N that enters an aquatic ecosystem remains biologically available, but for inorganic P, much can be adsorbed to soil and sediment particles. Some of this may become available in recipient ecosystems but some will not. Adsorbed P is much more likely to desorb and become biologically available in estuaries and coastal marine ecosystems than in lakes, due to higher ionic strength and more ions competing for sorption sites as salinity increases (Froelich 1988; Howarth et al. 1995).

Influence of watershed loads on N vs P limitation

For lakes, those that receive nutrient loads from their watershed with an N:P molar ratio for total nutrients above the relative N and P requirements of phytoplankton would be expected to be P-limited, unless biogeochemical processes within the lake decrease the availability of N relative to P, lowering the N:P ratio for available nutrients below that needed by the primary producers. Of course, processes within the lake can lower the ratio of available N:P, for instance if denitrification rates are greater than rates of N fixation. The rates of biogeochemical processes that control N and P availabilities vary as a function of lake size, water residence time, and other factors (Moss et al. 2013; They et al. 2017). Nonetheless, the watershed nutrient inputs set the stage for the relative availability of N and P that support primary production.

What are the relative requirements of N and P for algae and cyanobacteria? Redfield (1958) originally proposed that marine phytoplankton need N and P in a 16:1 molar ratio. There is some plasticity in this ratio for N:P needs across different primary producers, both in the oceans and in lakes (Hecky and Kilham 1988; Karl et al. 1997). Limnologists often use the N:P ratio of seston to infer whether a lake is more limited by N or P (Hecky and Kilham 1988; Dolman et al. 2016), although it is difficult to tease apart how much of the sestonic particulate N and P is in living phytoplankton vs detritus, including that of terrestrial origin. Nonetheless, the N:P ratio of seston in lakes, while generally > 16:1, is low compared to the loading ratios we show in Fig. 1d, ranging from 17:1 to 39:1 in a variety of lakes, with a mean of 23:1 for small lakes and 26:1 for large lakes (Sterner et al. 2008).

The high N:P loading ratios shown in Fig. 1d would be expected to lead toward a tendency for P limitation in the receiving waterbodies, particularly considering that much of the particle-adsorbed P may not become available in a lake. This is consistent with the widespread prevalence of P limitation of net primary productivity across many North American lakes (Vollenweider 1975). We want to emphasize, though, that lakes undoubtedly vary greatly in the N:P ratio of their nutrient inputs, and no individual lake receives the average N and P loads indicated by this regional approach based on NANI and NAPI. Combined with biogeochemical processes within lakes that can alter N and P availability, it should come as no surprise that many lakes are N limited. Lakes in catchments that are totally forested will have lower total loads of both N and P, with no inputs from agricultural sources; in regions where atmospheric N deposition is low, lakes in forested catchments are often oligotrophic and N limited or co-limited by N and P, but high deposition rates over time can shift the N:P stoichiometry, leading to P limitation (Bergström and Jansson 2006; Elser et al. 2009). Note further than even lakes in largely forested catchments may have inputs of both N and P related to the food and feed term of NANI and NAPI, for instance due to nutrient inputs from onsite wastewater treatment systems from nearby homes.

This same construct of the N:P ratio for external nutrient inputs setting the stage for limitation of primary production by N vs P also should apply to estuaries and coastal marine ecosystems: unless internal biogeochemical processes are sufficient to compensate, nutrient inputs from the watershed with an N:P ratio substantially above the needs of the phytoplankton would be expected to lead to P limitation. The N:P ratios shown in Fig. 1d are far above the Redfield ratio of 16:1 or the seston N:P ratios commonly found in marine ecosystems, which tend to be below 20:1 in coastal waters and productive areas and range between 15:1 and 28:1 in oligotrophic ocean systems (Karl et al. 1997; Sterner et al. 2008; Martiny et al. 2013). Why then is N limitation commonly observed in most coastal marine ecosystems, in sharp contrast to so many lakes (Howarth and Marino 2006)? Part of the answer undoubtedly lies in the differences between the relative influence of the biogeochemical processes that occur within coastal ecosystems compared to lakes. These differences include a lack of planktonic N fixation in saline estuaries (Howarth et al. 1988; Marino et al. 2002; Marino and Howarth 2016) and the desorption of P from particles as salinity increases in an estuary (Froelich 1988; Howarth et al. 1995). Another important part of the answer, although one less discussed in the literature, is the N:P ratio of external inputs to the ecosystem: many estuaries and other coastal marine systems receive a significant portion of their nutrient inputs from ocean sources, not just from their watersheds.

Nutrient inputs from ocean sources

In many cases, oceanic inputs to estuaries are not measured and are not included in nutrient budgets (as discussed by Nixon et al. 1996). However, just as salt from coastal ocean systems is mixed and advected up into an estuary, nutrients from these coastal waters are also transported. Nutrient inputs from the ocean can be substantial, particularly for P, as can be seen for coastal systems with datasets rich enough to make good estimates. For example, an estimated 24% of the P inputs to Chesapeake Bay and 36% of P inputs to the Potomac River estuary come from ocean sources, with virtually no input of N from the ocean to these systems (Boynton et al. 1995). For both the Changjiang River estuary (Li et al. 2011) and the heavily N-loaded West Falmouth Harbor (Hayn et al. 2014; Howarth et al. 2014), P inputs from the ocean source are sufficient to maintain N limitation. In both of these systems, P inputs were proportionally greater than N inputs relative to the Redfield ratio of 16:1.

Figure 3 shows molar N:P ratios for inorganic nutrients (as nitrate + nitrite and SRP) in the deep waters below the thermocline for the world’s oceans from the World Ocean Atlas 2018 (WOA18; Garcia et al. 2018). Values represent observations from the deepest bins within each 1° gridded cell. We present data on inorganic N and SRP, rather than total N and P as we did for the watershed loading estimates, simply because inorganic data are common and total N and P data rare in oceanographic studies and data bases. Note, though, that in deep ocean waters, concentrations of inorganic N and SRP are very high and greatly dominate both the total N and total P pools (Michaels et al. 1996). For most of the world’s oceans, bottom waters have molar N:P ratios near the Redfield ratio of 16:1. Mixing of these deep-ocean waters into the surface waters of the subtropical gyres provides N and P in a ratio very near the Redfield ratio, so it should not be surprising that these open ocean systems are nearly co-limited by N and P, as originally described by Redfield (1958). On the other hand, the N:P ratio in the bottom waters near the coast are often well below the Redfield ratio (Fig. 3). The flux of nutrients from these near-bottom waters which are depleted in N relative to P would thus help maintain the deficits of inorganic N in coastal waters that flow into estuaries. The advection of dissolved materials into an estuary from coastal waters with low N:P ratios will import relatively more P and strengthen the likelihood of persistent N-limitation.

Fig. 3
figure 3

Molar ratio of nitrate + nitrite to SRP in the bottom waters of the world’s oceans

We observe that the N:P ratio for inorganic nutrients of bottom waters often appears to be lower along coasts that have a substantial continental shelf, such as the East and Gulf coasts of the United States, the European coast of the North Sea, and the coast of China (Fig. 3). Much of the research on nutrient limitation in estuaries has occurred in these areas, where N limitation has been found to be so common. Interestingly, oligotrophic lagoons along the French Mediterranean coast are P limited (Souchu et al. 2010), and the nearby bottom waters in the Mediterranean Sea have a very high N:P ratio for inorganic nutrients compared to most ocean regions (Fig. 3).

To further examine the effects of continental shelf biogeochemical processes in structuring N and P inputs to estuaries, we plot the inorganic nutrient concentrations for surface waters on the continental shelves of the United States as a function of the width of the continental shelf at the sampling site. Since the 1° resolution of the WOA18 can conflate shelf and open ocean observations, we used the GLobal Ocean Data Analysis Project (GLODAP) database (Olsen et al. 2016) which offers highly quality-controlled observations from precise sampling locations. The data are for individual samples from cruises over many years and are snapshots in time, not longer-term mean values.

For both nitrate + nitrite (Fig. 4a) and SRP (Fig. 4b), concentrations are quite variable in locations where the continental shelves are narrow. For areas where the shelf is less than 30 km wide, nitrate + nitrite varies from 0 to 33 μmol kg−1 and SRP varies from 0 to 2.4 μmol kg−1. In areas where the continental shelves are wider, concentrations of nitrate + nitrite and of SRP are not as high as sometimes observed at sites on the narrower shelves. For sites from areas where the shelf is wider than 100 km, nitrate + nitrite concentrations are generally at or near 0 μmol kg−1 (Fig. 4a), while SRP concentrations are < 0.5 μmol kg−1 (Fig. 4b). These low concentrations for nitrate + nitrite and SRP reflect both uptake by phytoplankton and net exchanges with sediments. The shallower water depths of the continental shelves can encourage phytoplankton blooms, with concomitant drawdown of both nitrate + nitrite and SRP. Further, the sediments on continental shelves are the major site both for burial of P (Ruttenberg and Berner 1993; Howarth et al. 1995) and of denitrification (Falkowski 1997; Galloway et al. 2004) globally in the world’s oceans. Denitrification can have a substantial influence on the inorganic N availability in regions where the continental shelf is wide. For example, denitrification on the continental shelf along the east coast of the United States and Canada is estimated to be the major N sink there (Seitzinger and Giblin 1996; Fennel et al. 2006).

Fig. 4
figure 4

a Nitrate + nitrite concentration in surface waters for sites on the continental shelves of the United States plotted as a function of the width of the shelf at that site. East coast and west coast shown separately. b SRP concentration in surface waters for sites on the continental shelves of the United States plotted as a function of the width of the shelf at that site. East coast and west coast shown separately. c Molar ratio of nitrate + nitrite to SRP in surface waters for sites on the continental shelves of the United States plotted as a function of the width of the shelf at that site. East coast and west coast shown separately

The magnitude as well as the relative importance of deep-ocean nutrient inputs onto the continental shelves can vary from region to region. All along the east coast of the United States and Canada, the advection of deep-ocean, nutrient-rich waters is the largest supply of nutrients to the shelves (Nixon et al. 1996; Michaels et al. 1996; Howarth et al. 1996; Howarth 1998). For example, the continental shelf along the northeastern United States from Maine south through Delaware receives 75% of its N inputs from the deep-ocean source, 1.54 Tg N year−1 compared to only 0.27 Tg N year−1 from watershed sources and 0.21 Tg N year−1 from direct atmospheric deposition (Howarth 1998). In other regions, the deep-ocean nutrient source can be less important. The northern Gulf of Mexico receives only 5% of its nitrogen from the deep ocean, 0.14 Tg N year−1 compared to 2.1 Tg N year−1 from the Mississippi River basin and 0.28 Tg N year−1 from direct atmospheric deposition (Howarth 1998).

The input of nutrients from the ocean (continental shelf) source is only one factor that determines overall nutrient availability and limitation in an estuary, but in many cases, this ocean-sourced water can be substantial and may contain proportionately more SRP than nitrate + nitrite, resulting in low N:P ratios (Fig. 4c). This is an important driver that can counterbalance the higher N:P ratio inputs from watersheds, and which may therefore make N limitation more likely in many estuaries than in lakes. Based on this factor, we suspect that N limitation is more likely in estuaries near a wide continental shelf in close proximity to deep-ocean waters, such as the East Coast of the United States or the North Sea, than those near a narrow continental shelf such as the west coast of the United States.

We want to emphasize that the data in Figs. 3, 4a–c are for inorganic N and SRP, not total N and P. In deep ocean waters, organic N and P are small compared to inorganic N and P (Michaels et al. 1996), but this is less true in surface waters or on the continental shelves, where inorganic N and SRP concentrations are lower. In a qualitative sense, we believe the inorganic data well represent the patterns of relative availabilities of N and P as sources of nutrients to surface waters and estuaries. If dissolved inorganic N remains low compared to SRP in the waters entering an estuary, it suggests a disproportionately high input of available P. Nonetheless, we urge more emphasis on patterns and fluxes of total N and P in future studies, with routine measurement of total N and P by oceanographers working at continental shelf sites as well as by estuarine scientists, and more interaction of biogeochemists with physical oceanographers in estimating total N and P fluxes in ecosystem budgets.

Concluding thoughts and suggested future research

Whether an estuary or other coastal marine ecosystem is N or P limited will be determined by the balance of N and P inputs from land, of N and P inputs from ocean sources, and of biogeochemical processes within the ecosystem that alter the availability of N and P. We suggest that the nutrient inputs from the deep ocean onto the continental shelf and into the coastal ocean and estuaries may be much more important in the overall balance than generally recognized, and we encourage more effort in characterizing these inputs and including them in ecosystem-scale N and P budgets. Further, whether primary productivity in any specific estuary is more limited by N or P should be evaluated with consideration of coastal ocean inputs, specifically considering the N:P ratio of the coastal waters and how these vary among regions partially due to differences in the width of the continental shelves and the rate of advection of deep-ocean waters onto the shelves.

As noted above, nutrient budgets for coastal marine ecosystems typically do not include the ocean source (Nixon et al. 1996). Of some of the previous efforts that did include ocean inputs in the nutrient budgets of coastal systems, such as the LOICZ project in the 1990s and 2000s, only dissolved inorganic nutrients were considered and not total N and P (Crossland et al. 2005; Swaney et al. 2011). There is great value in examining patterns of inorganic N and P in estuaries and other coastal marine ecosystems. For example, it is from these patterns for inorganic nutrients that Li et al. (2011) were able to infer the importance of the ocean source of P for maintaining N limitation in the Changjiang River estuary. However, we suggest that future ecosystem-scale budgets concentrate on fluxes of total N and total P, as NANI and NAPI do for watershed fluxes. This facilitates the construction of mass balances and also acknowledges that organic forms of N and P can be recycled and become available to primary producers within the ecosystem on ecologically meaningful time scales.

The input of nutrients from coastal ocean waters into an estuary depends strongly on the physics of the estuary, as this controls the rate of upstream advection. All else being equal, the upstream advection is higher in systems that have a greater tidal amplitude, but the rate of freshwater discharge into the estuary is also important. We hypothesize that an estuary can be oligotrophic only in cases where not only are watershed inputs of nutrients low but the advection into the system from coastal waters is also low, as in the micro-tidal lagoons along the French Mediterranean coast (Souchu et al. 2010). The corollary is that estuaries such as those along the New England coast would have been naturally mesotrophic before European settlement (with concomitant increase in the watershed nutrient load) simply because of the large tidal inputs of nutrients into these systems from ocean sources. While the nutrient inputs to pristine lakes will almost always be very low, this is not the case for many estuaries because of the ocean-source nutrients. As a generality, we would expect primary production in pristine estuaries to be greater than in pristine lakes.

Our analysis suggests that the width of the continental shelf may be one important factor regulating inorganic N and P concentrations as well as the N:P ratio of inorganic nutrients in water advected into estuaries and shallow coastal systems. However, the rate of advection of deep ocean water onto the continental shelf, the water residence time on the shelf, and the rate of advection of water and nutrients across the shelf to an estuary all seem likely to be important drivers controlling the nutrient inputs to individual systems. We therefore urge more collaborative studies between coastal oceanographers and estuarine biogeochemists, to evaluate how these factors affect the nutrient inputs to an estuary. In this context, we again suggest that both patterns of inorganic nutrients and fluxes of total N and P be considered.

Finally, human-caused climate change has already altered circulation of the world’s oceans, including near the continents (Doney et al. 2012; Rahmstorf et al. 2015). We hypothesize that changes in these oceanographic dynamics may alter the ocean source of nutrients to estuaries, and we speculate that such changes may already be occurring. As an example from our long-term study site, the P concentrations in the water entering West Falmouth Harbor from the coastal waters of Buzzards Bay now appear to be lower than a decade or so ago, both for SRP and total P (Hayn et al., manuscript in preparation). The total N concentration in Buzzards Bay also appears to be lower now than observed a decade ago. This may reflect less advection of deep ocean water onto the continental shelf south of Cape Cod or a longer residence time of the water on this shelf, leading to lower concentrations of N and P.