Grand Canyon National Park

The argument that the threat posed by introduced species is overblown is often buttressed by the observation that native species sometimes also become invasive. An examination of the literature on plant invasions in the United States shows that six times more nonnative species have been termed invasive than native species, and that a member of the naturalized nonnative pool is 40 times more likely than a native species to be perceived as invasive. In the great majority of instances in which a native plant species is seen as invasive, the invasion is associated with an anthropogenic disturbance, especially changed fire or hydrological regime, livestock grazing, and changes wrought by an introduced species. These results suggest that natives are significantly less likely than nonnatives to be problematic for local ecosystems.

Convergent archeological, geological, palynological, dendrochronological, and radiometric data provide a paleoenvironmental record for the American Southwest at a level of detail and time resolution not previously achieved. Many prehistoric cultural and demographic changes on the Colorado Plateaus coincided with environmental fluctuations defined by precisely dated geoclimatic and bioclimatic indicators. These coincidences support the interpretation that socioeconomic changes and population displacements were commonly triggered by environmental stress.

Abstract Background Forest and nonforest ecosystems of the western United States are experiencing major transformations in response to land-use change, climate warming, and their interactive effects with wildland fire. Some ecosystems are transitioning to persistent alternative types, hereafter called “vegetation type conversion” (VTC). VTC is one of the most pressing management issues in the southwestern US, yet current strategies to intervene and address change often use trial-and-error approaches devised after the fact. To better understand how to manage VTC, we gathered managers, scientists, and practitioners from across the southwestern US to collect their experiences with VTC challenges, management responses, and outcomes. Results Participants in two workshops provided 11 descriptive case studies and 61 examples of VTC from their own field observations. These experiences demonstrate the extent and complexity of ecological reorganization across the region. High-severity fire was the predominant driver of VTC in semi-arid coniferous forests. By a large margin, these forests converted to shrubland, with fewer conversions to native or non-native herbaceous communities. Chaparral and sagebrush areas nearly always converted to non-native grasses through interactions among land use, climate, and fire. Management interventions in VTC areas most often attempted to reverse changes, although we found that these efforts cover only a small portion of high-severity burn areas undergoing VTC. Some areas incurred long (>10 years) observational periods prior to initiating interventions. Efforts to facilitate VTC were rare, but could cover large spatial areas. Conclusions Our findings underscore that type conversion is a common outcome of high-severity wildland fire in the southwestern US. Ecosystem managers are frontline observers of these far-reaching and potentially persistent changes, making their experiences valuable in further developing intervention strategies and research agendas. As its drivers increase with climate change, VTC appears increasingly likely in many ecological contexts and may require management paradigms to transition as well. Approaches to VTC potentially include developing new models of desired conditions, the use of experimentation by managers, and broader implementation of adaptive management strategies. Continuing to support and develop science-manager partnerships and peer learning groups will help to shape our response to ongoing rapid ecological transformations.

Artificial nightlight is increasingly recognized as an important environmental disturbance that influences the habitats and fitness of numerous species. However, its effects on wide‐ranging vertebrates and their interactions remain unclear. Light pollution has the potential to amplify land‐use change, and as such, answering the question of how this sensory stimulant affects behavior and habitat use of species valued for their ecological roles and economic impacts is critical for conservation and land‐use planning. Here, we combined satellite‐derived estimates of light pollution, with GPS‐data from cougars Puma concolor (n = 56), mule deer Odocoileus hemionus (n = 263) and locations of cougar‐killed deer (n = 1562 carcasses), to assess the effects of light exposure on mammal behavior and predator–prey relationships across wildland–urban gradients in the southwestern United States. Our results indicate that deer used the anthropogenic environments to access forage and were more active at night than their wildland conspecifics. Despite higher nightlight levels, cougars killed deer at the wildland–urban interface, but hunted them in the relatively darkest locations. Light had the greatest effect of all covariates on where cougars killed deer at the wildland–urban interface. Both species exhibited functional responses to light pollution at fine scales; individual cougars and deer with less light exposure increasingly avoided illuminated areas when exposed to greater radiance, whereas deer living in the wildland–urban interface selected elevated light levels. We conclude that integrating estimates of light pollution into ecological studies provides crucial insights into how the dynamic human footprint can alter animal behavior and ecosystem function across spatial scales.

