My research focuses on reconstructing past climate in southern Arizona using cave deposits called speleothems. However, this necessitates a broader perspective than simply a geochemical time series, and therefore, I also investigate modern cave systems using a combination of modeling and observational datasets. Finally, cave deposits are fundamentally non-renewable resources, and sampling for past climate reconstruction can be destructive, unlike other cave uses. My last investigation is focused on developing possible best practice recommendations for paleoclimate scientists and other cave stakeholders moving forward. We developed two new stalagmite records of past climate variability in southern Arizona over the past 7000 years. Past climate reconstruction from two caves (Cave of the Bells and Fort Huachuca Cave) highlights insolation control of southern Arizona hydroclimate from 7000-2000 years before present. Additionally, comparison between two stalagmites with different seasonal sensitivities uncovers a few eras of multi-decade long droughts in southern Arizona, which align with other regional reconstructions of past climates and elucidate forcings on Southwest paleoclimate as emergent from both external (insolation) and internal climate variability in the Pacific and Atlantic Ocean basins. Although the oxygen isotopic signal of cave calcite in speleothems is complex, agreement with these other records indicates that the speleothem records from these caves primarily record a climate signal.Modeling and monitoring of modern caves both helps us interpret paleoclimate records and enhances our understanding of cave systems in their own right. Modeling of Cave of the Bells dripwaters demonstrates the effect of storage and mixing on the dripwater oxygen isotope signal; non-climate processes can imprint on dripwater variability on multidecadal timescales. Monitoring shows that on very small spatial scales, every cave is different, and even sites within the same cave respond uniquely to surface climate. Most notably, calcite oxygen isotopic composition, used to reconstruct past climate, shows seasonal variability unrelated to dripwater and surface rainfall oxygen isotope variability. Substantial oxygen isotope disequilibrium is identified at numerous caves sites in southern Arizona, and this understanding aligns with a growing number of cave studies that demonstrate the long-held assumption of isotopic equilibrium in cave systems may not always be valid or that the way in which we define isotopic equilibrium insufficiently captures the variety of processes controlling the oxygen isotopic composition of speleothems. Overall, however, monitoring can identify stalagmites that are more sensitive to surface climate and less sensitive to these in-cave processes by identifying sites with dripwater variability responses to surface rainfall variability and sites that precipitate close to oxygen isotopic equilibrium. Finally, a major missing component in speleothem research is the fact that speleothems take thousands and sometimes hundreds of thousands of years to form. They are non-renewable resources on human timescales, and habitat for myriad microbes that have yet to be identified. Removal of speleothems for paleoclimate research is one of the only destructive uses of these deposits. With that in mind, I also analyze current methods of collecting speleothems and develop a framework based on two surveys of scientists and stakeholders to assist scientists and managers when evaluating potential methods of incorporating cave conservation into the speleothem sampling process. Thus, I approach caves from a variety of angles and timescales, from the past through the present to the future, illuminating caves as complex scientific and social systems.

Streams and groundwaters interact in distinctly different ways during flood versus base flow periods. Recent research in the Upper San Pedro River using isotopic and chemical data shows that (1) near-stream, or 'riparian,' groundwater recharged during high streamflow periods is a major contributor to streamflow for the rest of the year, and (2) the amount of riparian groundwater derived from this flood recharge can vary widely (10-90%) along the river. Riparian groundwater in gaining reaches is almost entirely basin groundwater, whereas losing reaches are dominated by prior streamflow.This description of streamflow gives rise to the questions of (1) how much flood recharge occurs at the river-scale, and (2) subsequently, what is the relative importance of flood recharge and basin groundwater in maintaining the hydrologic state of the riparian system. To address these questions, a coupled hydrologic-solute model was constructed for 45 km of the Upper San Pedro riparian system.

A commonly noted problem in the simulation of warm season convection in the North American monsoon region has been the inability of atmospheric models at the meso- scales (10 s to 100 s of kilometers) to simulate organized convection, principally mesoscale convective systems. With the use of convective parameterization, high precipitation biases in model simulations are typically observed over the peaks of mountain ranges. To address this issue, the Kain-Fritsch (KF) cumulus parameterization scheme has been modified with new diagnostic equations to compute the updraft velocity, the convective available potential energy closure assumption, and the convective trigger function. The scheme has been adapted for use in the Weather Research and Forecasting (WRF). A numerical weather prediction-type simulation is conducted for the North American Monsoon Experiment Intensive Observing Period 2 and a regional climate simulation is performed, by dynamically downscaling. In both of these applications, there are notable improvements in the WRF model-simulated precipitation due to the better representation of organized, propagating convection. The use of the modified KF scheme for atmospheric model simulations may provide a more computationally economical alternative to improve the representation of organized convection, as compared to convective-permitting simulations at the kilometer scale or a super-parameterization approach.

