Models, data, and graphical results

for draft publication titled "A data assimilation framework for improved pumping-induced subsidence estimation in data-scarce settings."

Authors: Jenny T. Soonthornrangsan¹, Femke C. Vossepoel¹, Mark Bakker²

¹Department of Geoscience & Engineering, Faculty of Civil Engineering and Geoscience, Delft University of Technology, Stevinweg 1, 2628 CN Delft, The Netherlands

²Department of Water Management, Faculty of Civil Engineering and Geoscience, Delft University of Technology, Stevinweg 1, 2628 CN Delft, The Netherlands

Current release: [

Libraries: Pastas: 1.3.0 Pandas: 2.0.3 Numpy: 1.26.3

Various python scripts are provided that create different graphical results.

• Figures_Paper2.py: Produces the figures shown in the main text of the paper

• ESMDA_BKK.py: Runs ESMDA algorithm using real observations for well nests in BKK

• ESMDA_Synthetic.py: Runs ESMDA algorithm using synthetic observations for well nests. Script also generates the synthetic observations.

• LS_BKK.py: Runs least squares algorithm using real observations for well nests in BKK

• LS_Synthetic.py: Runs least squares algorithm using the same synthetic observations from ESMDA

• Plot_ESMDA_LS_BKK.py: Plots the results of the ESMDA and LS runs for the real experiment for well nests in BKK

• Plot_ESMDA_LS_Synthetic.py: Plots the results of the ESMDA and LS runs for the synthetic experiments for well nests

figures\ : figures created from the scripts

figures_test\: figures used in publication to be compared to output figures from script. Should match.

models\: ESMDA/LS results created from ESMDA_BKK.py, ESMDA_Synthetic.py, LS_Synthetic.py, and LS_BKK.py script. Model outputs in Dpred files for each ESMDA step, parameters in Mprior files for each ESMDA step

inputs\: inputs needed for scripts

• Groundwater observations for each well nest. Note that in the paper, each well nest is referred to without "LC" (LC******.xlsx)

• Subsidence observations for each well nest (SurveyingLevels.xlsx)

• Parameters required for subsidence calcultion for each well nest e.g. Sske, Sskv, K, thickness (SUBParameters.xlsx)

• Initial estimate of steady state heads for each well nest (SS_Head_GWWellLocs.xlsx)

• Land surface elevation for each well nest (LandSurfElev_GWWellLocs.xlsx)

  • Coastal DEM 2.1: https://www.climatecentral.org/coastaldem-v2.1
  • Kulp, S. A., and B. H. Strauss. 2021. CoastalDEM v2. 1: A high-accuracy and high-resolution global coastal elevation model trained on ICESat-2 satellite lidar.

• Location, well depth, screen depth, available observation years for each well for each well nest (GroundwaterWellLocs.xls)

• Basin-wide pumping data 1954-2060 (BasinPumping.xlsx)

  • 500,000 m³/day scenario: EstTotalPump_54-60_Int50 sheet
  • 250,000 m³/day scenario: EstTotalPump_54-60_IntF25 sheet
  • 1,000,000 m³/day scenario: EstTotalPump_54-60_IntF100 sheet
  • 500,000 to 250,000 m³/day in 2040 scenario: EstTotalPump_54-60_IntF50_25 sheet
  • No pumping scenario: EstTotalPump_54-60_IntF0 sheet