In this work, we develop a software suite for studies of atmosphere – underlying SNOW- spaceborne optical receiver light TRANsmission calculations (SNOWTRAN) with applications for the solution of forward and inverse radiative transfer problems in polar regions. Assuming that aerosol load is extremely low, the proposed theory does not require the numerical procedures for the solution of the radiative transfer equation and is based on analytical equations for the spectral nadir reflectance and simple approximations for local optical properties of atmosphere and snow. The model developed is validated using EnMAP and PRISMA spaceborne imaging spectroscopy data close to the Concordia research station in Antarctica. A new fast technique for the determination of the snow grain size and assessment of the snowpack vertical inhomogeneity is then proposed and further demonstrated on EnMAP imagery over the Aviator Glacier and in the vicinity of the Concordia research station in Antarctica. The results revealed a large increase in precipitable water vapor at the Concordia research station in February 2023 linked to warming event, and a 4 times larger grain size at Aviator Glacier compared to Dome C.

Sentinel-2 provides the opportunity to map the snow cover at unprecedented spatial and temporal resolution at global scale. Here we calibrate and evaluate a simple empirical function to estimate the fractional snow cover (FSC) in open terrain using the normalized difference snow index (NDSI) from 20 m resolution Sentinel-2 images. The NDSI is computed from flat surface reflectances after masking cloud and snow-free areas. The NDSI-FSC function is calibrated using Pléiades very high resolution images and evaluated using independent datasets including SPOT 6/7 satellite images, time lapse camera photographs, terrestrial lidar scans and crowd-sourced in situ measurements. The calibration results show that the FSC can be represented with a sigmoid-shaped function 0.5×tanh(a×NDSI+b)+0.5 where a = 2.65 and b = -1.42 yielding a root mean square error of 25%. Similar RMSE are obtained with different evaluation datasets with a high topographic variability. With this function, we estimate that the confidence interval on the FSC retrievals is 38% at the 95% confidence level.

This note first summarizes the history of the manuscript "On a Continuum Model for Avalanche Flow and Its Simplified Variants" by Grigorian and Ostroumov―published in the Special Issue on snow avalanche dynamics of Geosciences―since the early 1990s and explains the guiding principles in editing it for publication. The changes are then detailed and some explanatory notes given for the benefit of readers who are not familiar with the early Russian work on snow avalanche dynamics. Finally, the editor's personal views as to why he still considers this paper of relevance for avalanche dynamics research today are presented in brief essays on key aspects of the paper, namely the role of simple and complex models in avalanche research and mitigation work, the status and possible applications of Grigorian's stress-limited friction law, and non-monotonicity of the dynamics of the Grigorian–Ostroumov model in the friction coefficient. A comparison of the erosion model proposed by those authors with two other models suggests to enhance it with an additional equation for the balance of tangential momentum across the shock front. A preliminary analysis indicates that continuous scouring entrainment is possible only in a restricted parameter range and that there is a second erosion regime with delayed entrainment.

Snowpack seasonality in the conterminous United States (U.S.) is explored using a daily,14 km horizontal resolution gridded snow water equivalent and snow depth reanalysis product. I2calculated seasonal snowpacks using two established methods: (1) the classic Sturm approach that3requires 60 days of snow cover with a peak depth >50 cm and (2) the snow seasonality metric (SSM)4that only requires 60 days of continuous snow cover. The latter approach yields continuous values5from -1 to +1, where -1 (+1) indicates an ephemeral (seasonal) snowpack. Both approaches identify6seasonal snowpacks in western mountains and the northernmost central and eastern U.S. By relaxing7the depth constraint and providing continuous values, the SSM identifies greater areas of seasonal8snowpacks compared to the Sturm method, particularly in the upper Midwest, New England, and the9Intermountain West. Ephemeral snowpacks are identified throughout lower elevation regions of the10western U.S. and across a broad swath centered near 35°N spanning the lee of the Rocky Mountains11to the Atlantic coast. Because it lacks a depth constraint, the SSM approach is sensitive to interannual12variability, indicating it may inform the location of shallow but long-duration snowpacks at risk13of transitioning to becoming ephemeral with climatic change. A case study in Oregon during an14extreme snow drought year highlights seasonal to ephemeral snowpack transitions.

The snow melt from the High Atlas is a critical water resource in Morocco. In spite of its importance, monitoring the spatio-temporal evolution of key snow cover properties like the snow water equivalent remains challenging due to the lack of in situ measurements at high elevation. Since 2015, the Sentinel-2 mission provides high spatial resolution images with a 5 day revisit time, which offers new opportunities to characterize snow cover distribution in mountain regions. Here we present a new data assimilation scheme to estimate the state of the snowpack without in situ data. The model was forced using MERRA-2 data and a particle filter was developed to dynamically reduce the biases in temperature and precipitation using Sentinel-2 observations of the snow cover area. The assimilation scheme was implemented using SnowModel, a distributed energy-balance snowpack model and tested in a pilot catchment in the High Atlas. The study period covers 2015-2016 snow season which corresponds to the first operational year of Sentinel-2A, therefore the full revisit capacity was not yet achieved. Yet, we show that the data assimilation led to a better agreement with independent observations of the snow height at an automatic weather station and the snow cover extent from MODIS. The performance of the data assimilation scheme should benefit from the continuous improvements in MERRA-2 reanalyses and the full revisit capacity of Sentinel-2.

The Indus basin is considered as the one with the highest dependence on snowmelt runoff in High Mountain Asia. The recent High Mountain Asia snow reanalysis enables to go beyond previous studies by evaluating both snowmelt and snow sublimation at the basin scale. Over 2000-2016, basin-average snowmelt is 101 11 Gt.a^-1 (121 ± 13 mm.a^-1), which represents about 25-30% of basin-average annual precipitation. Snow sublimation accounts for 11% of the mean annual snow ablation, but with a large spatial variability across the basin.

Cloud seeding is a weather modification technique for enhancing precipitation in arid and semi-arid regions, including Western U.S. However, designing an optimal cloud seeding operation based on comprehensive evaluation metrics, such as seeding agent dispersion and atmospheric conditions, has yet to be thoroughly explored for this region. This study investigated the impacts of cloud seeding, particularly utilizing silver iodide, on ice particle growth within clouds through numerical modeling. By leveraging the Snow Growth Model for Rimed Snowfall (SGMR), the microphysical processes involved in cloud seeding across five distinct events were simulated. The events were in the Lake Tahoe region, nestled within the Sierra Nevada Mountain ranges in the western U.S. This model was executed based on six primary variables, including cloud top height, cloud base height, cloud top temperature, cloud base temperature, liquid water content, and ice water content. This study incorporated datasets from the Modern-Era Retrospective Analysis for Research and Applications Version 2 and the Clouds and the Earth Radiant Energy System. Findings revealed the importance of ice nucleation, aggregation, diffusion, and riming processes and highlighted the effectiveness of cloud seeding in enhancing ice particle number concentration, ice water content, and snowfall rates. However, event-specific analyses indicated diverse precipitation responses to cloud seeding based on initial atmospheric conditions. The SGMR modeling hints at the importance of improving ice microphysical processes and provides insights for future cloud seeding research using regional and global climate models.

The optical signals detected on multiple satellite platforms over snow surfaces are determined by the optical properties of snow surface and atmosphere. The solution of both direct and inverse problems of an atmosphere – underlying snow system requires simple relationships between top-of-atmosphere (TOA) reflectance R and microphysical/optical characteristics of both snow and atmosphere. The task of this paper is to present a simple analytical relationship between the value of R as detected on a satellite with atmosphere/snow properties. Such a relationship can be established using a numerical solution of integro - differential radiative transfer equation (RTE) (Liou, 2022). However, this path is quite complicated and time consuming. The analytical solutions of RTE are needed for the solution of various applied atmospheric and snow optics problems (Cachorro et al., 2022; Mei et al., 2020, 2022; Kokhanovsky, 2021). This is the main driver of this work. To simplify the problem under study we consider the case of Antarctica, where both snow and atmosphere are almost free of pollutants. This work is focused on the simulation of the moderate spectral resolution TOA measurements (1nm or so) and the spectral range 400-1000nm.