1. Biological invasions threaten biodiversity, and understanding the factors that influence a community’s susceptibility to invasion informs both management of invasive species and conservation efforts towards promoting biodiversity. 2. In this study, we examined the native–exotic richness relationships (NERR) at two spatial scales and asked what variables mediate the relationship between native and exotic species richness in 1-m2 plots among old fields. We also tested whether the favourable environment and spatial heterogeneity mechanisms alter the relationship between native and exotic richness. Additionally, we identified biotic, abiotic and landscape factors that accounted for patterns of exotic species richness across spatial scales, from 50-m transects (i.e. local scales) to entire old fields (i.e. landscape scales). 3. We found that native and exotic richness were positively related at both local and landscape spatial scales, with the strength of their relationship greater at the landscape scale. In old-field communities with lower foliar cover, native and exotic plant richness were negatively related across 1-m2 plots, whereas in old fields with higher foliar cover, native and exotic plant richness were positively related across 1-m2 plots. Overall, predictors of exotic species richness depended on spatial scale: at local scales, native plant richness, light and landscape factors accounted for most of the variation, but biotic factors alone accounted for most of the variation at the landscape scale. 4. Our findings suggest that fields favourable for native species are also suitable for exotic species, and management efforts towards exotic richness should be tailored to spatial scale.

We studied the effects of water temperature, density as measured by density index (DI), and food deprivation (alternating cycles of daily feeding and fasting) on dorsal fin erosion of steelhead Oncorhynchus mykiss through a factorial experiment lasting 22 weeks. Dorsal fin index (DFI) correlated inversely with age, length, and weight of fish and was significantly affected by crowding (i.e., DI). Increasing DI (from 0.0 to 0.5) also improved feed conversions and slightly depressed survival, but DI had no significant effects on weight gain, condition (K), or carcass composition. Water temperature (10°C or 15°C) had a strong and significant effect on DFI; colder water improved DFI while slowing body growth. Low temperature treatments improved survival slightly and altered carcass composition by increasing moisture at the expense of fat and protein. The effects of density and temperature were also strongly interactive. Food deprivation slightly improved final DFI at low DIs, but greatly reduced the average size of fish affected by severe fin erosion. Food-deprived fish exhibited carcass composition trends similar to fish reared in 10°C water. Feed conversion improved significantly with alternating fasting and feeding cycles. Histological evaluation proved inconclusive for physical fin nipping evidence or for the presence of detrimental microorganisms.

Most studies assessing vegetation response following control of invasive Tamarix trees along southwestern U.S. rivers have been small in scale (e.g., river reach), or at a regional scale but with poor spatial-temporal replication, and most have not included testing the effects of a now widely used biological control. We monitored plant composition following Tamarix control along hydrologic, soil, and climatic gradients in 244 treated and 172 reference sites across six U.S. states. This represents the largest comprehensive assessment to date on the vegetation response to the four most common Tamarix control treatments. Biocontrol by a defoliating beetle (treatment 1) reduced the abundance of Tamarix less than active removal by mechanically using hand and chain-saws (2), heavy machinery (3) or burning (4). Tamarix abundance also decreased with lower temperatures, higher precipitation, and follow-up treatments for Tamarix resprouting. Native cover generally increased over time in active Tamarix removal sites, however, the increases observed were small and was not consistently increased by active revegetation. Overall, native cover was correlated to permanent stream flow, lower grazing pressure, lower soil salinity and temperatures, and higher precipitation. Species diversity also increased where Tamarix was removed. However, Tamarix treatments, especially those generating the highest disturbance (burning and heavy machinery), also often promoted secondary invasions of exotic forbs. The abundance of hydrophytic species was much lower in treated than in reference sites, suggesting that management of southwestern U.S. rivers has focused too much on weed control, overlooking restoration of fluvial processes that provide habitat for hydrophytic and floodplain vegetation. These results can help inform future management of Tamarix-infested rivers to restore hydrogeomorphic processes, increase native biodiversity and reduce abundance of noxious species.

Symptomatic hyponatremia has been reported in athletes participating in marathon and ultra-marathon events. We report four cases of symptomatic hyponatremia in hikers in the Grand Canyon. These are the first cases reported in recreational hikers in a wilderness setting. The overall fluid status of these patients is not well defined. In most instances this probably represents dilutional hyponatremia from sweat loss replaced with plain water. Recreational wilderness participants may require electrolyte replacement similar to endurance athletes.

