FVCOM_Petermann_Discharge_Input_Files

<h2>Overview</h2> <p>This repository contains all the input files that are required to run the subglacial discharge experiments for the Petermann Ice Shelf-fjord system using the Finite Volume Community Ocean Model (FVCOM). You can access the corresponding peer-reviewed, open-access research article here: <a href="https://doi.org/10.1038/s41467-025-59469-9" target="_blank" rel="noopener">Enhanced subglacial discharge amplifies Petermann Ice Shelf melting when ocean thermal forcing saturates</a></p> <h2>Pre-requisites</h2> <p>The open source code Finite Volume Community Ocean Model version 4.0 (FVCOM v4.0) augmented by an ice shelf (Zhou and Hattermann, 2020) and sea ice module (Prakash et al., 2022) is required to conduct the numerical experiments, and is made publicly available at <a href="https://doi.org/10.5281/zenodo.15084570" target="_blank" rel="noopener">https://doi.org/10.5281/zenodo.15084570</a>. A detailed description of our high-resolution (200 m) unstructured grid 3-D regional model setup centered on the Petermann ice shelf and fjord is available from Prakash et al. (2022) and Prakash et al. (2023).</p> <h2>File description</h2> <p>The model input files can be classified under 2 broad categories, namely&nbsp;</p> <ul> <li>static, and&nbsp;</li> <li>time-varying.</li> </ul> <p>Static : These include the following files</p> <ol> <li><em>zisf.dat</em> (realistic Petermann ice shelf draft derived from BedMachine v3 (Morlighem et al., 2017)),&nbsp;</li> <li><em>dep.dat</em> (smoothed bathymetry, which has been modified to provide an improved (with respect to BedMachine v3 (Morlighem et al., 2017)) representation of the sub-ice shelf topography),&nbsp;</li> <li><em>cor.dat</em> (Coriolis),</li> <li><em>sigma.dat </em>(terrain-following vertical coordinate system),</li> <li><em>grd.dat</em> (FVCOM Petermann Fjord regional model grid)</li> </ol> <p>Time-varying: These include the following files</p> <p>(a) Atmospheric Forcing</p> <ol> <li><em>ocnice_YYYY_wnd.nc</em> (wind forcing),&nbsp;</li> <li><em>ocnice_YYYY_fwf.nc</em> (freshwater flux), and</li> <li><em>ocnice_YYYY_hf.nc</em> (heat flux),&nbsp;</li> </ol> <p>where YYYY represents the forcing year, e.g. 2014.</p> <p>(b) Sea ice forcing</p> <ol> <li><em>ocnice_YYYY_icenudge.nc</em> (includes sea ice concentration and thickness, bulk ice salinity, and sea ice velocities),</li> </ol> <p>where YYYY represents the forcing year, e.g. 2014.</p> <p>(c) Ocean boundary conditions</p> <ol> <li><em>pf_l_bcorr_climaYY.nc</em> (bias-corrected lateral ocean boundary conditions which includes temperature, salinity, ocean velocities, and sea surface elevation)</li> </ol> <p>where YY represents the forcing year, e.g. 14 (for 2014).</p> <p>In addition to these two broad categories, a restart file (<em>*_restart_*.nc</em>) is included in each run directory, which is used to provide "realistic" (a.k.a. "hotstart") initial ocean conditions. The restart file in the directory&nbsp;<em>control_run_2014</em> includes a "realistic" initial (July 01, 2014 00:00:00 UTC) ocean condition derived from a 4-km pan-Arctic (A4) ROMS model (Hattermann et al., 2016). Thereafter, the FVCOM model solution on January 01, 2015 00:00:00 UTC is used to "hotstart" the 2015 calendar year. Lastly, the stable FVCOM model solution from January 01, 2016 00:00:00 UTC is used to "hotstart" the 2016 calendar year, as well as the discharge experiments. Note that the discharge experiments depart from the control experiment (<em>control_run_2016</em>) by introducing subglacial discharge of varying magnitudes. All other conditions are kept the same, in order to isolate and investigate the impact of subglacial discharge on the basal melting of the Petermann ice shelf. The namelist file (*.nml) for each run is also included in the corresponding directory.</p> <p>To minimize storage space, each run directory is populated with only the input files that are unique to that run.&nbsp;For e.g., <em>control_run_2015/2016</em> and the <em>discharge_experiments</em> do not include the static files from <em>control_run_2014</em>, however, they are required to conduct each of these runs. Likewise, <em>discharge_experiments</em> do not include the time-varying forcing files, since they are identical to the ones contained in <em>control_run_2016</em>. The <em>discharge_experiments</em> do however include the discharge forcings for the present, RCP 8.5, and median(present,RCP 8.5) scenarios. These are the riverdata (riverdata.nc (present), riverdata_expC.nc (median), and riverdata_expD.nc (RCP 8.5)) and the corresponding rivernamelist (riverNamelist.nml, riverNamelist_expC.nml, riverNamelist_expD.nml) files. Note that a rivernamelist is different from a namelist file, in that it only contains information about the subglacial discharge.</p> <p>Development of the smoothed model topography datasets for our unstructured grid (horizontal resolutions of 200 m, 2km, and 4km) regional model domain, including localized improvements made within the Petermann Fjord; the ocean nesting and surface forcing methodologies (models used, downscaling procedure etc.); bias-correction, among other technical details, are detailed in Prakash et al. (2022) and Prakash et al. (2023).</p>