Seasonal forecasting of snow resources at Alpine sites

Climate warming in mountain regions is resulting in glacier shrinking, seasonal snow cover reduction, and changes in the amount and seasonality of meltwater runoff, with consequences on water availability. Droughts are expected to become more severe in the future with economical and environmental losses both locally and downstream. Effective adaptation strategies involve multiple timescales, and seasonal forecasts can help in the optimization of the available snow and water resources with a lead time of several months. We developed a prototype to generate seasonal forecasts of snow depth and snow water equivalent with a starting date of 1 November and a lead time of 7 months, so up to 31 May of the following year. The prototype has been co-designed with end users in the field of water management, hydropower production and mountain ski tourism, meeting their needs in terms of indicators, time resolution of the forecasts and visualization of the forecast outputs. In this paper we present the modelling chain, based on the seasonal forecasts of the ECMWF and Meteo-France seasonal prediction systems, made available through the Copernicus Climate Change Service (C3S) Climate Data Store. Seasonal forecasts of precipitation, near-surface air temperature, radiative fluxes, wind and relative humidity are bias-corrected and downscaled to three sites in the Western Italian Alps and finally used as input for the physically based multi-layer snow model SNOWPACK. Precipitation is bias-corrected with a quantile mapping method using ERA5 reanalysis as a reference and then downscaled with the RainFARM stochastic procedure in order to allow an estimate of uncertainties due to the downscaling method. The impacts of precipitation bias adjustment and downscaling on the forecast skill have been investigated. The skill of the prototype in predicting the deviation of monthly snow depth with respect to the normal conditions from November to May in each season of the hindcast period 1995-2015 is demonstrated using both deterministic and probabilistic metrics. Forecast skills are determined with respect to a simple forecasting method based on the climatology, and station measurements are used as reference data. The prototype shows good skills at predicting the tercile category, i.e. snow depth below and above normal, in the winter (lead times: 2-3-4 months) and spring (lead times: 5-6-7 months) ahead: snow depth is predicted with higher accuracy (Brier skill score) and higher discrimination (area under the relative operating characteristics (ROC) curve skill score) with respect to a simple forecasting method based on the climatology. Ensemble mean monthly snow depth forecasts are significantly correlated with observations not only at short lead times of 1 and 2 months (November and December) but also at lead times of 5 and 6 months (March and April) when employing the ECMWFS5 forcing. Moreover the prototype shows skill at predicting extremely dry seasons, i.e. seasons with snow depth below the 10th percentile, while skills at predicting snow depth above the 90th percentile are model-, station- and score-dependent.The bias correction of precipitation forecasts is essential in the case of large biases in the global seasonal forecast system (MFS6) to reconstruct a realistic snow depth climatology; however, no remarkable differences are found among the skill scores when the precipitation input is bias-corrected, downscaled, or bias-corrected and downscaled, compared to the case in which raw data are employed, suggesting that skill scores are weakly sensitive to the treatment of the precipitation input

Seasonal forecasting of snow resources at Alpine sites

Climate warming in mountain regions is resulting in glacier shrinking, seasonal snow cover reduction, and changes in the amount and seasonality of meltwater runoff, with consequences on water availability. Droughts are expected to become more severe in the future with economical and environmental losses both locally and downstream. Effective adaptation strategies involve multiple timescales, and seasonal forecasts can help in the optimization of the available snow and water resources with a lead time of several months. We developed a prototype to generate seasonal forecasts of snow depth and snow water equivalent with a starting date of 1 November and a lead time of 7 months, so up to 31 May of the following year. The prototype has been co-designed with end users in the field of water management, hydropower production and mountain ski tourism, meeting their needs in terms of indicators, time resolution of the forecasts and visualization of the forecast outputs. In this paper we present the modelling chain, based on the seasonal forecasts of the ECMWF and Météo-France seasonal prediction systems, made available through the Copernicus Climate Change Service (C3S) Climate Data Store. Seasonal forecasts of precipitation, near-surface air temperature, radiative fluxes, wind and relative humidity are bias-corrected and downscaled to three sites in the Western Italian Alps and finally used as input for the physically based multi-layer snow model SNOWPACK. Precipitation is bias-corrected with a quantile mapping method using ERA5 reanalysis as a reference and then downscaled with the RainFARM stochastic procedure in order to allow an estimate of uncertainties due to the downscaling method. The impacts of precipitation bias adjustment and downscaling on the forecast skill have been investigated. The skill of the prototype in predicting the deviation of monthly snow depth with respect to the normal conditions from November to May in each season of the hindcast period 1995–2015 is demonstrated using both deterministic and probabilistic metrics. Forecast skills are determined with respect to a simple forecasting method based on the climatology, and station measurements are used as reference data. The prototype shows good skills at predicting the tercile category, i.e. snow depth below and above normal, in the winter (lead times: 2–3–4 months) and spring (lead times: 5–6–7 months) ahead: snow depth is predicted with higher accuracy (Brier skill score) and higher discrimination (area under the relative operating characteristics (ROC) curve skill score) with respect to a simple forecasting method based on the climatology. Ensemble mean monthly snow depth forecasts are significantly correlated with observations not only at short lead times of 1 and 2 months (November and December) but also at lead times of 5 and 6 months (March and April) when employing the ECMWFS5 forcing. Moreover the prototype shows skill at predicting extremely dry seasons, i.e. seasons with snow depth below the 10th percentile, while skills at predicting snow depth above the 90th percentile are model-, station- and score-dependent. The bias correction of precipitation forecasts is essential in the case of large biases in the global seasonal forecast system (MFS6) to reconstruct a realistic snow depth climatology; however, no remarkable differences are found among the skill scores when the precipitation input is bias-corrected, downscaled, or bias-corrected and downscaled, compared to the case in which raw data are employed, suggesting that skill scores are weakly sensitive to the treatment of the precipitation input.</p

Detecting barriers to transport: A review of different techniques

We review and discuss some different techniques for describing local dispersion properties in fluids. A recent Lagrangian diagnostics, based on the Finite Scale Lyapunov Exponent (FSLE), is presented and compared to the Finite Time Lyapunov Exponent (FTLE), and to the Okubo-Weiss (OW) and Hua-Klein (HK) criteria. We show that the OW and HK are a limiting case of the FTLE, and that the FSLE is the most efficient method for detecting the presence of cross-stream barriers. We illustrate our findings by considering two examples of geophysical interest: a kinematic meandering jet model, and Lagrangian tracers advected by stratospheric circulation.Comment: 15 pages, 9 figures, submitted to Physica

Linking Vegetation-Climate-Fire Relationships in Sub-Saharan Africa to Key Ecological Processes in Two Dynamic Global Vegetation Models

Africa is largely influenced by fires, which play an important ecological role influencing the distribution and structure of grassland, savanna and forest biomes. Here vegetation strongly interacts with climate and other environmental factors, such as herbivory and humans. Fire-enabled Dynamic Global Vegetation Models (DGVMs) display high uncertainty in predicting the distribution of current tropical biomes and the associated transitions, mainly due to the way they represent the main ecological processes and feedbacks related to water and fire. The aim of this study is to evaluate the outcomes of two state-of-the–art DGVMs, LPJ-GUESS and JSBACH, also currently used in two Earth System Models (ESMs), in order to assess which key ecological processes need to be included or improved to represent realistic interactions between vegetation cover, precipitation and fires in sub-Saharan Africa. To this end, we compare models and remote-sensing data, analyzing the relationships between tree and grass cover, mean annual rainfall, average rainfall seasonality and average fire intervals, using generalized linear models, and we compare the patterns of grasslands, savannas, and forests in sub-Saharan Africa. Our analysis suggests that LPJ-GUESS (with a simple fire-model and complex vegetation description) performs well in regions of low precipitation, while in humid and mesic areas the representation of the fire process should probably be improved to obtain more open savannas. JSBACH (with a complex fire-model and a simple vegetation description) can simulate a vegetation-fire feedback that can maintain open savannas at intermediate and high precipitation, although this feedback seems to have stronger effects than observed, while at low precipitation JSBACH needs improvements in the representation of tree-grass competition and drought effects. This comparative process-based analysis permits to highlight the main factors that determine the tropical vegetation distribution in models and observations in sub-Saharan Africa, suggesting possible improvements in DGVMs and, consequently, in ESM simulations for future projections. Given the need to use carbon storage in vegetation as a climate mitigation measure, these models represent a valuable tool to improve our understanding of the sustainability of vegetation carbon pools as a carbon sink and the vulnerability to disturbances such as fire

Precision Measurements of Stretching and Compression in Fluid Mixing

The mixing of an impurity into a flowing fluid is an important process in many areas of science, including geophysical processes, chemical reactors, and microfluidic devices. In some cases, for example periodic flows, the concepts of nonlinear dynamics provide a deep theoretical basis for understanding mixing. Unfortunately, the building blocks of this theory, i.e. the fixed points and invariant manifolds of the associated Poincare map, have remained inaccessible to direct experimental study, thus limiting the insight that could be obtained. Using precision measurements of tracer particle trajectories in a two-dimensional fluid flow producing chaotic mixing, we directly measure the time-dependent stretching and compression fields. These quantities, previously available only numerically, attain local maxima along lines coinciding with the stable and unstable manifolds, thus revealing the dynamical structures that control mixing. Contours or level sets of a passive impurity field are found to be aligned parallel to the lines of large compression (unstable manifolds) at each instant. This connection appears to persist as the onset of turbulence is approached.Comment: 5 pages, 5 figure

Aquifer recharge in the Piedmont Alpine zone: Historical trends and future scenarios

The spatial and temporal variability of air temperature, precipitation, actual evapotranspiration (AET) and their related water balance components, as well as their responses to anthropogenic climate change, provide fundamental information for an effective management of water resources and for a proactive involvement of users and stakeholders, in order to develop and apply adaptation and mitigation strategies at the local level. In this study, using an interdisciplinary research approach tailored to water management needs, we evaluate the past, present and future quantity of water potentially available for drinking supply in the water catchments feeding the about 2.3 million inhabitants of the Turin metropolitan area (the former Province of Turin, north-western Italy), considering climatologies at the quarterly and yearly timescales. Observed daily maximum surface air temperature and precipitation data from 1959 to 2017 were analysed to assess historical trends, their significance and the possible cross-correlations between the water balance components. Regional climate model (RCM) simulations from a small ensemble were analysed to provide mid-century projections of the difference between precipitation and AET for the area of interest in the future CMIP5 scenarios RCP4.5 (stabilization) and RCP8.5 (business as usual). Temporal and spatial variations in recharge were approximated with variations of drainage. The impact of irrigation, and of snowpack variability, on the latter was also assessed. The other terms of water balance were disregarded because they are affected by higher uncertainty. The analysis over the historical period indicated that the driest area of the study region displayed significant negative annual (and spring) trends of both precipitation and drainage. Results from field experiments were used to model irrigation, and we found that relatively wetter watersheds in the northern and in the southern parts behave differently, with a significant increase of AET due to irrigation. The analysis of future projections suggested almost stationary conditions for annual data. Regarding quarterly data, a slight decrease in summer drainage was found in three out of five models in both emission scenarios. The RCM ensemble exhibits a large spread in the representation of the future drainage trends. The large interannual variability of precipitation was also quantified and identified as a relevant risk factor for water management, expected to play a major role also in future decades

Sensitivity of snow models to the accuracy of meteorological forcings in mountain environments

