The main goal of this dissertation is to develop a seasonal water balance model for evaporation, runoff and water storage change based on observations from a large number of watersheds, and further to obtain a comprehensive understanding on the dominant physical controls on intra-annual water balance. Meanwhile, the method for estimating evaporation and water storage based on recession analysis is improved by quantifying the seasonal pattern of the partial contributing area and contributing storage to base flow during low flow seasons. A new method for quantifying seasonality is developed in this research. The difference between precipitation and soil water storage change, defined as effective precipitation, is considered as the available water. As an analog to climate aridity index, the ratio between monthly potential evaporation and effective precipitation is defined as a monthly aridity index. Water-limited or energy-limited months are defined based on the threshold of 1. Water-limited or energy-limited seasons are defined by aggregating water-limited or energy-limited months, respectively. Seasonal evaporation is modeled by extending the Budyko hypothesis, which is originally for mean annual water balance; while seasonal surface runoff and base flow are modeled by generalizing the proportionality hypothesis originating from the SCS curve number model for surface runoff at the event scale. The developed seasonal evaporation and runoff models are evaluated based on watersheds across the United States. For the extended Budyko model, 250 out of 277 study watersheds have a Nash-Sutcliff efficiency (NSE) higher than 0.5, and for the seasonal runoff model, 179 out of 203 study watersheds have a NSE higher than 0.5. Furthermore, the connection between the seasonal parameters of the developed model and a variety of physical factors in the study watersheds is investigated. For the extended Budyko model, vegetation is identified as an important physical factor that related to the seasonal model parameters. However, the relationship is only strong in water-limited seasons, due to the seasonality of the vegetation coverage. In the seasonal runoff model, the key controlling factors for wetting capacity and initial wetting are soil hydraulic conductivity and maximum rainfall intensity respectively. As for initial evaporation, vegetation is identified as the strongest controlling factor. Besides long-term climate, this research identifies the key controlling factors on seasonal water balance: the effects of soil water storage, vegetation, soil hydraulic conductivity, and storminess. The developed model is applied to the Chipola River watershed and the Apalachicola River basin in Florida for assessing potential climate change impact on the seasonal water balance. The developed model performance is compared with a physically-based distributed hydrologic model of the Soil Water Assessment Tool, showing a good performance for seasonal runoff, evaporation and storage change.

River export has a strong influence on the productivity of coastal waters. During storm events, rivers deliver disproportionate amounts of nutrients and organic matter to estuaries. Anthropogenic changes to the land use/cover (LULC) and water use also have a strong influence on the export of nutrients and organic matter to estuaries. This study specifically addressed the following questions: 1) How does river water chemistry vary across LULC patterns in the Mission and Aransas river watersheds? 2) How do fluxes of water, nutrients, and organic matter in the rivers vary between base flow and storm flow? 3) How do variations in nutrient/organic matter concentrations and stable isotope ratios of particulate organic matter (POM) in Copano Bay relate to river inputs? Water was collected from the Mission and Aransas rivers and Copano Bay from July, 2007 through November, 2008 and analyzed for concentrations of nitrate, ammonium, soluble reactive phosphorus (SRP), dissolved organic nitrogen, dissolved organic carbon, particulate organic nitrogen, particulate organic carbon (POC), and the stable C and N isotope ratios of the POM. The first half of the study period captured relatively wet conditions and the second half was relatively dry compared to long term climatology. Riverine export was calculated using the USGS LOADEST model. The percentage of annual constituent export during storms in 2007 was much greater than in 2008. Concentration-discharge relationships for inorganic nutrients varied between rivers, but concentrations were much higher in the Aransas River due to waste water contributions. Organic matter concentrations increased with flow in both rivers, but POM concentrations in the Aransas River were two fold higher due to large percentages of cultivated crop land. Values of [delta]¹³C-POC show a shift from autochthonous to allochthonous organic matter during storm events. Following storm events in Copano Bay, increases and quick draw down of nitrate and ammonium concentrations coupled with increases and slow draw down of SRP illustrate nitrogen limitation. Organic matter concentrations remained elevated for ~9 months following storm events. The [delta]¹³C-POC data show that increased concentrations were specifically related to increased autochthonous production. Linkages between LULC and nutrient loading to coastal waters are widely recognized, but patterns of nutrient delivery (i.e. timing, duration, and magnitude of watershed export) are often not considered. This study demonstrates the importance of sampling during storm events and defining system-specific discharge-concentration relationships for accurate watershed export estimation. This study also shows that storm inputs can support increased production for extended periods after events. Consideration of nutrient delivery patterns in addition to more traditional studies of LULC effects would support more effective management of coastal ecosystems in the future. / text

