Next Article in Journal
Application of HEC-HMS Model for Flow Simulation in the Lake Tana Basin: The Case of Gilgel Abay Catchment, Upper Blue Nile Basin, Ethiopia
Previous Article in Journal
Modeling of GRACE-Derived Groundwater Information in the Colorado River Basin
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Hydrologic Trends in the Upper Nueces River Basin of Texas—Implications for Water Resource Management and Ecological Health

1
Department of Chemistry, Geosciences, and Physics, Tarleton State University, Stephenville, TX 76402, USA
2
Department of Engineering and Computer Science, Tarleton State University, Stephenville, TX 76402, USA
3
Department of Biological Sciences, Tarleton State University, Stephenville, TX 76402, USA
4
Texas Institute for Applied Environmental Research, Tarleton State University, Stephenville, TX 76402, USA
*
Author to whom correspondence should be addressed.
Hydrology 2019, 6(1), 20; https://doi.org/10.3390/hydrology6010020
Submission received: 14 February 2019 / Revised: 4 March 2019 / Accepted: 5 March 2019 / Published: 8 March 2019

Abstract

:
Reliable water sources are central to human and environmental health. In south Texas, USA, the Nueces River Basin (NRB) directly or indirectly plays that important role for many counties. Several NRB stream segments are designated as ecologically significant because they serve crucial hydrologic, ecologic, and biologic functions. The hydrologically significant streams recharge the Edwards Aquifer, an essential water source for the region’s agricultural, industrial, and residential activities. Unfortunately, the semiarid to arid south Texas climate leads to large inter-annual precipitation variability which impacts streamflow, and as a consequence, the aquifer’s recharge. In this study, we used a suite of hydrologic metrics to evaluate the NRB’s hydroclimatic trends and assess their potential impacts on the watershed’s ecologically significant stream segments using precipitation and streamflow data from the National Climatic Data Center (NCDC) and Hydroclimatic Data Network (HCDN) respectively from 1970 to 2014. The results consistently showed statistically significant decreasing streamflow for certain low-flow indicators over various temporal scales, likely due to water rights diversions and minimal land use changes. This research could help decision-makers develop the necessary tools to manage water resources in south Texas, given the NRB’s significance as a source of water for domestic consumption and ecological health.

1. Introduction

Access to fresh water resources plays a critical role in sustainable development, particularly in water-scarce regions. Semi-arid areas such as south Texas are highly water-stressed due to limited annual precipitation, high evapotranspirative losses, and prolonged droughts [1,2] as well as competing water demands [3]. Annual rainfall in this region varies from about 20 inches to 40 inches but annual potential evapotranspiration exceeds precipitation by 2 to 5 times [1].
Future climate change effects are not expected to bring any respite to water woes in the region either. The authors of [4] projected a trend towards more dry days as well as warming temperatures in this region, likely leading to increased evapotranspirative losses. Conclusions drawn by the researchers of [5] indicated that there may be no significant change in the amount of precipitation over the next 100 years, but rather that the duration between successive precipitation events will increase in the region. Several recent studies, for example [6,7,8], have shown that increasing temperatures may lead to greater unpredictability in precipitation patterns, which will in turn impact streamflow. These studies also emphasized the increased propensity for extreme hydrologic events such as droughts and floods in the future, resulting in negative impacts on water resource management and on aquatic ecosystems, particularly with regards to environmental flows.
Furthermore, the authors of [4] reported that variations in streamflow in the United States are predominantly caused by precipitation. Elsewhere, as in southwestern Australia, streams changed from perennial to ephemeral flow because of sharp declines in precipitation during the mid to late 2000s [9]. Several studies, for instance [10,11] have used the concept of elasticity developed by [12] to evaluate the changes in streamflow in response to changes in hydrologic variables such as precipitation or environmental variables such as land use and land cover.
As such, the south Texas region is characterized by large inter-annual variability in precipitation which results in uncertain streamflows, rendering them an unreliable source of water supply [1]. For instance, highly irregular rainfall characteristics in the NRB result in brief periods of high flows interspersed between long spells of low or zero flows [13]. Approximately 2% of the south Texas region consists of perennial water bodies; despite the sparse occurrence of surface water and the sporadic nature of streams in the region, most urban areas (with the exception of San Antonio) rely on surface water supplies to meet their needs.
Upon the passing of Senate Bill 1 (1997), water planning in Texas became a regional process, which required the Texas Water Development Board [TWDB] to identify 16 regional water planning regions (RWPAs). Of these, the majority of the NRB covers Region L (South Central Texas), but also extends into Region J (Plateau), Region M (Rio Grande), and Region N (Coastal Bend) as shown in Figure 1. The European Union has also adopted a similar style of regional water planning with their River Basin Management Plan [14]. The NRB is also a significant source of water for Region N, particularly for the City of Corpus Christi. According to [15], two of the three reservoirs that serve this city, Choke Canyon Reservoir and Lake Corpus Christi [jointly referred to as CCR/LCC], are located on the NRB. As such, the State of Texas owns all surface water in the state and appropriation thereof is governed by water rights; for example, water stored in the CCR/LCC represents nearly 98% of the water rights in the basin. However, the reliability of these rights is contingent on hydrologic conditions, particularly the occurrence of droughts. Another study [16], found that while droughts are common in this region, their severity has progressively increased; meanwhile, annual inflows to the CCR/LCC system have decreased by 40% from the 1950s to the 1990s. An 83% increase in total water use in Region N over the next 50 years has also been projected, underlining the need for prudent planning of water supply and use.
The NRB is unique from a hydrologic perspective due to the complex interactions that occur between streamflows and the subsurface environment, which can significantly alter groundwater recharge and thus indirectly impact water availability. One of the most productive karst aquifer systems in the world, the Edwards Aquifer of Texas, serves the water needs of the south-central Texas region, including the city of San Antonio, and supports several endangered species such as the Texas blind salamander (Typhlomolge rathbuni). This system receives significant recharge contributions from above-ground sources such as the Nueces River [18,19,20]. About 60% of the drainage basin recharges the groundwater system with about 14 billion ft3 of water per year [21]. Roughly 85% of Edwards Aquifer’s recharge comes from the watershed while the remainder originates from direct precipitation and subsurface flows from neighboring aquifers such as the Trinity [22]. The Nueces River, in turn, receives contributions from surface runoff as a result of rainfall events as well as groundwater (in the form of springs). The authors of [13] emphasized that the upper portion of the NRB that overlies the Edwards Aquifer outcrop region is “of particular interest to water managers in south-central Texas because appreciable streamflow gains and losses are observed along various reaches in the area”. In addition to its hydrologic value to the Edwards Aquifer and the reservoir systems that supply the Coastal Bend region, parts of the Nueces River have also been deemed to have significant ecological value, which is determined by criteria such as biological function, hydrologic function, the presence of threatened, endangered, or unique communities [23].
While there is an abundance of literature on (a) short-term or prevailing hydrologic conditions and water quality in the basin in the form of United States Geologic Survey [USGS] scientific reports (e.g., [13,24,25]), (b) the environmental conditions of downstream sections of the basin such as the Nueces estuary and associated ecosystem (e.g., [26,27,28]) and (c) groundwater-river water interactions in the basin (e.g., [18,29,30]), few studies have focused on hydrologic trends and their impacts on the basin, particularly on the upstream-most segments designated as ecologically-significant. Additionally, literature on the impact of variations, if any, in precipitation and land use and land cover (LULC) on streamflow in this basin is lacking.
It is apparent from the preceding discussion that alterations in streamflow of the upper Nueces River and those of its associated tributaries such as the Frio and Sabinal can have adverse consequences on inflows in to the CCR/LCC systems. This in turn will affect water availability in the Coastal Bend region, critically-impact recharge to the Edwards aquifer, as well as detract from the ecological uniqueness of areas designated as such. The sporadic nature of the river as well as its history of recurrent droughts led the authors of [1] to label the climate of this region as ‘problematic’. The over-arching goal of this study is, therefore, to evaluate the variability in streamflow in the upper-most part of this basin [near the headwaters] both from a water-planning perspective as well as for stream health assessment reasons. Specifically, we investigate historical spatio-temporal trends in streamflow using a suite of hydrologic indicators to identify any congruity between these metrics and variables such as precipitation and evapotranspiration (ET) in streamflow. The period of 1970 to 2014 was chosen for this study as consistent and reliable data for both streamflow and climate variables are available for this timeframe. Gauges that are part of the USGS HCDN 2009 were used as they reflect on sites with minimal or no human interference or disturbance.

2. Methodology

2.1. Background Information