The paper presents the first retrievals of clean snow properties using spaceborne hyperspectral The Environmental Mapping and Analysis Program (EnMAP) observations. The location close to the Dome C in Antarctica has been selected. At this location the atmospheric effects except molecular light scattering and absorption are weak and the simplified atmospheric correction scheme can be applied. The ice grain size, snow specific surface area and snow spectral and broadband albedos have been retrieved using single view EnMAP measurements. In addition, we propose the technique to retrieve trace gas concentrations (e.g., water vapor, ozone) from EnMAP observations over the snow surfaces. Close correspondence of satellite and ground measured parameters has been found.

Radiation is the major driver of snowmelt, and hence its estimation is critically important. Net radiation reaching the forest floor is influenced by vegetation density. Previous studies in mid-latitude conifer forests have confirmed that net radiation decreases and then subsequently increases with increasing vegetation density, for clear sky conditions. This leads to existence of a net radiation minimum at an intermediate vegetation density. With increasing cloud cover, the minimum radiation shifts toward lower densities, sometimes resulting in a monotonically increasing radiation with vegetation density. The net radiation trend, however, is expected to change across sites, affecting the magnitude and timing of individual radiation components. This research explores the variability of net radiation on snow-covered forest floor for different vegetation densities along a latitudinal gradient. We especially investigate how the magnitude of minimum/maximum radiation and the corresponding vegetation density change with the site geographical location. To evaluate these, the net radiation is evaluated using the Forest Radiation Model at six different locations in predominantly white spruce (Picea glauca) canopy cover across North America, ranging from 45 to 66°N latitudes. Results show that the variation of net radiation with vegetation density considerably varies with latitude. In higher latitude forests, the magnitude of net radiation is generally smaller, and the minimum radiation is exhibited at relatively sparser vegetation densities, under clear sky conditions. For interspersed cloudy sky conditions, net radiation non-monotonically varies with latitude across the sites, depending on the seasonal sky cloudiness and air temperature. Latitudinal sensitivity of net radiation is lower on north-facing hillslopes than on south-facing sites.

Subseasonal-to-seasonal (S2S) weather forecasting has improved in recent years, thanks partly to better representation of physical variables in models. For instance, realistic initializations of snow and soil moisture in models yield enhanced predictability on S2S time scales. Snow depth and soil moisture also mediate month-to-month persistence of near-surface air temperature. Here the role of snow depth as predictor of temperature one month ahead in the Northern Hemisphere is probed via two causal pathways. Through the first pathway, snow depth anomalies in month 1 cause snow depth anomalies in month 2, which then cause temperature anomalies in month 2. This pathway represents the snow–albedo feedback, as well as cooling due to insulation, emissivity and heat loss. It is active from fall to summer, and its effect peaks in March/April in the midlatitudes and in May/June at high latitudes. A complementary second pathway, where snow depth anomalies in month 1 cause soil moisture anomalies in month 2, which then cause temperature anomalies in month 2 through soil moisture–temperature feedbacks, is only active in spring and summer. Its effect peaks later in the warm season than the effect of the first pathway. Geographically, snow depth mediates north of, and soil moisture south of, the areas with the highest temperature predictability from snow depth. These results indicate that the two pathways describe complementary physical mechanisms. The first pathway embodies month-to-month persistence of snow depth, and the second pathway represents melting of snow from one month to the next.

The dependence of the width of annual rings of woody plants on the characteristics of the snow cover is analyzed in various natural zones of the West Siberian Plain and adjacent mountainous areas: the maximum thickness and water reserve for the entire winter period and for individual months, the dates of disappearing, establishment, and duration of the occurrence of a stable snow cover. It has been shown that the role of the thickness and water content of snow cover for the radial growth of trees is differentiated by geographical location. On the plain, it intensifies in the forest-tundra and dry steppe. The response of radial growth to snow cover in the upper and lower parts of the forest belt is often opposite. Dates of establishment of stable snow cover are more important for tree growth compared to dates of disappearance. Dates of disappearance of stable snow cover are more significant in the southern regions than in the northern ones. The value of the duration of the period with stable snow cover for tree growth is higher in the southern regions.

This study investigates the spatial signatures of seasonal snow in Synthetic Aperture Radar (SAR) observations at different spatial scales and for different physiographic regions. Sentinel-1 C-band (SAR) backscattering coefficients (BSC) were analyzed in the Swiss Alps (SA), in high elevation forest and grasslands in Grand Mesa (GM), Colorado, and in North Dakota (ND) croplands. GM BSC exhibit 10dB sensitivity to wetness at small scales (~100 m) over homogeneous grassland. Sensitivity decreases to 5 dB in the presence of trees, and it is demonstrated that VH BSC sensitivity enables wet snow mapping below the tree-line. Area-variance scaling relationships show minima at ~100 m and 150-250 m respectively in barren and grasslands in SA and GM, increasing up to 1 km and longer in GM forests and ND agricultural fields. The spatial organization of BSC (as described by 1D-directional BSC wavelength spectra) exhibits multi-scaling behavior in the 100 -1,000 m range with a break at (180-360 m) that is also present in UAVSAR L-band measurements in GM. Spectral slopes in GM forested areas steepen during accumulation and flatten in the melting season with mirror behavior for grasslands reflecting changes in scattering mechanisms with snow depth and wetness, and vegetation mass and structure. Overall, this study reveals persistent patterns of SAR scattering variability spatially organized by land-cover, topography and regional winds with large inter-annual variability tied to precipitation. This dynamic scaling behavior emerges as an integral physical expression of snowpack variability that can be used to model sub-km scales and for downscaling applications.

Soil and Water Assessment Tool (SWAT) was used to simulate five glacierized river basins that are global in coverage and vary in climate. The river basins included the Narayani (Nepal), Vakhsh (Central Asia), Rhone (Switzerland), Mendoza (Central Andes, Argentina), and Central Dry Andes (Chile) with a total area of 85,000 km2. A modified SWAT snow algorithm was applied in order to consider spatial variation of associated snow melt /accumulation by elevation band across each subbasin. In the previous studies, melt rates varied as a function of elevation resulting from an air temperature gradient while the snow parameters were constant throughout the entire basin. A major improvement of the new snow algorithm is separating the glaciers from seasonal snow based on their characteristics. Two SWAT snow algorithms were evaluated in simulation of monthly runoff from glaciered watershed: 1) the snow parameters are lumped (i.e. constant throughout the entire basin) and 2) the snow parameters are spatially variable based on elevation band-subbasin (i.e. modified snow algorithm). Applying the distributed SWAT snow algorithm improved the model performance in simulation of monthly runoff with snow-glacial regime, so that mean RSR decreased to 0.49 from 0.55 and NSE increased to 0.75 from 0.69. Improvement of model performance was negligible in simulation of monthly runoff from the basins with monsoon runoff regime.

We present a simplified atmospheric correction algorithm for the snow/ice albedo retrieval using single view satellite measurements. The validation of the technique is performed using Ocean and Land Colour Instrument (OLCI) on board Copernicus Sentinel - 3 satellite and ground spectral or broadband albedo measurements from locations on the Greenland ice sheet and in the French Alps. Through comparison with independent ground observations, the technique is shown to perform accurately in a range of conditions from a 2100 m elevation mid-latitude location in the French Alps to a network of 15 locations across a 2390 m elevation range in seven regions across the Greenland ice sheet. Retrieved broadband albedo is accurate within 5% over a wide (0.5) broadband albedo range of the (N = 4,155) Greenland observations and with no apparent bias.