Extreme climate events are transforming plant communities in the desert Southwest of the United States. Abundant precipitation in 1998 associated with El Niño Southern Oscillation (ENSO) stimulated exceptional alien annual plant production in the Mojave Desert that fueled wildfires in 1999. Exacerbated by protracted drought, 80% of the burned Yucca brevifolia, a long-lived arborescent monocot, and 26% of unburned plants died at Joshua Tree National Park by 2004. Many burned plants <1 m tall died immediately, and survival of all but the tallest, oldest plants declined to the same low level by 2004. Postfire sprouting prolonged survival, but only at the wetter, high-elevation sites. During succeeding dry years, herbaceous plants were scarce, and individuals of Thomomys bottae (pocket gopher) gnawed the periderm and hollowed stems of Y. brevifolia causing many of them to topple. Thomomys bottae damage reduced plant survivorship at low-elevation, unburned sites and diminished survival of burned plants in all but the driest site, which already had low survival. Accentuated ENSO episodes and more frequent wildfires are expected for the desert Southwest and will likely shift Y. brevifolia population structure toward tall, old adults with fewer opportunities for plant recruitment, thus imperiling the persistence of this unique plant community.

Abstract Gully erosion of cultural sites in Grand Canyon National Park is an urgent management problem that has intensified in recent decades, potentially related to the effects of Glen Canyon Dam. We studied 25 gullies at nine sites in Grand Canyon over the 2002 monsoon–erosion season to better understand the geomorphology of the gully erosion and the effectiveness of erosion‐control structures (ECS) installed by the park under the direction of the Zuni Conservation Program. Field results indicate that Hortonian overland flow leads to concentrated flow in gullies and erosion focused at knickpoints along channels as well as at gully heads. Though groundcover type, soil shear strength and permeability vary systemat‐ically across catchments, gradient and, to a lesser degree, contributing drainage area seem to be the first‐order controls on gully extent, location of new knickpoints, and ECS damage. The installed ECS do reduce erosion relative to reaches without them and initial data suggest woody checkdams are preferable to rock linings, but maintenance is essential because damaged structures can exacerbate erosion. Topographic data from intensive field surveys and detailed photogrammetry provide slope–contributing area data for gully heads that have a trend consistent with previous empirical and theoretical formulations from a variety of landscapes. The same scaling holds below gully heads for knickpoint and ECS topographic data, with threshold coefficients the lowest for gully heads, slightly higher for knickpoints, and notably higher for damaged ECS. These topographic thresholds were used with 10‐cm digital elevation models to create simple predictive models for gully extent and structure damage. The model predictions accounted for the observed gullies but there are also many false‐positives. Purely topographical models are probably inadequate at this scale and application, but models that also parameterize the variable soil properties across sites would be useful for predicting erosion problems and ECS failure. Copyright © 2006 John Wiley &amp; Sons, Ltd.

We reviewed all North American radiotelemetry studies of within-home-range habitat selection by northern goshawks (Accipiter gentilis) and all studies relating territory occupancy and productivity to landscape habitat features. Goshawks selected habitats in the home range with structural characteristics of mature to old-growth forests, such as large trees and high canopy closure. We documented selection of these structures in a number of forest types, including ponderosa pine (Pinus ponderosa), mixed conifer, lodgepole pine (P. contorta), temperate rainforest, upland hardwood, and boreal forest, suggesting that goshawks are selecting forests for their structure rather than for species composition. Goshawks did not select stands with the greatest prey abundance. Selection for natural openings, edges, and stand diversity was inconclusive. Habitat selection patterns suggest that current goshawk management plans in the western United States may be inadequate.

that occurs during or up to 24 h after prolonged physical activity. It is reported to occur in individual physical activities or during organized endurance events conducted in environments in which medical care is limited and often not available, and patient evacuation to definitive care is often greatly delayed. Rapid recognition and appropriate treatment are essential in the severe form to increase the likelihood of a positive outcome. To mitigate the risk of EAH mismanagement, care providers in the prehospital and in hospital settings must differentiate from other causes that present with similar signs and symptoms. EAH most commonly has overlapping signs and symptoms with heat exhaustion and exertional heat stroke. Failure in this regard is a recognized cause of worsened morbidity and mortality. In an effort to produce best practice guidelines for EAH management, the Wilderness Medical Society convened an expert panel in May 2018. The panel was charged with updating the WMS Practice Guidelines for Treatment of Exercise-Associated Hyponatremia published in 2014 using evidence-based guidelines for the prevention, recognition, and treatment of EAH. Recommendations are made based on presenting with symptomatic EAH, particularly when point-of-care blood sodium testing is unavailable in the field. These recommendations are graded on the basis of the quality of supporting evidence and balanced between the benefits and risks/burdens for each parameter according to the methodology stipulated by the American College of Chest Physicians.

