Input files and model output for a schematized 1-D morphodynamic model of the Lower Rhine System (from Bonn in Germany to Vuren, Schoonhoven, and Keteldiep in the Netherlands) 

These files are related to the manuscript " Climate Change Impacts on Flow Partitioning in a River Bifurcation System " submitted to Geophysical Research Letters.

The Input files include

	- Model_cross_sections: cross-section data, with the main channel and floodplain width at each cross-section
	- initial mean bed elevation: Initial mean bed elevation of all the branches
	- initial bed composition: Initial bed and subsurface composition per grid cell formatted as input to SOBEK RE
	- internal and external morphological boundary conditions: Sediment flux at the upstream boundary and nodal point relations at the two bifurcation nodes formatted as input to SOBEK RE
	- Boundary condition_Upstream discharge: Hydrograph at the upstream boundary for the reference case and the scenarios
	- Boundary condition_Downstream water level: Water level at the downstream boundaries for the reference case and the scenarios (Please check SI to separate which boundary condition fits with which scenario)
	- GRAINP: Mixed size sediment input file for SOBEK 
	
	
Model results

	- flow_discharge_output.xlsx : Flow discharge at Lobith and the upstream part of the bifurcates for all the model runs
	- bed_level_output.xlsx: yearly averaged bed level for the entire reach for all the model runs
	- Dg_transport_output/DG_surface_output: 20-year averaged geometric mean grain size of the flux and bed surface over the entire reach for the reference case, hydrograph scenario Hn, and combined scenario Hn-SSP 585
	- transportvolume_output: 20-year averaged volume of total sediment flux over the entire reach for the reference case, hydrograph scenario Hn, and combined scenario Hn-SSP 585

In addition, the field data collected by Rijkswaterstaat on bed surface grain size in the Dutch Rhine in 2020 (2020_bed_surface_data.xlsx) is also included. These data have been used as information to derive the model input grain size fraction content.
	
***************Please use notepad or notepad++ to open the text files. 	
	
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
.txt file specs are explained below
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%% Model_cross_sections.txt %%%%%

description "#cross_section_name" MainChannelWidth(m) FloodplainWidth(m) SedimentTransportWidth(m)
{
          elevation_LowestLevel (m)     flowWidthLowestLevel (m)         totalWidthLowestLevel (m)
          .				.				 .	
	  .				.				 .
			
          	        
         elevation_HighestLevel (m)     flowWidthHighestLevel (m)         totalWidthHighestLevel (m)
}

.
.
.

crosssection "#profile_name" "#associated_cross_section_name_as_tagged_in_description" "branch" chainage_relative_to_branch(m) vertical_translation(m)
.
.


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%% initial bed composition.txt %%%%%

$GSINIT BRANCH num_branch AT chainage_relative_to_branch(m)
LAYNUM = 1 fraction_content_finest_fraction ..  .. .. ..  fraction_content_coarsest_fraction -------- %this is the layer closest to the surface, each layer is 0.5 m thick
.
.
.
LAYNUM = 20 fraction_content_finest_fraction ..  .. .. ..  fraction_content_coarsest_fraction -------- %this is the deepest layer



$GSLEVUNLA Branch num_branch AT chainage_relative_to_branch(m) ZZ
 ZZ=depth of the deepest layer per grid cell
 
 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%% internal and external morphological boundary conditions.txt %%%%%
%%%%% GRAINP.txt %%%%%

These text files contain explanations. 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%% #Location#scenario.txt in the 'Downstream water level per boundary' zip folder%%%%%

"YYYY/MM/DD; hh:mm:ss" water_level(m)

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%% #scenario#_QLobith.txt in the 'Boundary condition_Upstream discharge' zip folder%%%%%

"YYYY/MM/DD; hh:mm:ss" discharge(m^3/s)