What is a 'storm'?

[RondaStrauch]

• Define a 'storm' based on DHSVM modeling that reflects annual wettest soil conditions

Annual peak saturation extent - (3 hours peak in a day)
Thus, the time step is daily...the 3 hour max in a given day, and select maximum for each year
In the end we have a time series of maximum annual daily 3-hr max saturation extents.
We have this data for 1915 to 2010.

Pseudo-code:

1. Read in saturation extent.txt
2. Save Year, Month, Day, Hour, Percent
3. Select peak value per year
4. Print events in .csv: Year, Month, Day, Hour, Percent
5. Print events in textfile for DHSVM input: e.g. Map Layer DHSVM groundwater output #1
   5a. Count of dates
   5b. List of Dates with ordering and '1=' labels

Example:
Number of Maps 1 = 26 # Number of maps you would like to
Map Date 1 1 = 08/01/1987-00 # output for this variable
Map Date 2 1 = 08/01/1988-00 # Vary the first number from
Map Date 3 1 = 08/01/1989-00
Map Date 4 1 = 08/01/1990-00
Map Date 5 1 = 08/01/1991-00
Map Date 6 1 = 08/01/1992-00
Map Date 7 1 = 08/01/1993-00
Map Date 8 1 = 08/01/1994-00
Map Date 9 1 = 08/01/1995-00

[ChristinaB]

@NCristea
Please see the saturation_extent files in the code. We can make a new folder for the SaturatedStorm work.
Time, Percent Saturated Area
See this short version of the data for code testing

We need to run 1915-2010, and 2010-2099 for 20 other model folders.

[ChristinaB]

@NCristea Thanks for pushing your Notebook and looking into the parsing to final DHSVM format.
I will work on Precip and Streamflow plotting
What is a storm?
Option 1. Peak annual daily saturated area (%) (@NCristea code)
Option 2. Peak annual daily maximum streamflow at outlet (use suspended sediment model)
Option 3. Peak annual daily precipitation average/sqkm2 (use OGH)

[ChristinaB]

Here are sediment results to look at.
https://public.tableau.com/profile/christina.bandaragoda#!/vizhome/SaukFutureSediment/LowerSaukbcwackyIPSL

[ChristinaB]

Discussion Section - Options 1-3
For now we will stick with Option 1 and the original plan.

[ChristinaB]

@NCristea Small updated needed.

Your code prints:
Map date 1 01/08/2039-21
Map date 2 05/20/2040-21
Map date 3 11/03/2040-21
Map date 4 06/07/2042-15
Map date 5 12/10/2042-21

We need to add the ' 1 = ' in the middle
Map Date 1 1 = 08/01/1987-00 # output for this variable
Map Date 2 1 = 08/01/1988-00 # Vary the first number from
Map Date 3 1 = 08/01/1989-00
Map Date 4 1 = 08/01/1990-00

[ChristinaB]

For May - finalize code and add text to methods section in OneDrive draft paper

[ChristinaB]

[RondaStrauch]

• Reviewing the procedures we used to select storm dates.
• Discussed that we may want to limit historical dates 1969 to 2001 and exclude 2001-2011 because there is more effect of climate change trends in more recent years and thus, this would influence change detection with future time periods of mid and end of century.
• Create notebook edits that will have a box plot or histogram of individual model, include historic.
• Called distribution, not probability distribution, around a certain date because of concerns of stationarity.

[RondaStrauch]

• Decided to extract the max saturation extent dates from the rerun DHSVM model for the entire Skagit watershed rather than the Sauk sub watershed to ensure that we capture larger scale storms and overlap with the landslide domain.

• Christina is extracting the saturation extent files to provide new dates to Nicoleta that the storm code will used to create figures on saturation extent analytics (e.g., box plot).

• Nicoleta using the Skagit water table depth files derived from Sauk saturation extent to create workflow to produce figures (will need to redo workflow when we have Skagit watershed depth to groundwater)

@RondaStrauch - Send Nicoleta David Montgomery's figure of Olympic Peninsula with distributions by grid cell as an example.