DOI: 10.5281/zenodo.3958846
title: Coupled Hydrological and Thermal Models of Rockwall Permafrost

description: This dataset contains forcing data and selected output of coupled thermal and hydrological simulations applied to a high-elevated rockwall site (the Aiguille du Midi, 3842 m asl, Mont Blanc massif, France). All data are provided as .shp and .shx for display, as well as a .dbf file for quick reading. They are made of 5 columns, whose: « Node » is the Node ID, « X » is the position (in m) on the x axis, « Y » is the position (in m) on the y axis, « xINIT » is the calculated value for the parameter indicated in the file name(head, temperature, etc.) « Time » is the time step at which the value is calculated. The model output are gathered according to various cases studies of saturation and water flows. The model settings of the various cases studies are outlined in this file but more details about the mathematical approach and modeling settings and strategy are provided in the study to which the dataset belongs and which was submited for the first time to Journal of Geophysical Research: Earth Surface in July 2020. <strong><em>SaFl </em></strong>corresponds to a saturated with forced water flows case study by assuming a constant recharge and discharge. <strong><em>SaNF</em></strong> is a saturated case study with no water flows. <strong><em>uSFl</em></strong> corresponds to an unsaturated case study with forced water flows in selected fractures only. <strong><em>uSLF</em></strong> is unsaturated with limited water flows. For <strong><em>SaFl</em></strong> the following output are provided: Darcy flux (m.s<sup>-1</sup>) at various time step after of transient simulations (1550 AD, 2000 AD, 2015 AD, 2030 AD) such as displayed in Figures S2 and S3. Hydraulic heads (m) after initialization (0 AD) and at various time steps of the transient simulations (1850 AD, 2000 AD, 2015 AD, 2030 AD and 2100 AD) such as in Figure 6 and S3. The ice bulk volumetric fraction at various time step after of transient simulations (2000 AD, 2015 AD, 2030 AD) such as in Figure S3. Temperature (°C) after initialization (0 AD) and at various time steps of the transient simulations (1850 AD, 2000 AD, 2015 AD, 2030 AD and 2100 AD) such as in Figure 4 and S3. For <strong><em>SaNF </em></strong>the following output are provided: Temperature (°C) after initialization (0 AD) and for 1850 AD and 2100 AD such as in Figure 4. Hydraulic heads (m) after initialization (0 AD) and for 1850 AD and 2100 AD such as in Figure 6. For <strong><em>uSFl </em></strong>the following output are provided: Hydraulic heads (m) and temperature (°C) at 2100 AD such as in Figure 4 and 6. Saturation in 852 AD, 853 AD and 854 AD such as in Figure 5. For <strong><em>uSLF</em></strong> the following output are provided: Darcy flux (m.s<sup>-1</sup>) at 1550 AD and 2015 AD such as in Figure S2 Hydraulic heads (m) after initialization (0 AD) and at 1850 AD and 2100 AD) such as in Figure 6. Temperature (°C) after initialization (0 AD) and for 1850 AD and 2100 AD such as in Figure 4. In addition, the temperature and hydrological (hydraulic head changes for the unsaturated cases studies only) forcing data (« Input_data_boundary_conditions.csv ») are provided. This last file contains: <em>A to D</em>: the surface points extracted along the 4-m resolution DEM transect (Horizontal (X) position in m, Vertical (Y) position (elevation in m)) and MARST map (for the period 1961-1990: MARST<sub>init</sub>) as illustrated on Fig. S1, together with the adjusted MARST to run the model initialization simulations. <em>F to I</em>: the surface points taken on for forcingthe model at its upper boundary and in between which the forcing data were interpolated. <em>L to BH</em>: Data used for transient simulations with the time in year (<em>L</em>), the MARST anomaly applied to the adjusted MARST (<em>M</em>), the time in days (<em>N</em>), the MARST value applied at each surface node from 1 to 23 (<em>O </em>to <em>AK</em>), as well as the changes in head values applied at at each surface node from 1 to 23 (<em>AL </em>to <em>BH</em>) for the concerned simulations. More information can be made available by contacting Florence Magnin or Jean-Yves Josnin at <em>florence.magnin@univ-smb.fr</em><em> </em>or <em>jean-yves.josnin@univ-smb.fr</em>