Base flow recession curve during a dry period is a distinct hydrologic signature of a watershed. The base flow recession analysis for both streamflow and spring flow has been extensively studied in the literature. Studies have shown that the recession behaviors during the early stage and the late stage are different in many watersheds. However, research on the transition from early stage to late stage is limited and the hydrologic control on the transition is not completely understood. In this dissertation, a novel cumulative regression analysis method is developed to identify the transition flow objectively for individual recession events in the well-studied Panola Mountain Research Watershed in Georgia, USA. The streamflow at the watershed outlet is identified when the streamflow at the perennial stream head approaches zero, i.e., flowing streams contract to perennial streams. The identified transition flows are then compared with observed flows when the flowing stream contracts to the perennial stream head. As evidenced by a correlation coefficient of 0.90, these two characteristics of streamflow are found to be highly correlated, suggesting a fundamental linkage between the transition of base flow recession from early to late stages and the drying up of ephemeral streams. At the early stage, the contraction of ephemeral streams mostly controls the recession behavior. At the late stage, perennial streams dominate the flowing streams and groundwater hydraulics governs the recession behavior. The ephemeral stream densities vary from arid regions to humid regions. Therefore, the characteristics of transition flow across the climate gradients are also tested in 40 watersheds. It is found that climate, which is represented by climate aridity index, is the dominant controlling factor on transition flows from early to late recession stages. Transition flows and long-term average base flows are highly correlated with a correlation coefficient of 0.82. Long-term average base flow and the transition flow of recession are base flow characteristics at two temporal scales, i.e., the long-term scale and the event scale during a recession period. This is a signature of the co-evolution of climate, vegetation, soil, and topography at the watershed scale. The characteristics of early and late recession are applied for quantifying human impacts on streamflow in agricultural watersheds with extensive groundwater pumping for irrigation. A recession model is developed to incorporate the impacts of human activities (such as groundwater pumping) and climate variability (such as evapotranspiration) on base flow recession. Groundwater pumping is estimated based on the change of observed base flow recession in watersheds in the High Plains Aquifer. The estimated groundwater pumping rate is found consistent compared with the observed data of groundwater uses for irrigation. Besides streamflow recession analysis, this dissertation also presents a novel spring recession model for Silver Springs in Florida by incorporating groundwater head, spring pool altitude, and net recharge into the existing Torricelli model. The results show that the effective springshed area has continuously declined since 1988. The net recharge has declined since the 1970s with a significant drop in 2002. Subsequent to 2002, the net recharge increased modestly but not to the levels prior to the 1990s. The decreases in effective springshed area and net recharge caused by changes in hydroclimatic conditions including rainfall and temperature, along with groundwater withdrawals, contribute to the declined spring flow.

Turbulent boundary layers approximating those found on the NASA Orion Multi-Purpose Crew Vehicle (MPCV) thermal protection system during atmospheric reentry from the International Space Station have been studied by direct numerical simulation, with the ultimate goal of reducing aerothermodynamic heating prediction uncertainty. Simulations were performed using a new, well-verified, openly available Fourier/B-spline pseudospectral code called Suzerain equipped with a ``slow growth'' spatiotemporal homogenization approximation recently developed by Topalian et al. A first study aimed to reduce turbulence-driven heating prediction uncertainty by providing high-quality data suitable for calibrating Reynolds-averaged Navier--Stokes turbulence models to address the atypical boundary layer characteristics found in such reentry problems. The two data sets generated were Ma[approximate symbol] 0.9 and 1.15 homogenized boundary layers possessing Re[subscript theta, approximate symbol] 382 and 531, respectively. Edge-to-wall temperature ratios, T[subscript e]/T[subscript w], were close to 4.15 and wall blowing velocities, v[subscript w, superscript plus symbol]= v[subscript w]/u[subscript tau], were about 8 x 10-3 . The favorable pressure gradients had Pohlhausen parameters between 25 and 42. Skin frictions coefficients around 6 x10-3 and Nusselt numbers under 22 were observed. Near-wall vorticity fluctuations show qualitatively different profiles than observed by Spalart (J. Fluid Mech. 187 (1988)) or Guarini et al. (J. Fluid Mech. 414 (2000)). Small or negative displacement effects are evident. Uncertainty estimates and Favre-averaged equation budgets are provided. A second study aimed to reduce transition-driven uncertainty by determining where on the thermal protection system surface the boundary layer could sustain turbulence. Local boundary layer conditions were extracted from a laminar flow solution over the MPCV which included the bow shock, aerothermochemistry, heat shield surface curvature, and ablation. That information, as a function of leeward distance from the stagnation point, was approximated by Re[subscript theta], Ma[subscript e], [mathematical equation], v[subscript w, superscript plus sign], and T[subscript e]/T[subscript w] along with perfect gas assumptions. Homogenized turbulent boundary layers were initialized at those local conditions and evolved until either stationarity, implying the conditions could sustain turbulence, or relaminarization, implying the conditions could not. Fully turbulent fields relaminarized subject to conditions 4.134 m and 3.199 m leeward of the stagnation point. However, different initial conditions produced long-lived fluctuations at leeward position 2.299 m. Locations more than 1.389 m leeward of the stagnation point are predicted to sustain turbulence in this scenario. / text