Many ecological insights into the function of rivers and watersheds emerge from quantifying the flux of solutes or suspended materials in rivers. Numerous methods for flux estimation have been described, and each has its strengths and weaknesses. Currently, the largest practical challenges in flux estimation are to select among these methods and to implement or apply whichever method is chosen. To ease this process of method selection and application, we have written an R software package called loadflex that implements several of the most popular methods for flux estimation, including regressions, interpolations, and the special case of interpolation known as the period‐weighted approach. Our package also implements a lesser‐known and empirically promising approach called the “composite method,” to which we have added an algorithm for estimating prediction uncertainty. Here we describe the structure and key features of loadflex, with a special emphasis on the rationale and details of our composite method implementation. We then demonstrate the use of loadflex by fitting four different models to nitrate data from the Lamprey River in southeastern New Hampshire, where two large floods in 2006–2007 are hypothesized to have driven a long‐term shift in nitrate concentrations and fluxes from the watershed. The models each give believable estimates, and yet they yield different answers for whether and how the floods altered nitrate loads. In general, the best modeling approach for each new dataset will depend on the specific site and solute of interest, and researchers need to make an informed choice among the many possible models. Our package addresses this need by making it simple to apply and compare multiple load estimation models, ultimately allowing researchers to estimate riverine concentrations and fluxes with greater ease and accuracy.

Tropical cyclones play an increasingly important role in shaping ecosystems. Understanding and generalizing their responses is challenging because of meteorological variability among storms and its interaction with ecosystems. We present a research framework designed to compare tropical cyclone effects within and across ecosystems that: a) uses a disaggregating approach that measures the responses of individual ecosystem components, b) links the response of ecosystem components at fine temporal scales to meteorology and antecedent conditions, and c) examines responses of ecosystem using a resistance–resilience perspective by quantifying the magnitude of change and recovery time. We demonstrate the utility of the framework using three examples of ecosystem response: gross primary productivity, stream biogeochemical export, and organismal abundances. Finally, we present the case for a network of sentinel sites with consistent monitoring to measure and compare ecosystem responses to cyclones across the United States, which could help improve coastal ecosystem resilience.

Concentration-discharge relationships are a key tool for understanding the sources and transport of material from watersheds to fluvial networks. Storm events in particular provide insight into variability in the sources of solutes and sediment within watersheds, and the hydrologic pathways that connect hillslope to stream channel. Here we examine high-frequency sensor-based specific conductance and turbidity data from multiple storm events across two watersheds (Quebrada Sonadora and Rio Icacos) with different lithology in the Luquillo Mountains of Puerto Rico, a forested tropical ecosystem. Our analyses include Hurricane Maria, a category 5 hurricane. To analyze hysteresis, we used a recently developed set of metrics to describe and quantify storm events including the hysteresis index (HI), which describes the directionality of hysteresis loops, and the flushing index (FI), which can be used to infer whether the mobilization of material is source or transport limited. We also examine the role of antecedent discharge to predict hysteretic behavior during storms. Overall, specific conductance and turbidity showed contrasting responses to storms. The hysteretic behavior of specific conductance was similar across sites, displaying clockwise hysteresis and a negative FI indicating proximal sources of solutes and consistent source limitation. In contrast, the directionality of turbidity hysteresis was significantly different between watersheds, although both had strong flushing behavior indicative of transport limitation. Overall, models that included antecedent discharge did not perform any better than models with peak discharge alone, suggesting that the magnitude and trajectory of an individual event was the strongest driver of material flux and hysteretic behavior. Hurricane Maria produced unique hysteresis metrics within both watersheds, indicating a distinctive response to this major hydrological event. The similarity in response of specific conductance to storms suggests that solute sources and pathways are similar in the two watersheds. The divergence in behavior for turbidity suggests that sources and pathways of particulate matter vary between the two watersheds. The use of high-frequency sensor data allows the quantification of storm events while index-based metrics of hysteresis allow for the direct comparison of complex storm events across a heterogeneous landscape and variable flow conditions.

Tropical cyclones drive coastal ecosystem dynamics, and their frequency, intensity, and spatial distribution are predicted to shift with climate change. Patterns of resistance and resilience were synthesized for 4138 ecosystem time series from n = 26 storms occurring between 1985 and 2018 in the Northern Hemisphere to predict how coastal ecosystems will respond to future disturbance regimes. Data were grouped by ecosystems (fresh water, salt water, terrestrial, and wetland) and response categories (biogeochemistry, hydrography, mobile biota, sedentary fauna, and vascular plants). We observed a repeated pattern of trade-offs between resistance and resilience across analyses. These patterns are likely the outcomes of evolutionary adaptation, they conform to disturbance theories, and they indicate that consistent rules may govern ecosystem susceptibility to tropical cyclones.

Barriers within streams can affect riverine species' ability to access habitats and may reduce their population viability. Connectivity metrics attempt to quantify the impacts of barriers; however, little is known about their functioning when applied to dendritic habitats such as watersheds. Several graph-theoretic connectivity metrics were calculated on rivers originating in the Luquillo Mountains of Puerto Rico. These metrics were classified into two primary groups: metrics that count weighted paths through the stream network and metrics that predict the flow of organisms through a stream reach. Representative metrics from each of these categories were suggested to model the effects of dams and water intakes, respectively. Figure 6. a and b. These show Cloesness Centrality (CC) and Change in Coincidence Probability (DCP) values for the western branch of the Espritu Santo. This figure is available in colour online at wileyonlinelibrary.com/journal/rra CONNECTIVITY METRICS ON WATERSHEDS 263 Figure 7. a and b. These show Closeness Centrality (CC) and Change in Coincidence Probability (DCP) values for the western branch of the Icacos/Rio Blanco. This figure is available in colour online at wileyonlinelibrary.com/journal/rra U. MALVADKAR ET AL.