An Enkf-Based Scheme for Snow Multivariable Data Assimilation at an Alpine Site

Abstract The knowledge of snowpack dynamics is of critical importance to several real-time applications especially in mountain basins, such as agricultural production, water resource management, flood prevention, hydropower generation. Since simulations are affected by model biases and forcing data uncertainty, an increasing interest focuses on the assimilation of snow-related observations with the purpose of enhancing predictions on snowpack state. The study aims at investigating the effectiveness of snow multivariable data assimilation (DA) at an Alpine site. The system consists of a snow energy-balance model strengthened by a multivariable DA system. An Ensemble Kalman Filter (EnKF) scheme allows assimilating ground-based and remotely sensed snow observations in order to improve the model simulations. This research aims to investigate and discuss: (1) the limitations and constraints in implementing a multivariate EnKF scheme in the framework of snow modelling, and (2) its performance in consistently updating the snowpack state. The performance of the multivariable DA is shown for the study case of Torgnon station (Aosta Valley, Italy) in the period June 2012 - December 2013. The results of several experiments are discussed with the aim of analyzing system sensitivity to the DA frequency, the ensemble size, and the impact of assimilating different observations

Multivariate data assimilation in snow modelling at Alpine sites

The knowledge of snowpack dynamics is of critical importance to several real-time applications such as agricultural production, water resource management, flood prevention, hydropower generation, especially in mountain basins. Snowpack state can be estimated by models or from observations, even though both these sources of information are affected by several errors

Transfer of Snow Information across the Macro-to-Hillslope-Scale Gap Using a Physiographic Downscaling Approach: Implications for Hydrologic Modeling in Semiarid, Seasonally Snow-Dominated Watersheds

Snow and ice are substantial components of the global energy balance and hydrologic cycle. Seasonal snow covers an area of 47 million km2 at its average maximum extent, 98% of which occurs across the Northern Hemisphere. The earth’s radiation budget is largely controlled by the fraction of absorbed solar energy, a parameter that is dependent upon snow surface albedo. Mountain snowpacks act as natural reservoirs, storing large quantities of water throughout the winter until eventual release during the melting phase. Accurate characterization of snow-covered area (SCA) and snow water equivalent (SWE) in such terrain could substantially improve the estimation of timing and volume of melt water runoff. However, knowledge of these hydrologic states is limited in part by scarcely populated in situ observation networks and logistical constraints in field survey sampling. Thus, satellite remote sensing observations are often employed in conjunction with simulation models to improve the estimation of snowpack states and resultant fluxes. This study attempts to merge complementary datasets in order to predict spatially variable snow processes at high resolution in basins exhibiting complex terrain. Specifically, the goal is to provide a means to downscale existing remote sensing and snow modeling datasets using computationally efficient methods that utilize physiographic information regarding terrain and land cover. A linear combination model is proposed for downscaling fractional SCA from the Moderate Resolution Imaging Spectroradiometer (MODIS) instrument from its native resolution (500 m) to a hillslope-scale resolution (e.g., 10-30 m), preserving the predicted snow cover fraction at the basin scale. The model is calibrated to 30 m Landsat observations using elevation and incoming solar radiation indices for a study area in southwestern Idaho. Validation is performed with data not used during calibration. Results depict favorable model performance when comparing downscaled MODIS snow cover to Landsat binary observations. An “ideal” validation test is performed in which Landsat aggregate 500 m snow fraction informs the model with similarly positive results. The use of such an algorithm might benefit applications from flood forecasting to SWE reconstruction. In a snowmelt modeling application, the satellite-derived snow cover downscaling algorithm is applied as a binary mask to constrain spatial melt runoff data from the SNOw Data Assimilation System (SNODAS). Differential solar radiation, forest canopy, and snow albedo estimates are also used to further downscale the modeled melt. Comparison with available field lysimeter data show proper spatial disaggregation of modeled melt onto opposing hillslopes, though timing and magnitude issues exist. Implications for resolving snowmelt at hillslope scales are briefly discussed

Estimating The Spatial Distribution of Snow Water Equivalent Using in situ and Remote Sensing Observations