Research Article| January 01, 1996 Tributary debris fans and the late Holocene alluvial chronology of the Colorado River, eastern Grand Canyon, Arizona Richard Hereford; Richard Hereford 1U.S. Geological Survey, 2255 N. Gemini Drive, Flagstaff, Arizona 86001 Search for other works by this author on: GSW Google Scholar Kathryn S. Thompson; Kathryn S. Thompson 1U.S. Geological Survey, 2255 N. Gemini Drive, Flagstaff, Arizona 86001 Search for other works by this author on: GSW Google Scholar Kelly J. Burke; Kelly J. Burke 1U.S. Geological Survey, 2255 N. Gemini Drive, Flagstaff, Arizona 86001 Search for other works by this author on: GSW Google Scholar Helen C. Fairley Helen C. Fairley 2National Park Service, Grand Canyon National Park, P.O. Box 129, Grand Canyon, Arizona 86023-0129 Search for other works by this author on: GSW Google Scholar Author and Article Information Richard Hereford 1U.S. Geological Survey, 2255 N. Gemini Drive, Flagstaff, Arizona 86001 Kathryn S. Thompson 1U.S. Geological Survey, 2255 N. Gemini Drive, Flagstaff, Arizona 86001 Kelly J. Burke 1U.S. Geological Survey, 2255 N. Gemini Drive, Flagstaff, Arizona 86001 Helen C. Fairley 2National Park Service, Grand Canyon National Park, P.O. Box 129, Grand Canyon, Arizona 86023-0129 Publisher: Geological Society of America First Online: 01 Jun 2017 Online ISSN: 1943-2674 Print ISSN: 0016-7606 Geological Society of America GSA Bulletin (1996) 108 (1): 3–19. https://doi.org/10.1130/0016-7606(1996)108<0003:TDFATL>2.3.CO;2 Article history First Online: 01 Jun 2017 Cite View This Citation Add to Citation Manager Share Icon Share Facebook Twitter LinkedIn Email Permissions Search Site Citation Richard Hereford, Kathryn S. Thompson, Kelly J. Burke, Helen C. Fairley; Tributary debris fans and the late Holocene alluvial chronology of the Colorado River, eastern Grand Canyon, Arizona. GSA Bulletin 1996;; 108 (1): 3–19. doi: https://doi.org/10.1130/0016-7606(1996)108<0003:TDFATL>2.3.CO;2 Download citation file: Ris (Zotero) Refmanager EasyBib Bookends Mendeley Papers EndNote RefWorks BibTex toolbar search Search Dropdown Menu toolbar search search input Search input auto suggest filter your search All ContentBy SocietyGSA Bulletin Search Advanced Search Abstract Bouldery debris fans and sandy alluvial terraces of the Colorado River developed contemporaneously during the late Holocene at the mouths of nine major tributaries in eastern Grand Canyon. The age of the debris fans and alluvial terraces contributes to understanding river hydraulics and to the history of human activity along the river, which has been concentrated on these surfaces for at least two to three millennia. Poorly sorted, coarse-grained debris-flow deposits of several ages are interbedded with, overlie, or are overlapped by three terrace-forming alluviums. The alluvial deposits are of three age groups: the striped alluvium, deposited from before 770 b.c. to about a.d. 300; the alluvium of Pueblo II age deposited from about a.d. 700 to December 1900; and the alluvium of the upper mesquite terrace, deposited from about a.d. 1400 to 1880. Two elements define the geomorphology of a typical debris fan: the large, inactive surface of the fan and a smaller, entrenched, active debris-flow channel and fan that is about one-sixth the area of the inactive fan. The inactive fan is segmented into at least three surfaces with distinctive weathering characteristics. These surfaces are conformable with underlying debris-flow deposits that date from before 770 b.c. to around a.d. 660, a.d. 660 to before a.d. 1200, and from a.d. 1200 to slightly before 1890, respectively, based on late-19th-century photographs, radiocarbon and archaeologic dating of the three stratigraphically related alluviums, and radiocarbon dating of fine-grained debris-flow deposits. These debris flows aggraded the fans in at least three stages beginning about 2.8 ka, if not earlier in the late Holocene. Several main-stem floods eroded the margin of the segmented fans, reducing fan symmetry. The entrenched, active debris-flow channels contain deposits <100 yr old, which form debris fans at the mouth of the channel adjacent to the river. Early and middle Holocene debris-flow and alluvial deposits have not been recognized, as they were evidently not preserved adjacent to the river or are buried by younger deposits. This content is PDF only. Please click on the PDF icon to access. First Page Preview Close Modal You do not have access to this content, please speak to your institutional administrator if you feel you should have access.