A brief review of existing literature is presented in this section. The researchers of [21] assessed the strength of correlation between precipitation and climate indices such as the Pacific Decadal Oscillation [PDO] and El Nino Southern Oscillation [ENSO] in the Nueces, San Antonio, and Guadalupe River basins of south Texas; additionally, they also evaluated statistical trends in precipitation and streamflow between 1950 and 2009. They concluded that while precipitation and streamflow generally showed increasing trends during the period of interest, the correlation between these two variables was poor. Their inferences about hydrologic trends in the NRB are made from data recorded at a single USGS gauge located just upstream of the CCR/LCC system and do not reflect the most recent decades when significant warming was noted to have occurred in this region [6]. The authors of [31] simulated past (1986 to 2005) hydrologic conditions in the basin using a Variable Infiltration Capacity (VIC) model coupled with outputs from a suite of General Circulation Models [GCMs]. The VIC model was calibrated using streamflow data from two USGS HCDN gauging sites in the basin and then run with historical climate data output from four GCMs under the Representative Concentration Pathway [RCP] 4.5 over the period of 1966 to 2005. They showed that the performance of GCMs in capturing observed trends varied spatially but the GCMs were generally in agreement about the reduction in water availability at the outlet of the basin. Earlier studies by the authors of [32,33] suggested declining trends in runoff per unit rainfall in the Atascosa River Basin between 1935 and 1994. The Atascosa River is a tributary of the Frio River, which itself is a tributary of the Nueces River. Recently, the researchers of [7] conducted an extensive study of trends in hydrological drought indicators such as annual 7-day low flows [q7] and annual number of hydrological dry days across major US water resource regions; these areas are designated as USGS Hydrologic Unit Code Level-1 (HUC1) or 2-digit watersheds that encompass major river basins and associated watersheds. They used streamflow data from the USGS HCDN 2009 network [34] which consists of only those gauges that are minimally [or not] affected by human activities or influences. Their results indicated decreasing trends in q7 flows in ‘southeastern Texas’ but it must be noted that this is reflective on general trends across many HCDN sites in Texas and not confined to the NRB.

2.2. Site Description

The NRB study area, shown in Figure 2, is approximately 17,000 square miles in area and spans 20 south Texas counties delimited by Edwards County to the north, Webb County to the south, Maverick County to the west, and Atascoca County to the east [21]. The watershed is approximately 1600 ft. above mean sea level [35], and the Nueces River flows roughly 315 miles. to the Gulf of Mexico via Nueces Bay [31]. The NCDC classifies Texas into 10 climate divisions based on similarities in vegetation, weather, and climate. Nueces Basin’s main areas fall under climate divisions 6 - Edwards Plateau and 9 - South Texas Plains [36]. The region’s climate ranges from subtropical steppe to subtropical subhumid and is characterized by semi-arid to arid conditions with hot summers and dry winters [37]. Precipitation in the NRB region averages 20 inches to 30 inches annually [6].
As of 2009, the year the USGS reevaluated then reported the most recent list of unimpaired gauges in the HCDN [34], the NRB has 36 USGS-managed streamflow gauge sites [38] which continuously gather discharge data. Furthermore, the northern section of the NRB, which as previously mentioned is primarily covered by Region L, has 31 streams that are designated as ecologically significant [23]. As indicated in Figure 3, four of the six gauges that were used in the study are located on ecologically unique stream segments while all of the chosen gauge sites contribute to Edwards Aquifer recharge. Beyond the hydrologic function of the streams that is depicted in Figure 3, Table 1 further highlights the numerous invaluable roles the stream segments fulfill based on the TWDB’s ecologically-significant designation criteria [39].

2.3. Data Sources

The period 1970 to 2014 was used to determine the hydroclimatic variability for streamflow and precipitation in the study area. That particular timeframe was chosen because (a) it provided consistent and reliable data for both hydrological and climate variables, (b) the 45-year duration ensured statistical legitimacy for the study since authors of [40,41] indicated that at least a 15- to 25-year period is required to determine statistical validity for streamflow trend analysis, and (c) we aimed to extend the most recent time series for hydroclimatic analyses in the region done by [21] which went until the year 2009, the year the USGS made the HCDN revisions.
United States Geological Survey-verified continuous daily mean discharge data were acquired from the HCDN online database [42] in units of cubic feet per second (cfs). Using HCDN gauges ensured that minimal anthropogenic influence occurred thereby increasing the likelihood of attributing any streamflow variation to factors other than direct human influence. The work done by [43] documented a comprehensive USGS-sanctioned Geographic Information Systems (GIS) analysis of U.S. land-use trends from 1974 through 2012 as shown in Figure 4. The findings showed minimal land modification from under conservation to developed production in the northern NRB during the aforementioned time period, thereby fulfilling the main HCDN criteria. The northern NRB was also delineated into the subbasins that contribute to each stream gauge, also shown in Figure 4, to focus on specific areas of LULC change in relation to the locations of the stream gauges to estimate any likely effects on the HCDN sites if trends were later determined. The most notable land alterations, from conservation to developed, were relatively small and occurred during the 1974 to 1982 time period in areas that contribute to the 08196000 and 08198500 gauges. Additionally, the range of daily mean flows measured at each gauge over the 45-year study period is shown in Figure 5. Flow variations have been presented gauge-wise, from westernmost (gauge 08195000; left corner) to easternmost (gauge 08200000; right corner) in this figure. These violin plots are essentially a combination of box plots and kernel density plots that allow for visualization of the shape of the distribution and identification of possible mass clusters. The median of daily flows over the 45-year period have been shown using the diamond-plus symbol; the median flows for all of the selected gauges are less than 150 cfs. It can also be seen from Figure 5 that a large percent of the mass is concentrated on the low flows with gauge 08196000 showing the most pronounced ‘bulge’ at the low-flow ends. The large variability in flows, spanning four orders of magnitude is also evident from the figure.
Daily precipitation depths, measured in inches, were acquired from National Oceanic and Atmospheric Administration’s (NOAA) NCDC online database [44]. The choice of precipitation gauge stations used in the study was made entirely on the basis of availability of serially-continuous data. Some stations had near-complete records for the 45-year period of analysis but there were occasional gaps in the series. These gaps were filled in using spatial interpolation of synchronous observations from adjacent stations. Other stations did not have appreciable completeness or continuity, so data from multiple stations had to be consolidated in these cases. Therefore, the rationale here is the completeness of the precipitation data (without gaps) rather than the nearness of precipitation station to the discharge measurement location. For instance, Uvalde and Medina Counties had sufficient precipitation data while Kinney County did not, so gauges from nearby Zavala County was used as its substitute, as shown in Table 2.
Evapotranspiration is a major hydrological output in the water budget of south Texas where annual PET can exceed annual precipitation by 2 to 5 times [1]. Evapotranspiration data were adapted from [45] and evaluated for each selected gauge location. The author of [45] estimated the average change in actual ET trends over the contiguous United States for the period 1979 to 2015 using complementary relationship (CR) based ET estimation methods suggested by [46]. The research in [45] validated the ET estimates using precipitation and runoff data from the Parameter-Elevation Regressions on Independent Slopes Model (PRISM) and USGS Hydrologic Unit Code Level-6 (HUC6) or 12-digit watersheds respectively.

2.4. Hydroclimatic Data Analysis

To evaluate the relationships between the hydroclimatic variables, a range of hydrological assessment metrics specific for nonparametric datasets where applicable were used: modified Mann-Kendall test (MMK), annual minimum and 7-day low flow indicators, streamflow elasticity, and drought indices [Streamflow Drought Index (SDI) and Standardized Precipitation Index (SPI)]. Using a number of indicators was critical to determine potential relationships between hydroclimatic variables which could ultimately help indicate the extent to which variations in precipitation affected stream discharge. Consequently, trend analysis using MMK as well as hydrological sensitivity assessments between precipitation and streamflow were investigated along with streamflow elasticity and drought indices to understand the possible relationships. Low flows, which are necessary for maintaining certain essential ecological functions and overall stream health, were assessed using the previously stated low flow metrics and analyzed in conjunction with the trends and hydrological sensitivity to determine if potential streamflow variations were mainly associated with changes in precipitation.

2.4.1. Trend Analysis (MMK)