Mountain snowpack is one of the primary surface water sources for about one-sixth of the global population. More than 75% of the total runoff originates from mountain snowpacks in the Western U.S. Snowmelt water recharges reservoirs and aquifers gradually in the melting season, providing vital water supplies for urban and agricultural areas. Therefore, accurately monitoring the spatial and temporal distribution of mountain snowpack &ndash; often measured as snow water equivalent (SWE) &ndash; is crucial for effective water management. While existing SWE estimation approaches remain highly uncertain, particularly when applied over large mountainous regions, the remotely-sensed snow data provide new opportunities to better characterize the spatial distributions of mountain snowpack. This dissertation investigates the approaches that optimally blend satellite, airborne, and ground snow observations to improve (near) real-time SWE estimation over mountainous terrain. The second chapter of this dissertation evaluates the accuracy of existing SWE estimation models in Sierra Nevada California. Five large-scale SWE datasets at fine spatial resolutions (&lt;= 1000 m) are comprehensively validated and compared with the Airborne Snow Observatory (ASO) SWE data in the Tuolumne River Basin (2013-2017), and ground snow pillow and snow course SWE observations across the Sierra Nevada (2004-2014). These SWE datasets include REC-INT, REC-ParBal, a Sierra Nevada SWE reanalysis (REC-DA), and two operational SWE datasets from the Snow Data Assimilation System (SNODAS) and the National Water Model (NWM-SWE), respectively. The results show that the REC-DA overall provides the most accurate SWE estimates across the Sierra Nevada (R2 = 0.87, MAE = 66 mm, PBIAS = 8.3%), followed by the REC-ParBal (R2 = 0.73, MAE = 83 mm, PBIAS = -6.4%), which is the least biased SWE estimates. Generally, SNODAS (R2 = 0.66, MAE = 106 mm, PBIAS = 9.3%) and REC-INT (R2 = 0.61, MAE = 131 mm, PBIAS = -28.3%) exhibit comparable but lower accuracy than the earlier mentioned two datasets, while NWM-SWE (R2 = 0.49, MAE = 142 mm, PBIAS = -25.2%) shows the least accuracy among the five SWE datasets. Given that REC-DA is not applicable in real-time, in the third chapter, a SWE data-fusion framework is developed, which integrates the historical SWE patterns derived from REC-DA into a statistically-based linear regression model (LRM) to estimate SWE in real-time. To investigate the influence of satellite-observed daily mean fractional snow-covered area (DMFSCA) on SWE estimation accuracy, two LRMs are compared: a baseline regression model (LRM-baseline) in which physiographic data and historical SWE patterns are used as independent variables, and an FSCA-informed regression model (LRM-FSCA) in which the DMFSCA from MODIS satellite imagery is included as an additional independent variable. By incorporating DMFSCA, LRM-FSCA outperforms LRM-baseline with improved R2 from 0.54 to 0.60, and reduced PBIAS from 2.6% to 2.2% in snow pillow cross-validation. The improvement in LRM-FSCA&rsquo;s performance is more significant during snow accumulation periods than during the snowmelt seasons. Compared to the ASO SWE, the LRM-FSCA explains 85% of the variance on average, which is at least 21% higher than the operational SNODAS (R2 = 0.64) and NWM-SWE (R2 = 0.33) in comparison. In chapter 4, a SWE bias correction framework (SWE-BCF) is developed that incorporates the ASO SWE and machine learning (ML) algorithms to further improve LRM SWE estimates in real-time. The performance of a wide range of commonly used machine learning algorithms is examined in the SWE-BCF including Gaussian Process Regression (GPR), Support Vector Machine (SVM), Bayesian Regularized Neural Networks (BRNN), Random Forest (RF), and Gradient Boosting Machine (GBM). The results indicate that all ML algorithms are capable of improving LRM-SWE accuracy substantially. While no single model performs significantly better than others, GPR, overall, shows the best performance with a 20% (0.14) increase in mean R2 value, a 31% (51 mm) reduction in mean RMSE, and a 61% (18.0%) reduction in absolute PBIAS compared with the original LRM using ASO SWE data for model validation. RF shows the most robust and stable performance in SWE bias correction with a 10% (0.08) increase in median R2 and a 41% (50 mm) reduction in median RMSE compared with the original LRM.</p