Seasonal snow cover is a fundamental component of both the global energy budget and the water cycle. Its properties such as fractional snow cover or albedo are particularly affected by climate change. Several methods based on satellite data products are available to estimate these properties, each one with its pros and cons. This work presents a novel methodology that integrates three indexes applied to MODIS satellite data (Spectral Mixture Analysis (SMA), Normalized Difference Snow Index (NDSI) and Melt Area Detection Index (MADI)), to perform a spatio-temporal reconstruction of the fractional snow cover and albedo at 250 m spatial resolution in the Brunswick Peninsula, southwest Patagonia during the cold season (April-October) for the period 2000-2020. Three main steps are included: (1) the increase of the spatial resolution of MODIS (MOD09) data to 250 m using a spectral fusion technique; (2) the snow-cloud discrimination; (3) the daily spatio-temporal reconstruction of snow extent and its albedo signature with subpixel detection using the endmembers extraction and spectral mixture analysis. The results show a 98% agreement between MODIS snow detection and ground-based snow measurement at the automatic weather station Tres Morros (53.3174 °S, 71.2790 °W), with fractional snow cover values between 20% and 50%, showing a close relationship between snow and vegetation type. The number of snow days compiled from the MODIS data indicates a good performance (Pearson correlation of 0.9) compared with the number of skiing days at Cerro Mirador ski centre near Punta Arenas. Although the number of seasonal snow days show a significant increase trend of 0.54 days/year in Brunswick Peninsula during the 2000-2020 period a significant decreasing trend of -4.64 days/year was detected during 2010-2020, and also below the 400 m a.s.l. elevation, which is the most affected area. A reconstruction of the monthly mean temperature using the ERA5 Land reanalysis product shows a significant warming trend in May (0.068 ºC/year) and October (0.098 ºC/year) in the 2000-2020 period. Under the future emission scenario RCP8.5, the regional climate model RegCM4 predicts further warming during 2020-2050 of 0.059 ºC/year in July, 0.088 ºC/year in August, and 0.019 ºC/year in October, which will further reduce snow cover.

Snowpack loss in midlatitude mountains is ubiquitously projected by Earth system models, though the magnitudes, persistence and time horizons of decline vary. Using daily downscaled hydroclimate and snow projections we examine changes in snow seasonality across the U.S. Pacific Southwest region during a simulated severe 20-year dry spell in the 21st century (2051–2070) developed as part of the 4th California Climate Change Assessment to provide a "stress test" for water resources. Across California’s mountains, substantial declines (30–100% loss) in median peak annual snow water equivalent accompany changes in snow seasonality throughout the region compared to the historic period. We find 80% of historic seasonal snowpacks transition to ephemeral conditions. Subsetting empirical-statistical wildfire projections for California by snow seasonality transition regions indicates a two-to-four-fold increase in burned area, consistent with recent observations of high elevation wildfires following extended drought conditions. By analyzing six of the major California snow-fed river systems we demonstrate snowpack reductions and seasonality transitions result in concomitant declines in annual runoff (47-58% of historic values). The negative impacts to statewide water supply reliability by the projected dry spell will likely be magnified by changes in snowpack seasonality and increased wildfire activity.

We have proposed a simple algorithm to retrieve the total ozone column and snow properties (spectral albedo and effective light absorption path) using the high spatial resolution single – view MSI/S-2 measurements over Antarctica.

Accordingly, the present study is aimed at investigating quaternary climate changes in Binalod Heights. To identify glacial effects, Morphic indices, field evidence and effects, climatic evidence, and (laboratory) experimental analysis were employed. Determining the permanent snow line in the region was conducted using the Right Method and 65 cirques which are considered as much enriched feeding resources for the formation an ice cover in the region. The expansion of settlements in the region are lower than the permanent snow life is more accumulated than above the border of the permanent snow line. This issue indicates that refrigeration cells do not have the ability to create civil nuclear. Regarding quaternary climate changes and the gradual warming of the climate, the initial core of the City of Mashhad ranges from the center of the Kashf rood River to northeastern heights of Binalod. In addition, the existence of glacial cirques in heights as an important factor in feeding refrigerating conditions has been effective on the expansion of urbanization of Mashhad in the past time. Our new geomorphological mapping and landsystem reconstructions provide an important insight into the response of temperate Binalod glaciers to rapidly-warming climate

The role of Black Carbon (BC) as a contributor to glacial retreat is of particular interest to the scientific community and decision makers, due to its impact on snow albedo and glacier melt. In this study, a thermal-optical instrument (LAHM) was used to measure effective Black Carbon (eBC) in a series of surface snow samples collected from the Vallunaraju glacier, Cordillera Blanca, between April 2019 and May 2020. The time series obtained indicates a marked seasonal variability of eBC with maximum concentrations during the dry season and dramatic decrease during the wet season. The concentrations detected ranged between a minimum of 3.73 ng/g and 4.23 ng/g during the wet season and a maximum of 214.13 ng/g and 181.60 ng/g during the dry season, in the accumulation and ablation zone. Using SNICAR model, the reduction of albedo was estimated at 6.36% and 6.60% during the dry season and 0.68% and 0.95% during the wet season, which represents an average radiative forcing of 4.52 ± 1.84 W/m2 and 4.69 ± 1.59 W/m2 in the accumulation zone, and 0.49 ± 0.27 W/m2 and 0.68 ± 0.43 W/m2 in the ablation zone. The melting of snow due to the eBC translates into 80.18 ± 37.30 kg/m2 and 83.16 ± 32.75 kg/m2 during the dry season, and 7.91 ± 4.29 kg/m2 and 10.85 ± 6.62 kg/m2 during the wet season, in the accumulation and ablation zones, respectively. Finally, the HYSPLIT trajectory assessment shows that aerosols predominate in the Amazon rainforest, especially when forest fires are most abundant according to VIIRS images.

This document describes the theoretical basis of the algorithms used to determine properties of snow and ice from the measurements of the Ocean and Land Color Instrument (OLCI) onboard Sentinel-3 satellites within the Pre-operational Sentinel-3 snow and ice products (SICE) project: http://snow.geus.dk/. The code used for the SICE retrieval and its documentation can be found at https://github.com/GEUS-SICE/pySICE. The algorithms were developed after the work from Kokhanovsky et al. (2018, 2019, 2020).

Knowledge of the spatio-temporal occurrence of avalanche activity is critical for avalanche forecasting and hazard mapping. We present a near-real time automatic avalanche monitoring system that outputs detected avalanche polygons within roughly 10 min after Sentinel- 1 SAR data download. Our avalanche detection algorithm has an average probability of detection of 67.2 % with a false alarm rate averaging 45.9, with maximum POD's over 85 % and minimum FAR's of 24.9 % compared to manual detection of avalanches. The high variability in performance stems from the dynamic nature of snow in the Sentinel-1 data. After tuning parameters of the detection algorithm, we processed five years of Sentinel-1 images acquired over a 150 x 100 km large area in Northern Norway, with the best setup. Compared to a dataset of field-observed avalanches, 77.3 % were manually detectable. Using these manual detections as benchmark, the avalanche detection algorithm achieved an accuracy of 79 % with high POD's in cases of medium to large wet snow avalanches. For the first time, we can present a dataset of spatiotemporal avalanche activity over several winters from a large region. This unique dataset allows for research into the relationship between avalanche activity and triggering meteorological factors, mapping of avalanche prone areas and near-real time avalanche activity monitoring to assist public avalanche forecasting. Currently, the Norwegian Avalanche Warning Service is using our processing system for pre-operational use in three regions in Norway.

Snow-covered regions are the main source of reflection of incident shortwave radiation (ISR) on Earth’s surface. The deposition of light-absorbing particles (LAPs) on these regions increases the capacity of snow to absorb ISR and decreases surface snow albedo (SSA), which intensifies the radiative forcing leading to accelerated snowmelt and modifications of the hydrologic cycle. In this work we investigate changes in SSA and radiative forcing (RF) induced by LAPs in the Upper Aconcagua River Basin (Chilean central Andes) using remote sensing satellite data (MODIS), in-situ spectral SSA measurements, and the ISR (Chilean Solar Explorer platform) during the austral-winter months (May to August) for the 2004-2016 period. To estimate the changes in SSA and RF, we define two spectral ranges: i) an enclosed range (Ecr) between 841-876 nm, which isolates effects of Black-Carbon, an important LAP derived from anthropogenic activities, and ii) a broadband range (Bbr) between 300-2500 nm. Our results show that percent variations in SSA in the Ecr are higher than in the Bbr, regardless of the total amount of radiation received, which may be attributed to the presence of LAPs as these particles have a greater impact on SSA at wavelengths in the Ecr band than in the Bbr band.

