Reliable precipitation forcing is essential for calculating the water balance, seasonal snowpack, glacier mass balance, streamflow, and other hydrological variables. However, satellite precipitation is often the only forcing available to run hydrological models in data-scarce regions, compromising hydrological calculations when unreliable. The IMERG product estimates precipitation quasi-globally from a combination of passive microwave and infrared satellites, which are intercalibrated based on GPM’s DPR and GMI instruments. Current GPM-DPR radar algorithms have satisfactorily estimated rainfall, but a limited consideration of PSD, attenuation correction, and ground clutter have degraded snowfall estimation, especially in mountain regions. This study aims to improve satellite radar snowfall estimates for this situation. Nearly two years (between 2019 and 2022) of aloft precipitation concentration, surface hydrometeor size, number and fall velocity, and surface precipitation rate from a high elevation site in the Canadian Rockies and collocated GPM-DPR reflectivities were used to develop a new snowfall estimation algorithm. Snowfall estimates using the new algorithm and measured GPM-DPR reflectivities were compared to other GPM-DPR-based products, including CORRA, which is employed to intercalibrate IMERG. Snowfall rates estimated with measured Ka reflectivities, and from CORRA were compared to MRR-2 observations, and had correlation, bias, and RMSE of 0.58 and 0.07, 0.43 and -0.38 mm h-1, and 0.83 and 0.85 mm h-1, respectively. Predictions using measured Ka reflectivity suggest that enhanced satellite radar snowfall estimates can be achieved using a simple measured reflectivity algorithm. These improved snowfall estimates can be adopted to intercalibrate IMERG in cold mountain regions, thereby improving regional precipitation estimates.

Wildfires and heatwaves have recently affected the hydrological system in unprecedented ways due to climate change. In cold regions, these extremes cause rapid reductions in snow and ice albedo due to soot deposition and unseasonal melt. Snow and ice albedo dynamics control net shortwave radiation and the available energy for melt and runoff generation. Many albedo algorithms in hydrological models cannot accurately simulate albedo dynamics because they were developed or parameterised based on historical observations. Remotely sensed albedo data assimilation (DA) can potentially improve model performance by updating modelled albedo with observations. This study seeks to diagnose the effects of remotely sensed snow and ice albedo DA on the prediction of streamflow from glacierized basins during wildfires and heatwaves. Sentinel-2 20-m albedo estimates were assimilated into a glacio-hydrological model created using the Cold Regions Hydrological Modelling Platform (CRHM) in two Canadian Rockies glacierized basins, Athabasca Glacier Research Basin (AGRB) and Peyto Glacier Research Basin (PGRB). The study was conducted in 2018 (wildfires), 2019 (soot/algae), 2020 (normal), and 2021 (heatwaves). DA was employed to assimilate albedo into CRHM to simulate streamflow and was compared to a control run (CTRL) using off-the-shelf albedo parameters. Albedo DA benefited streamflow predictions during wildfires for both basins, with a KGE coefficient improvement of 0.18 and 0.20 in AGRB and PGRB, respectively. Four-year DA streamflow predictions were superior to CTRL in PGRB, but DA was slightly better in AGRB. DA was not beneficial to streamflow predictions during heatwaves. These results show that albedo DA can reveal otherwise unknown albedo and snowpack dynamics occurring in remote glacier accumulation zones that are not well simulated by model predictions alone. These findings corroborate the power of observational tools to incorporate near real-time information into hydrological models to better inform water managers of the streamflow response to wildfires and heatwaves.