## Jonathan Martin Frame

#### National Water Center;University of Alabama - Tuscaloosa

### Recent Activity

**ABSTRACT: **

This project contains preliminary 1D and 2D hydraulic models for a reach the Arroyo Seco (Los Angeles River tributary) upstream of the Brown Mountain Debris Barrier Dam. These models are **not** ready for use, they are **only** a starting point for further analysis.

**ABSTRACT: **

This repository holds supplemental figures for the paper "On Strictly Enforced Mass Conservation Constraints for Modeling the Rainfall-Runoff Process". These include distributions of event-based runoff ratios, scatter plots comparing event-based runoff ratios, initial (antecedent) flows, and rainfall totals, and hydrographs of every qualifying "event" that spans both Daymet and NLDAS forcing data.

Contents: 1) "hydrographs" is a directory with images showing hydrographs for every event in the record where the event criteria is met by both NLDAS and Daymet forcings coincidentally, 2) "basin_event_distributions_and_scatter_plots" is a directory with images showing distributions of runoff coefficients (runoff ratios) from 100 nearest-neighbour events; empirical membership values of potential runoff coefficients from nearest-neighbour storms based on Mahalabonis distances based on antecedent streamflow and rainfall total, and 3) "regional_q_q_plots" is a directory with images showing quantile-quantile plots of the percentile bins of the count of the events against the corresponding idealized percentile count. This is commonly referred to as a quantile-quantile plot. We then also plot the distance (and cumulative distance) of these quantiles from the one-to-one line.

**ABSTRACT: **

This repository holds supplemental figures for the paper "On Strictly Enforced Mass Conservation Constraints for Modeling the Rainfall-Runoff Process". These include distributions of event-based runoff ratios, scatter plots comparing event-based runoff ratios, initial (antecedent) flows, and rainfall totals, and hydrographs of every qualifying "event" that spans both Daymet and NLDAS forcing data.

**ABSTRACT: **

Runs from two papers exploring the use of mass conserving LSTM. Model results used in the papers are 1) model_outputs_for_analysis_extreme_events_paper.tar.gz, and 2) model_outputs_for_analysis_mass_balance_paper.tar.gz.

The models here are trained/calibrated on three different time periods. Standard Time Split (time split 1): test period(1989-1999) is the same period used by previous studies which allows us to confirm that the deep learning models (LSTM andMC-LSTM) trained for this project perform as expected relative to prior work. NWM Time Split (time split 2): The second test period (1995-2014) allows us to benchmark against the NWM-Rv2, which does not provide data prior to 1995. Return period split: The third test period (based on return periods) allows us to benchmark only on water years that contain streamflow events that are larger (per basin) than anything seen in the training data (<= 5-year return periods in training and > 5-year return periods in testing).

Also included are an ensemble of model runs for LSTM, MC-LSTM for the "standard" training period and two forcing products. These files are provided in the format "<model_type>_<forcing_product>_standard_training.tar.gz". Note that these individual ensemble member runs we used to produce the runs in the files "model_outputs_for_analysis_<*>_paper.tar.gz".

IMPORTANT NOTE: This python environment should be used to extract and load the data: https://github.com/jmframe/mclstm_2021_extrapolate/blob/main/python_environment.yml, as the pickle files serialized the data with specific versions of python libraries. Specifically, the pickle serialization was done with xarray=0.16.1.

Code to interpret these runs can be found here:

https://github.com/jmframe/mclstm_2021_extrapolate

https://github.com/jmframe/mclstm_2021_mass_balance

Papers are available here:

https://hess.copernicus.org/preprints/hess-2021-423/

**ABSTRACT: **

