This is the output from high-resolution model simulations of ocean conditions and melting beneath the floating part of Thwaites Glacier. The model is designed to study how these conditions change as the geometry of Thwaites Glacier evolved from 2011-2022. There is one simulation using the geometry from each year during this period, derived from satellite observations. The simulations are repeated for different ocean model forcing conditions, as described in the associated paper. PH was supported by the NERC/NSF Thwaites-MELT project (NE/S006656/1). ITGC contribution number 099. *******PLEASE BE ADVISED TO USE VERSION 2.0 DATA******* (Version 1 had the seabed bathymetry and ice shelf topography files incorrectly oriented.)

This is the output from high-resolution model simulations of ocean conditions and melting beneath the floating part of Thwaites Glacier. The model is designed to study how these conditions change as the geometry of Thwaites Glacier evolved from 2011-2022. There is one simulation using the geometry from each year during this period, derived from satellite observations. The simulations are repeated for different ocean model forcing conditions, as described in the associated paper. PH was supported by the NERC/NSF Thwaites-MELT project (NE/S006656/1). ITGC contribution number 099. *******PLEASE BE ADVISED TO USE VERSION 2.0 DATA******* Version 2 is available at https://doi.org/10.5285/473eb97c-63a8-4002-8b72-e7f07b2ab228. (Version 1 has the seabed bathymetry and ice shelf topography files incorrectly oriented.)

These are MITgcm ocean model outputs under Pine Island Glacier Ice Shelf. The simulations were designed to investigate the relative role of ocean conditions and ice shelf geometric changes between 2011 and 2021. Each set of runs contains one run for each year. BOTH-2011 to BOTH-2021 contain annual runs using the ocean and geometry from the corresponding year, OCEAN-2011 to OCEAN-2021 contain annual runs using the annual ocean conditions and 2012 ice shelf geometry and GEOM-2011 to GEOM-2021 contain annual runs using the 2012 ocean conditions and the annual ice shelf geometry. The ice shelf geometries used are derived from CryoSat-2 Digital Elevation Models (Lowery et al., 2025) and converted to ice shelf draft by assuming the ice is in hydrostatic equilibrium. The ocean boundary conditions are from observations from two moorings in Pine Island Bay. The output data contain fields of potential temperature, salinity, velocity and 2D variables from the boundary layer such as freshwater flux, thermal driving and sub-ice shelf velocity. The work took place as part of the Natural Environment Research Council (NERC) Satellite Data in Environmental Science (SENSE) Centre for Doctoral Training (grant no. NE/T00939X/1).