The impact of snow darkening by deposition of light absorbing aerosols (LAAs) on snow cover over the Himalaya-Tibetan-Plateau (HTP) and influence on the Asian monsoon are investigated using the NASA Goddard Earth Observing System Model Version 5 (GEOS-5). We find that during April-May-June, deposition of LAAs on snow leads to a reduction in surface albedo, initiating a sequence of feedback processes, starting with increased surface solar radiation, rapid snowmelt in HTP and warming of the surface and upper troposphere, followed by enhanced low-level southwesterlies and increased dust loading over the Himalayas-Indo-Gangetic Plain. The warming is amplified by increased dust aerosol heating, and subsequently amplified by latent heating from enhanced precipitation over the Himalaya foothills and northern India, via the Elevated Heat Pump (EHP) effect during June-July-August. The reduced snow cover in the HTP anchors the enhanced heating over the Tibetan Plateau and its southern slopes, in conjunction with an enhancement of the Tibetan Anticyclone, and the development of an anomalous Rossby wavetrain over East Asia, leading to weakening of the subtropical westerly jet, and northward displacement and intensification of the Mei-Yu rainbelt. Our results suggest that atmosphere-land heating by LAAs, particularly desert dust play a fundamental role in physical processes underpinning the snow-monsoon relationship proposed by Blandford more than a century ago.

In this work, we present two datasets for specific areas located in South Tyrol (Italy) and in Tyrol (Austria) that can be exploited to monitor and understand water resource dynamics in mountain regions. The idea is to provide the reader with information about the different sources of water supply over five defined test areas. The Snow Cover Fraction (SCF) and Soil Moisture Content (SMC) datasets are derived from machine learning algorithms based on remote sensing data. Both SCF and SMC products are characterized by a spatial resolution of 20 m and are provided for the period from October 2020 to May 2023 (SCF) and from October 2019 to September 2022 (SMC) respectively, each covering 3 seasons of interest, winter for SCF and spring-summer for SMC. For SCF maps, the validation with very high-resolution images shows high correlation coefficients of around 0.9. The SMC products were originally produced with an algorithm validated at global scale, but here, to obtain more insight in the specific alpine mountain environment, the values estimated from the maps are compared with ground measurements of automatic stations located at different altitudes and characterized by different aspects in the Val Mazia catchment in South Tyrol (Italy). In this case a MAE between 0.05 and 0.08 and an unbiased RMSE between 0.05 and 0.09 m3·m-3 were achieved. The datasets presented can be used as input for hydrological models as well as to hydrologically characterize the study alpine area starting from different sources of information.

In this work, backscattering signatures of snow-covered lake ice and sea ice from X- and Ku-band synthetic aperture radar (SAR) data are investigated. The SAR data were acquired with the ESA airborne SnowSAR sensor in winter 2012 over Lake Orajärvi in northern Finland and over landfast ice in the Bay of Bothnia of the Baltic Sea. Co-incident with the SnowSAR acquisitions in-situ snow and ice data were measured. In addition, time series of TerraSAR-X images and ice mass balance buoy data were acquired for the Lake Orajärvi in 2011-12. The main objective of our study was to investigate relationships between SAR backscattering signatures and snow depth over lake and sea ice, with the ultimate objective of assessing the feasibility of retrieval of snow characteristics using X- and Ku-band dual-polarisation (VV and VH) SAR over freshwater or sea ice. This study constitutes the first comprehensive survey of snow backscattering signatures at these two combined frequencies over both lake and sea ice. For lake ice, we show that X-band VH-polarized backscattering coefficient (σo) and the Ku-band VV/VH-ratio exhibited the highest sensitivity to the snow depth. For sea ice, the highest sensitivity to the snow depth was found from the Ku-band VV-polarised σo and the Ku-band VV/VH-ratio. However, the observed relations were relatively weak, indicating that at least for the prevailing snow conditions, obtaining reliable estimates of snow depth over lake and sea ice would be challenging using only X- and Ku-band backscattering information.

Snow parameters have traditionally been retrieved using discontinuous, multi-band sensors; however, continuous hyperspectral sensor are now being developed as an alternative. In this paper we investigate the performance of various sensor configurations using machine learning neural networks trained on a simulated dataset. Our results show improvements in accuracy of retrievals of snow grain size and impurity concentration for continuous hyperspectral channel configurations. Retrieval accuracy of snow albedo was found to be similar for all channel configurations.

Accurate measurement of solid precipitation (S) has a critical importance for proper understanding of the Earth’s hydrological cycle, validation of emerging technologies and weather prediction models, and developing parameterizations of severe weather elements such as visibility (Vis). However, measuring S is still a challenging problem mainly because of wind effects. The wind effects are normally mitigated by using a Double-Fence Automated Reference (DFAR) system to reduce the wind speed (Ug). To contribute towards addressing some of these issues we have analyzed data sets collected at a site located in Southern Ontario, Canada using several instruments. The instruments include two Geonor gauges, one placed inside a DFAR (SDFAR) and the other inside a double Alter shield (DASG), a Pluvio2 gauge inside a single Alter shield (SASP), a HotPlate, a PARSIVEL2 Disdrometer that measures S and fall velocity (V), and a FD12P senor that measures S and type and Vis. The results show that for the Ug observed in this study (Ug < 6 ms-1), both DASG and SASP have similar collection efficiency (CE) of near 70%. The transfer functions (TF) for DASG and SASP as a function of Ug and also Ug, and V have been derived. The TF for the DASG that includes both Ug and V showed better agreement with observation than just Ug alone. The S measured using all the other instruments were correlated well with SDFAR, but the PARSIVEL2 and FD12P overestimated and underestimated the snow amount respectively as compared the SDFAR. However, the HotPlate captured similar amount of S as the SDFAR. According to this study, the SDFAR showed good correlation with Vis.

Snow cover plays a critical role in the global energy and water cycles. Snow cover on the Tibetan Plateau (TP) provides vital water sources in western China and Himalayan regions in addition to its weather and climate significance. The massive high mountain topography of the TP is main conditions for the presence and persistence of snow cover on the plateau at the mid-low latitudes of the Northern Hemisphere (NH). However, how mountain topography controls snow cover distribution on the TP is largely remain unclear and the relationship is not well quantified. Here, the spatial distribution of snow cover and topographic controls on snow cover on the TP are examined based on snow cover frequency (SCF) derived from MODIS snow cover product (MOD10A2 v5) and digital elevation model (DEM). The results show that snow cover on the TP is spatially unevenly distributed and is characterized by rich snow and high SCF on the interior and surrounding high mountain ranges, and less snow and low SCF in inland basins and river valleys. Snow cover on the TP presents elevation dependence with the higher the altitude, the higher the SCF, the longer the snow cover duration and the more stable the intra-annual variation. Annual mean SCF below 3000 m above sea level (masl) is less than 4% and it reaches 77% above 6000 masl. The intra-annual snow cover variation below 4000 masl features a unimodal distribution, while above 4000 masl it presents a bimodal distribution. The mean minimum SCF below 6000 masl occurs in summer, while above 6000 masl it occurs in winter. Because of differences in solar radiation and water vapor sources, mean SCF generally increases with mountain slopes and it is the highest on the north-facing aspect, whereas the lowest is observed on the south-facing aspect.

We present the update of the Snow and Ice (SICE) property retrieval algorithm proposed initially by Kokhanovsky et al. (2019). The algorithm is based on the spectral measurements of Ocean and Land Color Instrument (OLCI) onboard Sentinel-3 satellites combined with the asymptotic radiative transfer theory valid for weakly absorbing turbid media. The main improvements include the introduction of a new atmospheric correction, retrieval of snow impurity load and properties, retrievals for partially snow-covered ground and also accounting for various thresholds to be used to assess the retrieval quality. The algorithm is available as python and Fortran packages at https://github.com/GEUS-SICE/pySICE. The technique can be applied to various optical sensors (satellite and ground-based) operated in the visible and near infrared regions of electromagnetic spectra.

