SNOWPACK simulations for Thwaites Cavity and Channel AMIGOS sites in West Antarctica

<p>Compressed zip file, containing the SNOWPACK model simulations used in the manuscript: Maclennan, M. L., Lenaerts, J. T. M., Shields, C. A., Hoffman, A. O., Wever, N., Thompson-Munson, M., Winters, A. C., Pettit, E. C., Scambos, T. A., and Wille, J. D.: <em>Climatology and Surface Impacts of Atmospheric Rivers on West Antarctica</em>, The Cryosphere Discuss. [preprint], https://doi.org/10.5194/tc-2022-101, in review, 2022.</p> <p>MeteoIO and SNOWPACK are software published under the GNU LGPLv3 license by the WSL Institute for Snow and Avalanche Research SLF, Davos, Switzerland at <a href="https://gitlabext.wsl.ch/snow-models">https://gitlabext.wsl.ch/snow-models</a>. The repository used to develop the versions of MeteoIO and SNOWPACK used in this study can be accessed at <a href="https://github.com/snowpack-model/snowpack">https://github.com/snowpack-model/snowpack</a> (<a href="https://doi.org/10.5281/zenodo.3891845">doi: 10.5281/zenodo.3891845</a>) with the exact version corresponding to commit 149c586 (doi: <a href="https://doi.org/10.5281/zenodo.7331794">10.5281/zenodo.7331794</a>).</p> <p>The workflow consists of two phases. First, using MERRA-2 reanalysis as forcing, a spinup is performed for the MERRA-2 grid point closest to the AMIGOS weather stations Channel (CHA) and Cavity (CAV).</p> <p>To perform the spinup:</p> <ol> <li>navigate to the <code>./simulations_spinup/</code> directory</li> <li>edit the following line in the file <code>spinup.rc</code>, to modify the path where the time_shift_script can be found: <code>time_shift_script=&quot;/path/to/snowpack/Source/snowpack/tools/timeshift_sno_files.sh&quot;</code> This script can be found in the SNOWPACK repository.</li> <li>Execute the two SNOWPACK spinup simulations specified below. Change the path to point to the SNOWPACK binary executable. The simulations are performed for the two sites, which use the same meteorological forcing, but have a different end date specified for the simulation (corresponding to the installation date of the AMIGOS weather station):</li> </ol> <pre><code class="language-bash">bash spinup.sh "/path/to/usr/bin/snowpack -c cfgfiles/THWAITES_CAV.ini -e 2020-01-10T01:00 &gt; log/THWAITES_SPINUP.log 2&gt;&amp;1" startover=1 bash spinup.sh "/path/to/usr/bin/snowpack -c cfgfiles/THWAITES_CHA.ini -e 2020-01-16T01:00 &gt; log/THWAITES_SPINUP.log 2&gt;&amp;1" startover=1</code></pre> <p>Brief overview of the files in the <code>simulations_spinup</code> directory:</p> <ul> <li><code>base.ini</code>: SNOWPACK settings shared by all simulations</li> <li><code>spinup.ini</code>: Specific settings used when the spinups are being performed</li> <li><code>final.ini</code>: Specific settings used when the final simulation is being performed, after the spinup has completed</li> <li><code>spinup.rc</code>: Specific spinup settings imported by the <code>spinup.sh</code> script</li> <li><code>spinup.sh</code>: bash script to perform spinups</li> <li><code>cfgfiles/THWAITES_CAV.ini</code>: Specific model settings for the Cavity site. This <code>*.ini</code> file is specified to SNOWPACK, and it directs to include other <code>*ini</code> files.</li> <li><code>cfgfiles/THWAITES_CHA.ini</code>: Specific model settings for the Channel site. This *.ini file is specified to SNOWPACK, and it directs to include other <code>*ini</code> files.</li> <li><code>current_snow/THWAITES.sno</code>: <code>*.sno</code> file used for the restarts within the spinup procedure</li> <li><code>log</code>: directory to store simulation logs. Note that in the commands above, we redirect <code>stdout</code> and <code>stderr</code> to files in this directory</li> <li><code>output</code>: directory to store model outputs. See SNOWPACK documentation for more details. <code>*.smet</code> and <code>*.pro</code> files can be visualized using <a href="https://niviz.org">NiViz</a></li> <li><code>smet</code>: directory to store the meteorological forcing data for the simulations</li> <li><code>snow_init/THWAITES.sno</code>: initial <code>*.sno</code> files used to start the spinup</li> </ul> <p><br> To perform the final simulations, covering the atmospheric river period:</p> <ol> <li>navigate to the <code>./simulations_final/</code> directory</li> <li>execute a snowpack simulation for each <code>*.ini</code> file in the <code>cfgfiles</code> directory:</li> </ol> <pre><code class="language-bash">for f in cfgfiles/*ini do /path/to/usr/bin/snowpack -c ${f} -e NOW done</code></pre> <p>The following simulations will be performed:</p> <table> <tbody> <tr> <td><strong>cfgfiles/*ini file:</strong></td> <td><strong>Label output files</strong></td> <td><strong>Mode</strong></td> <td><strong>AMIGOS AWS</strong></td> <td><strong>Parameters from AMIGOS AWS</strong></td> <td><strong>Reanalysis</strong></td> <td><strong>Parameters from reanalysis</strong></td> </tr> <tr> <td>THWAITES_CAV_ERA5.ini</td> <td>THWAITES_ERA5_HSDRIVEN_CAV_ERA5</td> <td>Snow height driven</td> <td>Cavity</td> <td>TA RH VW HS</td> <td>ERA5</td> <td>ILWR ISWR</td> </tr> <tr> <td>THWAITES_CAVHS_ERA5.ini</td> <td>THWAITES_ERA5_HSDRIVEN_CAVHS_ERA5</td> <td>Snow height driven</td> <td>Cavity</td> <td>HS</td> <td>ERA5</td> <td>TA RH VW ILWR ISWR</td> </tr> <tr> <td>THWAITES_CHA_ERA5.ini</td> <td>THWAITES_ERA5_HSDRIVEN_CHA_ERA5</td> <td>Snow height driven</td> <td>Channel</td> <td>TA RH VW HS</td> <td>ERA5</td> <td>ILWR ISWR</td> </tr> <tr> <td>THWAITES_CHAHS_ERA5.ini</td> <td>THWAITES_ERA5_HSDRIVEN_CHAHS_ERA5</td> <td>Snow height driven</td> <td>Channel</td> <td>HS</td> <td>ERA5</td> <td>TA RH VW ILWR ISWR</td> </tr> <tr> <td>THWAITES_REFCAV_ERA5.ini</td> <td>THWAITES_ERA5_PSUM_CAV_ERA5</td> <td>Precipitation driven&nbsp;</td> <td>Cavity</td> <td>&nbsp;</td> <td>ERA5</td> <td>PSUM TA RH VW ILWR ISWR</td> </tr> <tr> <td>THWAITES_REFCHA_ERA5.ini</td> <td>THWAITES_ERA5_PSUM_CHA_ERA5</td> <td>Precipitation driven&nbsp;</td> <td>Channel</td> <td>&nbsp;</td> <td>ERA5</td> <td>PSUM TA RH VW ILWR ISWR</td> </tr> <tr> <td>THWAITES_CAV_MERRA2.ini</td> <td>THWAITES_MERRA2_HSDRIVEN_CAV_MERRA2</td> <td>Snow height driven</td> <td>Cavity</td> <td>TA RH VW HS</td> <td>MERRA-2</td> <td>ILWR ISWR</td> </tr> <tr> <td>THWAITES_CAVHS_MERRA2.ini</td> <td>THWAITES_MERRA2_HSDRIVEN_CAVHS_MERRA2</td> <td>Snow height driven</td> <td>Cavity</td> <td>HS</td> <td>MERRA-2</td> <td>TA RH VW ILWR ISWR</td> </tr> <tr> <td>THWAITES_CHA_MERRA2.ini</td> <td>THWAITES_MERRA2_HSDRIVEN_CHA_MERRA2</td> <td>Snow height driven</td> <td>Channel</td> <td>TA RH VW HS</td> <td>MERRA-2</td> <td>ILWR ISWR</td> </tr> <tr> <td>THWAITES_CHAHS_MERRA2.ini</td> <td>THWAITES_MERRA2_HSDRIVEN_CHAHS_MERRA2</td> <td>Snow height driven</td> <td>Channel</td> <td>HS</td> <td>MERRA-2</td> <td>TA RH VW ILWR ISWR</td> </tr> <tr> <td>THWAITES_REFCHA_MERRA2.ini</td> <td>THWAITES_MERRA2_PSUM_CAV_MERRA2</td> <td>Precipitation driven&nbsp;</td> <td>Cavity</td> <td>&nbsp;</td> <td>MERRA-2</td> <td>PSUM TA RH VW ILWR ISWR</td> </tr> <tr> <td>THWAITES_REFCAV_MERRA2.ini</td> <td>THWAITES_MERRA2_PSUM_CHA_MERRA2</td> <td>Precipitation driven&nbsp;</td> <td>Channel</td> <td>&nbsp;</td> <td>MERRA-2</td> <td>PSUM TA RH VW ILWR ISWR</td> </tr> </tbody> </table> <p>&nbsp;</p> <p>Brief overview of the files in the <code>simulations_final</code> directory:</p> <ul> <li><code>base.ini</code>: SNOWPACK settings shared by all simulations (identical to <code>base.ini</code> in simulations_spinup directory</li> <li><code>cfgfiles</code>: Specific model settings for each simulation as listed in the table above.</li> <li><code>current_snow</code>: The initial *.sno file containing the initial firn conditions</li> <li><code>final.ini</code>: Settings used by all simulations in this directory</li> <li><code>output</code>: directory to store model outputs. See SNOWPACK documentation for more details. <code>*.smet</code> and <code>*.pro</code> files can be visualized using <a href="https://niviz.org">NiViz</a></li> <li><code>prepare_snowpack.sh</code>: bash script setting up the final simulations after the simulations in simulations_spinup completed. Among others, the code to mark the installation reference layers for the snow depth measurements is included in this script.</li> <li><code>smet</code>: ERA5 and MERRA-2 forcing data, as well as snow height for Cavity and Channel (<code>cavity_hs.smet</code> and <code>channel_hs.smet</code>, respectively) as well as the other meteorological parameters (<code>cavity_meteo.smet</code> and <code>channel_meteo.smet</code>, respectively).</li> </ul> <p>&nbsp;</p> <p><br> Note:<br> Please refer to <a href="https://doi.org/10.15784/601552">doi: 10.15784/601552</a> and <a href="https://doi.org/10.15784/601549">doi: 10.15784/601549</a> for the AMIGOS weather station data</p>