Peak flow observations for each of the CAMELS basins from the USGS National Water Information System (https://nwis.waterdata.usgs.gov/usa/nwis/peak; accessed 10 June 2021).

We used all available annual peak flows for each of the CAMELS basins and fit these values to the Pearson Type III distribution with log transformation using the method of moments, as described in U.S. Interagency Committee on Water Data,Bulletin 17b (IACWD, 1982). The probability density function is: f(x|\tau, \alpha, \beta) = \frac{(\frac{x-\tau}{\beta})^{\alpha-1}exp(-\frac{x-\tau}{\beta})}{|\beta|\gamma(\alpha)}, The CAMELS basins are suitable for this method under the assumptions of flood flows that are not appreciably altered by reservoir regulation, watershed changes or where the possibility of unusual events, such as dam failures.

We used Matlab code from Mathworks File Exchange to fit the peak annual flow events to the distribution to obtain return period estimates for each basin (Burkey, 2009). These return period calculations can also be done with free and open source software available from the USGS (https://water.usgs.gov/software/PeakFQ/; accessed 10 June 2021). We classified the water year of each basin (basin-year) according to the return period of its observed peak annual discharge

### Contact

(Log in to send email) |

All | 0 |

Collection | 0 |

Resource | 0 |

App Connector | 0 |

**Resource**

**ABSTRACT: **

Annual peak flows for the CAMELS basins; downloaded from USGS National WaterInformation System (https://nwis.waterdata.usgs.gov/usa/nwis/peak; accessed 10 June 2021)

**Resource**

##### Created: June 14, 2021, 12:55 a.m.

##### Authors: Frame, Jonathan

**ABSTRACT: **

Peak flow observations for each of the CAMELS basins from the USGS National Water Information System (https://nwis.waterdata.usgs.gov/usa/nwis/peak; accessed 10 June 2021).

We used all available annual peak flows for each of the CAMELS basins and fit these values to the Pearson Type III distribution with log transformation using the method of moments, as described in U.S. Interagency Committee on Water Data,Bulletin 17b (IACWD, 1982). The probability density function is: f(x|\tau, \alpha, \beta) = \frac{(\frac{x-\tau}{\beta})^{\alpha-1}exp(-\frac{x-\tau}{\beta})}{|\beta|\gamma(\alpha)}, The CAMELS basins are suitable for this method under the assumptions of flood flows that are not appreciably altered by reservoir regulation, watershed changes or where the possibility of unusual events, such as dam failures.

We used Matlab code from Mathworks File Exchange to fit the peak annual flow events to the distribution to obtain return period estimates for each basin (Burkey, 2009). These return period calculations can also be done with free and open source software available from the USGS (https://water.usgs.gov/software/PeakFQ/; accessed 10 June 2021). We classified the water year of each basin (basin-year) according to the return period of its observed peak annual discharge

**Resource**

**ABSTRACT: **

Runs from two papers exploring the use of mass conserving LSTM. Model results used in the papers are 1) model_outputs_for_analysis_extreme_events_paper.tar.gz, and 2) model_outputs_for_analysis_mass_balance_paper.tar.gz.

The models here are trained/calibrated on three different time periods. Standard Time Split (time split 1): test period(1989-1999) is the same period used by previous studies which allows us to confirm that the deep learning models (LSTM andMC-LSTM) trained for this project perform as expected relative to prior work. NWM Time Split (time split 2): The second test period (1995-2014) allows us to benchmark against the NWM-Rv2, which does not provide data prior to 1995. Return period split: The third test period (based on return periods) allows us to benchmark only on water years that contain streamflow events that are larger (per basin) than anything seen in the training data (<= 5-year return periods in training and > 5-year return periods in testing).

Also included are an ensemble of model runs for LSTM, MC-LSTM for the "standard" training period and two forcing products. These files are provided in the format "<model_type>_<forcing_product>_standard_training.tar.gz". Note that these individual ensemble member runs we used to produce the runs in the files "model_outputs_for_analysis_<*>_paper.tar.gz".

IMPORTANT NOTE: This python environment should be used to extract and load the data: https://github.com/jmframe/mclstm_2021_extrapolate/blob/main/python_environment.yml, as the pickle files serialized the data with specific versions of python libraries. Specifically, the pickle serialization was done with xarray=0.16.1.

Code to interpret these runs can be found here:

https://github.com/jmframe/mclstm_2021_extrapolate

https://github.com/jmframe/mclstm_2021_mass_balance

Papers are available here:

https://hess.copernicus.org/preprints/hess-2021-423/

**Resource**

##### Created: July 24, 2022, 8:01 p.m.

##### Authors: Frame, Jonathan Martin

**ABSTRACT: **

This repository holds supplemental figures for the paper "On Strictly Enforced Mass Conservation Constraints for Modeling the Rainfall-Runoff Process". These include distributions of event-based runoff ratios, scatter plots comparing event-based runoff ratios, initial (antecedent) flows, and rainfall totals, and hydrographs of every qualifying "event" that spans both Daymet and NLDAS forcing data.

**Resource**

##### Created: July 26, 2022, 6:05 a.m.

##### Authors: Frame, Jonathan Martin

**ABSTRACT: **

This repository holds supplemental figures for the paper "On Strictly Enforced Mass Conservation Constraints for Modeling the Rainfall-Runoff Process". These include distributions of event-based runoff ratios, scatter plots comparing event-based runoff ratios, initial (antecedent) flows, and rainfall totals, and hydrographs of every qualifying "event" that spans both Daymet and NLDAS forcing data.

Contents: 1) "hydrographs" is a directory with images showing hydrographs for every event in the record where the event criteria is met by both NLDAS and Daymet forcings coincidentally, 2) "basin_event_distributions_and_scatter_plots" is a directory with images showing distributions of runoff coefficients (runoff ratios) from 100 nearest-neighbour events; empirical membership values of potential runoff coefficients from nearest-neighbour storms based on Mahalabonis distances based on antecedent streamflow and rainfall total, and 3) "regional_q_q_plots" is a directory with images showing quantile-quantile plots of the percentile bins of the count of the events against the corresponding idealized percentile count. This is commonly referred to as a quantile-quantile plot. We then also plot the distance (and cumulative distance) of these quantiles from the one-to-one line.

**Resource**

##### Created: June 26, 2023, 6:59 p.m.

##### Authors: Frame, Jonathan Martin

**ABSTRACT: **

This project contains preliminary 1D and 2D hydraulic models for a reach the Arroyo Seco (Los Angeles River tributary) upstream of the Brown Mountain Debris Barrier Dam. These models are **not** ready for use, they are **only** a starting point for further analysis.