BISICLES ice-sheet model simulations for the Amundsen Sea Embayment. Full details of the model set-up and ensemble design are described in the attached manuscript which has been accepted for publication in Journal of Glaciology. In brief, a 213-member ensemble of simulations was created by varying four different model parameters. The parameters are the u0 value in a regularised Coulomb friction law, the rate of imposed thinning of floating ice (∂h/∂t(Ωf)), and scaling factors for sliding and viscosity coefficients (C and ϕ) between 0.9 and 1.1. We attach a summary text file of results, as well as NetCDF files of simulated variables land ice thickness and u and v components of velocity. ASE2050_bisicles.csv contains annual (2007 to 2050, columns 5 to 48) sea level equivalent (mm) mass losses of ice from the Pine Island and Thwaites Glacier catchment basins. The parameters, given in columns 1 to 4, respectively, are the u0 (m/a), the rate of imposed thinning of floating ice (m/a), and the scaling factors for sliding and viscosity coefficients. The NetCDF files in ASE_BISICLES.tar.gz contain annual (2007 to 2050) simulated output variables for the Amundsen Sea region at a spatial resolution of 1 km, with one file per ensemble member. The variables follow the ISMIP6 naming protocol: (https://www.climate-cryosphere.org/wiki/index.php?title=ISMIP6-Projections-Antarctica#A2.3_Model_output_variables_and_README_file). We include state variables lithk, uvelmean, and vvelmean. Each file is named according to the variable, the simulation parameters, and the resultant 2050 SLE value of ice loss (mm). For example, lithk_ASE_BISICLES.uj_20.dhfdt_5.C_0.90.phi_0.90.slr_43.06.nc is the land ice thickness data for simulation u0=20 m/a, ∂h/∂t(Ωf) = 5 m/a, C scaled by 0.9, ϕ scaled by 0.9, and a final SLE of 43.06 mm.