Land and ocean ecosystems are strongly connected and mutually interactive. As climate changes and other anthropogenic stressors intensify, the complex pathways that link these systems will strengthen or weaken in ways that are currently beyond reliable prediction. In this review we offer a framework of land–ocean couplings and their role in shaping marine ecosystems in coastal temperate rainforest (CTR) ecoregions, where high freshwater and materials flux result in particularly strong land–ocean connections. Using the largest contiguous expanse of CTR on Earth—the Northeast Pacific CTR (NPCTR)—as a case study, we integrate current understanding of the spatial and temporal scales of interacting processes across the land–ocean continuum, and examine how these processes structure and are defining features of marine ecosystems from nearshore to offshore domains. We look ahead to the potential effects of climate and other anthropogenic changes on the coupled land–ocean meta-ecosystem. Finally, we review key data gaps and provide research recommendations for an integrated, transdisciplinary approach with the intent to guide future evaluations of and management recommendations for ongoing impacts to marine ecosystems of the NPCTR and other CTRs globally. In the light of extreme events including heatwaves, fire, and flooding, which are occurring almost annually, this integrative agenda is not only necessary but urgent.

Riverine export of carbon (C) is an important part of the global C cycle; however, most riverine C budgets focus on individual forms of C and fail to comprehensively measure both organic and inorganic C species in concert. To address this knowledge gap, we conducted high frequency sampling of multiple C forms, including dissolved organic C (DOC), inorganic carbon (as alkalinity), particulate organic C (POC), coarse particulate organic C (CPOC), and invertebrate biomass C across the main run-off season in a predominantly rain-fed watershed in Southeast Alaska. Streamwater concentrations were used to model daily watershed C export from May through October. Concentration and modeled yield data indicated that DOC was the primary form of riverine C export (8708 kg C/km2), except during low flow periods when alkalinity (3125 kg C/km2) was the dominant form of C export. Relative to DOC and alkalinity, export of particulate organic C (POC: 992 kg C/km2; CPOC: 313 kg C/km2) and invertebrates (40 kg C/km2) was small, but these forms of organic matter could disproportionately impact downstream food webs because of their higher quality, assessed via C to nitrogen ratios. These seasonal and flow driven changes to C form and export likely provide subsidies to downstream and nearshore ecosystems such that predicted shifts in regional hydroclimate could substantially impact C transfer and incorporation into aquatic food webs.

Large-magnitude snow avalanches pose a hazard to humans and infrastructure worldwide. Analyzing the spatiotemporal behavior of avalanches and the contributory climate factors is important for understanding historical variability in climate-avalanche relationships as well as improving avalanche forecasting. This study uses established dendrochronological methods to develop long-term regional avalanche chronologies for three different climate types: high-latitude maritime climate of southeast Alaska, intermountain climate of the northern Rocky Mountains, and continental climate of Colorado. In the maritime study area, we collected 434 cross sections throughout six avalanche paths near Juneau, Alaska. This resulted in 2706 identified avalanche growth disturbances between year 1720 and 2018 Common Era (CE), which allowed us to reconstruct 82 years with large magnitude avalanche activity across three sub-regions. By combining this tree-ring derived avalanche dataset with a suite of climate and atmospheric variables and applying a generalized linear model to fit a binomial regression, we found February and March precipitation and the Oceanic Niño Index (ONI) were significant predictors of large magnitude avalanche activity in the southeast Alaska study area. In the intermountain climate study area, tree-rings from 647 trees exhibited 2134 avalanche-related growth disturbances in the northern Rocky Mountains of northwest Montana from 1867 to 2019. The data show that the amount of snowpack across the northern Rocky Mountain region is directly related to avalanche probability. Coincident with warming and regional snowpack reductions, a decline of ~ 14% (~ 2% per decade) in overall large magnitude avalanche probability is apparent through the period 1950–2017 CE. In the continental climate of Colorado, we sampled 24 avalanche paths throughout the state and collected 1188 total samples with 4135 identified growth disturbances from 1698 to 2019. Preliminary results suggest years with large magnitude avalanche activity across the sub-regions of this study area are generally characterized by stormy winters with above average snowpack development but that early and late winter temperature and precipitation also play an important role in large avalanche activity. Characterizing historical climate-avalanche relationships across different climate types provides a broad baseline for understanding potential future changes in avalanche activity. Overall, this work helps forecasters and planners better understand the influence of climate on large magnitude avalanche frequency, and how potential changes in avalanche character and occurrence will affect their operations in the context of a warming climate.