Stretching along the border of North Dakota and Minnesota, The Red River Valley (RRV) of the North has the highest frequency of reported blizzards within the contiguous United States. Despite the numerous impacts these events have, few systematic studies exist discussing the meteorological properties of blizzards. As a result, forecasting these events and lesser blowing snow events is an ongoing forecast challenge. This study presents a climatology of atmospheric patterns associated with RRV blizzards for the winter seasons of 1979-1980 to 2017-2018. Patterns were identified using subjective and objective techniques using meteorological fields from the North American Regional Reanalysis (NARR). The RRV experiences on average, 2.6 events per year. Blizzard frequency is bimodal with peaks occurring in December and March. The events can largely be typed into four meteorological categories dependent on the forcing that drives the blizzard: Alberta Clippers, Arctic Fronts, Colorado Lows, and Hybrids. Objective classification of these blizzards using a competitive neural network known as the Self-Organizing Map (SOM) demonstrates that gross segregation of the events can be achieved with a small (8-class) map. This implies that objective analysis techniques can be used to identify these events in weather and climate model output that may aid future forecasting and risk assessment projects.

On 5-7 April 2018 a landfalling atmospheric river resulted in widespread heavy precipitation in the Sierra Nevada of California and Nevada. Observed snow levels during this event were among the highest snow levels recorded since observations began in 2002 and exceeded 2.75 km for 31 hours in the northern Sierra Nevada and 3.75 km for 12 hours in the southern Sierra Nevada. The anomalously high snow levels and over 80 mm of precipitation caused flooding, debris flows, and wet snow avalanches in the upper elevations of the Sierra Nevada. The origin of this atmospheric river was super typhoon Jelawat, whose moisture remnants were entrained and maintained by an extratropical cyclone in the northeast Pacific. This event was notable due to its April occurrence, as six other typhoon remnants that caused heavy precipitation with high snow levels (mean = 2.92 km) in the northern Sierra Nevada all occurred during October.

Communication pathways are an underestimated source of synthetic particles in the environment. This study investigated the impact of traffic volume on microplastics (MPs) and tire wear particles (TWPs) pollution in road snow. An examination was conducted in a medium-sized city situated in northeastern Poland, known for being one of the cleanest regions in the country. MPs and TWPs were found at all 53 sites, regardless of the intensity of traffic. The average concentration for all samples was 354.72 pcs/L. Statistically significant differences were found between the average values of particles concentration on low, medium, and heavy traffic roads amounted 62.32 pcs/L, 335.97 pcs/L and 792.76 pcs/L respectively. Within all three studied groups of roads, MPs and TWPs with the smallest size, ranging from 50 to 200 μm, were prevalent. In all studied groups of roads, 4 analyzed shapes of particles were found, with irregular fragments being the most abundant form (89.23%). The most frequently recorded color among the collected samples was black (99.85%), and the least frequently recorded was blue, constituting only 0.01%. This study suggests that snow cover on the roads may act like temporary storage of pollutants during winter in the temperate climate zone and after thaw can significantly increase the concentration of MPs and TWPs in surface waters. Possible measures to decrease the release of MPs and TWPs into the environment in the city may include reducing traffic volume and speed, implementing street sweeping, utilizing filtration chambers, and installing stormwater bioretention systems or settling ponds.

With advances in the development of autonomous vehicles (AVs), more attention has been paid to the effects caused by adverse weather conditions on them. It is well known that the performance of self-driving vehicles is reduced when they are exposed to stressors that impair visibility or generate water or snow accumulation on sensor surfaces. This paper proposes a model to quantify weather precipitation, such as rain and snow, perceived by moving vehicles based on outdoor data. The modelling covers a wide range of parameters, such as varying wind direction and realistic particle size distributions. The model allows the calculation of precipitation intensity on inclined surfaces of different orientations and on a circular driving path. The modelling results were partially validated against direct measurements carried out by a test vehicle. The model outputs showed strong correlation with the experimental data for both rain and snow. Mitigation strategies for heavy precipitation on vehicles can be developed and correlations between precipitation rate and accumulation level can be traced using the presented analytical model. Dimensional Analysis of the problem has highlighted the critical parameters that can help the design of future experiments. The obtained results highlight the importance of the angle of the sensing surface on the perceived precipitation level. The proposed model is used to analyze optimal orientations for minimization of the precipitation flux, which can help to determine the positioning of sensors on the surface of autonomous vehicles.

Snow cover products are primarily derived from the Moderate-resolution Imaging Spectrometer (MODIS) and Advanced Very High-Resolution Radiometer (AVHRR) datasets. MODIS achieves both snow/non-snow discrimination and snow cover fractional retrieval, while early AVHRR-based snow cover products only focused on snow/non-snow discrimination. The AVHRR Climate Data Record (AVHRR -CDR) provides a nearly 40-year global dataset that has a potential to fill a gap in long term snow cover fractional monitoring. Our study selects the Qinghai-Tibet Plateau as the experimental area, utilizing AVHRR-CDR surface reflectance data (AVH09C1) and calibrating with MODIS snow product MOD10A1. The snow cover percentage retrieval from the AVHRR dataset is performed using Surface reflectance at 0.64μm (SR1) and Surface reflectance at 0.86μm (SR2), along with a simulated Normalized Difference Snow Index (NDSI) model. Also, in order to detect the effects of land cover type and topography on snow inversion, we tested the accuracy of the algorithm with and without these influences, respectively (vanilla algorithm & improved algorithm). The accuracy of the AVHRR snow cover percentage data product is evaluated using MOD10A1、 ground snow-depth measurements and ERA5. Results indicate that the logic model based on NDSI has the best fitting effect, with R-square and RMSE values of 0.83 and 0.10, respectively. Meanwhile, the accuracy was improved after taking into account the effects of land cover type and topography. The model is validated using MOD10A1 snow-covered areas, showing snow cover area differences of less than 4% across 6 temporal phases. Improved algorithm results in better consistency with MOD10A1 than vanilla algorithm. Moreover, the RMSE reaches greater levels when the elevation is below 2,000 meters or above 6,000 meters, and is lower when the slope is between 16° and 20°. Using ground snow-depth measurements as ground truth, the multi-year recall rates are mostly above 0.7, with an average recall rate of 0.81. Results also show a high degree of consistency with ERA5. The validation results demonstrate that the AVHRR snow cover percentage remote sensing product proposed in this study exhibits high accuracy in the Tibetan Plateau region, also demonstrating that land cover type and topographic factors are important to the algorithm. Our study lays the foundation for the global snow cover percentage product based on AVHRR-CDR，furthermore lays a basic work for generating a long-term AVHRR-MODIS fractional snow cover dataset.

Snow crab replaced groundfish fisheries as the primary focus of the Newfoundland and Labrador (NL) fishing industry in the mid-1990s. Over the past three decades, management philosophies in this male-only fishery have shifted from promoting industry expansion to promoting industry rationalization to implementing Precautionary Approach (PA) management. Until the mid-2010s, there were regionalized management strategies characterized by higher exploitation rates in northern than southern sub-stock units along the NL marine shelves. However, in the late 2010s, exploitation rates were permitted to elevate to high levels across the entire stock range in association with a large resource decline and evidence of biological harm through fishing emerged in chronically depleted areas. In 2019, when stock and fishery productivity were near historic lows, a multi-indicator PA system was informally introduced into management of the resource. This event coincided with anticipated improvements in stock and fishery performance. This paper examines the extent to which the multi-indicator PA management system, focused on promotion of both biological protection and maximization of fishing efficiency, may have contributed to recent improvements in stock and fishery status. A suite of indicators from areas implementing PA guidance into management in 2019 prior to formal adoption of the PA system in 2023 suggest the system has been highly beneficial in promoting rapid recovery, improving stock and fishery performance metrics to levels matching or exceeding historic levels under similar conditions. We discuss capacity of the system to better safeguard biological aspects of resource and fisheries productivity moving forward under a scenario of an expectant resource decline.

