Title: Dataset: A numerical simulation of the ocean, sea ice and ice shelves in the Amundsen Sea (Antarctica) over the period 2006-2022 and its associated code and input files
A three-dimensional numerical model of the Amundsen Sea (Antarctica) was used to simulate the period Jan.2006-Mar.2022 under consistent atmospheric/oceanic forcings, bathymetry/ice shelf topography, and model equations/parameters. The model is an implementation of the Regional Ocean Modeling System (ROMS, https://www.myroms.org/) with extensions for sea ice (Budgell 2005) and ice shelves (Dinniman et al. 2011). It simulates the ocean hydrography and circulation, sea ice thermodynamics and dynamics, and the basal melt of the ice shelves, with a uniform horizontal mesh of 1.5km and 20 topography-following vertical levels. Forcings include the ERA5 reanalysis (3-hourly), 10 tidal constituents from CATS 2008, and ocean/sea ice conditions at the edges of the model domain taken from the 5km-resolution circumpolar model of Dinniman et al. 2020 and from daily SSM/I satellite images. The model outputs are divided into nine directories each containing two years worth of model results (run661-669) in the NetCDF format. Each directory contains: daily-averaged model fields (roms_avg_xxxx.nc), instantaneous snapshots every 3 hours for select fields (roms_qck_xxxx.nc), and instantaneous snapshots every 30 days (roms_his_xxxx.nc). All the metadata information necessary for the interpretation of the model outputs (dimensions, units, etc) is included inside the NetCDF files. The NetCDF files follow the CF conventions and can be opened with various software that are open source and freely available over the Internet. In addition to the model outputs, this archive includes the computer code as well as the input files necessary for reproducing the model outputs of this archive.  more » « less
Award ID(s):
1941292
PAR ID:
10532664
Author(s) / Creator(s):
Publisher / Repository:
William & Mary. Virginia Institute of Marine Science
Date Published:
Edition / Version:
1.0
Subject(s) / Keyword(s):
Amundsen Sea Antarctica model ocean sea ice ice shelves continental shelf polynyas ROMS polar Southern Ocean
Format(s):
Medium: X Other: NetCDF
Location:
Amundsen Sea, Antarctica
Institution:
Virginia Institute of Marine Science (VIMS), William & Mary
Sponsoring Org:
National Science Foundation
More Like this
  1. Satellite images from Antarctica reveal important changes in the coastal icescape (fast-ice, icebergs and ice shelves) but these yearly changes and their impacts on the coastal circulation and ice shelf basal melt rates are not represented in the Earth System Models used to project future sea level rise. The impacts of these yearly icescape changes are thus investigated using a high-resolution regional ocean-ice shelves-sea ice coupled model of the Amundsen Sea (Antarctica). A set of nine semi-idealized experiments were designed to highlight the impacts of (a) the collapse of the Thwaites Glacier Tongue, (b) the disappearance of the Bear Ridge Iceberg Chain and tabular iceberg B22, and (c) presence/absence of a fast-ice cover between Thwaites and Pine Island ice shelves, in both cold and warm background hydrological conditions. The dataset features the results of the nine experiments and reveals changes in sea ice concentrations, coastal oceanic circulation and oceanic heat supply to the ice shelf cavities, ice shelf basal melt rates, hydrological conditions, and fluxes of heat/freshwater at the sea surface. These model results are archived in self-documented NetCDF files with the appropriate metadata for each variable. The dataset includes a 'readme file' providing an overview of the archive as well as additional information regarding the model results. 
    more » « less
  2. Earliest snowmelt estimation dates calculated for the year 2003 are provided using sea ice brightness temperatures from AMSR-E (Cavalieri et al., 2014) and DMSP SSM/I-SSMIS (Meier et al., 2019), as well as simulated sea ice brightness temperatures from the CESM2 JRA-55 (Danabasoglu et al., 2020; Kobayashi et al., 2015; Tsujino et al., 2018), which were created using the Arctic Ocean Observation Operator (ARC3O; Burgard et al, 2020a,b). Scripts and README files are provided for preparing the model data to act as input to ARC3O. 
    more » « less
  3. {"Abstract":["The original "CESM2 Tuned Sea Ice Albedo" experiments were run to investigate the impact of sea ice snow albedo tuning on Arctic mean state and global responses. These simulations were run by David Bailey, and the results were written up in Kay et al. 2022 [https://doi.org/10.1029/2021MS002679]. These were configured the same as the standard CMIP6 runs, but with adjusted r_snw=1.5 and dt_mlt=1.0 in the sea ice model (cice) namelist. (For local disk access, these data are all in the main directory and are not separated by case name.)Original experimentsControl run: b.e21.B1850.f09_g17.CMIP6-piControl.001_branch2 * Data are generally available years 811-1350 * These runs are fully coupled * Monthly Ocean, Sea Ice, and Atmosphere variables are availableHistorical runs: b.e21.BHIST.f09_g17.CMIP6-historical.01? * There are four ensemble members * Ensemble number 012 started in 1850 and ran to 2015 * Ensemble numbers 013, 014, 015 start in 1920, branched from ensemble 012 with a pertlim, and go to 2015 * These runs are fully coupled * Monthly Ocean, Sea Ice, and Atmosphere variables are availableFuture runs: b.e21.BSSP370cmip6.f09_g17.CMIP6-SSP3-7.0.? * The first three SSP runs (012,013,014) were configured the same as the CMIP6 scenario runs, but there was a bug in the future forcing * New runs were then done with the updated forcing (112, 113, 114, 115, 116) * The runs use standard CMIP6 SSP370 standard forcing, start in 2015 and go through 2100 * These runs are fully coupled * Monthly Ocean, Sea Ice, and Atmosphere variables are availableThere are also some experiments that are not fully coupled experiments: * e.e21.E1850.f09_g17.CMIP6-piControl.* * d.e21.f09_g17.{001,002,003,004,005}Additional ExperimentsThe additional CESM2 Tuned Sea Ice Albedo experiments were run to extend the dataset and include more frequent (daily) output as well as some modifications to the external forcing. This was done as part of a NSF Navigating the New Arctic (NNA) project (NSF #1928119): NNA Track 1: Collaborative Research: Maritime transportation in a changing Arctic: Navigating climate and sea ice uncertainties. These simulations were run by Alice DuVivier, and the data will be used in various studies. (For local disk access, these data are all in the main directory and are separated by case name.)Historical runs: b.e21.BHISTsmbb.f09_g17.CMIP6-historical.1? * There are five ensemble members: 118, 119, 120, 121, 122 * These ensembles start in 1990 and run through 2015 * These use smoothed biomass burning forcing (smbb), which uses the smoothed biomass burning forcing starting the 1990s. Information about the smoothed biomass burning impacts can be found in DeRepentigny et al. 2022 [https://doi.org/10.1126/sciadv.abo2405] * These runs are fully coupled * Daily and Monthly Ocean, Sea Ice, and Atmosphere variables are available * Restart files are also available upon request to co-authorsFuture runs: b.e21.BSSP370cmip6.f09_g17.CMIP6-SSP3-7.0.2? * There are five ensemble members: 218, 219, 220, 221, 222 * The runs use standard CMIP6 SSP370 standard forcing, start in 2015 and go through 2100. They are continued from the SMBB historical runs (above) * These runs are fully coupled * Daily and Monthly Ocean, Sea Ice, and Atmosphere variables are available * Restart files are also available upon request to co-authors * These runs were done on the Derecho supercomputer (instead of Cheyenne, which was used for all the above experiments)"]} 
    more » « less
  4. This dataset contains the daily Arctic sea ice area (SIA) and sea ice extent (SIE) data for all CMIP6 models and the historical period based on the NOAA/NSIDC Climate Data Record (CDR) created for Heuzé and Jahn, The first ice-free day in the Arctic Ocean could occur before 2030, accepted, Nature Communications. This is a derived dataset based on publicly available underlying data: - For the CMIP6 data, the SIA and SIE data included here is based on the daily siconc and siconca CMIP6 model output freely available on the CMIP6 data portals (https://pcmdi.llnl.gov/CMIP6/). These pan-Arctic daily SIA and SIE were calculated north of 30N, on each model's native grid, using each models grid area data (areacello or areacella). SIA was defined as sea ice concentration multiplied by the grid cell area and summed over all grid cells. SIE was defined as the sum of the grid cell area for all grid cells where the sea ice concentration was larger than 0.15. All processed SIA and SIE data is included in this dataset, even if the model was later excluded from the analysis for one reason or another (see Heuzé and Jahn 2024, Methods section). All data included has the same number of days as the underlying model. The historical data spans 1980-2014 and can be found in the CMIP6_historical_data.zip file, and the scenario data spans 2015 to the end of the 21st century simulation, for multiple scenarios (SSPs), and can be found in CMIP6_ssp_data.zip. Files are provided as .zip files to make it easy to download all data at once, as the SIA and SIE data is saved in one file per model and ensemble member, and for the scenario simulations, also per ssp. - For the NOAA/NSIDC Climate Data Record (CDR), the SIA and SIE data included here is based on the NOAA/NSIDC Climate Data Record of Passive Microwave Sea Ice Concentration, Version 4, doi:10.7265/efmz-2t65, Meier et al 2021. The sea ice concentration is multiplied by the grid size of each grid box, for this data, 25x25 kilometers (km) = 625 kilometers squared (km2), and then summed over the full domain. In doing that, we include the interpolated data in the pole hole as included in the sea ice concentration data, but exclude all land/coastal grid points (i.e., values > 2.5 in the underlying data). As the filename indicates, we removed all leap year data from this data (dropped every Feb 29th) so that all years have 365 days. Note that while the file name says this data is for 19790101 to 20231231, it does indeed include 1978 as first year (so 1978-01-01-2023-12-31), with daily data starting on 1978-10-25 (nan before then). We did not change the name of the data file to still allow all archived scripts using this datafile to run. Scripts that work on this data associated with Heuzé and Jahn (2024) can be found at: https://zenodo.org/records/14008665, doi:10.5281/zenodo.14006059 References: Meier, W. N., F. Fetterer, A. K. Windnagel, and S. Stewart. 2021. NOAA/NSIDC Climate Data Record of Passive Microwave Sea Ice Concentration, Version 4. Boulder, Colorado, USA. NSIDC: National Snow and Ice Data Center https://doi.org/10.7265/efmz-2t65 
    more » « less
  5. Abstract. Free-drift estimates of sea ice motion are necessary to produce a seamless observational record combining buoy and satellite-derived sea ice motionvectors. We develop a new parameterization for the free drift of sea ice based on wind forcing, wind turning angle, sea ice state variables(thickness and concentration), and estimates of the ocean currents. Given the fact that the spatial distribution of the wind–ice–ocean transfercoefficient has a similar structure to that of the spatial distribution of sea ice thickness, we take the standard free-drift equation and introducea wind–ice–ocean transfer coefficient that scales linearly with ice thickness. Results show a mean bias error of −0.5 cm s−1(low-speed bias) and a root-mean-square error of 5.1 cm s−1, considering daily buoy drift data as truth. This represents a 35 %reduction of the error on drift speed compared to the free-drift estimates used in the Polar Pathfinder dataset (Tschudi et al., 2019b). Thethickness-dependent transfer coefficient provides an improved seasonality and long-term trend of the sea ice drift speed, with a minimum (maximum)drift speed in May (October), compared to July (January) for the constant transfer coefficient parameterizations which simply follow the peak inmean surface wind stresses. Over the 1979–2019 period, the trend in sea ice drift in this new model is +0.45 cm s−1 per decadecompared with +0.39 cm s−1 per decade from the buoy observations, whereas there is essentially no trend in a free-driftparameterization with a constant transfer coefficient (−0.09 cm s−1 per decade) or the Polar Pathfinder free-drift input data(−0.01 cm s−1 per decade). The optimal wind turning angle obtained from a least-squares fitting is 25∘, resulting in a meanerror and a root-mean-square error of +3 and 42∘ on the direction of the drift, respectively. The ocean current estimates obtained from theminimization procedure resolve key large-scale features such as the Beaufort Gyre and Transpolar Drift Stream and are in good agreement with oceanstate estimates from the ECCO, GLORYS, and PIOMAS ice–ocean reanalyses, as well as geostrophic currents from dynamical ocean topography, with aroot-mean-square difference of 2.4, 2.9, 2.6, and 3.8 cm s−1, respectively. Finally, a repeat of the analysis on two sub-sections of thetime series (pre- and post-2000) clearly shows the acceleration of the Beaufort Gyre (particularly along the Alaskan coastline) and an expansion ofthe gyre in the post-2000s, concurrent with a thinning of the sea ice cover and the observed acceleration of the ice drift speed and oceancurrents. This new dataset is publicly available for complementing merged observation-based sea ice drift datasets that include satellite and buoydrift records. 
    more » « less