Drainage areas (or watersheds) were calculated for each stream segment in the National Hydrography Dataset High Resolution Plus Beta (NHD HR Plus Beta) using non-end junction points and the provided hydrologically enforced Digital Elevation Model (DEM). The DEM was used to create a flow direction raster. The pour points used were determined from junction points at each stream intersection excluding stream ends (i.e. channel heads are not included). Watershed calculations were automated using the ArcGIS "watershed" tool and then edited manually for accuracy. In western San Diego County (non-desert) there were 39,980 drainage areas calculated. The percent impervious in each watershed was calculated using the National Land Cover Dataset 2019 Impervious surface.