Two ends of the fire regime spectrum are a frequent low-intensity fire regime and an infrequent high-intensity fire regime, but intermediate fire regimes combine high- and low-severity fire over space and time. We used fire-scar and tree-age data to reconstruct fire regime attributes of mixed-conifer and aspen forests in the North Rim area of Grand Canyon National Park, with a goal of estimating patch sizes of historical high-severity fire and comparing them with modern patch sizes. We used three methods based on (i) aspen groves, (ii) even-aged stands, and (iii) inverse distance weighting, to estimate occurrence and patch sizes of historical high-severity fire. Evidence of high-severity fire was common in the 1800s, and high-severity fire years were associated with drought. High-severity fire patch sizes likely ranged from 10 −1 to 10 2 ha. However, the forest is quite young, and we cannot rule out a historical large high-severity fire that could have reinitiated much of the 1400 ha study area. Fire scars, which are indicative of low-severity fire, were also common. Historical fire was likely heterogeneous across the landscape. Maintaining heterogeneity of fire severity, size, and frequency would promote heterogeneity of forest structure and composition and resilience to future disturbances.

Recovery of imperiled fishes can be achieved through suppression of invasives, but outcomes may vary with environmental conditions. We studied the response of imperiled desert fishes to an invasive brown (Salmo trutta) and rainbow trout (Oncorhynchus mykiss) suppression program in a Colorado River tributary, with natural flow and longitudinal variation in thermal characteristics. We investigated trends in fish populations related to suppression and tested hypotheses about the impacts of salmonid densities, hydrologic variation, and spatial–thermal gradients on the distribution and abundance of native fish species using zero-inflated generalized linear mixed effects models. Between 2012 and 2018, salmonids declined 89%, and native fishes increased dramatically (∼480%) once trout suppression surpassed ∼60%. Temperature and trout density were consistently retained in the top models predicting the abundance and distribution of native fishes. The greatest increases occurred in warmer reaches and in years with spring flooding. Surprisingly, given the evolution of native fishes in disturbance-prone systems, intense, monsoon-driven flooding limited native fish recruitment. Applied concertedly, invasive species suppression and efforts to mimic natural flow and thermal regimes may allow rapid and widespread native fish recovery.

Landscape genetic studies based on neutral genetic markers have contributed to our understanding of the influence of landscape composition and configuration on gene flow and genetic variation. However, the potential for species to adapt to changing landscapes will depend on how natural selection influences adaptive genetic variation. We demonstrate how landscape resistance models can be combined with genetic simulations incorporating natural selection to explore how the spread of adaptive variation is affected by landscape characteristics, using desert bighorn sheep (Ovis canadensis nelsoni) in three differing regions of the southwestern United States as an example. We conducted genetic sampling and least-cost path modeling to optimize landscape resistance models independently for each region, and then simulated the spread of an adaptive allele favored by selection across each region. Optimized landscape resistance models differed between regions with respect to landscape variables included and their relationships to resistance, but the slope of terrain and the presence of water barriers and major roads had the greatest impacts on gene flow. Genetic simulations showed that differences among landscapes strongly influenced spread of adaptive genetic variation, with faster spread (1) in landscapes with more continuously distributed habitat and (2) when a pre-existing allele (i.e., standing genetic variation) rather than a novel allele (i.e., mutation) served as the source of adaptive genetic variation. The combination of landscape resistance models and genetic simulations has broad conservation applications and can facilitate comparisons of adaptive potential within and between landscapes.

