skip to main content


The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 11:00 PM ET on Friday, April 12 until 2:00 AM ET on Saturday, April 13 due to maintenance. We apologize for the inconvenience.

Title: Snow Phenology and Hydrologic Timing in the Yukon River Basin, AK, USA
The Yukon River basin encompasses over 832,000 km2 of boreal Arctic Alaska and northwest Canada, providing a major transportation corridor and multiple natural resources to regional communities. The river seasonal hydrology is defined by a long winter frozen season and a snowmelt-driven spring flood pulse. Capabilities for accurate monitoring and forecasting of the annual spring freshet and river ice breakup (RIB) in the Yukon and other northern rivers is limited, but critical for understanding hydrologic processes related to snow, and for assessing flood-related risks to regional communities. We developed a regional snow phenology record using satellite passive microwave remote sensing to elucidate interactions between the timing of upland snowmelt and the downstream spring flood pulse and RIB in the Yukon. The seasonal snow metrics included annual Main Melt Onset Date (MMOD), Snowoff (SO) and Snowmelt Duration (SMD) derived from multifrequency (18.7 and 36.5 GHz) daily brightness temperatures and a physically-based Gradient Ratio Polarization (GRP) retrieval algorithm. The resulting snow phenology record extends over a 29-year period (1988–2016) with 6.25 km grid resolution. The MMOD retrievals showed good agreement with similar snow metrics derived from in situ weather station measurements of snowpack water equivalence (r = 0.48, bias = −3.63 days) and surface air temperatures (r = 0.69, bias = 1 day). The MMOD and SO impact on the spring freshet was investigated by comparing areal quantiles of the remotely sensed snow metrics with measured streamflow quantiles over selected sub-basins. The SO 50% quantile showed the strongest (p < 0.1) correspondence with the measured spring flood pulse at Stevens Village (r = 0.71) and Pilot (r = 0.63) river gaging stations, representing two major Yukon sub-basins. MMOD quantiles indicating 20% and 50% of a catchment under active snowmelt corresponded favorably with downstream RIB (r = 0.61) from 19 river observation stations spanning a range of Yukon sub-basins; these results also revealed a 14–27 day lag between MMOD and subsequent RIB. Together, the satellite based MMOD and SO metrics show potential value for regional monitoring and forecasting of the spring flood pulse and RIB timing in the Yukon and other boreal Arctic basins.  more » « less
Award ID(s):
Author(s) / Creator(s):
; ; ; ;
Date Published:
Journal Name:
Remote Sensing
Page Range / eLocation ID:
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Arctic landscapes are warming and becoming wetter due to changes in precipitation and the timing of snowmelt which consequently alters seasonal runoff and river discharge patterns. These changes in hydrology lead to increased mobilization and transport of terrestrial dissolved organic matter (DOM) to Arctic coastal seas where significant impacts on biogeochemical cycling can occur. Here, we present measurements of dissolved organic carbon (DOC) and chromophoric DOM (CDOM) in the Yukon River-to-Bering Sea system and two river plumes on the Alaska North Slope which flow into the Beaufort Sea. Our sampling characterized optical and biogeochemical properties of DOM during high and low river discharge periods for the Yukon River-Bering Sea system. The average DOC concentration at the multiple Yukon River mouths ranged from a high of 10.36 mg C L -1 during the ascending limb of the 2019 freshet (late May), 6.4 mg C L -1 during the descending limb of the 2019 freshet (late June), and a low of 3.86 mg C L -1 during low river discharge in August 2018. CDOM absorption coefficient at 412 nm ( a CDOM (412)) averaged 8.23 m -1 , 5.07 m -1 , and 1.9 m -1 , respectively. Several approaches to model DOC concentration based on its relationship with CDOM properties demonstrated cross-system seasonal and spatial robustness for these Arctic coastal systems despite spanning an order of magnitude decrease in DOC concentration from the lower Yukon River to the Northern Bering Sea as well as the North Slope systems. “Snapshot” fluxes of DOC and CDOM across the Yukon River Delta to Norton Sound were calculated from our measurements and modeled water fluxes forced with upstream USGS river gauge data. Our findings suggest that during high river flow, DOM reaches the delta largely unaltered by inputs or physical and biogeochemical processing and that the transformations of Yukon River DOM largely occur in the plume. However, during low summer discharge, multiple processes including local precipitation events, microbial decomposition, photochemistry, and likely others can alter the DOM properties within the lower Yukon River and Delta prior to flowing into Norton Sound. 
    more » « less
  2. null (Ed.)
    Among all the natural hazards throughout the world, floods occur most frequently. They occur in high latitude regions, such as: 82% of the area of North America; most of Russia; Norway, Finland, and Sweden in North Europe; China and Japan in Asia. River flooding due to ice jams may happen during the spring breakup season. The Northeast and North Central region, and some areas of the western United States, are especially harmed by floods due to ice jams and snowmelt. In this study, observations from operational satellites are used to map and monitor floods due to ice jams and snowmelt. For a coarse-to-moderate resolution sensor on board the operational satellites, like the Visible Infrared Imaging Radiometer Suite (VIIRS) on board the National Polar-orbiting Partnership (NPP) and the Joint Polar Satellite System (JPSS) series, and the Advanced Baseline Imager (ABI) on board the GOES-R series, a pixel is usually composed of a mix of water and land. Water fraction can provide more information and can be estimated through mixed-pixel decomposition. The flood map can be derived from the water fraction difference after and before flooding. In high latitude areas, while conventional observations are usually sparse, multiple observations can be available from polar-orbiting satellites during a single day, and river forecasters can observe ice movement, snowmelt status and flood water evolution from satellite-based flood maps, which is very helpful in ice jam determination and flood prediction. The high temporal resolution of geostationary satellite imagery, like that of the ABI, can provide the greatest extent of flood signals, and multi-day composite flood products from higher spatial resolution imagery, such as VIIRS, can pinpoint areas of interest to uncover more details. One unique feature of our JPSS and GOES-R flood products is that they include not only normal flood type, but also a special flood type as the supra-snow/ice flood, and moreover, snow and ice masks. Following the demonstrations in this study, it is expected that the JPSS and GOES-R flood products, with ice and snow information, can allow dynamic monitoring and prediction of floods due to ice jams and snowmelt for wide-end users. 
    more » « less
  3. Abstract

    Arctic rivers are sensitive to climate and environmental change, but the biogeochemical response remains poorly understood. Monthly size‐fractionated dissolved organic matter (DOM) samples from the lower Yukon River were characterized using UV–visible, fluorescence, and Fourier transform‐infrared (FT‐IR) spectroscopy techniques. The EEM‐PARAFAC analysis revealed three major fluorescent DOM components, including two humic‐like components (C480and C400) and one protein‐like component (C310), with their relative importance following the order of C480 ≥ C400 > C310in the high‐molecular‐weight DOM (1 kDa–0.4 μm) and C400 > C480 > C310in the low‐molecular‐weight DOM pool (< 1 kDa). Transformation in DOM and change in sources were manifested in major fluorescent components and optical properties, including biological index (BIX), humification index (HIX), spectral slope (S275–295) and specific UV absorbance at 254 nm (SUVA254). These changes occurred within different DOM size‐fractions and among ice‐covered, spring freshet, and open seasons. Joint analysis of EEM and FT‐IR spectra using a data fusion technique showed that humic‐like DOM is mostly associated with C─H, C═C, and C─O bonds, while protein‐like DOM is correlated more with C─N and N─H related structures. DOM aromaticity and the ratios of HIX to BIX and protein‐like to humic‐like components may be used as a compelling proxy to measure change in source waters and to infer permafrost dynamics. Our results provide insight into the seasonal variation in DOM composition for different size‐fractions in the lower Yukon River, and a baseline dataset against which future changes can be understood in the context of arctic basin biogeochemical cycling.

    more » « less
  4. Forecasting the timing and magnitude of snowmelt and runoff is critical to managing mountain water resources. Warming temperatures are increasing the rain–snow transition elevation and are limiting the forecasting skill of statistical models relating historical snow water equivalent to streamflow. While physically based methods are available, they require accurate estimations of the spatial and temporal distribution of meteorological variables in complex terrain. Across many mountainous areas, measurements of precipitation and other meteorological variables are limited to a few reference stations and are not adequate to resolve the complex interactions between topography and atmospheric flow. In this paper, we evaluate the ability of the Weather Research and Forecasting (WRF) Model to approximate the inputs required for a physics-based snow model, iSnobal, instead of using meteorological measurements, for the Boise River Basin (BRB) in Idaho, United States. An iSnobal simulation using station data from 40 locations in and around the BRB resulted in an average root-mean-square error (RMSE) of 4.5 mm compared with 12 SNOTEL measurements. Applying WRF forcings alone was associated with an RMSE of 10.5 mm, while including a simple bias correction to the WRF outputs of temperature and precipitation reduced the RMSE to 6.5 mm. The results highlight the utility of using WRF outputs as input to snowmelt models, as all required input variables are spatiotemporally complete. This will have important benefits in areas with sparse measurement networks and will aid snowmelt and runoff forecasting in mountainous basins.

    more » « less
  5. Abstract

    The springtime transition to regional‐scale onset of photosynthesis and net ecosystem carbon uptake in boreal and tundra ecosystems are linked to the soil freeze–thaw state. We present evidence from diagnostic and inversion models constrained by satellite fluorescence and airborneCO2from 2012 to 2014 indicating the timing and magnitude of spring carbon uptake in Alaska correlates with landscape thaw and ecoregion. Landscape thaw in boreal forests typically occurs in late April (DOY111 ± 7) with a 29 ± 6 day lag until photosynthetic onset. North Slope tundra thaws 3 weeks later (DOY133 ± 5) but experiences only a 20 ± 5 day lag until photosynthetic onset. These time lag differences reflect efficient cold season adaptation in tundra shrub and the longer dehardening period for boreal evergreens. Despite the short transition from thaw to photosynthetic onset in tundra, synchrony of tundra respiration with snow melt and landscape thaw delays the transition from net carbon loss (at photosynthetic onset) to net uptake by 13 ± 7 days, thus reducing the tundra net carbon uptake period. Two globalCO2inversions using aCASAGFEDmodel prior estimate earlier northern high latitude net carbon uptake compared to our regional inversion, which we attribute to (i) early photosynthetic‐onset model prior bias, (ii) inverse method (scaling factor + optimization window), and (iii) sparsity of available AlaskanCO2observations. Another global inversion with zero prior estimates the same timing for net carbon uptake as the regional model but smaller seasonal amplitude. The analysis of Alaskan eddy covariance observations confirms regional scale findings for tundra, but indicates that photosynthesis and net carbon uptake occur up to 1 month earlier in evergreens than captured by models orCO2inversions, with better correlation to above‐freezing air temperature than date of primary thaw. Further collection and analysis of boreal evergreen species over multiple years and at additional subarctic flux towers are critically needed.

    more » « less