The North American monsoon (NAM) is the principal driver of summer severe weather in the Southwest U.S. With sufficient atmospheric instability and moisture, monsoon convection initiates during daytime in the mountains and later may organize, principally into mesoscale convective systems (MCSs). Most monsoon-related severe weather occurs in association with organized convection, including microbursts, dust storms, flash flooding and lightning. The overarching theme of this dissertation research is to investigate simulation of monsoon severe weather due to organized convection within the use of regional atmospheric modeling. A commonly used cumulus parameterization scheme has been modified to better account for dynamic pressure effects, resulting in an improved representation of a simulated MCS during the North American monsoon experiment and the climatology of warm season precipitation in a long-term regional climate model simulation. The effect of urbanization on organized convection occurring in Phoenix is evaluated in model sensitivity experiments using an urban canopy model (UCM) and urban land cover compared to pre-settlement natural desert land cover. The presence of vegetation and irrigation makes Phoenix a "heat sink" in comparison to its surrounding desert, and as a result the modeled precipitation in response to urbanization decreases within the Phoenix urban area and increase on its periphery. Finally, analysis of how monsoon severe weather is changing in association with observed global climate change is considered within the context of a series of retrospectively simulated severe weather events during the period 1948-2010 in a numerical weather prediction paradigm. The individual severe weather events are identified by favorable thermodynamic conditions of instability and atmospheric moisture (precipitable water). Changes in precipitation extremes are evaluated with extreme value statistics. During the last several decades, there has been intensification of organized convective precipitation, but these events occur with less frequency. A more favorable thermodynamic environment for monsoon thunderstorms is the driver of these changes, which is consistent with the broader notion that anthropogenic climate change is presently intensifying weather extremes worldwide.

The North American Monsoon System (NAMS) is a significant weather and climate phenomenon that brings critical rainfall to the southwestern United States and northwestern Mexico. As a result of the North American Monsoon Experiment, and research efforts surrounding the field campaign, the understanding of the NAMS has increased considerably over the last 15 years. In addition questions concerning potential flash flood causing mechanisms of the NAMS have not been thoroughly investigated. This dissertation is comprised of two papers that collectively address the aspects of the literary understanding of the NAMS as we know it today and conduct an investigation into the complex interactions between various weather systems that may influence the NAMS. In the first paper, a review of the major research of the NAMS literature since the last comprehensive review 15 years ago is conducted. The results of his review are assessed for where our understanding has been improved and where future research needs to be guided for purposes of the second paper. Based upon the results from the literature review, the second paper focuses on identification of inverted troughs and gulf surges based upon lower- and mid-level atmospheric parameters for purposes of assessing the impacts on National Weather Service Storm Report flash flood dates. This research contributes to the synthesis of the current knowledge of the NAMS in general and to the specific regional impacts that do occur during periods of heavy precipitation over the NAMS region for purposes of improving meteorological predictability of flash flooding. The results can (1) gauge our understanding of the NAMS literature to date and (2) improve meteorological forecasts through the recognition of synoptic and sub-synoptic patterns related to the NAMS that are most likely to cause flash floods.

This dissertation aims to better understand how various climate modeling approaches affect the fidelity of the North American Monsoon (NAM), as well as the sensitivity of the future state of the NAM under a global warming scenario. Here, we improved over current fully-coupled general circulation models (GCM), which struggle to fully resolve the controlling dynamics responsible for the development and maintenance of the NAM. To accomplish this, we dynamically downscaled a GCM with a regional climate model (RCM). The advantage here being a higher model resolution that improves the representation of processes on scales beyond that which GCMs can resolve. However, as all RCM applications are subject to the transference of biases inherent to the parent GCM, this study developed and evaluated a process to reduce these biases. Pertaining to both precipitation and the various controlling dynamics of the NAM, we found simulations driven by these bias-corrected forcing conditions performed moderately better across a 32-year historical climatology than simulations driven by the original GCM data. Current GCM consensus suggests future tropospheric warming associated with increased radiative forcing as greenhouse gas concentrations increase will suppress the NAM convective environment through greater atmospheric stability. This mechanism yields later onset dates and a generally drier season, but a slight increase to the intensity during July-August. After comparing downscaled simulations forced with original and corrected forcing conditions, we argue that the role of unresolved GCM surface features such as changes to the Gulf of California evaporation lead to a more convective environment. Even when downscaling the original GCM data with known biases, the inclusion of these surface features altered and in some cases reversed GCM trends throughout the southwest United States. This reversal towards a wetter NAM is further magnified in future bias-corrected simulations, which suggest (1) fewer average number of dry days by the end of the 21st century (2) onset occurring up to two to three weeks earlier than the historical average, and (3) more extreme daily precipitation values. However, consistent across each GCM and RCM model is the increase in inter-annual variability, suggesting greater susceptibility to drought conditions in the future.