Topographic measures are frequently used in a variety of landscape ecology applications, in their simplest form as elevation, slope, and aspect, but increasingly more complex measures are being employed. We explore terrain metric similarity with changes in scale, both grain and extent, and examine how selecting the best measures is sensitive to changes in application. There are three types of topographic measures: 1) those that relate to orientation for approximating solar input, 2) those that capture variability in terrain configuration, and 3) those that provide metrics about landform features. Many biodiversity hotspots and predators have been found to coincide with areas of complexity, yet most complexity measures cannot differentiate between terrain steepness and uneven and broken terrain. Currently characterizing terrain in landscape-level analyses can be challenging, especially at coarser spatial resolutions but developing methods that improve landscape-level assessments include multivariate approaches and the use of neighborhood statistics. Some measures are sensitive to the spatial grain of calculation, the physiography of the landscape, and the scale of application. We show which measures have the potential to be multi-collinear, and illustrate with a case study how the selection of the best measures can change depending on the question at hand using mountain lion (Puma concolor) occurrence data. The case study showed a combination of infrequently employed metrics, such as view-shed analysis and focal statistics, outperform more commonly employed singular metrics. The use of focal statistics as a measure of topographic complexity shows promise for improving how mountain lions use terrain features.

Abstract Restoration treatments were tested on the South Rim (ponderosa pine) and North Rim (mixed conifer) of Grand Canyon National Park, never‐harvested forests where the historically frequent surface regime was interrupted in the late 19th century. Treatments were designed for “minimal impact” by limiting the size of trees to be thinned and minimizing mechanical equipment. Treatments included (1) thinning of small trees (diameter ≤ 12.7 cm) and prescribed fire; (2) thinning of small trees located close to large old trees and prescribed fire; (3) prescribed fire only; and (4) control. On the South Rim, density declined by 23–45% but basal area only by 9–14%. On the North Rim, density was reduced 32–68% and basal area declined 18–31%. Declines of 5–8% were observed in these same variables in the controls. Surface fuels were significantly reduced in all burned treatments (70–84% reduction in forest floor, 66–76% reduction in woody debris). However, canopy cover was nearly unchanged, and canopy fuel loads were not significantly reduced, although canopy base heights increased. The experiment accomplished the goal of minimally impacting the forest ecosystem but the effects of small‐tree thinning plus burning were nearly indistinguishable from those of burning alone. All the treated units were left with modest gain in resistance to severe wildfires but in conditions still far removed from prefire‐exclusion conditions.

The Trail of Time exhibition under construction at Grand Canyon National Park is the world's largest geoscience exhibition at one of the world's grandest geologic landscapes. It is a 2-km-long interpretive walking timeline trail that leverages Grand Canyon vistas and rocks to guide visitors to ponder, explore, and understand the magnitude of geologic time and the stories encoded by Grand Canyon rock layers and landscapes. As one of a new generation of geoscience education exhibits, the Trail of Time targets multiple cognitive and affective levels with accurate content, active geoscience inquiry and interpretation, and place-based cultural integration. It developed as an outgrowth of sustained geoscience research funded by the National Science Foundation, with scientists as the conceivers and coordinators of the project. It benefits from a high level of synergy with the National Park Service interpretation division, as well as extensive on-site and off-site evaluation of pedagogic effectiveness in the outdoor informal science environment. The Trail of Time will impact many of the five million annual visitors to the National Park. Associated cognitive research on public understanding of “deep time” offers opportunities to inform more effective geoscience pedagogy for informal and formal educational settings.

Abstract The Kaibab Plateau and Grand Canyon National Park in the USA contain both shallow and deep karst systems, which interact in ways that are not well known, although recent studies have allowed better interpretations of this unique system. Detailed characterization of sinkholes and their distribution on the surface using geographical information system and LiDAR data can be used to relate the infiltration points to the overall hydrogeological system. Flow paths through the deep regional geological structure were delineated using non-toxic fluorescent dyes. The flow characteristics of the coupled aquifer system were evaluated using hydrograph recession curve analysis via discharge data from Roaring Springs, the sole source of the water supply for the Grand Canyon National Park. The interactions between these coupled surface and deep karst systems are complex and challenging to understand. Although the surface karst behaves in much the same way as karst in other similar regions, the deep karst has a base flow recession coefficient an order of magnitude lower than many other karst aquifers throughout the world. Dye trace analysis reveals rapid, conduit-dominated flow that demonstrates fracture connectivity along faults between the surface and deep karst. An understanding of this coupled karst system will better inform aquifer management and research in other complex karst systems.