skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


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. 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
  4. 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
  5. The High-resolution (6 kilometer (km)) Ice-Ocean Modeling and Assimilation System (HIOMAS) is used to simulate the evolution of sea ice for the Arctic Ocean and adjacent areas, including the Barents Sea, Norwegian Sea, Greenland Sea, Baffin Bay, and waters along Northwest Passage, over the period 2010 to 2069. The hindcast and future forcing over the period is from one of the Coupled Model Intercomparison Project Phase 6 (CMIP6) models, the CNRM-CM6-1-HR global climate model (GCM) run at the National Center for Meteorological Research, Météo-France and CNRS Laboratory (CNRM). Monthly mean sea ice thickness (meters (m)) is provided over 2010 to 2069 in NETCDF file format, with model grid information such as latitudes and longitudes of model grid cells included. I have archived future projection of monthly mean sea ice thickness files over 2010 to 2069 in https://pscfiles.apl.uw.edu/zhang/HIOMAS_6km/. The files are in netcdf format, which are created by running HIOMAS using the future projection forcing of the CNRM-CM6-1-HR GCM run conducted at the National Center for Meteorological Research, Météo-France and CNRS Laboratory (CNRM). Martin interpolated the forcing onto the new, expanded HIOMAS grid. There is a readme.txt file. 
    more » « less