This dissertation aims to identify the main changes in monsoon activity observed over the American monsoon and Amazonian regions during the last decades and the possible links between such changes. To address this, several observational and reanalysis datasets were used. The results suggest the occurrence of two regime types of the North American monsoon during 1948-2009: two dry regimes during 1948-1959 and 1990-2009 and one wet regime during 1960-1989. The occurrence of such regimes is modulated by the Atlantic Multidecadal Oscillation. However, the two dry regimes have different causes. In particular, the more recent dry regime is mainly due to both an anomalous westward expansion of the North Atlantic Subtropical High and a northward displacement of the subtropical jet stream over the United States. The former enhances the low-level flow from the Gulf of Mexico to the Great Plains and weakens moisture transport to Mexico and the southwestern US. In addition to such a weakening of the North American monsoon during the last two decades, this research shows that the American monsoon systems have shortened after 1978 due to a trend toward earlier retreats of the North American monsoon and delayed onsets of the southern Amazon wet season. These changes produce a longer transition season between both monsoon systems. Whether these changes are caused by a common factor or they are the consequence of independent and unrelated causes was not clear previously. The results discussed here indicate that the observed changes in the American monsoons are partially a consequence of the westward expansion of the North Atlantic surface high observed since 1978. Such a westward expansion enhances the activity of easterly waves over the southern Caribbean Sea and northern South America, producing a dominant easterly flow over the region, which in turn prevents the reversal of the cross-equatorial flow necessary to transport moisture to the southern Amazon and the South American monsoon domain and contributes to its delayed onset. This investigation provides evidence that the shortening and weakening of the American monsoons and the lengthening of the transition season between them are associated with the same large-scale forcing, which may be caused by anthropogenic influence. / text

Western science and Indigenous knowledge understand Southwest climate and the North American monsoon from different cultural perspectives. However, scant literature exists relating to climate and Indigenous communities in the Southwest. On the contrary, substantial climate research has occurred with Arctic Indigenous communities; however, a general aspiration among communities is Indigenous-led climate research and education. This requires more Native scientists and culturally responsive climate science curricula. Southwest Indigenous communities are primed to do this. This dissertation examines 1) the current scientific understanding of the North American monsoon, 2) the state of climate research in Indigenous communities, and 3) the development of culturally responsive climate science curricula. The first paper synthesizes the current scientific understanding of the monsoon and its interannual variability. Pacific Ocean-based teleconnections, such as ENSO-PDO combined indices do add skill in early-season monsoon forecasting. However, general circulation models continue to deal with computational-spatial resolution limitations challenging their application in future climate change projections of the monsoon. The second paper focuses on climate-related research in Indigenous communities in the Arctic and the Southwest to highlight lessons-learned. Climate researchers working with Native communities must exercise cultural considerations for Indigenous relationships with the climate and Indigenous protocols for acquiring and disseminating knowledge. Furthermore, increasing the number of Native students in science and Native scientists are ways to improve climate-related research in Indigenous communities. The third paper is a participatory action research approach to develop a culturally responsive climate science curriculum for Tohono O'odham high school and college students. This project worked with a community advisory board as well as Tohono O'odham Community College instructors and student interns. Pre-assessment surveys were given to community members learn of the most relevant weather and climate topics. The curriculum was developed incorporating local, culturally relevant topics. Climate workshops were offered in the communities using activities developed for the curriculum. Workshop evaluations were positive; however, they also addressed the need for more culturally relevant examples. The overlapping theme for these dissertation papers is cultural understanding for climate research and education in Indigenous communities toward a means for Indigenous-led climate research/education within their own communities.

