This paper is the outcome of a community initiative to identify major unsolved scientific problems in hydrology motivated by a need for stronger harmonisation of research efforts. The procedure involved a public consultation through on-line media, followed by two workshops through which a large number of potential science questions were collated, prioritised, and synthesised. In spite of the diversity of the participants (230 scientists in total), the process revealed much about community priorities and the state of our science: a preference for continuity in research questions rather than radical departures or redirections from past and current work. Questions remain focussed on process-based understanding of hydrological variability and causality at all space and time scales. Increased attention to environmental change drives a new emphasis on understanding how change propagates across interfaces within the hydrological system and across disciplinary boundaries. In particular, the expansion of the human footprint raises a new set of questions related to human interactions with nature and water cycle feedbacks in the context of complex water management problems. We hope that this reflection and synthesis of the 23 unsolved problems in hydrology will help guide research efforts for some years to come.

A deadly cascade: A catastrophic landslide in Uttarakhand state in India on February 2021 damaged two hydropower plants, and more than 200 people were killed or are missing. Shugar et al. describe the cascade of events that led to this disaster. A massive rock and ice avalanche roared down a Himalayan valley, turning into a deadly debris flow upstream from the first of the two hydropower plants. The sequence of events highlights the increasing risk in the Himalayas caused by increased warming and development. Science , abh4455, this issue p. 300

The Theia Snow collection routinely provides high-resolution maps of the snow-covered area from Sentinel-2 and Landsat-8 observations. The collection covers selected areas worldwide, including the main mountain regions in western Europe (e.g. Alps, Pyrenees) and the High Atlas in Morocco. Each product of the Theia Snow collection contains four classes: snow, no snow, cloud and no data. We present the algorithm to generate the snow products and provide an evaluation of the accuracy of Sentinel-2 snow products using in situ snow depth measurements, higher-resolution snow maps and visual control. The results suggest that the snow is accurately detected in the Theia snow collection and that the snow detection is more accurate than the Sen2Cor outputs (ESA level 2 product). An issue that should be addressed in a future release is the occurrence of false snow detection in some large clouds. The snow maps are currently produced and freely distributed on average 5 d after the image acquisition as raster and vector files via the Theia portal (https://doi.org/10.24400/329360/F7Q52MNK).

The seasonal snow in the Pyrenees is critical for hydropower production, crop irrigation and tourism in France, Spain and Andorra. Complementary to in situ observations, satellite remote sensing is useful to monitor the effect of climate on the snow dynamics. The MODIS daily snow products (Terra/MOD10A1 and Aqua/MYD10A1) are widely used to generate snow cover climatologies, yet it is preferable to assess their accuracies prior to their use. Here, we use both in situ snow observations and remote sensing data to evaluate the MODIS snow products in the Pyrenees. First, we compare the MODIS products to in situ snow depth (SD) and snow water equivalent (SWE) measurements. We estimate the values of the SWE and SD best detection thresholds to 40 mm water equivalent (w.e.) and 150 mm, respectively, for both MOD10A1 and MYD10A1. κ coefficients are within 0.74 and 0.92 depending on the product and the variable for these thresholds. However, we also find a seasonal trend in the optimal SWE and SD thresholds, reflecting the hysteresis in the relationship between the depth of the snowpack (or SWE) and its extent within a MODIS pixel. Then, a set of Landsat images is used to validate MOD10A1 and MYD10A1 for 157 dates between 2002 and 2010. The resulting accuracies are 97% (κ = 0.85) for MOD10A1 and 96% (κ = 0.81) for MYD10A1, which indicates a good agreement between both data sets. The effect of vegetation on the results is analyzed by filtering the forested areas using a land cover map. As expected, the accuracies decrease over the forests but the agreement remains acceptable (MOD10A1: 96%, κ = 0.77; MYD10A1: 95%, κ = 0.67). We conclude that MODIS snow products have a sufficient accuracy for hydroclimate studies at the scale of the Pyrenees range. Using a gap-filling algorithm we generate a consistent snow cover climatology, which allows us to compute the mean monthly snow cover duration per elevation band and aspect classes. There is snow on the ground at least 50% of the time above 1600 m between December and April. We finally analyze the snow patterns for the atypical winter 2011–2012. Snow cover duration anomalies reveal a deficient snowpack on the Spanish side of the Pyrenees, which seems to have caused a drop in the national hydropower production.