The MMK is a reliable and popular nonparametric trend analysis tool that was used to detect the hydroclimatic trends. Numerous researchers, for example [21,47,48,49,50,51,52,53], have used various forms of the Mann-Kendall test. The authors of [50] as well as [51] indicated that for accurate hydroclimatic trend analyses the input data must be serially independent to prevent serial correlation or autocorrelation. Serial correlation occurs when errors in a particular timespan transfer into subsequent time periods [54], which could cause an overestimation or underestimation of trends [50,51]. Therefore, the MMK, which was derived by [55] and designed to account for serial correlation by modifying the variance of the original Mann-Kendall test [51], was chosen.
The original Mann-Kendall test statistic, which is shown in Equation 1, was developed by [56,57]. A detailed explanation along with the appropriate equations for the original trend analysis method Equations (1) and (2) where n, Xi, Xj, and D are sample size, consecutive data values, and the difference between successive values respectively were presented by authors of [10,47,51].
S =   i = 1 n 1 j = i + 1 n s g n ( X j   X i )
s g n ( X j X i ) =   { 1   i f   D > 0 0   i f   D = 0 1   i f   D < 0
The aforementioned authors, [10,47,51], also provided the test statistic Equation (3) and equations for the modified version of the Mann-Kendall test as shown in Equation (3) through Equation (7) where V(S)* is the modified variance from V(S) in Equation (5) and rk is the lag-k autocorrelation coefficient Equation (7):
Z =   { S 1 V ( S ) *   f o r   S > 0 0   f o r   S = 0 S + 1 V ( S ) *   f o r   S < 0
V ( S ) *   = V ( S ) n n *
V ( S ) = n ( n 1 ) ( 2 n + S ) 18
n n * = 1 +   2 n ( n 1 ) ( n 2 )   i = 1 n 1 ( n i ) ( n i 1 ) ( n i 2 ) r i
r k = 1 n k i = 1 n k ( X i X ¯ ) ( X i + k X ¯ ) 1 n i = 1 n ( X i X ¯ ) 2
The null hypothesis (H0) for the MMK is “no monotonic trends exist in the time series” [10], and it was rejected at a ≤ 0.01 significance level, while the alternate hypothesis (HA) is “there are monotonic trends in the time series”. The hydroclimatic trends were determined using a combination of their Z values and corresponding P values. Increasing and decreasing trends were represented by positive and negative Z values respectively [10,47,48].
Hydrological trend analysis was performed for minimum, median, and maximum stream discharge and precipitation over monthly, seasonal, and annual temporal scales. Annual assessments were done on a calendar year scale. June to August represented summer months, September to November was fall, December to February signified winter, and March to May denoted spring. Trend assessments were done for the periods of 1970 to 2014. Additionally, we performed trend analyses for the period 1994 to 2014. Trends in mean annual temperatures in South Texas between 1930 and 2001 were studied by [1]. These authors fit regression equations for temperature based on data collected from 16 stations across the study region, which generally overlaps the NRB. Their findings suggest an upward trend in temperature around the mid-1990s, thus serving as the rationale behind distinguishing the 1994 to 2014 period for separate analysis. The results from both sets of trend analysis were later compared against each other. The “fume” package in R [58] was used to perform the MMK.

2.4.2. Low Flow Metrics (Annual Minimum and 7-day Low Flows)

The annual 7-day low flow of a stream is the lowest average flow that occurs within a consecutive 7-day period during any given year while the annual minimum flow measures the lowest discharge for each year. Low flow plays an essential role in the ecological balance of riverine systems. Such flow is critical for maintaining environmental flows, water quality, biodiversity, natural migration trends, ecosystem integrity, and overall stream health. The south Texas region is known to experience harsh summer conditions, unpredictable precipitation patterns, and widely variable flow, as exemplified in Figure 6, where representative examples were chosen. Low flow assessments for the selected gauges in the NRB were therefore necessary because of the important hydrological, ecological, and various other functions the stream segments provide to the downstream regions due to their ecologically significant designations. The low flow metrics analysis was done using the “FlowScreen” package in R, which was developed by [59], while the MMK was used to determine the trends from the low flow indicators.

2.4.3. Streamflow-Precipitation Elasticity

Streamflow elasticity is a measure of mean streamflow variation in proportion to changes in another driving variable such as precipitation [11,60]. The analysis was used because it quantifies stream discharge sensitivity to environmental changes [11,60]. The author of [12] first presented the streamflow elasticity concept for precipitation variations using the expression shown in Equation (8), where Q and P are stream discharge and precipitation respectively [11]. The researchers of [60] refined the concept to standardize results using median as well as mean streamflow and precipitation values as shown in Equation (9). In Equation (9), P ¯ ,   Q ¯ , Q t ,   and   P t denote mean annual streamflow and precipitation as well as yearly stream discharge and precipitation respectively [10,11,60]. To perform the elasticity calculations, the stream discharge and precipitation units were both standardized to depths in inches.
ϵ p ( P ,   Q ) = d Q / Q d P / P = d Q d P   P Q
ϵ p = m e d i a n ( Q t Q ¯ P t P ¯   Q ¯ P ¯ )

2.4.4. Drought Indices (SDI and SPI)

The SDI and SPI were evaluated to assess the NRB’s hydrologic extremes. The developers of the SDI, [61], along with [62], noted that the indices require minimal amounts of data and computational effort, which make them highly advantageous over earlier indices developed for determining water availability and variability patterns. The SDI creators, [61], indicated that to compute the SPI a times series of monthly precipitation, which is denoted by Pi,j where i is the hydrological year and j is the month in the hydrological year, is required to produce the series as shown in Equation 10. From the series depicted by Equation (10), a sequence where Rj,k represents the depth of cumulative precipitation for the k-th reference period, which consists of overlapping time intervals of 3, 6, 9, and 12 months in the i-th hydrological year, and reference periods k is derived. In the SPI calculation shown in Equation (11), the mean and standard deviation of the cumulative precipitation depths are represented by R ¯ k and SR,k respectively for the k-th reference period for statistical estimations over an extensive timespan.
R i , k = j = 1 3 k P i , j   i   = 1 , 2 , ,   j = 1 , 2 , , 12 ,   k = 1 ,   2 ,   3 ,   4
S P I i , k =   R i , k R ¯ k S R , k   i = 1 , 2 ,   k = 1 , 2 , 3 , 4
Additionally, during drought indices computation the usually skewed probability distributions of hydroclimatic data require transformation to a normal distribution, which is typically done using a Gamma distribution [62]. The computational representations for the logarithmic transformation are shown in Equations (12) and (13) where ln(Ri,k), w ¯ k , and Sw,k are the natural logarithms of cumulative precipitation, mean, and standard deviation respectively over a lengthy timeframe.
S P I i , k = w i , k w ¯ k S w , k   i = 1 , 2 ,   k = 1 , 2 , 3 , 4
w i , k = ln ( R i , k )   i   = 1 , 2 , ,   k = 1 ,   2 ,   3 ,   4
Converting the time series to normal distributions also ensures that they have a mean of zero and standard deviation of 1, which allows for easy interpretation of the results. Index values between −1 and +1 represent average conditions, <−1 indicates varying degrees of drought, and >+1 shows wetter than average conditions. Furthermore, the SPI and SDI computations use the same methodology [61] but different univariate hydroclimatic input data. Precipitation and streamflow data are required for the SPI and SDI respectively. When both indices are analyzed together, the SDI lags the SPI. The drought indices were calculated using a 12-month time-scale as droughts usually take multiple seasons or even years to manifest in hydrologic variables.

3. Results

The results for each hydrological assessment metric are presented sequentially, from the westernmost gauge, 08190500, to easternmost gauge, 08200000, in Kinney and Medina Counties respectively. Additionally, the notations Qmin, Qmed, and Qmax correspondingly represent minimum, median, and maximum discharge.

3.1. MMK Trend Analyses – 1970 to 2014

Few statistically significant (p ≤ 0.01) precipitation trends were detected for the entire study period. Significant decreases in annual maximum precipitation were only identified in Uvalde County during the spring and fall seasons as well as the months of April. Contrastingly, many statistically significant annual and seasonal stream discharge trends were found for streamflow, as indicated in Table 3. An overwhelmingly greater number of decreasing trends were identified in the Qmin and Qmed flows when compared to Qmax. Decreasing flows were most often detected in Qmin and Qmed in the NRB’s central to easternmost gauges (08196000, 08195000, 08198500, and 08200000) found in Uvalde and Medina Counties for various temporal periods. Except for a few instances, the previously listed gauges simultaneously experienced decreasing trends over the annual time scale, fall and winter periods, as well as most of the months that corresponded to the said seasons (September, October, November, December, and January) for both Qmin and Qmed. The gauges that are located in the western part of the basin (08190500 and 08190000) collectively experienced the least number of statistically significant trends for any flow level. No trends were detected in the westernmost gauge, 08190500, over any of the temporal scales during the 1970 to 2014 timespan, while very few trends were identified for gauge 08190000 throughout the same period. Statistically significant decreases in Qmax were only found in gauges 08195000 (winter, fall, January, and October) as well as 08200000 (August, September, October, and November). July was the only month in which no statistically significant trends were identified.

3.2. MMK Trend Analyses – 1994 to 2014

No statistically significant precipitation trends were found for the last 20 years of the study, but significant decreasing streamflow trends were detected for annual and seasonal periods, as shown in Table 4. The 1994 to 2014 period had more significant declines in Qmax in comparison to the full 45-year assessment. Statistically significant decreases in Qmin, Qmed, and Qmax were more often detected in the westernmost to central gauges, which again differed from the 1970 to 2014 timespan, where the central to easternmost gauges largely experienced a greater number of declines. Decreases in Qmed and Qmax were also identified in each gauge except those on the eastern end of the basin in the fall and winter seasons, while only gauges 08196000 and 08200000 did not have significant declines during winter in Qmin. Few statistically significant trends were found for the annual stream discharge (08196000, 08195000, and 08200000 for Qmed as well as 08190500 and 081900000 for Qmax), the spring seasons (08190500, 08195000, and 08198500 in Qmin), and summer periods (08200000 in Qmin only). The months that collectively experienced the most streamflow declines in Qmin, Qmed, and Qmax were October, November, and December. However, the months of April to September which practically covered the entire spring and summer seasons together had the least decreases in streamflow in Qmin, Qmed, and Qmax. No statistically significant trends were found for the months of May and August.

3.3. Low Flow Indicators

Statistically significant decreasing trends were detected in both the annual minimum and 7-day low flows for two central gauges (08196000 and 08195000) and the easternmost gauge (08200000). Contrastingly, no significant trends were found for the westernmost gauges (08190500 and 08190000) as well as one central gauge (08198500) for the low flow metrics analyses.

3.4. Streamflow-Precipitation Elasticity

Table 5 shows the streamflow-precipitation elasticity that was calculated for each gauge site. The elasticity of gauges 08190500 and 08198500 with Zavala and Uvalde Counties respectively were close to 1:1 ratios. In both cases, the hydroclimatic sensitivity essentially indicated that a 1% increase in annual precipitation would have theoretically resulted in an almost equivalent increase in annual stream discharge. For the remaining gauges however, there was less sensitivity between the hydroclimatic variables since the elasticity coefficients were less than 1. The results signaled that annual precipitation fluctuations had a smaller impact on yearly streamflow. In theory, a 1% variation in annual precipitation would have produced a 0.47%, 0.60%, 0.41%, and 0.72% change in annual streamflow for gauges 08190000, 08196000, 08195000, and 08200000 respectively.

3.5. Drought Indices (SDI and SPI)

The SPI and SDI were evaluated to assess NRB’s hydroclimatic variables for the 45-year study period for anomalous hydrologic occurrences such as droughts in the SDI, which is defined as having results that are <−1 [62]. The changes in precipitation were later followed by similar variations in streamflow. In each corresponding stream discharge and precipitation dataset, there were multiple fluctuations between wetter and drier conditions as well as droughts in varying degrees of severity. As shown in Figure 7, normal to wetter conditions dominated in the first 30 years for both indices, with few recorded instances of droughts. When droughts occurred during the first 30 years, they were mild, most often in the SPI, and had short durations. In the last 15 years of each dataset however, downward trends were evident for both drought indices. The SDI and SPI also revealed that each stream segment and its corresponding county underwent multiple periods of abnormally dry conditions during the last 15 years of the 45-year study period. In each case, more severe drought conditions were first manifested in the SPI due to fluctuations in precipitation then lagged by less harsh droughts in the SDI.

3.6. ET Trends

Table 6 shows the ET results for each gauge location that was extracted from [45], who estimated the average changes in ET over the conterminous United States for the period 1979 to 2015. The effect of the ET changes ranged from positive to negative values, which indicated increases and decreases in ET respectively, across the western to eastern areas of the watershed. The data showed that the westernmost gauges (08190500 and 08190000) were most affected by ET. The central gauge sites (08196000 and 08195000) were practically unaffected by changes in ET while the eastern ones (08198500 and 08200000) were the least impacted.

4. Discussion

This study used a collection of hydrologic indicators to robustly investigate the hydroclimatic relationships of ecologically significant streams in the upper regions of the NRB and the potential impacts of hydrologic alterations for the period 1970 to 2014. The novel aspects of this study are that we (a) have focused on the upper-NRB where the ecologically significant segments lie and where no literature regarding the impact of hydrologic changes on these segments exists and (b) use a diverse suite of hydrological indicators to capture the response and sensitivity of the river to changes in external drivers. We have also updated the most recent time series used in hydroclimatic analyses of the NRB by the authors of [21] from 2009 to 2014. Additionally, we make recommendations to policy-makers and natural resource managers on holistic and sustainable approaches to water resource planning in the region.
The results showed occurrences of statistically significant streamflow declines in the selected gauges, but precipitation in the corresponding counties overwhelmingly showed no trends. The MMK assessments detected significant decreasing stream discharge trends in the central to easternmost gauges, mainly in Qmin and Qmed for the 1970 to 2014 period, as well as more uniform declines over each flow level (Qmin, Qmed, and Qmax) for the westernmost and central gauges during the most recent 20 years of the study. The decreases in stream discharge were consistent with earlier studies of the spatiotemporal hydroclimatic relationships in the basin in [32,33], which concluded that runoff per unit precipitation in the Atascosa River Basin, a subbasin in the upper NRB, from 1935 and 1994 showed declining trends.
The droughts that were experienced in the NRB, particularly ones during last 15 years, were first manifested in the SPI followed by the SDI. The SPI, a precipitation based index, detected the trends in drought early on; whereas, the SDI, which is based on streamflow, also identified the drought trends but with an expected lag or delay. The lag between the two drought indices suggested that a relationship exists between precipitation and streamflow. The degree to which the hydroclimatic variables are related was shown to be moderate to high by the streamflow-precipitation elasticity. However, because statistically significant decreases in precipitation were only detected in Uvalde County during spring, fall, and the months of April for the 45-year study period, it was clear that precipitation was not the only reason for stream discharge declines. The authors of [63] indicated that precipitation that eventually becomes streamflow generally encounters a lag period because of factors such as basin size and antecedent soil moisture. Consequently, the precipitation reductions in the fall seasons coupled with extended lag times could have led to streamflow declines in the successive fall and winter periods. Plus, with the statistically significant declines in spring precipitation, decreased flows in Qmin and Qmed were mainly detected in gauges in Uvalde County during that season. Summer and its corresponding months virtually showed very few statistically significant decreasing trends in streamflow for both 45- and 20-year time periods. Texas’ semi-arid to arid climate experiences hot summers [37] and variable precipitation [1]. During summer months, many Texas streams suffer from reduced and in some cases zero flows. Therefore, the typical summer conditions could possibly explain why practically no streamflow trends were detected during that time period. However, due to the demonstrable lack of statistically significant precipitation trends, it was apparent that other reasons contributed to the declining streamflows.
The authors of [63] noted that stream discharge trends sometimes do not reflect changes in precipitation because of a combination of climatological variables. The dataset from [45] assessed changes in actual ET trends from 1979 to 2015, and it was used as a surrogate for temperature variation for each gauge location. The effects of ET on streamflow, however, were not conclusive for either the annual 45- or 20-year periods. The assessments from [45] indicated that changes in ET trends were positive in the western areas of the basin, which suggests that ET had a greater effect on the westernmost gauges (08190500 and 08190000). The central to easternmost gauges (08196000, 08195000, 081985000, and 08200000) all experienced little to no impact from ET as the [45] dataset showed negligible changes in ET trends in those areas. The ET results were not indicative of the likely locations for decreasing discharge in the 45-year period as trend analyses showed that streamflow declines overwhelmingly occurred in the central to easternmost gauges. When the changes in ET trends were compared with the annual decreases seen in the last 20 years, only in Qmax did statistically significant declines often occur in the western portion of the basin, where changes in ET trends were also positive. In the other flow levels, no significant trends were found for Qmin while Qmed only experienced decreases in the locations that span the central to eastern gauges. Similar to the 45-year period, the trend analyses in Qmed for the 20-year timespan conflicted with the ET findings which showed that the central to eastern areas had the least ET effect. The general discrepancy between the streamflow and changing ET trends was likely because the ET assessments were not basin-specific as noted by the researchers of [64]. The authors suggested that hydroclimatic variables should be analyzed at the local rather than regional scale for the most accurate analysis since trend effects on extensive scales essentially normalize the final result which may not account for the subtle variances on a localized level. Consequently, other environmental factors were considered.
As previously mentioned, [43] documented U.S. LULC trends from 1974 through 2012 using GIS analysis. Different land uses either impede or facilitate stream discharge. Therefore, investigating temporal changes in land use from being under conservation to being developed with time was necessary to indicate if the NRB region was affected by manmade influences that would have altered the basin’s hydrological characteristics, hence affecting the area’s HCDN status. The analysis indicated that the NRB had undergone LULC alterations, particularly in its southern end, as shown in Figure 4. More importantly however, the analysis revealed that the upper-most basin regions, which contribute to the gauges in our study, remained virtually unchanged during the 1974 to 2012 timespan with only marginal modifications to undeveloped areas thereby fulfilling the HCDN’s minimal anthropogenic disturbance criteria. As Figure 4 also highlights, the most modified regions fed the 08195000 and 08198500 gauges, both in Uvalde County, and those changes occurred sometime within the period of 1974 to 1982.
Vegetative cover, and by extension LULC, directly impact different areas of the hydrological cycle such as soil moisture, ET, runoff, and infiltration [48,65] as well as streamflow. Greater vegetative cover buffers precipitation and ultimately decreases runoff, and as a consequence, stream discharge [66]. Some brush species also impact certain hydrologic processes as they alter spring flows and stream discharge that is interconnected with aquifer systems [67]. As a result, land use changes in the NRB, while minimal, could have still impacted the hydrological dynamics in the region. Furthermore, studies by [68,69] explain that land use changes can lead to declines in stream discharge in areas that are susceptible to low flows. The minimal LULC modifications could therefore signal an increase in human consumption along the course of the different stream segments, potentially leading to decreased streamflows as less water arrives to the HCDN gauges.
The authors of [49,51,63] all agree that when stream discharge trends show results that are inconsistent with climatological variation, anthropogenic influences such as water management might be a reason for the changes. As mentioned previously, Texas’ state-owned water is managed using water rights, which are granted to select users via water rights permits. The Texas Watermaster Program enforces water permits, but only in a handful of watersheds such as the Rio Grande and Brazos does the body manage the entire river basin. In the case of the NRB, some regions are under the control of the South Texas Watermaster [70], but the upper sections, which were part of our study, are not. Therefore, water rights owners exceeding their water use quotas could be another possible reason why the declines were detected by the gauges. Additionally, many water permit holders have access to the ecologically significant stream segments that were studied. Their water usage activities are concentrated upstream of the gauges and primarily involve diversion points, which are authorized locations on the channel of a stream where water can be legally removed. After closer review of the non-Watermaster records, it was confirmed that many permit owners, particularly in Uvalde and Medina Counties, are found in areas without Watermaster supervision. The impact of watershed management practices on water quality trends in rivers in the United States has also been recently pointed out by [71].
Furthermore, water used by permit holders in the non-Watermaster areas are supposed to be annually self-reported [72]. When archived reports for the previously named counties were evaluated from [72], some permittees did not report their usage during different time periods, which is mandatory even if they did not utilize their designated quotas [73]. Therefore, should permit holders who operate in non-Watermaster areas subvert the honor system and not comply with their permits, they could potentially overuse and underreport the amount of water that they were assigned. Years of sustained overexploitation could have been reflected in the declining stream discharge, particularly during low streamflow periods, which could explain why the low flow metrics showed instances of statistically significant decreases. Water rights are prioritized based on seniority among permit owners using priority dates, so permittees with earlier water rights are given greater preference during lower flow periods [74]. Consequently, users in non-Watermaster areas could withdraw more than their permits allow, yet go unnoticed unless more senior water rights holders are impacted to the point that they file complaints for investigation.

5. Conclusions

In conclusion, the study used a collection of hydrologic metrics (MMK trend analyses, low flow indicators, streamflow-precipitation elasticity, as well as drought indices) on six USGS gauges in the upper NRB to robustly investigate the area’s hydroclimatic relationships and the potential impacts of those associations to south Texas. The hydroclimatic analyses of the chosen stream gauges showed the relationship between the studied hydrological variables in the upper NRB existed. The assessments also detected significant decreasing stream discharge trends in a number of gauges for both the 1970 to 2014 period and during the last 20 years of the study. Although the causes of the declining streamflow trends could not be totally attributed to rainfall fluctuations, a combination of other climatological, environmental, and water management factors were assessed as probable contributors to reductions in flow. The decreasing trends should be of concern, however. The declining stream discharge could see negative impacts on the stream segments that were designated as ecologically significant, particularly for their hydrologic, biological, and ecological functions to the NRB, and by extension, south Texas. More importantly, the fact that decreasing streamflow was detected in all of the gauges at some point or the other over the study period should be of interest to the region’s water and natural resource managers as the declines may have implications for current as well as future water-planning and aquatic health. It must also be noted that our findings will likely have broader impacts beyond the NRB considering that other parts of the state have stream segments designated as ecologically-significant. As indicated in [75], these include segments in Regions C, H, J, K, and L, shown in Figure 1 for spatial reference.
Continued declining stream discharge trends could directly impact the hydrologic functions of the NRB and potentially lead to negative implications for the economic, agricultural, and industrial sectors for south Texas, as the region’s surface and groundwater supply heavily depends on the drainage basin. Additionally, if the streamflow that traverses the NRB down to Nueces Bay continues to decline, changes in water chemistry, disruptions to nutrient movement via nutrient spiraling, biological functions, biodiversity, natural organismal migration patterns, reproduction habits of aquatic life, as well as the overall ecological health of the region could be impaired. Further decreases in streamflow could spell even more danger for the already vulnerable endemic, threatened, and endangered populations and ecosystems of the Blanco blind salamander (Eurycea robusta), Comal blind salamander (Eurycea tridentifera), Plateau shiner (Cyprinella lepida), and Blanco River Springs salamander (Eurycea pterophila), for example, which are highlighted in [23].
As it specifically relates to the Nueces Bay, reduced stream discharge could be detrimental to its aquatic species, and by extension, the fisheries industry in that marine environment. The researchers of [31] noted that the Nueces Bay relies on the flows from the NRB to dilute the saline water to make the coastal location habitable for brackish water species. With decreasing stream discharge and rising global sea levels because of changes in climate [76] as well as thermal expansion due to global warming [77], the ability to buffer the saline conditions in the Nueces Bay would become more difficult and possibly make the coastal environment inhospitable to coastal marine life.
Furthermore, riverine plants, specifically riparian vegetation and macrophytes, play vital roles in freshwater environments and overall stream health such as, for maintaining water quality, biodiversity, spawning and breeding grounds, aquatic habitats, as well as bank stabilization [78]. With sustained periods of unusual flow patterns such as reduced stream discharge, the diversity of riparian plant species changes from more to less water-tolerant [79] thereby impacting the important functions riverine plants provide to lotic environments and emphasizing the detrimental effects declining flows could cause.
As seen in the results and preceding content, the effects of declining stream discharge could be detrimental to numerous aspects of ecological health as well as the hydrological balance of the NRB. Therefore, it is imperative that water managers and resource planners consider strategies to mitigate the potential impacts of such effects. The possible consequences of the distribution of water rights should be evaluated and the mechanisms to emphasize the significance of truthful reporting of water use quotas, particularly since the upper regions of the NRB is in a non-Watermaster area, should be considered.

Author Contributions

Each author played a substantial role in the research effort. The data curation, formal analysis, and investigation were completed by E.D.T. Conceptualization, methodology, project administration, and supervision were done by K.V. The responsibility of completing the original draft and visualization fell on E.D.T. and K.V., while the project’s validation was handled by V.C., N.K. and K.V. Each participant provided resources at various points of the research as well as reviewed and edited the document.

Funding

This research received no external funding.

Acknowledgments

We would like to acknowledge the insightful comments and suggestions provided by the reviewers that helped enhance the quality of this work. Tarleton State University’s College of Science and Technology, particularly The Center for Environmental Studies, also deserves special mention for sponsoring E.D.T.’s graduate research.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Norwine, J.; John, K. The Changing Climate of South Texas, 1900-2100: Problems and Prospects, Impacts and Implications; CREST-RESSACA Texas A & M University-Kingsville: Kingsville, TX, USA, 2007. [Google Scholar]
  2. Sturdivant, A.W.; Rogers, C.S.; Rister, M.E.; Lacewell, R.D.; Norris, J.W.; Leal, J.; Garza, J.A.; Adams, J. Economic costs of desalination in south Texas: A case study. J. Contemp. Water Res. Educ. 2007, 137, 21–39. [Google Scholar] [CrossRef]
  3. Uddameri, V.; Honnungar, V.; Hernandez, E.A. Assessment of groundwater water quality in central and southern Gulf Coast Aquifer, TX using principal component analysis. Environ. Earth Sci. 2014, 71, 2653–2671. [Google Scholar] [CrossRef]
  4. Clark, J.S.; Melillo, J.; Mohan, J.; Salk, C. The seasonal timing of warming that controls the onset of the growing season. Glob. Change Biol. 2014, 20, 1136–1145. [Google Scholar] [CrossRef] [PubMed]
  5. Tebaldi, C.; Hayhoe, K.; Arblaster, J.M.; Meehl, G.A. Going to the extremes. Clim. Chan. 2006, 79, 185–211. [Google Scholar] [CrossRef]
  6. Venkataraman, K.; Tummuri, S.; Medina, A.; Perry, J. 21st century drought outlook for major climate divisions of Texas based on CMIP5 multimodel ensemble: Implications for water resource management. J. Hydrol. 2016, 534, 300–316. [Google Scholar] [CrossRef] [Green Version]
  7. Poshtiri, M.P.; Pal, I. Patterns of hydrological drought indicators in major U.S. river basins. Clim. Chan. 2015, 134, 549–563. [Google Scholar] [CrossRef]
  8. Hernandez, E.A.; Uddameri, V. Standardized precipitation evaporation index (SPEI)-based drought assessment in semi-arid south Texas. Environ. Earth Sci. 2014, 71, 2491–2501. [Google Scholar] [CrossRef]
  9. Petrone, K.C.; Hughes, J.D.; Van Niel, T.G.; Silberstein, R.P. Streamflow decline in southwestern Australia, 1950–2008. Geophys. Res. Lett. 2010, 37, 11. [Google Scholar] [CrossRef]
  10. Kibria, K.N.; Ahiablame, L.; Hay, C.; Djira, G. Streamflow trends and responses to climate variability and land cover change in South Dakota. Hydrology 2016, 3, 2. [Google Scholar] [CrossRef]
  11. Kim, B.S.; Hong, S.J.; Lee, H.D. The potential effects of climate change on streamflow in Rivers Basin of Korea using precipitation elasticity. Environ. Eng. Res. 2013, 18, 9–20. [Google Scholar] [CrossRef]
  12. Schaake, J.C. From climate to flow. In Climate Change and U.S. Water Resources; Waggoner, P.E., Ed.; John Wiley and Sons Inc.: New York, NY, USA, 1990; pp. 177–206. [Google Scholar]
  13. Banta, J.R.; Lambert, R.B.; Slattery, R.N.; Ockerman, D.J. Streamflow Gain and Loss and Water Quality in the Upper Nueces River Basin, South-Central Texas, 2008–10; U.S. Geological Survey Scientific Investigations Report 2012–5181; U.S. Geological Survey: Reston, VA, USA, 2012; 40p. Available online: http://pubs.usgs.gov/sir/2012/5181 (accessed on 15 April 2017).
  14. Skoulikaris, C.; Zafirakou, A. River Basin Management Plans as a tool for sustainable transboundary river basins’ management. Environ. Sci. Pollut. Res. 2019. [Google Scholar] [CrossRef] [PubMed]
  15. Coastal Bend Regional Water Planning Group; Henningson, Durham and Richardson Engineering, Inc. Coastal Bend Regional Water Planning Area Region N Regional Water Plan Volume II Water Management Strategies; Texas Water Development Board: Austin, TX, USA, 2010. Available online: http://www.twdb.texas.gov/waterplanning/rwp/plans/2011/N/Region_N_2011_RWPV2.pdf (accessed on 20 June 2017).
  16. Coastal Bend Regional Water Planning Group; Henningson, Durham and Richardson Engineering, Inc. Coastal Bend Regional Water Planning Area Region N; Texas Water Development Board: Austin, TX, USA, 2015. Available online: https://www.twdb.texas.gov/waterplanning/rwp/plans/2016/N/Region_N_2016_RWP.pdf (accessed on 20 June 2017).
  17. Texas A&M AgriLife. Water in Texas. 2018. Available online: https://agrilife.org/siwi/water_in_texas/ (accessed on 7 November 2018).
  18. Kromann, J. Surface Water Recharge in Karst: Edwards-Trinity Aquifers-Nueces River System; Texas Scholar Works; University of Texas Libraries: Austin, TX, USA, 2015; Available online: https://repositories.lib.utexas.edu/handle/2152/44401 (accessed on 18 June 2017).
  19. Mace, R.E.; Anaya, R. Estimate of Recharge to the Edwards (Balcones Fault Zone) and Edwards–Trinity (Plateau) Aquifers in Kinney County, Texas; Aquifers of the Edwards Plateau, Texas Water Development Board Report 360; Texas Water Development Board: Austin, TX, USA, 2004; pp. 345–366. Available online: https://www.twdb.texas.gov/publications/reports/numbered_reports/doc/R360/Ch18.pdf (accessed on 21 March 2017).
  20. Perkins, R.B.; Vaugh, S.K.; Rodriguez, A.D. Pilot recharge models of the Edwards Aquifer, Texas, 2002. Available online: http://infohouse.p2ric.org/ref/33/32661.pdf (accessed on 21 March 2017).
  21. Joseph, J.F.; Falcon, H.E.; Sharif, H.O. Hydrologic trends and correlations in south Texas River basins: 1950–2009. J. Hydrol. Eng. 2013, 18, 1653–1662. [Google Scholar] [CrossRef]
  22. Hamilton, J.M.; Johnson, S.; Esquilin, R.; Thompson, E.L.; Wiatrek, A.; Luevano, G.; Gregory, D.; Burgoon, C.; Mireles, J.; Gloyd, R.; et al. Edwards Aquifer Authority Hydrologic Data Report for 2003; Edwards Aquifer Authority: San Antonio, TX, USA, 2004; Available online: https://www.edwardsaquifer.org/science-and-maps/research-and-scientific-reports/hydrologic-data-reports (accessed on 29 April 2017).
  23. Norris, C.W.; Moulton, D.W.; El-Hage, A.; Bradsby, D. Ecologically Significant River and Stream Segments of Region L (South Central) Regional Water Planning Area; Water Resources Technical Series—WRTS-2005-01; Texas Parks and Wildlife: Austin, TX, USA, 2005. Available online: https://tpwd.texas.gov/publications/pwdpubs/media/pwd_rp_v3400_1163.pdf (accessed on 21 March 2017).
  24. Ockerman, D.J.; Fernandez, C.J. Hydrologic Conditions and Water Quality of Precipitation and Storm Runoff for Two Agricultural Areas of the Oso Creek Watershed, Nueces County, Texas, 2005–08; United States Geological Survey Scientific Investigations Report 2010–5136; United States Geological Survey: Reston, VA, USA, 2010. Available online: http://pubs.usgs.gov/sir/2010/5136/ (accessed on 21 March 2017).
  25. Ockerman, D.J.; Heitmuller, F.T. Simulation of Streamflow and Suspended Sediment Concentrations and Loads in the Lower Nueces River Watershed, Downstream from Lake Corpus Christi to the Nueces Estuary, South Texas, 1958–2008; United States Geological Survey Scientific Investigations Report 2010–5194; United States Geological Survey: Reston, VA, USA, 2010. Available online: https://pubs.usgs.gov/sir/2010/5194/pdf/sir2010-5194.pdf (accessed on 21 March 2017).
  26. Adams, J.B. A review of methods and frameworks used to determine the environmental water requirements of estuaries. Hydrol. Sci. J. 2014, 59, 451–465. [Google Scholar] [CrossRef] [Green Version]
  27. Adams, J.S.; Tunnell, J.W. Rincon Bayou salinity monitoring project. Coastal Bend Bays & Estuaries Program Final Report - Publication No. CBBEP-66. 2010. Available online: http://cbbep.org/publications/virtuallibrary/0921report.pdf (accessed on 20 June 2017).
  28. Montagna, P.A.; Hutchinson, L.M.; Scholz, D.; Palmer, T.; Arismendez, S.; Yoskowitz, D.; Nueces Estuary Ecosystem Management Initiative. Coastal Bend Bays & Estuaries Program - Publication No. CBBEP-72. 2011. Available online: http://www.cbbep.org/publications/virtuallibrary/1018.pdf (accessed on 21 March 2017).
  29. Murgulet, D.; Murgulet, N.; Spalt, N.; Douglas, A.; Hay, R.G. Impact of hydrological alterations on river-groundwater exchange in a semi-arid area: Nueces River, Texas. Sci. Total Environ. 2016, 572, 595–607. [Google Scholar] [CrossRef] [PubMed]
  30. Gary, M.; Kromann, J. Evaluation of Riparian Water Flux Patterns in the Upper Nueces, Sabinal, Frio, and Dry Frio Rivers, Texas in Relation to the Control of Arundo Donax; Edwards Aquifer Authority Report; Edwards Aquifer Authority: San Antonio, TX, USA, 2013; Available online: http://www.pullkillplant.org/S/pdfs/Arundofinal_10-13.pdf (accessed on 29 April 2017).
  31. Ji, J.; Sinha, T. Effects of climate change on hydrology of the Nueces River Basin in South Texas. In Proceedings of the 2016 ASABE Annual International Meeting, Orlando, FL, USA, 17–20 July 2016; Available online: https://elibrary.asabe.org/azdez.asp?AID=47279&t=2 (accessed on 2 April 2017).
  32. Henningson, Durham and Richardson Engineering, Inc.; Geraghty & Miller, Inc. Regional Water Supply Study Nueces River Basin. Volume II – Technical Report; Texas Water Development Board: Austin, TX, USA, 1991. Available online: https://www.twdb.texas.gov/publications/reports/contracted_reports/doc/90483772b.pdf (accessed on 18 June 2017).
  33. Henningson, Durham and Richardson Engineering, Inc. Water Availability in the Nueces River Basin. 1999. Available online: https://twdl-ir.tdl.org/handle/10850/1360 (accessed on 18 June 2017).
  34. Lins, H.F. USGS Hydro-Climatic Data Network 2009 (HCDN-2009); US Geological Survey Fact Sheet 2012-3047; U.S. Geological Survey: Reston, VA, USA, 2012; p. 4. Available online: https://pubs.usgs.gov/fs/2012/3047/pdf/fs2012-3047.pdf (accessed on 25 June 2017).
  35. United States Geological Survey. South-Central Texas Study Unit. 2013. Available online: https://pubs.usgs.gov/of//1999/ofr99243/SCTXstudyunit.htm (accessed on 10 March 2017).
  36. Texas Water Development Board. Chapter 4 Climate of Texas. 2012. Available online: http://www.twdb.texas.gov/publications/state_water_plan/2012/04.pdf (accessed on 21 March 2017).
  37. Texas Water Development Board. Climatic Atlas of Texas. 1983. Available online: http://www.twdb.texas.gov/publications/reports/limited_printing/doc/LP192.pdf (accessed on 21 March 2017).
  38. United States Geological Survey. Current conditions for Texas: Streamflow. 2018. Available online: https://waterdata.usgs.gov/tx/nwis/current/?type=flow&group_key=basin_cd (accessed on 27 December 2018).
  39. Texas Water Development Board. Water Planning Data for Region L (South Central Texas). 2012. Available online: http://tpwd.texas.gov/landwater/water/conservation/water_resources/water_quantity/sigsegs/regionl.phtml (accessed on 25 March 2017).
  40. Burn, D.H.; Elnur, M.A.H. Detection of hydrologic trends and variability. J. Hydrol. 2002, 255, 107–122. [Google Scholar] [CrossRef]
  41. Mather, J.R. Using computed stream flow in watershed analysis. J. Am. Water Resour. Assoc. 1981, 17, 474–482. [Google Scholar] [CrossRef]
  42. United States Geological Survey. Explore Real-Time Water Data Using New Products from USGS TXWSC. 2016. Available online: https://pubs.usgs.gov/of//1999/ofr99243/SCTXstudyunit.htm (accessed on 19 June 2016).
  43. Falcone, J.A. U.S. Conterminous Wall-to-Wall Anthropogenic Land Use Trends (NWALT), 1974–2012; United States Geological Survey Data Series 948; United States Geological Survey: Reston, VA, USA, 2015; 33p. Available online: https://pubs.usgs.gov/ds/0948/ds948.pdf (accessed on 15 February 2018).
  44. National Oceanic and Atmospheric Administration. Climate Data Online Search. 2017. Available online: https://www.ncdc.noaa.gov/cdo-web/search?datasetid=GHCND (accessed on 23 June 2017).
  45. Szilagyi, J. A calibration-free, robust estimation of monthly land surface evapotranspiration rates for continental-scale hydrology. Hydrol. Res. 2017. [Google Scholar] [CrossRef]
  46. McMahon, T.A.; Peel, M.C.; Lowe, L.; Srikanthan, R.; McVicar, T.R. Estimating actual, potential, reference crop and pan evaporation using standard meteorological data: A pragmatic synthesis. Hydrol. Earth Syst. Sci. 2013, 17, 1331–1363. [Google Scholar] [CrossRef]
  47. Barua, S.; Muttil, N.; Ng, A.W.M.; Perera, B.J.C. Precipitation trend and its implications for water resource management within the Yarra River catchment, Australia. Hydrol. Process. 2013, 27, 1727–1738. [Google Scholar] [CrossRef]
  48. Gu, C.; Mu, X.; Zhao, G.; Gao, P.; Sun, W.; Yu, Q. Changes in stream flow and their relationships with climatic variations and anthropogenic activities in the Poyang Lake Basin, China. Water 2016, 8, 564. [Google Scholar] [CrossRef]
  49. Kahya, E.; Kalayci, S. Trend analysis of streamflow in Turkey. J. Hydrol. 2003, 289, 128–144. [Google Scholar] [CrossRef]
  50. Kam, J.; Sheffield, J. Changes in the low flow regime over the eastern United States (1962–2011): Variability, trends, and attributions. Clim. Change 2015, 135, 639–653. [Google Scholar] [CrossRef]
  51. Kumar, S.; Merwade, V.; Kam, J.; Thurner, K. Streamflow trends in Indiana: Effects of long term persistence, precipitation, and subsurface drains. J. Hydrol. 2009, 374, 171–183. [Google Scholar] [CrossRef]
  52. Lins, H.F.; Slack, J.R. Streamflow trends in the United States. Geophys. Res. Lett. 1999, 26, 227–230. [Google Scholar] [CrossRef] [Green Version]
  53. Novotny, E.V.; Stefan, H.G. Stream flow in Minnesota: Indicator of climate change. J. Hydrol. 2006, 334, 319–333. [Google Scholar] [CrossRef]
  54. Williams, R. Serial Correlation [Optional; Very Brief Overview]. 2015. Available online: https://www3.nd.edu/~rwilliam/stats2/l26.pdf (accessed on 21 March 2017).
  55. Hamed, K.H.; Rao, A.R. A modified Mann-Kendall trend test for autocorrelated data. J. Hydrol. 1998, 204, 182–196. [Google Scholar] [CrossRef]
  56. Mann, H.B. Nonparametric tests against trend. Econom. J. Econom. Soc. 1945, 13, 245–259. [Google Scholar] [CrossRef]
  57. Kendall, M.G. Rank Correlation Methods; Charles Griffin & Co: London, UK, 1975. [Google Scholar]
  58. Santander Meteorology Group. Package ‘Fume’. 2013. Available online: http://www2.uaem.mx/r-mirror/web/packages/fume/fume.pdf (accessed on 21 March 2017).
  59. Dierauer, J.; Whitfield, P. Package ‘FlowScreen’. 2017. Available online: https://cran.r-project.org/web/packages/FlowScreen/FlowScreen.pdf (accessed on 6 June 2017).
  60. Sankarasubramanian, A.; Vogel, R.M.; Limbrunner, J.F. Climate elasticity of streamflow in the United States. Water Resour. Res. 2001, 37, 1771–1781. [Google Scholar] [CrossRef] [Green Version]
  61. Nalbantis, I.; Tsakiris, G. Assessment of hydrological drought revisited. Water Resour. Manag. 2009, 23, 881–897. [Google Scholar] [CrossRef]
  62. Tabari, H.; Nikbakht, J.; Talaee, P.H. Hydrological drought assessment in northwestern Iran based on streamflow drought index (SDI). Water Resour. Manag. 2012, 27, 137–151. [Google Scholar] [CrossRef]
  63. Lettenmaier, D.P.; Wood, E.F.; Wallis, J.R. Hydro-climatological trends in the continental United States, 1948-88. J. Clim. 1994, 7, 586–607. [Google Scholar] [CrossRef]
  64. Sharma, R.H.; Shakya, N.M. Hydrological changes and its impact on water resources of Bagmati watershed, Nepal. J. Hydrol. 2006, 327, 315–322. [Google Scholar] [CrossRef]
  65. Rientjes, T.H.M.; Haile, A.T.; Kebede, E.; Mannaerts, C.M.M.; Habib, E.; Steenhuis, T.S. Changes in land cover, rainfall and stream flow in Upper Gilgel Abbay catchment, Blue Nile Basin – Ethiopia. Hydrol. Earth Syst. Sci. 2011, 15, 1979–1989. [Google Scholar] [CrossRef]
  66. Armson, D.; Stringer, P.; Ennos, A.R. The effect of street trees and amenity grass on urban surface water runoff in Manchester, UK. Urban For. Urban Green. 2013, 12, 282–286. [Google Scholar] [CrossRef]
  67. Jones, C.A.; Gregory, L. Effects of Brush Management on Water Resources; Texas Water Resources Institute Technical Report – 338; Texas A&M AgriLife: College Station, TX, USA, 2008; Available online: http://riogrande.tamu.edu/reports/2008/tr338.pdf (accessed on 18 February 2018).
  68. Chakilu, G.G.; Moges, M.A. Assessing the land use/cover dynamics and its impact on the low flow of Gumara Watershed, Upper Blue Nile Basin, Ethiopia. Hydrol. Curr. Res. 2017, 7, 268. [Google Scholar] [CrossRef]
  69. Kim, J.; Choi, J.; Choi, C.; Park, S. Impacts of changes in climate and land use/land cover under IPCC RCP scenarios on streamflow in the Hoeya River Basin, Korea. Sci. Total Environ. 2013, 452, 181–195. [Google Scholar] [CrossRef] [PubMed]
  70. Texas Commission on Environmental Quality. South Texas Watermaster Program. 2018. Available online: https://www.tceq.texas.gov/permitting/water_rights/wmaster/stwr (accessed on 7 February 2018).
  71. Murphy, J.; Sprague, L. Water-quality trends in US rivers: Exploring effects from streamflow trends and changes in watershed management. Sci. Total. Environ. 2019, 656, 645–658. [Google Scholar] [CrossRef]
  72. Texas Commission on Environmental Quality. Water Rights and Water use Data. 2017. Available online: https://www.tceq.texas.gov/permitting/water_rights/wr-permitting/wrwud (accessed on 7 February 2018).
  73. Texas Commission on Environmental Quality. Water Rights Permits: Annual Water Use Reporting. 2018. Available online: https://www.tceq.texas.gov/permitting/water_rights/wr-permitting/water-right-permits-annual-water-use-report (accessed on 7 February 2018).
  74. Texas Commission on Environmental Quality. Am I Regulated? Water Rights in Texas. 2017. Available online: https://www.tceq.texas.gov/permitting/water_rights/wr-permitting/wr_amiregulated.html (accessed on 7 February 2018).
  75. Texas Parks & Wildlife. Ecologically Significant Stream Segments. Available online: https://tpwd.texas.gov/landwater/water/conservation/water_resources/water_quantity/sigsegs/ (accessed on 25 February 2019).
  76. Church, J.A.; White, N.J. A 20th century acceleration in global sea-level rise. Geophys. Res. Lett. 2006, 33, 1. [Google Scholar] [CrossRef]
  77. Wigley, T.M.; Raper, S.C.B. Thermal expansion of sea water associated with global warming. Nature 1987, 330, 127. [Google Scholar] [CrossRef]
  78. Texas Commission on Environmental Quality; Texas Parks and Wildlife Development; Texas Water Development Board. Texas Instream Flow Studies: Technical Overview Report 369. 2008. Available online: http://www.twdb.texas.gov/publications/reports/numbered_reports/doc/R369_InstreamFlows.pdf (accessed on 21 March 2017).
  79. Haney, J.A.; Turner, D.S.; Springer, A.E.; Stromberg, J.C.; Stevens, L.E.; Pearthree, P.A.; Supplee, V. Ecological Implications of Verde River Flows. Verde River Basin Partnership, Arizona Water Institute, and The Nature Conservancy Report. 2008. Available online: http://azconservation.org/dl/TNCAZ_VerdeRiver_Ecological_Flows.pdf (accessed on 3 June 2017).
Figure 1. Texas location in the United States and Regional Water Planning Areas in Texas [17].
Figure 1. Texas location in the United States and Regional Water Planning Areas in Texas [17].
Hydrology 06 00020 g001
Figure 2. Study area: Nueces River Basin in TX, USA.
Figure 2. Study area: Nueces River Basin in TX, USA.
Hydrology 06 00020 g002
Figure 3. Hydroclimatic Data Network gauge locations (see Table 1 for gauge information).
Figure 3. Hydroclimatic Data Network gauge locations (see Table 1 for gauge information).
Hydrology 06 00020 g003
Figure 4. Land use changes from 1974 to 2012.
Figure 4. Land use changes from 1974 to 2012.
Hydrology 06 00020 g004
Figure 5. Streamflow variability from 1970–2014 (median values for this time period are shown using the diamond plus symbol).
Figure 5. Streamflow variability from 1970–2014 (median values for this time period are shown using the diamond plus symbol).
Hydrology 06 00020 g005
Figure 6. Streamflow from the most western (top), central (middle), and eastern (bottom) gauges in the NRB.
Figure 6. Streamflow from the most western (top), central (middle), and eastern (bottom) gauges in the NRB.
Hydrology 06 00020 g006
Figure 7. Streamflow Drought Index (Red) and Standardized Precipitation Index (Blue).
Figure 7. Streamflow Drought Index (Red) and Standardized Precipitation Index (Blue).
Hydrology 06 00020 g007
Table 1. Functions of the significant streams in the study area.
Table 1. Functions of the significant streams in the study area.
IdentifierGauge IDStreamEcologically Significant Streams Criteria
BFHSRCHQTS
A08190500West Nueces River
B08190000Nueces River
C08195000Frio River
D08198500Sabinal River
1 BC – Biological function; HC – Hydrologic function; RC – Riparian conservation areas; HQ – High water quality, exceptional aquatic life, and high aesthetic value; TS – Threatened or endangered species, unique communities.
Table 2. Consolidated precipitation stations.
Table 2. Consolidated precipitation stations.
CountySubstituted forStation NameStation IDLatitudeLongitudeCorresponding Gauge
ZavalaKinneyLa Pryor TX USGHCND:USC0041492028.950000−99.83333008190500
Crystal City TX USGHCND:USC0041216028.683330−99.833330
Chaparrosa Ranch TX USGHCND:US1TXZV000128.885800−99.996900
McNally TX USGHCND:US1TXZV000828.839700−99.924200
UvaldeN/ASabinal TX USGHCND:USC0041787329.333330−99.48333008190000
Uvalde Tx USGHCND:USC0041926529.216670−99.76667008195000
Utopia Tx USGHCND:USC0041926029.616670−99.51667008196000
Uvalde 3 SW Tx USGHCND:USC0041926829.216670−99.75000008198500
MedinaN/ARiomedina TX USGHCND:USC0041762829.466670−98.86667008200000
Lytle 3 W USGHCND:USC0041545429.233330−98.800000
Hondo TX USGHCND:USC0041425429.336500−99.138600
Natalia TX USGHCND:USC0041620529.200000−98.866670
Table 3. Annual and seasonal NRB streamflow trends from 1970 to 2014 ([NT] no significant trend; [↓] statistically significant decreasing trend).
Table 3. Annual and seasonal NRB streamflow trends from 1970 to 2014 ([NT] no significant trend; [↓] statistically significant decreasing trend).
Qmin
Gauge IDAnnualSpringSummerFallWinter
08190500NTNTNTNTNT
08190000NTNTNT
08196000NT
08195000NT
08198500NTNT
08200000NTNT
Qmed
Gauge IDAnnualSpringSummerFallWinter
08190500NTNTNTNTNT
08190000NTNTNT
08196000NT
08195000NT
08198500NTNTNT
08200000NTNT
Qmax
Gauge IDAnnualSpringSummerFallWinter
08190500NTNTNTNTNT
08190000NTNTNTNTNT
08196000NTNTNTNTNT
08195000NTNTNT
08198500NTNTNTNTNT
08200000NTNTNTNTNT
Table 4. Annual and seasonal NRB streamflow trends from 1994 to 2014 ([NT] no significant trend; [↓] statistically significant decreasing trend).
Table 4. Annual and seasonal NRB streamflow trends from 1994 to 2014 ([NT] no significant trend; [↓] statistically significant decreasing trend).
Gauge IDQmin
AnnualSpringSummerFallWinter
08190500NTNTNT
08190000NTNTNTNT
08196000NTNTNTNTNT
08195000NTNT
08198500NTNTNT
08200000NTNTNTNT
Gauge IDQmed
AnnualSpringSummerFallWinter
08190500NTNTNT
08190000NTNTNT
08196000NTNT
08195000NTNT
08198500NTNTNTNT
08200000NTNTNTNT
Gauge IDQmax
AnnualSpringSummerFallWinter
08190500NTNT
08190000NTNT
08196000NTNTNT
08195000NTNTNT
08198500NTNTNTNT
08200000NTNTNTNTNT
Table 5. Streamflow precipitation elasticity for each gauge.
Table 5. Streamflow precipitation elasticity for each gauge.
Gauge IDGauge NameElasticityCounties
08190500W Nueces Rv nr Brackettville, TX0.99Zavala
08190000Nueces Rv at Laguna, TX0.47Uvalde
08196000Dry Frio Rv nr Reagan Wells, TX0.60
08195000Frio Rv at Concan, TX0.41
08198500Sabinal Rv at Sabinal, TX0.95
08200000Hondo Ck nr Tarpley, TX0.72Medina
Table 6. Upper NRB change in actual ET from 1979 to 2015 adapted from Szilagyi (2017).
Table 6. Upper NRB change in actual ET from 1979 to 2015 adapted from Szilagyi (2017).
Gauge IDGauge NameΔ Actual ET (in/year)Counties
08190500W Nueces Rv nr Brackettville, TX0.048Zavala
08190000Nueces Rv at Laguna, TX0.030Uvalde
08196000Dry Frio Rv nr Reagan Wells, TX0.003
08195000Frio Rv at Concan, TX0.001
08198500Sabinal Rv at Sabinal, TX−0.005
08200000Hondo Ck nr Tarpley, TX−0.024Medina

Share and Cite

MDPI and ACS Style

Thomas, E.D.; Venkataraman, K.; Chraibi, V.; Kannan, N. Hydrologic Trends in the Upper Nueces River Basin of Texas—Implications for Water Resource Management and Ecological Health. Hydrology 2019, 6, 20. https://doi.org/10.3390/hydrology6010020

AMA Style

Thomas ED, Venkataraman K, Chraibi V, Kannan N. Hydrologic Trends in the Upper Nueces River Basin of Texas—Implications for Water Resource Management and Ecological Health. Hydrology. 2019; 6(1):20. https://doi.org/10.3390/hydrology6010020

Chicago/Turabian Style

Thomas, E. Dave, Kartik Venkataraman, Victoria Chraibi, and Narayanan Kannan. 2019. "Hydrologic Trends in the Upper Nueces River Basin of Texas—Implications for Water Resource Management and Ecological Health" Hydrology 6, no. 1: 20. https://doi.org/10.3390/hydrology6010020

APA Style

Thomas, E. D., Venkataraman, K., Chraibi, V., & Kannan, N. (2019). Hydrologic Trends in the Upper Nueces River Basin of Texas—Implications for Water Resource Management and Ecological Health. Hydrology, 6(1), 20. https://doi.org/10.3390/hydrology6010020

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop