Pub Date : 2021-11-01DOI: 10.2113/gseegeosci.27.4.375
Dennis Staley, J. Lancaster, Alan J. Gallegos, T. Wasklewicz
{"title":"Introduction to Special Issue on Slope Stability in Memory of Jerome (Jerry) De Graff: Part 1","authors":"Dennis Staley, J. Lancaster, Alan J. Gallegos, T. Wasklewicz","doi":"10.2113/gseegeosci.27.4.375","DOIUrl":"https://doi.org/10.2113/gseegeosci.27.4.375","url":null,"abstract":"","PeriodicalId":138906,"journal":{"name":"Environmental and Engineering Geoscience","volume":"98 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2021-11-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"134276394","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Landslide runout has traditionally been quantified by the height-to-length ratio, H/L, which, in many cases, is strongly influenced by the slope of the runout path. In this study, we propose an alternative mobility measure, the unitless Runout Number, measured as the landslide length divided by the square root of the landslide area, which characterizes landslide shape in terms of elongation. We used a database of 158 landslides of varying runout distances from locations in northern California, Oregon, and Washington state to compare the two runout measurement methods and explore their predictability using parameters that can be measured or estimated using geographic information systems. The Runout Number better describes the overall runout for several landslide and slope geometries. The two mobility measures show very little correlation to each other, indicating that the two parameters describe different landslide mobility mechanisms. When compared to predictive parameters shown by prior research to relate to landslide runout, the two runout measurement methods show different correlations. H/L correlates more strongly to initial slope angle, upslope contributing area, landslide area, and grain size distribution (percent clay, silt, total fines, and sand). The Runout Number correlates more strongly to planimetric curvature, upslope contributing area normalized by landslide area, and percent sand. Although these correlations are not necessarily strong enough for prediction, they indicate the validity of both runout measurement methods and the benefit of including both numbers when characterizing landslide mobility.
{"title":"Runout Number: A New Metric for Landslide Runout Characterization","authors":"Cory S. Wallace, P. Santi","doi":"10.2113/eeg-d-20-00144","DOIUrl":"https://doi.org/10.2113/eeg-d-20-00144","url":null,"abstract":"\u0000 Landslide runout has traditionally been quantified by the height-to-length ratio, H/L, which, in many cases, is strongly influenced by the slope of the runout path. In this study, we propose an alternative mobility measure, the unitless Runout Number, measured as the landslide length divided by the square root of the landslide area, which characterizes landslide shape in terms of elongation. We used a database of 158 landslides of varying runout distances from locations in northern California, Oregon, and Washington state to compare the two runout measurement methods and explore their predictability using parameters that can be measured or estimated using geographic information systems. The Runout Number better describes the overall runout for several landslide and slope geometries. The two mobility measures show very little correlation to each other, indicating that the two parameters describe different landslide mobility mechanisms. When compared to predictive parameters shown by prior research to relate to landslide runout, the two runout measurement methods show different correlations. H/L correlates more strongly to initial slope angle, upslope contributing area, landslide area, and grain size distribution (percent clay, silt, total fines, and sand). The Runout Number correlates more strongly to planimetric curvature, upslope contributing area normalized by landslide area, and percent sand. Although these correlations are not necessarily strong enough for prediction, they indicate the validity of both runout measurement methods and the benefit of including both numbers when characterizing landslide mobility.","PeriodicalId":138906,"journal":{"name":"Environmental and Engineering Geoscience","volume":"18 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2021-11-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"133716832","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
C. Pluhar, K. R. Ford, G. Stock, J. Stone, S. Zimmerman
Yosemite National Park, California, is one of the best-documented sites of historical rockfalls and other rock slope failures; however, past work shows that this record does not capture the infrequent largest occurrences, prehistoric events orders of magnitude larger than the largest historic ones. These large prehistoric events are evident as voluminous bouldery landslide deposits, permitting volume and age quantification to better understand local volume–frequency relationships, potential triggering mechanisms, and the hazard such events might pose. The Tiltill rockslide in northern Yosemite is one such example, consisting of 2.1 × 106 m3 ± 1.6 × 106 m3 of talus (1.5 × 106 m3 original volume of rock mass) that slid across the floor of Tiltill Valley, partially damming Tiltill Creek to create a seasonal pond that drains through and around the rockslide mass. This volume and the rockslide's effective coefficient of friction, 0.47, place it near the boundary between long-runout landslides and ordinary Coulomb failure. Although the rockslide superficially appears to consist of two separate lobes, statistically indistinguishable 10Be exposure dates from eight samples indicate a single event that occurred at 13.0 ± 0.8 ka. The age of the Tiltill rockslide and its relatively low elevation compared to equilibrium line altitudes at this place and time make glacial debutressing a highly unlikely triggering mechanism. Seismic shaking associated with fault rupture along the eastern Sierra Nevada is shown to be a plausible but unverified trigger.
{"title":"Magnitude and Timing of the Tiltill Rockslide in Yosemite National Park, California","authors":"C. Pluhar, K. R. Ford, G. Stock, J. Stone, S. Zimmerman","doi":"10.2113/eeg-d-21-00033","DOIUrl":"https://doi.org/10.2113/eeg-d-21-00033","url":null,"abstract":"\u0000 Yosemite National Park, California, is one of the best-documented sites of historical rockfalls and other rock slope failures; however, past work shows that this record does not capture the infrequent largest occurrences, prehistoric events orders of magnitude larger than the largest historic ones. These large prehistoric events are evident as voluminous bouldery landslide deposits, permitting volume and age quantification to better understand local volume–frequency relationships, potential triggering mechanisms, and the hazard such events might pose. The Tiltill rockslide in northern Yosemite is one such example, consisting of 2.1 × 106 m3 ± 1.6 × 106 m3 of talus (1.5 × 106 m3 original volume of rock mass) that slid across the floor of Tiltill Valley, partially damming Tiltill Creek to create a seasonal pond that drains through and around the rockslide mass. This volume and the rockslide's effective coefficient of friction, 0.47, place it near the boundary between long-runout landslides and ordinary Coulomb failure. Although the rockslide superficially appears to consist of two separate lobes, statistically indistinguishable 10Be exposure dates from eight samples indicate a single event that occurred at 13.0 ± 0.8 ka. The age of the Tiltill rockslide and its relatively low elevation compared to equilibrium line altitudes at this place and time make glacial debutressing a highly unlikely triggering mechanism. Seismic shaking associated with fault rupture along the eastern Sierra Nevada is shown to be a plausible but unverified trigger.","PeriodicalId":138906,"journal":{"name":"Environmental and Engineering Geoscience","volume":"1 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2021-11-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"130999393","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2021-08-01DOI: 10.2113/gseegeosci.27.3.373
I. Pope
{"title":"A Hero On Mount St. Helens: The Life and Legacy of David A. Johnston","authors":"I. Pope","doi":"10.2113/gseegeosci.27.3.373","DOIUrl":"https://doi.org/10.2113/gseegeosci.27.3.373","url":null,"abstract":"","PeriodicalId":138906,"journal":{"name":"Environmental and Engineering Geoscience","volume":"42 4 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2021-08-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"123496413","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2021-08-01DOI: 10.2113/gseegeosci.27.3.371
R. V. Arsdale
{"title":"Timefulness: How Thinking Like a Geologist Can Help Save the World","authors":"R. V. Arsdale","doi":"10.2113/gseegeosci.27.3.371","DOIUrl":"https://doi.org/10.2113/gseegeosci.27.3.371","url":null,"abstract":"","PeriodicalId":138906,"journal":{"name":"Environmental and Engineering Geoscience","volume":"28 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2021-08-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"116051264","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Substantial advances have been achieved in various aspects of debris-flow hazard assessments over the past decade. These advances include sophisticated ways to date previous events, two- and three-dimensional runout models including multi-phase flows and debris entrainment options, and applications of extreme value statistics to assemble frequency–magnitude analyses. Pertinent questions have remained the same: How often, how big, how fast, how deep, how intense, and how far? Similarly, although major life loss attributable to debris flows can often, but not always, be avoided in developed nations, debris flows remain one of the principal geophysical killers in mountainous terrains. Substantial differences in debris-flow hazard persist between nations. Some rely on a design magnitude associated with a specific return period; others use relationships between intensity and frequency; and some allow for, but do not mandate, in-depth quantitative risk assessments. Differences exist in the management of debris-flow risks, from highly sophisticated and nation-wide applied protocols to retroaction in which catastrophic debris flows occur before they are considered for mitigation. Two factors conspire to challenge future generations of debris-flow researchers, practitioners, and decision makers: Population growth and climate change, which are increasingly manifested by augmenting hydroclimatic extremes. While researchers will undoubtedly finesse future remote sensing, dating, and runout techniques and models, practitioners will need to focus on translating those advances into practical cost-efficient tools and integrating those tools into long-term debris-flow risk management.
{"title":"Debris-Flow Hazard Assessments: A Practitioner's View","authors":"M. Jakob","doi":"10.2113/EEG-D-20-00110","DOIUrl":"https://doi.org/10.2113/EEG-D-20-00110","url":null,"abstract":"\u0000 Substantial advances have been achieved in various aspects of debris-flow hazard assessments over the past decade. These advances include sophisticated ways to date previous events, two- and three-dimensional runout models including multi-phase flows and debris entrainment options, and applications of extreme value statistics to assemble frequency–magnitude analyses. Pertinent questions have remained the same: How often, how big, how fast, how deep, how intense, and how far? Similarly, although major life loss attributable to debris flows can often, but not always, be avoided in developed nations, debris flows remain one of the principal geophysical killers in mountainous terrains. Substantial differences in debris-flow hazard persist between nations. Some rely on a design magnitude associated with a specific return period; others use relationships between intensity and frequency; and some allow for, but do not mandate, in-depth quantitative risk assessments. Differences exist in the management of debris-flow risks, from highly sophisticated and nation-wide applied protocols to retroaction in which catastrophic debris flows occur before they are considered for mitigation. Two factors conspire to challenge future generations of debris-flow researchers, practitioners, and decision makers: Population growth and climate change, which are increasingly manifested by augmenting hydroclimatic extremes. While researchers will undoubtedly finesse future remote sensing, dating, and runout techniques and models, practitioners will need to focus on translating those advances into practical cost-efficient tools and integrating those tools into long-term debris-flow risk management.","PeriodicalId":138906,"journal":{"name":"Environmental and Engineering Geoscience","volume":"60 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2021-03-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"116583368","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
The objective of this study was to evaluate the factors that contribute to the high frequency of landslides in the Kope Formation and the overlying colluvial soil present in the Cincinnati area, southwestern Ohio. The Kope Formation consists of approximately 80 percent shale inter-bedded with 20 percent limestone. The colluvium that forms from the weathering of the shale bedrock consists of a low-plasticity clay. Based on field observations, LiDAR data, and information gathered from city and county agencies, we created a landslide inventory map for the Cincinnati area, identifying 842 landslides. From the inventory map, we selected 10 landslides that included seven rotational and three translational slides for detailed investigations. Representative samples were collected from the landslide sites for determining natural water content, Atterberg limits, grain size distribution, shear strength parameters, and slake durability index. For the translational landslides, strength parameters were determined along the contact between the bedrock and the overlying colluvium. The results of the study indicate that multiple factors contribute to landslide susceptibility of the Kope Formation and the overlying colluvium, including low shear strength of the colluvial soil, development of porewater pressure within the slope, human activity such as loading the top or cutting the toe of a slope, low to very low durability of the bedrock that allows rapid disintegration of the bedrock and accumulation of colluvial soil, undercutting of the slope toe by stream water, and steepness of the slopes.
{"title":"Factors Contributing to Landslide Susceptibility of the Kope Formation, Cincinnati, Ohio","authors":"M. Glassmeyer, A. Shakoor","doi":"10.2113/EEG-D-20-00077","DOIUrl":"https://doi.org/10.2113/EEG-D-20-00077","url":null,"abstract":"\u0000 The objective of this study was to evaluate the factors that contribute to the high frequency of landslides in the Kope Formation and the overlying colluvial soil present in the Cincinnati area, southwestern Ohio. The Kope Formation consists of approximately 80 percent shale inter-bedded with 20 percent limestone. The colluvium that forms from the weathering of the shale bedrock consists of a low-plasticity clay. Based on field observations, LiDAR data, and information gathered from city and county agencies, we created a landslide inventory map for the Cincinnati area, identifying 842 landslides. From the inventory map, we selected 10 landslides that included seven rotational and three translational slides for detailed investigations. Representative samples were collected from the landslide sites for determining natural water content, Atterberg limits, grain size distribution, shear strength parameters, and slake durability index. For the translational landslides, strength parameters were determined along the contact between the bedrock and the overlying colluvium. The results of the study indicate that multiple factors contribute to landslide susceptibility of the Kope Formation and the overlying colluvium, including low shear strength of the colluvial soil, development of porewater pressure within the slope, human activity such as loading the top or cutting the toe of a slope, low to very low durability of the bedrock that allows rapid disintegration of the bedrock and accumulation of colluvial soil, undercutting of the slope toe by stream water, and steepness of the slopes.","PeriodicalId":138906,"journal":{"name":"Environmental and Engineering Geoscience","volume":"25 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2021-03-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"126601191","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
L. McGuire, F. Rengers, N. Oakley, J. Kean, D. Staley, Hui Tang, Marian Orla‐Barile, A. Youberg
The extreme heat from wildfire alters soil properties and incinerates vegetation, leading to changes in infiltration capacity, ground cover, soil erodibility, and rainfall interception. These changes promote elevated rates of runoff and sediment transport that increase the likelihood of runoff-generated debris flows. Debris flows are most common in the year immediately following wildfire, but temporal changes in the likelihood and magnitude of debris flows following wildfire are not well constrained. In this study, we combine measurements of soil-hydraulic properties with vegetation survey data and numerical modeling to understand how debris-flow threats are likely to change in steep, burned watersheds during the first 3 years of recovery. We focus on documenting recovery following the 2016 Fish Fire in the San Gabriel Mountains, California, and demonstrate how a numerical model can be used to predict temporal changes in debris-flow properties and initiation thresholds. Numerical modeling suggests that the 15-minute intensity-duration (ID) threshold for debris flows in post-fire year 1 can vary from 15 to 30 mm/hr, depending on how rainfall is temporally distributed within a storm. Simulations further demonstrate that expected debris-flow volumes would be reduced by more than a factor of three following 1 year of recovery and that the 15-minute rainfall ID threshold would increase from 15 to 30 mm/hr to greater than 60 mm/hr by post-fire year 3. These results provide constraints on debris-flow thresholds within the San Gabriel Mountains and highlight the importance of considering local rainfall characteristics when using numerical models to assess debris-flow and flood potential.
{"title":"Time Since Burning and Rainfall Characteristics Impact Post-Fire Debris-Flow Initiation and Magnitude","authors":"L. McGuire, F. Rengers, N. Oakley, J. Kean, D. Staley, Hui Tang, Marian Orla‐Barile, A. Youberg","doi":"10.2113/EEG-D-20-00029","DOIUrl":"https://doi.org/10.2113/EEG-D-20-00029","url":null,"abstract":"\u0000 The extreme heat from wildfire alters soil properties and incinerates vegetation, leading to changes in infiltration capacity, ground cover, soil erodibility, and rainfall interception. These changes promote elevated rates of runoff and sediment transport that increase the likelihood of runoff-generated debris flows. Debris flows are most common in the year immediately following wildfire, but temporal changes in the likelihood and magnitude of debris flows following wildfire are not well constrained. In this study, we combine measurements of soil-hydraulic properties with vegetation survey data and numerical modeling to understand how debris-flow threats are likely to change in steep, burned watersheds during the first 3 years of recovery. We focus on documenting recovery following the 2016 Fish Fire in the San Gabriel Mountains, California, and demonstrate how a numerical model can be used to predict temporal changes in debris-flow properties and initiation thresholds. Numerical modeling suggests that the 15-minute intensity-duration (ID) threshold for debris flows in post-fire year 1 can vary from 15 to 30 mm/hr, depending on how rainfall is temporally distributed within a storm. Simulations further demonstrate that expected debris-flow volumes would be reduced by more than a factor of three following 1 year of recovery and that the 15-minute rainfall ID threshold would increase from 15 to 30 mm/hr to greater than 60 mm/hr by post-fire year 3. These results provide constraints on debris-flow thresholds within the San Gabriel Mountains and highlight the importance of considering local rainfall characteristics when using numerical models to assess debris-flow and flood potential.","PeriodicalId":138906,"journal":{"name":"Environmental and Engineering Geoscience","volume":"5 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2021-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"133509071","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Debris flows involve the high-speed downslope motion of rocks, soil, and water. Their high flow velocity and high potential for impact loading make them one of the most hazardous types of gravitational mass flows. This study focused on the roles of particle size grading and degree of fluid saturation on impact behavior of fluid-saturated granular flows on a model rigid barrier in a small-scale flume. The use of a transparent debris-flow model and plane laser-induced fluorescence allowed the motion of particles and fluid within the medium to be examined and tracked using image processing. In this study, experiments were conducted on flows consisting of two uniform and one well-graded particle size gradings at three different fluid contents. The evolution of the velocity profiles, impact load, bed normal pressure, and fluid pore pressure for the different flows were measured and analyzed in order to gain a quantitative comparison of their behavior before, during, and after impact.
{"title":"Experimental Investigation on the Impact Dynamics of Saturated Granular Flows on Rigid Barriers","authors":"N. Sanvitale, E. Bowman, M. Cabrera","doi":"10.2113/EEG-D-20-00033","DOIUrl":"https://doi.org/10.2113/EEG-D-20-00033","url":null,"abstract":"\u0000 Debris flows involve the high-speed downslope motion of rocks, soil, and water. Their high flow velocity and high potential for impact loading make them one of the most hazardous types of gravitational mass flows. This study focused on the roles of particle size grading and degree of fluid saturation on impact behavior of fluid-saturated granular flows on a model rigid barrier in a small-scale flume. The use of a transparent debris-flow model and plane laser-induced fluorescence allowed the motion of particles and fluid within the medium to be examined and tracked using image processing. In this study, experiments were conducted on flows consisting of two uniform and one well-graded particle size gradings at three different fluid contents. The evolution of the velocity profiles, impact load, bed normal pressure, and fluid pore pressure for the different flows were measured and analyzed in order to gain a quantitative comparison of their behavior before, during, and after impact.","PeriodicalId":138906,"journal":{"name":"Environmental and Engineering Geoscience","volume":"11 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2021-01-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"116896427","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Stable O and H isotope data distinguish three sources for base flow in five reaches of the San Pedro River: (A) base flow and sub-flow from upstream reaches of the river; (B) bank storage derived from summer monsoon floodwater; and (C) water from the mountainous flanks of the river catchment. A and C support base flow in the sub-basin upstream of Sierra Vista. A, B, and C combine to support base flow near St. David. Source C in this area is ancient deep-basin groundwater. Source C dominates in Cascabel near Benson Narrows, with downstream additions from A. In Cascabel near Gamez Road, sources A and C combined to support base flow that had disappeared by 2019. Near Redington, source C appears to have operated through a limestone aquifer vulnerable to short-term drought. Groundwater sub-basins separated by impermeable sills in the riverbed are evolving into hydrologically separate sub-basins as base flow across the sills decreases. The decrease in base flow partly reflects regional long-term drought, which has been exacerbated by pumping. Additional groundwater demand from urban growth upstream of Benson is likely to cause further decline of base flow near St. David and Sierra Vista.
{"title":"Sources of Perennial Water Supporting Critical Ecosystems, San Pedro Valley, Arizona","authors":"C. Eastoe","doi":"10.2113/eeg-d-20-00040","DOIUrl":"https://doi.org/10.2113/eeg-d-20-00040","url":null,"abstract":"\u0000 Stable O and H isotope data distinguish three sources for base flow in five reaches of the San Pedro River: (A) base flow and sub-flow from upstream reaches of the river; (B) bank storage derived from summer monsoon floodwater; and (C) water from the mountainous flanks of the river catchment. A and C support base flow in the sub-basin upstream of Sierra Vista. A, B, and C combine to support base flow near St. David. Source C in this area is ancient deep-basin groundwater. Source C dominates in Cascabel near Benson Narrows, with downstream additions from A. In Cascabel near Gamez Road, sources A and C combined to support base flow that had disappeared by 2019. Near Redington, source C appears to have operated through a limestone aquifer vulnerable to short-term drought. Groundwater sub-basins separated by impermeable sills in the riverbed are evolving into hydrologically separate sub-basins as base flow across the sills decreases. The decrease in base flow partly reflects regional long-term drought, which has been exacerbated by pumping. Additional groundwater demand from urban growth upstream of Benson is likely to cause further decline of base flow near St. David and Sierra Vista.","PeriodicalId":138906,"journal":{"name":"Environmental and Engineering Geoscience","volume":"65 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2020-11-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"131973637","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}