DOI: 10.5281/zenodo.4879293
title: Dataset - Warming-induced monsoon precipitation phase change intensifies glacier mass loss in the southeastern Tibetan Plateau

description: Materials and data results needed to reproduce the findings of the study published in (PNAS) Proceedings of the National Academy of Sciences of the United States of America: "<em>Warming-induced monsoon precipitation phase change intensifies glacier mass loss in the southeastern Tibetan Plateau"</em>. A. Jouberton, T. E. Shaw, E. Miles, M. McCarthy, S. Fugger, S. Ren, A. Dehecq, W. Yang and F. Pellicciotti It includes the meteorological forcing time-series, an exhaustive list of the model parameters, the outputs of TOPKAPI-ETH, and Matlab scripts allowing to reproduce the figures and compute the numbers given in the main manuscript as well as in the Supplementary Information. --------------- <strong>Contents</strong> : Folder : "Matlab_scripts"<br> '<strong>Climate_import.m</strong>' : Organizes meteorological forcing and generates Figure S9<br> <strong> 'TOPKAPI_result_import.m' :</strong> Imports TOPKAPI's reference run outputs and prepares them for analysis<br> <strong> 'Experiment_analysis.m' : </strong>Analyses the results of the forcing experiments, generates Figure 4 and Figure S25<br> <strong> 'Main_text_results.m' : </strong>Analyses the results of TOPKAPI's reference runs, generates Figure 1D, FIgure 2 and Figure 3<br> <strong> 'TOPKAPI_validation.m' : </strong>Compares TOPKAPI's reference run results with several validation datasets, generates the figures and performance metrics of the model calibration and validation procedure.<br> <strong> 'Parlung_albedo_regional_analysis.m'</strong>: Computes the mean glacier albedo per elevation band for each glacier within the Southeastern Tibean Plateau and compares it to the albedo of Parlung No.4 glacier.<br> <strong> 'Parlung_GMB_regional_analysis.m'</strong>: Computes the mean glacier mass balance per elevation band for each glacier within the Southeastern Tibean Plateau and compares it to the glacier mass balance of Parlung No.4 glacier.<br> <strong> 'Precipitation_phase_sensitivity_analysis.m'</strong>: Performs a sensitivity analysis on the simulated monsoon snowfall ratio per elevation band and on the attribution of glacier mass loss to precipitation<br> phase change using Monte Carlo simulations.<br> <strong> 'TOPKAPI_MODIS_validation.m'</strong>: Compares the snow cover at Parlung No.4 catchment simulated by TOPKAPI-ETH and observed by MODIS, generates Figure S19. Folder : "Remote_sensing" : Sub-Folder: 'Hugonnet' = Glacier mass balance averaged over 2000-2020 covering the Southeastern Tibetan Plateau, 100m resolution, derived from Hugonnet et al. 2021<br> Sub-Folder: 'MODIS' = contains the snow cover at Parlung No.4 derived from the daily product MOD10A1 version 61, for the period 2000-2018<br> Sub-Folder: 'Regional_glacier_albedo' = contains the annual glacier surface albedo from 2000 to 2020, covering the Southeastern Tibetan Plateau, 500m resolution.<br> Sub-Folder: 'Shapefiles' = contains the Parlung No.4 glacier outlines in 1974 and from the RGI 6.0<br> <strong>'ASTER_Nyainqentanglha_15m_utm.tif'</strong> = ASTER Digital elevation model at 15m resolution covering the Southeastern Tibetan Plateau<br> <strong> 'parlung_mask_1974.mat' </strong>= Parlung No.4 glacier mask as a matlab file<br> <strong> 'dh_ASTER_SRTM_30m.tif' </strong>= Mean elevation change rate from 2000 to 2016 at Parlung No.4 catchment.<br> <strong> 'Geodetic_map.mat'</strong> = Elevation change maps for the periods 1974-2000 and 1974-2014, as a matlab file<br> <strong> 'GMB_geodetic.mat' </strong>= Geodetic mass balance (glacier-wide mean and profile per elevation band) used in Figure S13<br> <strong> 'parlung_30m_catchment_mask.tif'</strong> = Parlung No.4 catchment mask<br> <strong>'parlung_1974_30m_dem.tif' </strong>= DEM of Parlung No.4 catchment, 30 m resolution<br> <strong> 'parlung_1974_30m_gla.tif'</strong> = Parlung No. 4 glacier mask, 30m resolution<br> <strong>'parlung_1974_30m_glah.tif' </strong>= Reconstructed ice thickness of 1975 for Parlung No.4 glacier<br> <strong>'parlung_1974_2000_diff_24m.tif' </strong>= Elevation change from DEM differencing at Parlung No.4 catchment for 1974-2000<br> <strong> 'parlung_1974_2014_diff_24m.tif' </strong>= Elevation change from DEM differencing at Parlung No.4 catchment for 1974-2014<br> <strong> 'Parlung_1974_bedrock_dem_30m.tif' </strong>= Bedrock surface digital elevation model of the catchment, 30m spatial resolution<br> <strong>'RGI_KangriKarpo_100m_utm_id.tif' </strong>= Glacier mask covering the Kangri Karpo mountain region, 100m resolution, with glacier IDs in the attribute table<br> <strong> 'RGI_Nyainqentanglha_100m_utm_id.tif' </strong>= = Glacier mask covering the Southeastern Tibetan Plateau, 100m resolution, with glacier IDs in the attribute table Folder : "TOPKAPI_forcing" :<br> <strong> CCT_AWS4600_extended.csv : </strong>Hourly cloud cover transmissivity from 1975 to 2018 reconstructed at AWSoff location <br> <strong> Climate.mat : </strong>Organizes meteorological forcings, output from the matlab script '<strong>Climate_import.m</strong>'<br> <strong> LR_AWS4600_extended.csv :</strong> Hourly temperature lapse-rates from 1975 to 2018 reconstructed at AWSoff location <br> <strong> Precipitation_AWS4600_extended.csv : </strong>Hourly precipitation from 1975 to 2018 reconstructed at AWSoff location <br> <strong> Ta_AWS4600_extended.csv :</strong> Hourly air temperature from 1975 to 2018 reconstructed at AWSoff location<br> Sub-Folder: 'National_meteorological_stations' = Contains the daily air temperature and precipitation measured at the national meteorological stations of Bomi, Zayu, Zuogong and Basu<br> Sub-Folder: 'Reference_run_inputs' = Contains the input files necessary to run TOPKAPI-ETH to obtain the outputs from which the results of this study are based on. Folder : "TOPKAPI_output":<br> <strong> </strong> Sub-Folder : "Forcing experiment" = organized TOPKAPI outputs from the forcing experiment<br> Sub-Folder :" Reference_run_outputs" = raw TOPKAPI outputs from the reference run (catchment average, spatial and grid cells)<br> Sub-Folder : "Reference_run_results" = organized TOPKAPI outputs from the reference run<br> Sub-Folder : "Snow_ice_cover" = contains TOPKAPI-ETH derived snow cover maps (daily map outputs)<br> Sub-Folder : "Regional_analysis" =<br> 'Alb'= Table containing the mean glacier albedo (2000-2020) per normalized elevation band, for each glacier in the SETP (RGI 6.0)<br> 'GMB'= Table containing the mean glacier mass balance (2000-2020) per normalized elevation band, for each glacier in the SETP (RGI 6.0)<br> 'Hypso_xxm' = Table containing the percentage of glacier area per normalized elevation band, for each glacier in the SETP (RGI 6.0), resolution of 100/500m<br> 'NormEl_100m' = Table containing the elevation per normalized elevation band, for each glacier in the SETP (RGI 6.0), resolution of 100/500m<br> Sub-Folder : "Semi_distributed_outputs" = Precipitation phase and amounts resulting from TOPKAPI-ETH simulation per elevation band, for the reference run and for the Monte Carlo sensitivity analysis Folder : "Validation_data"<br> <strong> 'topkapi.out_reference_discharge2016' </strong>= <strong> </strong>raw TOPKAPI outputs run in 2016 with AWSoff air temperature<br> <strong> </strong><strong> 'master_file_parlung.mat' </strong>=<strong> </strong>matlab structure containing AWS measurements, necessary for running <strong>'TOPKAPI_validation.m'</strong><br> <strong> 'Qdigit.mat' </strong>= Discharge measured at the Parlung No.4 glacier outlet, from Li et al., (2016)<br> <strong> 'Parlung_Q_1970.mat'</strong> = 'Discharge time-series used to run TOPKAPI-ETH (goes back to 1975, but filled with 0 when no measurements are available) In order to run the Matlab scripts, it is recommended to download all folders and gather them into the same folder. Any request about data or questions on how to run the Matlab scripts can be asked to the author of the paper (at achille.jouberton@wsl.ch).