Désagrégation de l'humidité du sol issue des produits satellitaires micro-ondes passives et exploration de son utilisation pour l'amélioration de la modélisation et la prévision hydrologique

De plus en plus de produits satellitaires en micro-ondes passives sont disponibles. Cependant, leur large résolution spatiale (25-50 km) n’en font pas un outil adéquat pour des applications hydrologiques à une échelle locale telles que la modélisation et la prévision hydrologiques. Dans de nombreuses études, une désagrégation d’échelle de l’humidité du sol des produits satellites micro-ondes est faite puis validée avec des mesures in-situ. Toutefois, l’utilisation de ces données issues d’une désagrégation d’échelle n’a pas encore été pleinement étudiée pour des applications en hydrologie. Ainsi, l’objectif de cette thèse est de proposer une méthode de désagrégation d’échelle de l’humidité du sol issue de données satellitaires en micro-ondes passives (Satellite Passive Microwave Active and Passive - SMAP) à différentes résolutions spatiales afin d’évaluer leur apport sur l’amélioration potentielle des modélisations et prévisions hydrologiques. À partir d’un modèle de forêt aléatoire, une désagrégation d’échelle de l’humidité du sol de SMAP l’amène de 36-km de résolution initialement à des produits finaux à 9-, 3- et 1-km de résolution. Les prédicteurs utilisés sont à haute résolution spatiale et de sources différentes telles que Sentinel-1A, MODIS et SRTM. L'humidité du sol issue de cette désagrégation d’échelle est ensuite assimilée dans un modèle hydrologique distribué à base physique pour tenter d’améliorer les sorties de débit. Ces expériences sont menées sur les bassins versants des rivières Susquehanna (de grande taille) et Upper-Susquehanna (en comparaison de petite taille), tous deux situés aux États-Unis. De plus, le modèle assimile aussi des données d’humidité du sol en profondeur issue d’une extrapolation verticale des données SMAP. Par ailleurs, les données d’humidité du sol SMAP et les mesures in-situ sont combinées par la technique de fusion conditionnelle. Ce produit de fusion SMAP/in-situ est assimilé dans le modèle hydrologique pour tenter d’améliorer la prévision hydrologique sur le bassin versant Au Saumon situé au Québec. Les résultats montrent que l'utilisation de l’humidité du sol à fine résolution spatiale issue de la désagrégation d’échelle améliore la représentation de la variabilité spatiale de l’humidité du sol. En effet, le produit à 1- km de résolution fournit plus de détails que les produits à 3- et 9-km ou que le produit SMAP de base à 36-km de résolution. De même, l’utilisation du produit de fusion SMAP/ in-situ améliore la qualité et la représentation spatiale de l’humidité du sol. Sur le bassin versant Susquehanna, la modélisation hydrologique s’améliore avec l’assimilation du produit de désagrégation d’échelle à 9-km, sans avoir recours à des résolutions plus fines. En revanche, sur le bassin versant Upper-Susquehanna, c’est le produit avec la résolution spatiale la plus fine à 1- km qui offre les meilleurs résultats de modélisation hydrologique. L’assimilation de l’humidité du sol en profondeur issue de l’extrapolation verticale des données SMAP n’améliore que peu la qualité du modèle hydrologique. Par contre, l’assimilation du produit de fusion SMAP/in-situ sur le bassin versant Au Saumon améliore la qualité de la prévision du débit, même si celle-ci n’est pas très significative.Abstract: The availability of satellite passive microwave soil moisture is increasing, yet its spatial resolution (i.e., 25-50 km) is too coarse to use for local scale hydrological applications such as streamflow simulation and forecasting. Many studies have attempted to downscale satellite passive microwave soil moisture products for their validation with in-situ soil moisture measurements. However, their use for hydrological applications has not yet been fully explored. Thus, the objective of this thesis is to downscale the satellite passive microwave soil moisture (i.e., Satellite Microwave Active and Passive - SMAP) to a range of spatial resolutions and explore its value in improving streamflow simulation and forecasting. The random forest machine learning technique was used to downscale the SMAP soil moisture from 36-km to 9-, 3- and 1-km spatial resolutions. A combination of host of high-resolution predictors derived from different sources including Sentinel-1A, MODIS and SRTM were used for downscaling. The downscaled SMAP soil moisture was then assimilated into a physically-based distributed hydrological model for improving streamflow simulation for Susquehanna (larger in size) and Upper Susquehanna (relatively smaller in size) watersheds, located in the United States. In addition, the vertically extrapolated SMAP soil moisture was assimilated into the model. On the other hand, the SMAP and in-situ soil moisture were merged using the conditional merging technique and the merged SMAP/in-situ soil moisture was then assimilated into the model to improve streamflow forecast over the au Saumon watershed. The results show that the downscaling improved the spatial variability of soil moisture. Indeed, the 1-km downscaled SMAP soil moisture presented a higher spatial detail of soil moisture than the 3-, 9- or original resolution (36-km) SMAP product. Similarly, the merging of SMAP and in-situ soil moisture improved the accuracy as well as spatial representation soil moisture. Interestingly, the assimilation of the 9-km downscaled SMAP soil moisture significantly improved the accuracy of streamflow simulation for the Susquehanna watershed without the need of going to higher spatial resolution, whereas for the Upper Susquehanna watershed the 1-km downscaled SMAP showed better results than the coarser resolutions. The assimilation of vertically extrapolated SMAP soil moisture only slightly further improved the accuracy of the streamflow simulation. On the other hand, the assimilation of merged SMAP/in-situ soil moisture for the au Saumon watershed improved the accuracy of streamflow forecast, yet the improvement was not that significant. Overall, this study demonstrated the potential of satellite passive microwave soil moisture for streamflow simulation and forecasting