Dissolved inorganic carbon (DIC) represents an important but poorly constrained form of lateral carbon flux to the oceans. With high precipitation rates, large glaciers, and dense temperate rainforest, Southeast Alaska plays a critical role in the transport of carbon to the Gulf of Alaska (GOA). Previous estimates of DIC flux across the Southeast Alaska Drainage Basin (SEAKDB) are poorly constrained in space and time. Our goal was to incorporate recent measurements of DIC concentrations with previous measurements from the U.S. Geological Survey in order to model the spatial and temporal patterns of riverine DIC transport from SEAK to the GOA. We aggregated DIC concentration measurements from 1957 to 2020 and associated measurements of mean daily discharge. We then constructed load estimation models to generate concentration predictions across 24 watersheds. By spatially matching measurements of DIC with SEAKDB watersheds, we extrapolated concentration predictions across 2,455 watersheds encompassing approximately 190,000 km2. Models were aggregated according to two factors, the presence of karst and the discharge regime. Finally, monthly flux predictions were generated for each watershed using predicted concentrations and runoff estimates from the Distributed Climate Water Balance Model. Mean annual DIC flux from the SEAKDB was 2.36 Tg C with an average yield of 12.52 g C m−2. Both karst presence and flow regimes modified DIC flux and speciation across coastal marine areas. The high resolution of DIC flux estimates will provide useful inputs for describing seasonal C dynamics, and further refines our understanding of C budgets in the Pacific temperate rainforest and the surrounding marine environment.

Snow avalanches are a natural hazard in mountainous areas worldwide with severe impacts that include fatalities, damage to infrastructure, disruption to commerce, and landscape disturbance. Understanding long-term avalanche frequency patterns, and associated climate and weather influences, improves our understanding of how climate change may affect avalanche activity. We used dendrochronological techniques to evaluate the historical frequency of large magnitude avalanches (LMAs) in the high-latitude climate of southeast Alaska, United States. We collected 434 cross sections throughout six avalanche paths near Juneau, Alaska. This resulted in 2706 identified avalanche growth disturbances between 1720 and 2018, which allowed us to reconstruct 82 years with LMA activity across three sub-regions. By combining this tree-ring-derived avalanche data set with a suite of climate and atmospheric variables and applying a generalized linear model to fit a binomial regression, we found that February and March precipitation and the Oceanic Niño Index (ONI) were significant predictors of LMA activity in the study area. Specifically, LMA activity occurred during winters with substantial February and March precipitation and neutral or negative (cold) ONI values, while years not characterized by LMAs occur more frequently during warm winters (positive ONI values). Our examination of the climate-avalanche relationship in southeast Alaska sheds light on important climate variables and physical processes associated with LMA years. These results can be used to inform long-term infrastructure planning and avalanche mitigation operations in an urban area, such as Juneau, where critical infrastructure is subject to substantial avalanche hazard.

Biospheric particulate organic carbon (POCbio) burial and rock petrogenic particulate organic carbon (POCpetro) oxidation are opposing long-term controls on the global carbon cycle, sequestering and releasing carbon, respectively. Here, we examine how watershed glacierization impacts the POC source by assessing the concentration and isotopic composition (δ13C and Δ14C) of POC exported from four watersheds with 0%–49% glacier coverage across a melt season in Southeast Alaska. We used two mixing models (age-weight percent and dual carbon isotope) to calculate concentrations of POCbio and POCpetro within the bulk POC pool. The fraction POCpetro contribution was highest in the heavily glacierized watershed (age-weight percent: 0.39 ± 0.05; dual isotope: 0.42 (0.37–0.47)), demonstrating a glacial source of POCpetro to fjords. POCpetro was mobilized via glacier melt and subglacial flow, while POCbio was largely flushed from the non-glacierized landscape by rain. Flow normalized POCbio concentrations exceeded POCpetro concentrations for all streams, but surprisingly were highest in the heavily glacierized watershed (mean: 0.70 mgL−1; range 0.16–1.41 mgL−1), suggesting that glacier rivers can contribute substantial POCbio to coastal waters. Further, the most heavily glacierized watershed had the highest sediment concentration (207 mgL−1; 7–708 mgL−1), and thus may facilitate long-term POCbio protection via sediment burial in glacier-dominated fjords. Our results suggest that continuing glacial retreat will decrease POC concentrations and increase POCbio:POCpetro exported from currently glacierized watersheds. Glacier retreat may thus decrease carbon storage in marine sediments and provide a positive feedback mechanism to climate change that is sensitive to future changes in POCpetro oxidation.