In the southwestern United States, the North American monsoon (NAM) is the main driver of severe weather and accounts for nearly half the annual precipitation. How the monsoon has behaved in the past and how it will change in the future is a question of major importance for natural resource management and infrastructural planning. In this dissertation, I present the results of three studies that have investigated North American monsoon variability and change from the perspective of paleoclimate records, future climate change projections, and simulation of the low-frequency variability with the longest retrospective atmospheric reanalysis. In the first study, a monsoon-sensitive network of tree-ring chronologies is evaluated within its ability to reproduce NAM variability during the past four centuries. Matrix methods are used to detect the low-frequency spatiotemporal variability. The treering chronologies can reasonable characterizes the dominant modes of NAM climate variability. The monsoon tree-ring network is able to reproduce the interannual variability of cool and warm season precipitation, in a manner similar to the period of the instrumental record. Earlywood and latewood adjusted chronologies reveal low frequency climate variability at decadal and longer timescales that is beyond the ability of the instrumental record to temporally well resolve. This low-frequency climate variability seems to be part of a much larger cycle that coincides with the occurrence of multiyear persistent droughts. In the second study, we consider the modes of natural climate variability identified in the previous study to objectively assess the degree of physical uncertainty in climate change projections for NAM from Regional Climate Models (RCMs) used in the North American Regional Climate Change Assessment Program (NARCCAP). Climate change projection models are evaluated mainly on their ability to represent warm season driven by quasi-stationary Rossby wave trains and El Niño Southern Oscillation – Pacific Decadal Variability (ENSO-PDV). It is concluded that use of the NARCCAP model ensemble mean for NAM climate projections is probably not suitable. NARCCAP RCMs are largely a slave to their driving global models and their error in the specification of large-scale atmospheric circulation. Only one out of eight NARCCAP RCMs has a reasonable representation of the seasonal cycle of monsoon precipitation and ENSOdriven interannual variability in both the 20th and 21st centuries. No decadal variability was observed in any of the NARCCAP RCMs. In the third study, the low-frequency drought signal found with tree-ring chronologies is further explored within the framework of a regional climate modeling. Version 2 of the Twentieth-Century Reanalysis (DD-20CR) is dynamically downscaled over a contiguous U.S.-Mexico domain. Statistic analysis of the DD-20CR suggests that the low-frequency drought signal in the Southwest is driven by atmospheric circulation changes on global to continental scales that affect precipitation in Central American as well. DD-20CR reproduces the spatial patterns of precipitation associated with climate variability at decadal and longer timescales in a manner that compares well with observational records and tree-ring chronologies. Low-frequency climate variability is therefore likely responsible for the multiyear persistent droughts in the last four centuries, as independently evaluated from the tree-ring monsoon-sensitive network.

The aim of the research undertaken in this dissertation was to use medium-range to seasonal precipitation forecasts for hydrologic applications for catchments in the core North American Monsoon (NAM) region. To this end, it was necessary to develop a better understanding of the physical and statistical relationships between runoff processes and the temporal statistics of rainfall. To achieve this goal, development of statistically downscaled estimates of warm season precipitation over the core region of the North American Monsoon Experiment (NAME) were developed. Currently, NAM precipitation is poorly predicted on local and regional scales by Global Circulation Models (GCMs). The downscaling technique used here, the K-Nearest Neighbor (KNN) model, combines information from retrospective GCM forecasts with simultaneous historical observations to infer statistical relationships between the low-resolution GCM fields and the locally-observed precipitation records. The stochastic nature of monsoon rainfall presents significant challenges for downscaling efforts and, therefore, necessitate a regionalization and an ensemble or probabilistic-based approach to quantitative precipitation forecasting. It was found that regionalization of the precipitation climatology prior to downscaling using KNN offered significant advantages in terms of improved skill scores.Selected output variables from retrospective ensemble runs of the National Centers for Environmental Predictions medium-range forecast (MRF) model were fed into the KNN downscaling model. The quality of the downscaled precipitation forecasts was evaluated in terms of a standard suite of ensemble verification metrics. This study represents the first time the KNN model has been successfully applied within a warm season convective climate regime and shown to produce skillful and reliable ensemble forecasts of daily precipitation out to a lead time of four to six days, depending on the forecast month.Knowledge of the behavior of the regional hydrologic systems in NAM was transferred into a modeling framework aimed at improving intra-seasonal hydrologic predictions. To this end, a robust lumped-parameter computational model of intermediate conceptual complexity was calibrated and applied to generate streamflow in three unregulated test basins in the core region of the NAM. The modeled response to different time-accumulated KNN-generated precipitation forcing was investigated. Although the model had some difficulty in accurately simulating hydrologic fluxes on the basis of Hortonian runoff principles only, the preliminary results achieved from this study are encouraging. The primary and most novel finding from this study is an improved predictability of the NAM system using state-of-the-art ensemble forecasting systems. Additionally, this research significantly enhanced the utility of the MRF ensemble forecasts and made them reliable for regional hydrologic applications. Finally, monthly streamflow simulations (from an ensemble-based approach) have been demonstrated. Estimated ensemble forecasts provide quantitative estimates of uncertainty associated with our model forecasts.