Snow models are usually evaluated at sites providing high-quality meteorological data, so that the uncertainty in the meteorological input data can be neglected when assessing model performances. However, high-quality input data are rarely available in mountain areas and, in practical applications, the meteorological forcing used to drive snow models is typically derived from spatial interpolation of the available in situ data or from reanalyses, whose accuracy can be considerably lower. In order to fully characterize the performances of a snow model, the model sensitivity to errors in the input data should be quantified. In this study we test the ability of six snow models to reproduce snow water equivalent, snow density and snow depth when they are forced by meteorological input data with gradually lower accuracy. The SNOWPACK, GEOTOP, HTESSEL, UTOPIA, SMASH and S3M snow models are forced, first, with high-quality measurements performed at the experimental site of Torgnon, located at 2160ma.s.l. in the Italian Alps (control run). Then, the models are forced by data at gradually lower temporal and/or spatial resolution, obtained by (i) sampling the original Torgnon 30 min time series at 3, 6, and 12 h, (ii) spatially interpolating neighbouring in situ station measurements and (iii) extracting information from GLDAS, ERA5 and ERA-Interim reanalyses at the grid point closest to the Torgnon site. Since the selected models are characterized by different degrees of complexity, from highly sophisticated multi-layer snow models to simple, empirical, single-layer snow schemes, we also discuss the results of these experiments in relation to the model complexity. The results show that, when forced by accurate 30 min resolution weather station data, the single-layer, intermediatecomplexity snow models HTESSEL and UTOPIA provide similar skills to the more sophisticated multi-layer model SNOWPACK, and these three models show better agreement with observations and more robust performances over different seasons compared to the lower-complexity models SMASH and S3M. All models forced by 3-hourly data provide similar skills to the control run, while the use of 6- A nd 12-hourly temporal resolution forcings may lead to a reduction in model performances if the incoming shortwave radiation is not properly represented. The SMASH model generally shows low sensitivity to the temporal degradation of the input data. Spatially interpolated data from neighbouring stations and reanalyses are found to be adequate forcings, provided that temperature and precipitation variables are not affected by large biases over the considered period. However, a simple bias-adjustment technique applied to ERA-Interim temperatures allowed all models to achieve similar performances to the control run. Regardless of their complexity, all models show weaknesses in the representation of the snow density

Thermal niche predicts recent changes in range size for bird species

Species’ distributions are strongly affected by climate, and climate change is affecting species and populations. Thermal niches are widely used as proxies for estimating thermal sensitivity of species, and have been frequently related to community composition, population trends and latitudinal/elevational shifts in distribution. To our knowledge, no work has yet explored the relationship between thermal niche and change in range size (changes in the number of occupied spatial units over time) in birds. In this study, we related a 30 yr change in range size to species thermal index (STI: average temperature at occurrence sites) and to other factors (i.e. birds’ associated habitats, body mass, hunting status) potentially affecting bird populations/range size. We analysed trends of breeding bird range in Italy for a suite of poorly studied cold-adapted animals potentially sensitive to global warming, and for a related group of control species taxonomically similar and with comparable mass but mainly occurring at lower/warmer sites. We found a strong positive correlation between change in range size and STI, confirming that recent climatic warming has favoured species of warmer climates and adversely affected species occupying colder areas. A model including STI and birds’ associated habitats was not so strongly supported, with forest species performing better than alpine open habitat and agricultural ones. In line with previous works highlighting effects of recent climate change on community composition, species’ population trends and poleward/upward distributional shifts, we found STI to be the most important predictor of change in range size variation in breeding birds

Asian Monsoon and Elevated Heat Pump Mechanism in the CMCC coupled aerosol-climate model simulations

A coupled aerosol-atmosphere-ocean-seaice model is used to analyse the relationship between aerosol and the Asian summer monsoon. In this analysis the elevated heat pump hypothesis and the solar dimming effect associ- ated with aerosol loading are verified and are found to be consistent with our simulations. When increased aerosol loading is found on the Himalayas slopes in the pre-monsoon period (April-May), an intensification in early mon- soon rainfall over India is obtained. An increase in rainfall and cloudiness during the early monsoon has a cooling effect on the land surface. Here it is verified that this cooling is produced also through the solar dimming effect by the presence of more dust from the deserts brought by an increased westerly flow in early monsoon season. A subsequent reduction in monsoon rainfall over India is found, with a beginning of this decrease in northern India. As these results, obtained with a fully coupled model, reproduce a reasonably realistic pattern, it is possible to consider absorbing aerosols as a possible source of seasonal predictability of the Asian summer monsoon over the Indian subcontinent