Watershed classification has long been a key tool in the hydrological sciences, but few studies have been extended to biogeochemistry. We developed a combined hydro-biogeochemical classification for watersheds draining to the coastal margin of the Northeast Pacific coastal temperate rainforest (1,443,062 km2), including 2,695 small coastal rivers (SCR) and 10 large continental watersheds. We used cluster analysis to group SCR watersheds into 12 types, based on watershed properties. The most important variables for distinguishing SCR watershed types were evapotranspiration, slope, snowfall, and total precipitation. We used both streamflow and dissolved organic carbon (DOC) measurements from rivers (n = 104 and 90 watersheds respectively) to validate the classification. Watershed types corresponded with broad differences in streamflow regime, mean annual runoff, DOC seasonality, and mean DOC concentration. These links between watershed type and river conditions enabled the first region-wide empirical characterization of river hydrobiogeochemistry at the land-sea margin, spanning extensive ungauged and unsampled areas. We found very high annual runoff (mean > 3,000 mm, n = 10) in three watershed types totaling 59,024 km2 and ranging from heavily glacierized mountain watersheds with high flow in summer to a rain-fed mountain watershed type with high flow in fall-winter. DOC hotspots (mean > 4 mg L−1, n = 14) were found in three other watershed types (48,557 km2) with perhumid rainforest climates and less-mountainous topography. We described four patterns of DOC seasonality linked to watershed hydrology, with fall-flushing being widespread. Hydro-biogeochemical watershed classification may be useful for other complex regions with sparse observation networks.

Glacier retreat poses risks and benefits for species of cultural and economic importance. One example is Pacific salmon (Oncorhynchus spp.), supporting subsistence harvests, and commercial and recreational fisheries worth billions of dollars annually. Although decreases in summer streamflow and warming freshwater is reducing salmon habitat quality in parts of their range, glacier retreat is creating new streams and lakes that salmon can colonize. However, potential gains in future salmon habitat associated with glacier loss have yet to be quantified across the range of Pacific salmon. Here we project future gains in Pacific salmon freshwater habitat by linking a model of glacier mass change for 315 glaciers, forced by five different Global Climate Models, with a simple model of salmon stream habitat potential throughout the Pacific Mountain ranges of western North America. We project that by the year 2100 glacier retreat will create 6,146 (±1,619) km of new streams accessible for colonization by Pacific salmon, of which 1,930 (±569) km have the potential to be used for spawning and juvenile rearing, representing 0 to 27% gains within the 18 sub-regions we studied. These findings can inform proactive management and conservation of Pacific salmon in this era of rapid climate change.

With a unique biogeophysical signature relative to other freshwater sources, meltwater from glaciers plays a crucial role in the hydrological and ecological regime of high latitude coastal areas. Today, as glaciers worldwide exhibit persistent negative mass balance, glacier runoff is changing in both magnitude and timing, with potential downstream impacts on infrastructure, ecosystems, and ecosystem resources. However, runoff trends may be difficult to detect in coastal systems with large precipitation variability. Here, we use the coupled energy balance and water routing model SnowModel-HydroFlow to examine changes in timing and magnitude of runoff from the western Juneau Icefield in Southeast Alaska between 1980 and 2016. We find that under sustained glacier mass loss (−0.57 ± 0.12 m w. e. a−1), several hydrological variables related to runoff show increasing trends. This includes annual and spring glacier ice melt volumes (+10% and +16% decade−1) which, because of higher proportions of precipitation, translate to smaller increases in glacier runoff (+3% and +7% decade−1) and total watershed runoff (+1.4% and +3% decade−1). These results suggest that the western Juneau Icefield watersheds are still in an increasing glacier runoff period prior to reaching “peak water.” In terms of timing, we find that maximum glacier ice melt is occurring earlier (2.5 days decade−1), indicating a change in the source and quality of freshwater being delivered downstream in the early summer. Our findings highlight that even in maritime climates with large precipitation variability, high latitude coastal watersheds are experiencing hydrological regime change driven by ongoing glacier mass loss.