This dataset contains numerical simulation results of the wave fields in the Davis Sea from end of December 2019 to start of February 2020. Hindcasts were obtained through the third-generation spectral wave model WAVEWATCH-III (hereafter WW3).
A high resolution Davis Sea regional grid (resolution 0.1 degree, 60-80E longitude, 70-60S latitude) was nested into global grid domain (resolution 0.5 degree, 80S-80N latitude). The global model is forced with 0.5 degree sea ice concentration and 10m-wind fields from ECMWF's ERA5 reanalysis. The Davis sea model is forced with 0.1 degree 10m-wind fields from ECMWF's archived forecasts, and high-resolution (3.125km) AMSR2 satellite data for sea ice concentration (Beitsch et al., 2013 updated). Ice-induced wave attenuation is parameterized following Sutherland et al. (2019, doi:10.1016/j.apor.2019.03.023) whilst the break-up of sea ice is parameterized as 'broken' or 'unbroken' based on the break-up parameter of Voermans et al.(2020, doi:10.5194/tc-14-4265-2020). The numerical simulations have been calibrated using the buoy-observations of Voermans (2022, dataset, doi:10.26179/cdmx-n995). Sensitivity of the simulations to sea ice properties was tested and all results are provided in the dataset.
The data tree:
* global: model outputs for the global domain
- ncfield: gridded wave and ice data for this domain in netCDF-4 format
- nests: binary data used by WW3 for boundary conditions for the Davis Sea grid
- restarts: binary data used by WW3 for restarting this domain
* davis_sea: model outputs for the Davis Sea domain
- ncfield: gridded wave and ice data for this domain in netCDF-4 format
- ncpoint: spectral wave data for a few points in the Davis Sea in netCDF-4 format
- nctrack: spectral wave data following the wave buoys of Voermans et al (2022) in the Davis Sea in netCDF-4 format
- restarts: binary data used by WW3 for restarting this domain
- IHOT: binary text field of broken and unbroken ice for restarting this domain
File naming convention (by example):
ww3.20200101_20200103_M3D_IHOT_H0P0325_A0P01_YY9P0_SS0P1_HH0P55.nc
* 20200101_20200103 identifies the datespan of the simulation in YYYYMMDD format
* A0P01 refers to the attenuation coefficient of the model (where P stands for 'point'), in this case, A=0.01
* YY is the Young's Modulus timed 10^9, here, Y-9.0e9 Pa
* SS is the ice strength 'sigma' times 10^6, here sigma=0.1e6
* HH is the ice thickness, here h=0.55 m
* H0P0325 is proportional to the epsilon calibration coefficient (H=0.5ice_thicknessepsilon).
* M3D refers to the 3rd instantiation of the model
* IHOT refers to hot start using the ice breakup field from the previous week.
ww3.*_M3D_IHOT_H0P065_A0P05_YY6P0_SS0P55.nc is considered the baseline file (note, this simulation only covers the first two weeks of the study period).
Reference: Beitsch, A., Kaleschke, L. and Kern, S. (2013). "AMSR2 ASI 3.125 km Sea Ice Concentration Data, V0.1", Institute of Oceanography, University of Hamburg, Germany, digital media
