AK K64G Demo Quick Start Guide 2c: Hydrology (soil hydraulic parameters via texture and or OM)

NGEE-Arctic: AK-K64G simulation example 2 - Parameterization

This example 2C runs a single-point, extracted from half-degree Global data, ELM v2 simulation, within which NGEE Kougarok Mile 64 study site located, Using user-provided soil physical properties, PCT_SAND, PCT_CLAY, ORGANIC, to allow ELM derive its soil hydraulic (and thermal) parameters, namely Clapp-Hornberger constants.

(1) It will start from an initial restart file from full-run's spinup stage

(2) It will modify 'surfdata.nc', land surface properties for ELM to read, in a temporary work directory.

(3) Re-run ELM, using newly-generated 'surfdata.nc'.

(4) Check soil moisture or temperature output in jupyter notebook.

REPEAT (2)-(4), for purpose to obtain a reasonbly well soil moisture simulation for site if observed.

I. Parametering ELM Demo Example C

NOTE, assuming that

(1) you've already done with AK-K64G Demo Quick Start Guide 2. one year simulation of transient run.

AND,

(2)you should have TWO (2) containers, one for running ELM model (but stopped), another for visualizing outputs in jupyter notebook which is active in your internet brower. Your laptop layout should be in some ways like in Parametering ELM Demo Example A

**STEP 1. Prior to modification, inspect soil moisture simulation in Jupyter Notebook, such as following:

• Chose H2OSOI, volume fraction of soil water (for all soil layers), from drop-down box.

• AND, then Clicking open menu Run, and clicking Run Selected Cell and All Below.

STEP 2. Modifying ELM soil physical properties, PCT_SAND, PCT_CLAY, ORGANIC, to adjust hydraulic parameters

• if not yet, Click Jupyter Notebook for visualizing's Left Panel Folder brower so that it back to folder work/kougarok, like following:

• Click open file, surfdata.txt, like following:

NOTE that surfdata.txt is an ASCII file converted from surfdata.nc which actually read by ELM, and editable in juypter notebook text editor.

• Scrolling down the file surfdata.txt to Line 27-35 (shown in shade), where are the header information of soil physical property PCT_SAND, PCT_CLAY, ORGANIC. Pay attention to their units and dimensions.
  

• Further Scrolling down to Line 342 - 376.

NOTE that,

(1) The 3 data blocks, for 3 properties respectively, have 10 values each. Those 10 values are for 10 soil layers (of one grid cell).

(2) By default, the bottom depth (meters) of layers, from ground surface, are:

1.7512817916255204E-002
4.5091787175931458E-002
9.0561820418344652E-002
0.16552923140455322
0.28912959650683373
0.49291214751726548
0.82889277396569816
1.3828311793343830
2.2961212109234443
3.8018819123227208     

• For an example, modify first layer value of 81.9990416468083 of the ORGANIC,

To, 129.9990416468083, a value nearly representing peat (130 kg m-3 at maximum in ELM).

NOTE that Jupyter Notebook text editor will automatically save this modified surfdata.txt file. For sure, you may click open menu File and chose save text to save your change(s).

AND, this modified surfdata.txt will be converted back to overwrite surfdata.nc when running OLMT example script in STEP 3 next.

STEP 3. Re-Run the model of 1 year transient simulation using OLMT

• In one of Terminal (for running ELM docker container), Type the following command, OR, copy/paste it:

docker run -t -i --hostname=docker --user=modeluser -v elmdata:/inputdata/ -v elmoutput:/output/ -v elmwork:/tmp/ fasstsimulation/elm-builds:elm_v2-for-ngee /bin/bash /tools/OLMT/examples/site_fullrun_docker_transient_nyears.sh --transient_years=1

OR, by pressing UP-ARROW key one time (to get the same command as used before), as following:

STEP 4. Visualizing new ELM output variable, H2OSOI, TSOI, or others.

• AFTER STEP 3 (re-run) finished, Click-open UPPER Tab Plot_ELM_variable.ipynb file again. And move starting cell to [1].

• THEN, for an example, if click that restart the kernel & run all cells button, , it will run script in all-steps at once, and will plot **GPP** by default, for individual PFT(s) if exist(s).

• Scrolling up to move your mouse pointer to cell [4], like

• By clicking the dropdown arrow, it will show whole list of variables from ELM outputs, and Choose one, such as H2OSOI.

• AND, then Clicking open menu Run, and clicking Run Selected Cell and All Below

THEN, it will plot like following:

YOU may notice what have changed, comparing to that in STEP 1, especially first layer (in blue line).

YOU may check other variable, such as TSOI (soil temperature by layers), simulated by ELM.

REPEAT STEP 2 - 4, by changing those values, either by tuning or by real field data, to obtain reasonably well soil moisture, matrical potential, or soil temperaure for site.

FINALLY, save modified surfdata.nc into your file systems, like following (back to Docker App)