Seasonal snowpacks, characterized by their snow water equivalent (SWE), play a major role in the hydrological cycle, snow melt contributions to floods and the subsequent availability of water resources downstream. Accurately estimating SWE and understanding its spatial and temporal variations presents a considerable challenge, particularly within mountainous regions characterized by complex terrain and limited observational data. Seeking to enhance the performance of the widely used Soil and Water Assessment Tool (SWAT), we report a new approach characterising snowpack behaviour incorporating both modelled and remotely sensed derived SWE calibration data. We focus on the Chenab River Basin (CRB) a headwater catchment of the Indus Basin, globally significant in terms of human inhabitants and intensifying flood risk due to climate change. We conducted a thorough assessment of five satellite-derived and reanalysis-based precipitation datasets: ERA5-Land, CMORPH, TRMM, APHRODITE, and CPC UPP. This reveals significant levels of uncertainty in global precipitation products when compared to reference data from observed stations as well as in the resulting simulated streamflow from the SWAT model. Subsequently, we expanded the scope of the SWAT model to encompass the spatial and temporal simulation of SWE. This was achieved by incorporating information from remotely sensed and modelled SWE products, manually adjusting snow parameters in R-SWAT for both the main basin and at sub-basin scales. Integrating SWE from reference snow products into the calibration process, alongside streamflow data, substantially enhanced modelling accuracy to simulate SWE compared to the conventional auto-calibration and single-variable approaches reliant solely on streamflow data. This approach results in considerable improvement in SWE predictions and to some extent in streamflow simulation in catchments dominated by snow. This research highlights the potential of remote sensing and modelled SWE parameterisation in the absence of in-situ snowpack data in high-altitude environments. An improved understanding of SWE behaviour is vital for predicting hydrological responses spanning hazards to water resources in the populous downstream regions of the Indus Basin, especially in the face of climate change.

The land-surface developments of the European Centre for Medium-range Weather Forecasts (ECMWF) are based on the Carbon-Hydrology Tiled Scheme for Surface Exchanges over Land (CHTESSEL) and form an integral part of the Integrated Forecasting System (IFS), supporting a wide range of global weather, climate and environmental applications. In order to structure, coordinate and focus future developments and benefit from international collaboration in new areas, a flexible system named ECLand which would facilitates modular extensions to support numerical weather prediction (NWP) and society-relevant operational services, e.g. Copernicus, is presented . This paper introduces recent examples of novel ECLand developments on (i) vegetation, (ii) snow, (iii) soil, (iv) open water/lake (v) river/inundation, and (vi) urban areas. The developments are evaluated separately with long-range, atmosphere-forced surface offline simulations, and coupled land-atmosphere-ocean experiments. This illustrates the benchmark criteria for assessing both, process fidelity with regards to land surface fluxes and reservoirs of the water-energy-carbon exchange on the one hand, and on the other hand the requirements of ECMWF’s NWP, climate and atmospheric composition monitoring services using an Earth system assimilation prediction framework.

Oscillations in global modes of variability (MoV) form global teleconnections that affect regional climate variability and modify the potential for severe and damaging weather conditions. Understanding the link between certain MoVs and regional climate can improve the ability to more accurately predict environmental conditions that impact human life and health. In this study, we explore the connection between different MoVs, including the Arctic Oscillation (AO), Eurasian teleconnection, Indian Ocean Dipole (IOD), North Atlantic Oscillation (NAO), and El Niño Southern Oscillation (Nino34), with winter and summer climate in the High Mountain Asia (HMA) region, including geopotential height at 250 hPa (z250), 2-m air temperature (T2M), total precipitation (PRECTOT), and fractional snow cover area (fSCA). Relationships are explored for the same monthly period between the MoVs and the climate variables, but also using a lagged correlation analysis to investigate whether any relationship exists at different time lags. We find that T2M has a negative correlation with the Eurasian teleconnection in the Inner Tibetan Plateau and Central China in both winter and summer and a positive correlation in Western China in summer. PRECTOT has a positive correlation with all MoV in most regions in winter, especially with the IOD, and a negative correlation in summer, especially with the Eurasian teleconnection. Snow cover in winter is positively correlated with most indices throughout many regions in HMA, likely due to wintertime precipitation also being positively correlated with most indices. Generally, the AO and NAO show similar correlation patterns with all climate variables, especially in the winter, possibly due to their oscillations being so similar. Furthermore, the AO and NAO are shown to be less significant in explaining the variation in HMA climate compared to other MoVs such as the Eurasian teleconnection. Overall, our results pinpoint different time-windows and specific regions within HMA that exhibit high correlation between climate and MoVs, which might offer additional predictability of the MoVs as well as of climate and weather patterns in HMA and throughout the globe.

This article explores the cross-cultural significance of three iconic fairy tales: Cinderella, Snow White, and Sleeping Beauty. Using a comparative analysis approach, the study investigates recurring themes, motifs, and symbolism found in these tales across different cultures and time periods. The research draws on primary and secondary sources such as literary works and scholarly articles to uncover the underlying messages that resonate with audiences worldwide. The historical and cultural contexts of the tales are also examined to provide a deeper understanding of their significance in shaping cultural identity and values. The findings of the research demonstrate the enduring legacy of these tales and their continued relevance in contemporary society. The study highlights how the tales have evolved and adapted over time to reflect changing cultural values and beliefs. The interpretation of the findings offers insights into the role of storytelling in shaping cultural identity and values, and the ways in which these stories continue to inspire and captivate audiences worldwide. Overall, this research provides a comprehensive and nuanced understanding of the significance of these three fairy tales across cultures and time. By illuminating the similarities and differences in these tales, the study offers a fresh perspective on the power of storytelling and its ability to transcend cultural boundaries. The research contributes to the fields of folklore, cultural studies, and cross-cultural communication, inspiring further exploration into the captivating world of fairy tales.

In winter sports, the equipment often comes into contact with snow or ice, and this contact generates a force that resists motion. In some sports, such as cross-country skiing, this resistive force can significantly affect the outcome of a race, as a small reduction in this force can give an athlete an advantage. Researchers have examined the contact between skis and snow in detail, and to fully understand this friction, the entire ski must be studied at various scales. At the macro scale, the entire geometry of the ski is considered and the apparent contact between the ski and the snow is considered and at the micro scale the contact between the snow and the ski base ski-base textures. In the present work, a method for characterising contact between the ski-base texture and virtual snow will be presented. Six different ski-base textures will be considered. Five of them are stone-ground ski bases, and three of them have linear longitudinal textures with a varying number of lines and peak-to-valley height, and the other two are factory-ground “universal” ski bases. The sixth ski base has been fabricated by a steel-scraping procedure. In general, the results show that a ski base texture with a higher Spk-value has less real contact area, and that the mutual differences can be large for surfaces with similar Sa-values. The average interfacial separation is, in general, correlated with the Sa-value, where a “rougher” surface exhibits a larger average interfacial separation. The results for the reciprocal average interfacial separation, which is related to the Couette type of viscous friction, were in line with the general consensus that a “rougher” texture performs better at high speed than a “smoother” one, and it was found that a texture with high Sa and Spk values resulted in a low reciprocal average interfacial separation and consequently low viscous friction. The reciprocal average interfacial separation was found to increase with increasing real contact area, indicating a correlation between the real area of contact and the Couette part of viscous friction.

Analysis of the hourly Canadian Prairie data for the past 60 years has transformed our quantitative understanding of land-atmosphere-cloud coupling. The key reason is that trained observers made hourly estimates of opaque cloud fraction that obscures the sun, moon or stars, following the same protocol for 60 years at all stations. These 24 daily estimates of opaque cloud data are of sufficient quality that they can be calibrated against Baseline Surface Radiation Network data to give the climatology of the daily short-wave, longwave and total cloud forcing (SWCF, LWCF and CF). This key radiative forcing has not been available previously for climate datasets. Net cloud radiative forcing reverses sign from negative in the warm season to positive in the cold season, when reflective snow reduces the negative SWCF below the positive LWCF. This in turn leads to a large climate discontinuity with snow cover, with a systematic cooling of 10°C or more with snow cover. In addition, snow cover transforms the coupling between cloud cover and the diurnal range of temperature. In the warm season, maximum temperature increases with decreasing cloud, while minimum temperature barely changes; while in the cold season with snow cover, maximum temperature decreases with decreasing cloud and minimum temperature decreases even more. In the warm season, the diurnal ranges of temperature, relative humidity, equivalent potential temperature and the pressure height of the lifting condensation level are all tightly coupled to opaque cloud cover. Given over 600 station-years of hourly data, we are able to extract, perhaps for the first time, the coupling between cloud forcing and the warm season imbalance of the diurnal cycle; which changes monotonically from a warming and drying under clear skies to a cooling and moistening under cloudy skies with precipitation. Because we have the daily cloud radiative forci, which is large, we are able to show that the memory of water storage anomalies, from precipitation and the snowpack, goes back many months. The spring climatology shows the memory of snowfall back through the entire winter, and the memory in summer goes back to the months of snowmelt. Lagged precipitation anomalies modify the thermodynamic coupling of the diurnal cycle to the cloud forcing, and shift the diurnal cycle of mixing ratio which has a double peak. The seasonal extraction of the surface total water storage is a large damping of the interannual variability of precipitation anomalies in the growing season. The large land-use change from summer fallow to intensive cropping, which peaked in the early 1990s, has led to a coupled climate response that has cooled and moistened the growing season, lowering cloud-base, increasing equivalent potential temperature, and increasing precipitation. We show a simplified energy balance of the Prairies during the growing season and its dependence on reflective cloud.

Snow is a vital environmental parameter and dynamically responsive to climate change, particularly in mountainous regions. Snow cover can be monitored at variable spatial scales using Earth Observation (EO) data. Long-lasting remote sensing missions enable the generation of multi-decadal time series and thus the detection of long-term trends. However, there have been few attempts to use these to model future snow cover dynamics. In this study, we therefore explore the potential of such time series to forecast the Snow Line Elevation (SLE) in the European Alps. We generate monthly SLE time series from the entire Landsat archive (1985-2021) in 43 Alpine catchments. Positive long-term SLE change rates are detected, with the highest rates (5-8 m/y) in the Western and Central Alps. We utilize this SLE dataset to implement and evaluate seven uni-variate time series modeling and forecasting approaches. The best results were achieved by Random Forests, with a Nash-Sutcliffe efficiency (NSE) of 0.79 and a Mean Absolut Error (MAE) of 258 m, Telescope (0.76, 268 m), and seasonal ARIMA (0.75, 270 m). Since the model performance varies strongly with the input data, we developed a Combined forecast based on the best performing methods in each catchment. This approach was then used to forecast the SLE for the years 2022-2029. In the majority of the catchments the shift of the forecast median SLE level retained the sign of the long-term trend. In cases where a deviating SLE dynamic is forecast a discussion based on the unique properties of the catchment and past SLE dynamics is required. In the future, we expect major improvements in our SLE forecasting efforts by including external predictor variables in a multi-variate modeling approach.

Falling snow alters its own microwave signatures when it begins to accumulate on the ground, making retrieval of snowfall challenging. This paper investigates the effects of snow-cover depth and cloud liquid water content on microwave signatures of terrestrial snowfall using reanalysis data and multi-annual observations by the Global Precipitation Measurement (GPM) core satellite with particular emphasis on the 89 and 166 GHz channels. It is found that over shallow snow cover (snow water equivalent (SWE) ≤ 100 kg m^-2) and low values of cloud liquid water path (LWP 100–150 g m^-2), the scattering of light snowfall (intensities ≤ 0.5 mm h⁻¹) is detectable only at frequency 166 GHz, while for higher snowfall rates, the signal can also be detected at 89 GHz. However, when SWE exceeds 200 kg m^-2 and the LWP is greater than 100–150 g m^-2, the emission from the increased liquid water content in snowing clouds becomes the only surrogate microwave signal of snowfall that is stronger at frequency 89 than 166 GHz. The results also reveal that over high latitudes above 60°N where the SWE is greater than 200 kg m^-2 and LWP is lower than 100–150 g m^-2, the snowfall microwave signal could not be detected with GPM without considering a priori data about SWE and LWP. Our findings provide quantitative insights for improving retrieval of snowfall in particular over snow-covered terrain.

Quinacridone (QA) and 2,9-dimethylquinacridone (DQA) are synthetic substances suitable as a hard, abrasion-resistant, self-organizing gliding layer on ice, snow, and water. For sustainable use, a large number of parameters must be considered to demonstrate that these non-biogenic substances and their by-products and degradation products are harmless to humans and the environment in the quantities released. For this task, available experimental data are used and supplemented for all tautomers by numerous relevant physical, chemical, toxicological and ecotoxicological estimated values based on various Quantitative Structure Activity Relationship (QSAR) methods. On the one hand, the low solubility of QA and DQA leads to stable gliding layers and thus low abrasion and uptake by plants, animals, and humans. On the other hand, the four hydrogen bond forming functional groups per molecule allow nanoparticle decomposition and enzymatic degradation in natural environments. All available data justify a sustainable use of QA and DQA as a gliding layer. The assessment of the toxicological properties is complemented with an investigation of the size and morphology of DQA particles as well as field tests indicating an excellent performance as a gliding layer on snow.

Various products of the Integrated Multisatellite Retrievals for GPM (IMERG) and passive mi-crowave (PMW) sensors are assessed with respect to near-surface wet-bulb temperature (Tw), precipitation intensity, and surface type (i.e., with and without snow and ice on the surface) over the CONUS and using Stage-IV product as reference precipitation. IMERG products include precipitation estimates from infrared (IR), combined PMW, and their combination. PMW products generally have higher skills than IR over snow- and ice-free surfaces. Over snow- and ice-covered surfaces (1) PMW products (except AMSR-2) show a higher correlation coefficient than IR, (2) IR and PMW precipitation products tend to overestimate precipitation, but at colder temperatures (e.g., Tw<-10oC) PMW products tend to underestimate and IR product continues to show large overestimations, and (3) PMW sensors show higher overall skill in detecting precipitation oc-currence, but not necessarily at very cold Tw. The results suggest that the current approach of IMERG (i.e., replacing PMW with IR precipitation estimates over snow- and ice-surfaces) may need to be revised.

Alpine forest gaps can distribute snowfall, solar radiation and rainfall, thus inducing a heterogeneous hydrothermal microenvironment between the inside and outside areas of forest gaps. Additionally, the characteristics of the heterogeneous microenvironment could vary greatly across the gap location properties during winter and the growing season. To determine the response of total phenol loss (TPL) from the litter to alpine forest gap disturbance during decomposition, we conducted a field litterbag experiment within a representative fir (Abies faxoniana Rehd.) forest based on the gap location properties. The TPL and abundances of fungi and bacteria from two typical shrub species (willow, Salix paraplesia Schneid., and bamboo, Fargesia nitida (Mitford) Keng f.) were measured during the following periods over two years: snow formation (SF), snow cover (SC) snow melting (ST), the early growing season (EG) and the later growing season (LG). At the end of the study, we found that the snow cover depth, frequencies of the freeze-thaw cycle and the fungal copy g-1 to bacterial copy g-1 ratio had significant effects on the litter TPL. The abundances of fungi and bacteria decreased from the gap center to the closed canopy during the two SF, SC, ST and LG periods and reversed during the two EG periods. The TPL closely followed the same trend as the microbial abundance during the first year of incubation. In addition, both species had larger TPLs in the gap center during the first winter, first year and entire two years. These findings suggest that alpine forest gap formation accelerates litter TPL and plays a dual role during specific critical periods by distributing abiotic and biotic factors directly and indirectly. In conclusion, reduced snow cover depth and duration during winter warming under current climate change scenarios or as gaps vanish may slow litter TPL in alpine biomes.

The Sentinel Application Platform (SNAP) architecture facilitates Earth Observation data processing (http://step.esa.int/main/toolboxes/snap/). In this work we present results from a new Snow Processor for SNAP. We also describe physical principles behind the developed snow property retrieval technique based on the analysis of Ocean and Land Colour Instrument (OLCI) onboard Sentinel-3A/B measurements over clean and polluted snow fields. Using OLCI spectral reflectance measurements in the range 400-1020nm, we derive important snow properties such as spectral and broadband albedo, snow specific surface area, snow extent and grain size on the spatial grid of 300m. The algorithm also incorporates cloud screening and atmospheric correction procedures over snow surfaces. We present validation results using ground measurements from Antarctica, the Greenland ice sheet and the French Alps. We find the spectral albedo retrieved with accuracy of better than 3% on average, making our retrievals sufficient for a variety of applications. Broadband albedo is retrieved with the average accuracy of about 5% over snow. Therefore, the uncertainties of satellite retrievals are close to experimental errors of ground measurements. The retrieved surface grain size shows good agreement with ground observations. Snow specific surface area observations are also consistent with our OLCI retrievals. We present snow albedo and grain size mapping over the inland ice sheet of Greenland for areas including dry snow, melted/melting snow and impurity rich bare ice. The algorithm can be applied to OLCI Sentinel-3 measurements providing an opportunity for creation of long – term snow property records essential for climate monitoring and data assimilation studies - especially in the Arctic region, where we face rapid environmental changes including reduction of snow/ice extent and, therefore, planetary albedo.

From June 18 to 19, 2013, the Ésera river in the Pyrenees, Northern Spain, caused widespread damage due to flooding as a result of torrential rains and sustained snowmelt. We estimate the contribution of snow melt to total discharge applying a snow energy balance to the catchment. Precipitation is derived from sparse local measurements and the WRF-ARW model over three nested domains, down to a grid cell size of 2 km. Temperature profiles, precipitation and precipitation gradient are well simulated, although with a possible displacement regarding the observations. Snowpack melting was correctly reproduced and verified in three instrumented sites, and according to satellite images. We found that the hydrological simulations agree well with measured discharge. Snowmelt represented 33% of total runoff during the main flood event and 23% at peak flow. The snow energy balance model indicates that most of the energy for snow melt during the day of maximum precipitation came from turbulent fluxes. This approach forecast correctly peak flow and discharge during normal conditions at least 24h in advance and could give an early warning of the extreme event 2.5 days before.

Loess was deposited in the Negev during the Pleistocene, but such sediments seem to be missing for the Holocene. This could be due to erosion unless structures such as ruins offered protection. We studied soils developed on archaeological hilltop ruins in the Negev and the Petra region and compared them with local soils, paleosols, geological outcrops, and current dust. The ruin soils in both regions were found to consist of similarly complex mixtures of local and remote sediment sources. They differ from sediments deposited during current dust storms. This seems due to fixation processes: average accretion rates are estimated to ~0.14 mm/a, suggesting that only ~3% of the current dust that can be trapped with dry marble dust collectors is stored in the soils. Vegetation, biocrusts, and/or clast pavements associated with vesicular layers seem to act as sediment-fixing agents. As well, climate might play a role: rain, and in particular one snowstorm in the Petra region brought a high amount of sediment that was more similar to the ruin soils. Wet deposition and snow might catalyze dust deposition and enhance fixation by fostering vegetation and crust formation. Frequent snow during the Pleistocene might be one explanation of enhanced loess deposition.

Monitoring snowmelt dynamics in mountains is crucial to understand water releases downstream. Sentinel-1 (S-1) synthetic aperture radar (SAR) has become one of the most widely used techniques to achieve this aim due to its high frequency of acquisitions and all-weather capability. This work aims to understand the possibilities of S-1 SAR imagery to capture snowmelt dynamics and related changes in streamflow response in semiarid mountains. The results proved that S-1 SAR imagery was able not only to capture the final spring melting but also all melting cycles that commonly appear throughout the year in these types of environments. The general change detection approach to identify wet snow was adapted for these regions using as reference the average S-1 SAR image from the previous summer, and a threshold of -3.00 dB. In addition, four different type of melting runoff onsets depending on physical snow condition were identified. When translating that at the catchment scale, distributed melting runoff onset maps were defined to better understand the spatiotemporal evolution of melting dynamics. Finally, a linear connection between melting dynamics and streamflow was found for long-lasting melting cycles, with a determination coefficient (R2) ranging from 0.62 to 0.83 and an average delay between the melting onset and streamflow peak of about 21 days.

The aims of the review were to elucidate the spatial variation in sea ice algae primary production rates and biomasses (Chl a) in the Canadian Arctic – Greenland region, characterized by comparable physical settings. A database was compiled from 30 studies of sea ice algae primary production rates, biomasses (Chl a), snow and ice thicknesses, ice types, nutrient (Si(OH)4, PO4, (NO3+NO2), and (NH4) concentrations in ice and below ice from the region. Production rates were significantly high (463 mg C m-2 d-1) in Resolute Bay and in Northern Baffin Bay (317 mg C m-2 d-1) both in Canadian Arctic, as compared to a rate of 0.2 mg C m-2 d-1 in northeast Greenland. The biomasses reached 340 mg Chl a m-2 in Resolute Bay in comparison to a 0.02 mg Chl a m-2 in southwest Greenland. Primary production at other Canadian and Greenland sites were comparable but sea ice Chl a was higher (15.0 ± 13.4 mg Chl a m-2) at Canadian sites as compared to Greenland (0.8 ± 0.5 mg Chl a m-2). Resolute and Northern Baffin Bays production rates were significantly higher when compared to other Arctic Ocean sites outside the studied region. The review concludes that high production rates and biomasses in Resolute and Northern Baffin Bay’s were related to inflow and mixing of nutrient rich waters of Pacific origin. A conceptual model with drivers and inhibitors of sea ice algae primary production is projected, and the database compiles a dataset of published data for further studies.

Implementation of the German Climate change Strategy in the Federal State of Saxony-Anhalt is discussed in this paper. It shares the requirement and importance of sustainable development. An overview of strategy, The DAS Indicator System is provided with results of a portion of work being done for the ministry of agriculture by Deutscher Wetterdienst (DWD). Applicability of the indicator system is also shown by evaluation of results for specific indicators from 1961-2015.

Forest recreation can be successfully conducted for the purpose of psychological relaxation, as has been proven in previous scientific studies. During the winter in many countries, when snow cover occurs frequently, forest recreation (walking, relaxation, photography, etc.) is common. Nevertheless, whether forest therapy conducted in a forest environment with a snow cover will also have a positive effect on psychological indicators remains unknown. Furthermore, male subjects frequently participate in forest therapy experiments, whereas females are rarely involved. Thus, in this study, the effectuality of forest recreation during winter and with snow cover was tested on 32 young females. For these reasons, the experiment involved 15-minute periods of relaxation in a forest environment or in an urban environment, in addition to a pre-test under indoor conditions. Four psychological questionnaires (POMS, PANAS, ROS, SVS) were administered to participants before and after interventions. Results showed that participants’ levels of negative mood, as measured by different aspects of the POMS questionnaire (tension-anxiety, anger-hostility, depression-dejection, confusion, fatigue), decreased after exposure to the forest environment. In contrast, both tension-anxiety and anger-hostility increased in the urban environment. The indicator of negative affect from the PANAS questionnaire also increased after exposure to the urban environment, whereas the indicator of positive affect based on PANAS was higher in the forest environment than in the urban environment. Restorativeness and subjective vitality exhibited higher values after exposure to the forest environment in comparison to those from the control and pre-test. The changes in these indicators demonstrates that forest recreation in the snow during winter can significantly increase psychological relaxation in young females, as well as showing that recreation can be successfully conducted under these winter conditions.