Community Review of Southern Ocean Satellite Data Needs

This review represents the Southern Ocean community’s satellite data needs for the coming decade. Developed through widespread engagement, and incorporating perspectives from a range of stakeholders (both research and operational), it is designed as an important community-driven strategy paper that provides the rationale and information required for future planning and investment. The Southern Ocean is vast but globally connected, and the communities that require satellite-derived data in the region are diverse. This review includes many observable variables, including sea-ice properties, sea-surface temperature, sea-surface height, atmospheric parameters, marine biology (both micro and macro) and related activities, terrestrial cryospheric connections, sea-surface salinity, and a discussion of coincident and in situ data collection. Recommendations include commitment to data continuity, increase in particular capabilities (sensor types, spatial, temporal), improvements in dissemination of data/products/uncertainties, and innovation in calibration/validation capabilities. Full recommendations are detailed by variable as well as summarized. This review provides a starting point for scientists to understand more about Southern Ocean processes and their global roles, for funders to understand the desires of the community, for commercial operators to safely conduct their activities in the Southern Ocean, and for space agencies to gain greater impact from Southern Ocean-related acquisitions and missions.The authors acknowledge the Climate at the Cryosphere program and the Southern Ocean Observing System for initiating this community effort, WCRP, SCAR, and SCOR for endorsing the effort, and CliC, SOOS, and SCAR for supporting authors’ travel for collaboration on the review. Jamie Shutler’s time on this review was funded by the European Space Agency project OceanFlux Greenhouse Gases Evolution (Contract number 4000112091/14/I-LG)

Summaries of the Sixth Annual JPL Airborne Earth Science Workshop

This publication contains the summaries for the Sixth Annual JPL Airborne Earth Science Workshop, held in Pasadena, California, on March 4-8, 1996. The main workshop is divided into two smaller workshops as follows: (1) The Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) workshop, on March 4-6. The summaries for this workshop appear in Volume 1; (2) The Airborne Synthetic Aperture Radar (AIRSAR) workshop, on March 6-8. The summaries for this workshop appear in Volume 2