Forested watersheds have often been managed for flood mitigation. Studies have shown that forests have the potential to minimize peak flows during storm events, yet the relationship between forests and flooding is inexact. Forest roads, usually found in managed systems, can potentially magnify the effects of forest harvesting on water yields. A distributed hydrologic model (DHSVM) was calibrated for a 760 ha watershed in the Blue Ridge Mountains of North Carolina. The impacts of forest road density were evaluated by running the model using uniform input parameters but changing road densities. Road densities tested were 0.5, 1.0, 3.0, 4.3, 6.0 and 12.0 km km-2. Results indicate that increases in road density increased average streamflows at densities ≥ 4.3 km km-2. During small storm events, discharge was impacted at densities ≥ 6.0 km km-2 and streamflows were impacted during large rainfall events ≥ 3.0 km km-2 road densities. These findings indicate that forest roads can influence water yields and additional management efforts may be needed that can slow the water yield from forest roads. / Master of Science

Hydrology is the study of water. Hydrology tracks various attributes of water such as its quality and movement. As a tool Hydrology allows researchers to investigate topics such as the impacts of wildfires, logging, and commercial development. With perfect and complete data collection researchers could answer these questions with complete certainty. However, due to cost and potential sources of error this is impractical. As such researchers rely on simulations. The Distributed Hydrology Soil Vegetation Model(also referenced to as DHSVM) is a scientific mathematical model to numerically represent watersheds. Hydrology, as with all fields, continues to produce large amounts of data from researchers. As the stores of data increase the scientific models that process them require occasional improvements to better handle processing the masses of information. This paper investigates DHSVM as a serial C program. The paper implements and analyzes various high performance computing advancements to the original code base. Specifically this paper investigates compiler optimization, implementing par- allel computing with OpenMP, and adding distributed computing with OpenMPI. DHSVM was also tuned to run many instances on California Polytechnic State Uni- visity, San Luis Obispo’s high performance computer cluster. These additions to DHSVM help speed-up the results returned to researches, and improves DHSVM’s ability to be used with uncertainty analysis methods. iv This paper was able to improve the performance of DHSVM 2 times with serial and compiler optimization. In addition to the serial and compiler optimizations this paper found that OpenMP provided a noticeable speed up on hardware, that also scaled as the hardware improved. The pareallel optimization doubled DHSVM’s speed again on commodity hardware. Finally it was found that OpenMPI was best used for running multiple instances of DHSVM. All combined this paper was able to improve the performance of DHSVM by 4.4 times per instance, and allow it to run multiple instances on computing clusters.

Forest harvesting has been shown to cause various changes in water quantity and water quality parameters, highlighting the need for comprehensive forest practice rules. Studies show a myriad of impacts to ecosystems as a result of watershed level changes, such as forest harvesting. Being able to better understand the impact that forest harvesting can have on stream temperature is especially critical in locations where federally threatened or endangered fish species are located. The overall goal of this research project is to assess responses in stream temperature to various riparian and forest harvest treatments in a maritime, mountainous environment. The results of this study aim to inform decision makers with additional information pertaining to the effects of forest harvest on water temperature. Modeling is done as a part of the third Caspar Creek Paired Experimental Watershed study. Located in Mendocino County, the site provides a place for California researchers and decision makers to learn about the cumulative watershed effects of forest management operations on peak flows, sediment production, anadromous fish, macro-invertebrate communities, nutrient cycling and more. Historic data was used to calibrate the Distributed Hydrology Soil Vegetation Model (DHSVM) and River Basin Model (RBM) to measured stream temperatures in the South Fork of Caspar Creek (SFC) for hydrologic years 2010-2016. Critical summer time periods, when temperatures are highest and flows are low, are the primary concern for this work. The key modeling scenarios evaluated were (1) varying percentages of Watercourse and Lake Protection Zones (WLPZ) canopy cover, (2) the 2018-2019 SFC forest harvest and (3) an experimental design converting dominant riparian vegetation along 300-yard stream reaches. Modeling results showed that stream temperatures begin to rise above third-growth conditions when canopy cover is reduced to 25% and 0% retention levels. Larger increases in Maximum Weekly Maximum Temperature (MWMT) values, compared to Maximum Weekly Average Temperature (MWAT) values, were seen across all scenarios. There was essentially no difference between altering buffer areas along only class I streams, compared to along all stream classes. At the 0% canopy retention, MWMT values consistently rose above recommended thermal limits for Coho salmon (Oncorhynchus kisutch) and state regulations prohibiting more than a 5 degree F increase in waters. Clearcutting the entire watershed produced less of an effect than simulations clearing on only the riparian area, suggesting that groundwater inflows act to mitigate stream temperature rises in the SFC. The 2018-2019 harvest showed a relatively consistent increase in MWAT values (avg. 0.11 degree C) and more varied increases in MWMT values (avg. 0.32 degree C). Simulations converting dominant riparian vegetation by clearing could not be considered conclusive due to sensitivity analyses suggesting potentially unrealistic tracking of downstream temperatures. Additional sensitivity analyses suggest that tree height and the monthly extinction coefficient (a function of Leaf Area Index) are most influential on stream temperature changes in SFC. This is consistent with other modeling studies and suggests stream temperature management focus on tall, dense buffers as opposed to wider buffer widths.

Cette étude a été menée dans le bassin versant du lac Tana situé dans le haut bassin du Nil Bleu en Ethiopie, où les enjeux liés à l'exploitation des ressources en eau sont critiques et la disponibilité de données pertinentes est faible. L'étude est divisée en trois parties principales qui comprennent l'évaluation de la qualité et l'analyse des données existantes, l'évaluation de l'applicabilité d'un modèle hydrologique à base physique et le calcul des bilans en eau et des matières en suspension du lac Tana. Les couvertures spatiales et temporelles des variables hydro-météorologiques ont été évaluées à partir des méthodes recommandées dans la littérature. L'analyse des jeux de données hydro-météorologiques comprend le contrôle de la qualité des données, et la détermination de leurs caractéristiques spatiales et temporelles. L'applicabilité du modèle hydrologique à base physique a été évaluée sur deux bassins versants. Cet exercice de modélisation a montré l'intérêt des méthodes de désagrégation temporelle pour générer des données météorologiques à des échelles de temps plus courtes. Compte tenu de la qualité des données utilisées, la performance du modèle dans deux des bassins versants a été jugée satisfaisante. Le bilan en eau du lac Tana a été calculé au pas temps mensuel ce qui a impliqué l'estimation des différents termes. Un modèle dynamique et conceptuel a été utilisé pour estimer les contributions des bassins versants pour lesquels les données hydrométriques sont manquantes. Les résultats du calcul sont généralement bons et les niveaux observés du lac sont reproduits de manière satisfaisante. L'estimation du bilan des matières en suspension pour le lac a nécessité le développement et l'utilisation des courbes débit-concentration et un modèle régional. Les résultats présentés dans ce travail sont utiles pour la planification et la gestion des ressources en eau dans cette région d'Ethiopie.

Declines in glacier area and volume are widespread. These changes will have important hydrologic consequences since glaciers store tremendous amounts of fresh water and buffer seasonally low flows in many densely populated regions. In this thesis I focus on a region that is hydrologically vulnerable to glacier change, namely the Cordillera Blanca, Peru. I present three manuscripts that focus on measuring glacier area change, modeling the effect of this area change on the hydrology of one watershed, and isotopic sampling to elucidate hydrologic processes in this watershed and the entire Cordillera Blanca. In the first manuscript, I describe a methodology for mapping glaciers using satellite imagery. Satellite data, in conjunction with automated glacier mapping methods, are being used more frequently to map changes in glacier size. In contrast to the majority of studies using automated methods, I correct satellite images for atmospheric effects. Mapping glaciers with atmospherically-corrected satellite images resulted in an approximately 5% increase in glacier area, relative to glaciers mapped with non-atmospherically-corrected images. I also applied a consistent threshold that was validated using high-resolution satellite imagery. This helps to reduce error associated with change analysis. For the entire Cordillera Blanca, I calculated a 25% decrease in glacier area from 1987 to 2010. The rate of glacier area loss has increased significantly based on the most recent estimates. In the second manuscript, I use a physically-based, hydrologic model, the Distributed Hydrology Soil Vegetation Model (DHSVM) with a newly-coupled dynamic glacier model to simulate stream discharge and glacier change in the Llanganuco watershed of the Cordillera Blanca. I also examined statistical trends associated with historical records of temperature, precipitation, and discharge. I observed significant positive trends in annual temperature, but no trends in precipitation or discharge despite a 25% reduction in glacier area in this watershed over the same time. The model setup process and the results of sensitivity analyses are described. Of the input parameters I examined, I found that the model was particularly sensitive to changes in albedo and precipitation. Based on established efficiency criteria, the newly-coupled model did a decent job of simulating historical stream discharge and glacier area during 10 year calibration and validation periods. However, due to the lack of additional validation data and an inability to quantify uncertainty associated with model output, the model is not yet ready to be used for predicting future discharge based on different climate projections. In the third manuscript I describe the knowledge gained about hydrologic processes from isotopic sampling in the Llanganuco watershed, as well as other watersheds of the Cordillera Blanca. Thirty water samples from Llanganuco were collected in July 2011 and measured for stable isotopes of water, δ¹⁸O and δ²H. I first calculated the isotopic lapse rate, or the relationship between isotopic values and elevation. Lapse rates from this watershed are slightly more positive than global averages. This observation is best explained by the influence of glaciers. I also calculated the strength of the relationship between isotopic values and percent glacier cover. For Llanganuco, glacier cover is a better predictor of isotopic value than elevation. Based on examination of the same relationships at larger scales in the Cordillera Blanca, this relationship appears to be persistent at a regional scale. Finally, I used a simple two-component mixing model to estimate the relative contributions of glacier meltwater and groundwater in the Llanganuco watershed. Glacier meltwater made up approximately three-fourths of surface water that exited the watershed during this two week period in July, 2011. The importance of glacier meltwater is clearly demonstrated using stable isotopes, but further, more detailed monthly sampling is necessary to accurately determine annual and dry season streamflow contributions from glacier meltwater and groundwater. / Graduation date: 2013