{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Lab 5-1: Calculating the Water Balance for the area draining to Cedar Lake, Washington\n", "---" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "We are going to use observational data to test the conservation of mass in a watershed. Monthly climate data are from the Cedar Lake station at 1560 ft elevation, at a latitude of 47.25 N and longitude of 121.44 W: [monthly climate from 7/1/1898 to 6/10/2016](https://wrcc.dri.edu/cgi-bin/cliMAIN.pl?wa1233)). Streamflow is from USGS gauge 12116400 for years 2001-2019. The gauge is listed as having an area of 217 $km^2$. Monthly ET is from climatological [maps](http://www.cpc.ncep.noaa.gov/soilmst/e.shtml). Pan evaporation is from the [Western Regional Climate Center WRCC](https://wrcc.dri.edu/Climate/comp_table_show.php?stype=pan_evap_avg). The Maple Leaf Reservoir site in Seattle, at a lower elevation, operated 1941 to 1960, and the Lake Kachess site, at a higher elevation, operated 1931-1977. Note the different units and time periods. This is often what you get from observational data. Your job is to make sense of it. " ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import pandas as pd\n", "import scipy.stats as stats\n", "import matplotlib.pyplot as plt\n", "%matplotlib inline" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Load the csv file with these various measurements from the Cedar River Watershed or nearby vicinity. " ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", " | month | \n", "Avg_max_T_degF | \n", "Avg_min_T_degF | \n", "Avg_total_P_in | \n", "Avg_streamflow_cfs | \n", "Avg_ET_mm_per_month | \n", "Pan_evap_MapleLeaf_inches_per_month | \n", "Pan_evap_LakeKachess_inches_per_month | \n", "
---|---|---|---|---|---|---|---|---|
0 | \n", "1 | \n", "40.0 | \n", "30.1 | \n", "13.52 | \n", "218 | \n", "3 | \n", "0.61 | \n", "0.00 | \n", "
1 | \n", "2 | \n", "43.8 | \n", "31.5 | \n", "10.06 | \n", "141 | \n", "7 | \n", "0.82 | \n", "0.00 | \n", "
2 | \n", "3 | \n", "47.6 | \n", "33.0 | \n", "10.36 | \n", "132 | \n", "17 | \n", "1.80 | \n", "0.00 | \n", "
3 | \n", "4 | \n", "53.7 | \n", "36.2 | \n", "8.01 | \n", "130 | \n", "30 | \n", "3.26 | \n", "2.37 | \n", "
4 | \n", "5 | \n", "60.8 | \n", "41.4 | \n", "6.52 | \n", "142 | \n", "50 | \n", "4.64 | \n", "4.82 | \n", "