Utility of Gridded Rainfall for IHACRES Daily River Flow Predictions in Southern California Watersheds1

Abstract: Studies to regionalize conceptual hydrologic models generally require rainfall and river flow data from multiple watersheds. Besides the considerable time (cost) to assemble and process rainfall data for many watersheds, investigators often need to choose from a number of candidate gauges, subjectively weighing the relative importance of proximity and elevation to select a representative rainfall dataset. The Unified Raingauge Dataset (URD) is a gridded daily rainfall dataset that covers the conterminous United States at 0.25 × 0.25 degrees spatial resolution and is available from 1948 to present. The objective of this study was to determine whether uncertainty in daily river flow predictions using the conceptual hydrologic model IHACRES in small to moderate size watersheds (50‐400 km²) in southern California would increase if URD gridded rainfall data were used in place of single rain gauge data to calibrate the model. Rain gauge data were obtained from the gauge nearest the watershed centroid and the gauge closest in elevation to the watershed mean elevation. Results from 20 randomly selected watersheds indicated that IHACRES calibration performance was similar using rainfall data from the URD grids and rain gauge data. There was some evidence of greater uncertainties associated with the URD calibrations in areas where topography may affect rainfall amounts. In contrast to the URD data, monthly gridded data produced by the Parameter‐Elevation Regressions on Independent Slopes Model (PRISM) includes adjustments for elevation and produces gridded values at a finer spatial resolution (4 km²). A limited test on two watersheds demonstrated that scaling the URD daily rainfall estimates to match the PRISM monthly values may improve rainfall estimates and model simulation performance.     

Assessment of the TRMM 3B42 Precipitation Product in Southern Brazil

The increase in the use of satellite‐derived precipitation products generated by different methods and algorithms emphasizes the need for a deeper analysis of their quality and accuracy. Using the contingency table method, we evaluated the accuracy of versions 6 and 7 of the Tropical Rainfall Measuring Mission Precipitation (TRMM) 3B42 product in southern Brazil by comparing daily precipitation over 13 years (V6 was tested for historical context). The interpolated data from 25 rain gauges were compared with both versions of TRMM. The V7 product tended to produce a slight increase in PC (proportion correct). V7 also showed a slight increase in the correlation coefficient (CC) and a significant increase in the H (hit rate) and CSI (critical success) indexes. However, the upgraded version shows an undesirable increase in the false alarm ratio. When the rainfall volumes were compared, V6 clearly underestimated the total rainfall over the entire period, but the V7 product slightly overestimated the cumulative volume (11%) which still represented a more reliable estimate than from V6. Furthermore, the main improvement in V7 was a large increase in the quantitative recognition of extreme precipitation events: V6 detected only 1% of the daily rainfalls above 60 mm, whereas V7 detected 57% of the events.     

Using SWAT to Model Streamflow in Two River Basins With Ground and Satellite Precipitation Data1

Abstract: Both ground rain gauge and remotely sensed precipitation (Next Generation Weather Radar – NEXRAD Stage III) data have been used to support spatially distributed hydrological modeling. This study is unique in that it utilizes and compares the performance of National Weather Service (NWS) rain gauge, NEXRAD Stage III, and Tropical Rainfall Measurement Mission (TRMM) 3B42 (Version 6) data for the hydrological modeling of the Middle Nueces River Watershed in South Texas and Middle Rio Grande Watershed in South Texas and northern Mexico. The hydrologic model chosen for this study is the Soil and Water Assessment Tool (SWAT), which is a comprehensive, physical‐based tool that models watershed hydrology and water quality within stream reaches. Minor adjustments to selected model parameters were applied to make parameter values more realistic based on results from previous studies. In both watersheds, NEXRAD Stage III data yields results with low mass balance error between simulated and actual streamflow (±13%) and high monthly Nash‐Sutcliffe efficiency coefficients (NS > 0.60) for both calibration (July 1, 2003 to December 31, 2006) and validation (2007) periods. In the Middle Rio Grande Watershed NEXRAD Stage III data also yield robust daily results (time averaged over a three‐day period) with NS values of (0.60‐0.88). TRMM 3B42 data generate simulations for the Middle Rio Grande Watershed of variable qualtiy (MBE = +13 to ?16%; NS = 0.38‐0.94; RMSE = 0.07‐0.65), but greatly overestimates streamflow during the calibration period in the Middle Nueces Watershed. During the calibration period use of NWS rain gauge data does not generate acceptable simulations in both watersheds. Significantly, our study is the first to successfully demonstrate the utility of satellite‐estimated precipitation (TRMM 3B42) in supporting hydrologic modeling with SWAT; thereby, potentially extending the realm (between 50°N and 50°S) where remotely sensed precipitation data can support hydrologic modeling outside of regions that have modern, ground‐based radar networks (i.e., much of the third world).     

Evaluation of Gauge‐Radar Merging Methods Using a Semi‐Distributed Hydrological Model in the Upper Thames River Basin,Canada

Gauge‐radar merging methods combine rainfall estimates from rain gauges and radar to capitalize on the strengths of the individual instruments. The performance of four well‐known gauge‐radar merging methods, including mean field bias correction, Brandes spatial adjustment, local bias correction using kriging, and conditional merging, are examined using Environment Canada radar and the Upper Thames River Basin in southwestern Ontario, Canada, as a case study. The analysis assesses the effect of gauge‐radar merging methods on: (1) the accuracy of predicted rainfall accumulations; and (2) the accuracy of predicted streamflows using a semi‐distributed hydrological model. In addition, several influencing factors (i.e., gauge density, storm type, basin type, proximity to the radar tower, and time‐step of adjustment) are analyzed to determine their effect on the performance of the rainfall estimation techniques. Confirming results of previous studies, the merging methods provide an increase in the accuracy of both rainfall accumulation estimations and predicted streamflows. The results also indicate specific factors such as gauge density, rainfall intensity, and time‐step of adjustment can reduce the accuracy of merging methods and play a key role in the examination of its use for operational purposes. Results provide guidance for hydrologists and engineers assessing how best to apply corrected radar products to improve rainfall estimation and hydrological modeling accuracy.     

Physically Based,Hydrologic Model Results Based on Three Precipitation Products1

Abstract: The main objective of the study is to examine the accuracy of and differences among simulated streamflows driven by rainfall estimates from a network of 22 rain gauges spread over a 2,170 km² watershed, NEXRAD Stage III radar data, and Tropical Rainfall Measuring Mission (TRMM) 3B42 satellite data. The Gridded Surface Subsurface Hydrologic Analysis (GSSHA), a physically based, distributed parameter, grid‐structured, hydrologic model, was used to simulate the June‐2002 flooding event in the Upper Guadalupe River watershed in south central Texas. There were significant differences between the rainfall fields estimated by the three types of measurement technologies. These differences resulted in even larger differences in the simulated hydrologic response of the watershed. In general, simulations driven by radar rainfall yielded better results than those driven by satellite or rain‐gauge estimates. This study also presents an overview of effects of land cover changes on runoff and stream discharge. The results demonstrate that, for major rainfall events similar to the 2002 event, the effect of urbanization on the watershed in the past two decades would not have made any significant effect on the hydrologic response. The effect of urbanization on the hydrologic response increases as the size of the rainfall event decreases.     

Long‐Term High‐Resolution Radar Rainfall Fields for Urban Hydrology

Accurate records of high‐resolution rainfall fields are essential in urban hydrology, and are lacking in many areas. We develop a high‐resolution (15 min, 1 km²) radar rainfall data set for Charlotte, North Carolina during the 2001‐2010 period using the Hydro‐NEXRAD system with radar reflectivity from the National Weather Service Weather Surveillance Radar 1988 Doppler weather radar located in Greer, South Carolina. A dense network of 71 rain gages is used for estimating and correcting radar rainfall biases. Radar rainfall estimates with daily mean field bias (MFB) correction accurately capture the spatial and temporal structure of extreme rainfall, but bias correction at finer timescales can improve cold‐season and tropical cyclone rainfall estimates. Approximately 25 rain gages are sufficient to estimate daily MFB over an area of at least 2,500 km², suggesting that robust bias correction is feasible in many urban areas. Conditional (rain‐rate dependent) bias can be removed, but at the expense of other performance criteria such as mean square error. Hydro‐NEXRAD radar rainfall estimates are also compared with the coarser resolution (hourly, 16 km²) Stage IV operational rainfall product. Stage IV is adequate for flood water balance studies but is insufficient for applications such as urban flood modeling, in which the temporal and spatial scales of relevant hydrologic processes are short. We recommend the increased use of high‐resolution radar rainfall fields in urban hydrology.     

EVALUATION OF NEXRAD STAGE III PRECIPITATION DATA OVER A SEMIARID REGION1

This study examines NEXRAD Stage III product (hourly, cell size 4 km by 4 km) for its ability in estimating precipitation in central New Mexico, a semiarid area. A comparison between Stage III and a network of gauge precipitation estimates during 1995 to 2001 indicates that Stage III (1) overestimates the hourly conditional mean (CM) precipitation by 33 percent in the monsoon season and 55 percent in the nonmonsoon season; (2) overestimates the hourly CM precipitation for concurrent radar‐gauge pairs (nonzero value) by 13 percent in the monsoon season and 6 percent in the nonmonsoon season; (3) overestimates the seasonal precipitation accumulation by 11 to 88 percent in monsoon season and underestimates by 18 to 89 percent in the nonmonsoon season; and (4) either overestimates annual precipitation accumulation up to 28.2 percent or underestimates it up to 11.9 percent. A truncation of 57 to 72 percent of the total rainfall hours is observed in the Stage III data in the nonmonsoon season, which may be the main cause for both the underestimation of the radar rainfall accumulation and the lower conditional probability of radar rainfall detection in the nonmonsoon season. The study results indicate that the truncation caused loss of small rainfall amounts (events) is not effectively corrected by the real‐time rain gauge calibration that can adjust the rainfall rates but cannot recover the truncated small rainfall events. However, the truncation error in the monsoon season may be suppressed due to the larger rainfall rate and/or combined effect of overestimates by bright band and hail contaminations, virga, advection, etc. In general, improvement in NEXRAD performance since the monsoon season in 1998 is observed, which is consistent with the systematic improvement in the NEXRAD network.     

Rainfall Interception Based on Indirect Methods: A Case Study in Temperate Forests in Oaxaca,Mexico

Rainfall interception represents the amount of water trapped in natural cover that is not drained directly to the ground. Intercepted rainfall may evaporate after a rain event, making it one of the main drivers of water balance and hydrologic regionalization. This process can be affected by factors such as climate, altitude, vegetation type, and topography. Here is a simple method of calculating rainfall interception in temperate forests using in Santa Maria Yavesia, Oaxaca, and Mexico as an illustrative study area. We used two rain gauges to measure net precipitation (N_p) under the canopy at each study site and one gauge outside the canopy to obtain gross precipitation (G_p). Throughfall (T_h) was indirectly measured using hemispherical photographs. Rainfall interception was obtained through a combination T_h and G_p and N_p. The mean rainfall interception was 50.6% in the Abies forests, 23%–40% in the coniferous‐mixed forests, and 27.4% in the broad‐leaved forests. We classified rainfall events by intensity to determine the effect of canopy structure and precipitation and found that 75% of the events were weak events, 24% were moderate events, and 1% were strong events. In addition, we found that rainfall interception was lower when the intensity of precipitation was higher. Our method can be replicated in different ecosystems worldwide as a tool for assessing the influence of rainfall interception in terms of ecological services.     

Water Resources Modeling of the Ganges‐Brahmaputra‐Meghna River Basins Using Satellite Remote Sensing Data1

Nishat, Bushra and S.M. Mahbubur Rahman, 2009. Water Resources Modeling of the Ganges‐Brahmaputra‐Meghna River Basins Using Satellite Remote Sensing Data. Journal of the American Water Resources Association (JAWRA) 45(6):1313‐1327. Abstract: Large‐scale water resources modeling can provide useful insights on future water availability scenarios for downstream nations in anticipation of proposed upstream water resources projects in large international river basins (IRBs). However, model set up can be challenging due to the large amounts of data requirement on both static states (soils, vegetation, topography, drainage network, etc.) and dynamic variables (rainfall, streamflow, soil moisture, evapotranspiration, etc.) over the basin from multiple nations and data collection agencies. Under such circumstances, satellite remote sensing provides a more pragmatic and convenient alternative because of the vantage of space and easy availability from a single data platform. In this paper, we demonstrate a modeling effort to set up a water resources management model, MIKE BASIN, over the Ganges, Brahmaputra, and Meghna (GBM) river basins. The model is set up with the objective of providing Bangladesh, the lowermost riparian nation in the GBM basins, a framework for assessing proposed water diversion scenarios in the upstream transboundary regions of India and deriving quantitative impacts on water availability. Using an array of satellite remote sensing data on topography, vegetation, and rainfall from the transboundary regions, we demonstrate that it is possible to calibrate MIKE BASIN to a satisfactory level and predict streamflow in the Ganges and Brahmaputra rivers at the entry points of Bangladesh at relevant scales of water resources management. Simulated runoff for the Ganges and Brahmaputra rivers follow the trends in the rated discharge for the calibration period. However, monthly flow volume differs from the actual rated flow by (?) 8% to (+) 20% in the Ganges basin, by (?) 15 to (+) 12% in the Brahmaputra basin, and by (?) 15 to (+) 19% in the Meghna basin. Our large‐scale modeling initiative is generic enough for other downstream nations in IRBs to adopt for their own modeling needs.     

Satellite Rainfall Estimates Over South America – Possible Applicability to the Water Management of Large Watersheds1

Pereira Filho, Augusto J., Richard E. Carbone, John E. Janowiak, Phillip Arkin, Robert Joyce, Ricardo Hallak, and Camila G.M. Ramos, 2010. Satellite Rainfall Estimates Over South America – Possible Applicability to the Water Management of Large Watersheds. Journal of the American Water Resources Association (JAWRA) 46(2):344-360. DOI: 10.1111/j.1752-1688.2009.00406.x Abstract: This work analyzes high-resolution precipitation data from satellite-derived rainfall estimates over South America, especially over the Amazon Basin. The goal is to examine whether satellite-derived precipitation estimates can be used in hydrology and in the management of larger watersheds of South America. High spatial-temporal resolution precipitation estimates obtained with the CMORPH method serve this purpose while providing an additional hydrometeorological perspective on the convective regime over South America and its predictability. CMORPH rainfall estimates at 8-km spatial resolution for 2003 and 2004 were compared with available rain gauge measurements at daily, monthly, and yearly accumulation time scales. The results show the correlation between satellite-derived and gauge-measured precipitation increases with accumulation period from daily to monthly, especially during the rainy season. Time-longitude diagrams of CMORPH hourly rainfall show the genesis, strength, longevity, and phase speed of convective systems. Hourly rainfall analyses indicate that convection over the Amazon region is often more organized than previously thought, thus inferring that basin scale predictions of rainfall for hydrological and water management purposes have the potential to become more skillful. Flow estimates based on CMORPH and the rain gauge network are compared to long-term observed average flow. The results suggest this satellite-based rainfall estimation technique has considerable utility. Other statistics for monthly accumulations also suggest CMORPH can be an important source of rainfall information at smaller spatial scales where in situ observations are lacking.     

Validation of Satellite Precipitation Adjustment Methodology From Seven Basins in the Continental United States1

Tobin, Kenneth J. and Marvin E. Bennett, 2012. Validation of Satellite Precipitation Adjustment Methodology From Seven Basins in the Continental United States. Journal of the American Water Resources Association (JAWRA) 48(2): 221‐234. DOI: 10.1111/j.1752‐1688.2011.00604.x Abstract: The precipitation science community has expressed concern regarding the ability of satellite‐based precipitation products to accurately capture rainfall values over land. There has been some work that has focused on addressing the deficiencies of satellite precipitation products, particularly on the adjustment of bias. This article outlines a methodology that adjusts satellite products utilizing ground‐based precipitation data. The approach is not a simple bias adjustment, but is a three‐step process that transforms a satellite product based on a ground‐based precipitation product (NEXRAD‐derived Multisensor Precipitation Estimator [MPE] product or rain‐gauge data). The developed methodology was successfully applied to seven moderate‐to‐large sized watersheds from continental United States (CONUS) and northern Mexico over a spectrum of climatic regimes ranging from dry to humid settings. Methodology validation is based on comparison of observed and simulated streamflow generated with SWAT (Soil and Water Assessment Tool) model using unadjusted and adjusted precipitation products as input. Streamflow comparison is based on mass balance error and Nash‐Sutcliffe efficiency coefficient. Finally, the contribution of how adjustment to correct misses, false alarms, and bias impacts adjusted datasets and the potential impact that the adjustment methodology can have on hydrological applications such as water resource monitoring and flood prediction are explored.     

ACCOUNTING FOR DISCREPANCIES IN PAN EVAPORATION CALCULATIONS1

ABSTRACT: Four experiments were made to document and account for differences in evaporation data that were calculated using pans equipped with float-activated recorders and pans with hook gauge/rain gauge instrumentation. Paired in-pan comparisons indicated that evaporation differences were not due to the technique of measuring water level within the pan. Also, the recorder float-lag did not account for the differences. By sampling rainfall events, it was found that evaporation pans and standard (8 in. orifice) rain gauges record significantly different amounts of rain, which results in differences in calculated evaporation on rainy days. Monitoring networks with evaporation pans should have uniform instrumentation that accurately records rainfall into the pans for consistent results.     

SPATIAL ANALYSIS FOR MONTHLY RAINFALL IN SOUTH FLORIDA1

ABSTRACT: Several methods have been developed to interpolate point rainfall data and integrate areal rainfall data from any network of stations. From previous studies, it can be concluded that models for spatial analysis of rainfall are dependent on topography, area of analysis, type of rainfall, and density of gauging network. The purpose of this study is to evaluate a set of six appropriate models for point and areal rainfall estimations over a 4000 square mile area in South Florida. In this study, a case of developing spatial continuity model for monthly rainfall from a database that had various lengths of records and missing data is documented. The spatial correlation and variogram models for monthly rainfall were developed. Six methods of spatial interpolation were applied and the results validated with historical observations. The results of the study indicate that the multiquadric, kriging, and optimal interpolation schemes are the best three methods for interpolation of monthly rainfall within the study area. The optimal and kriging methods have the advantage of providing estimates of the error of interpolation. The optimal interpolation method uses the spatial correlation function and the kriging method uses the variogram function. The two spatial functions are related. Either of the two methods provide good estimates of monthly point and areal rainfall in the study area.     

El Niño‐Southern Oscillation,rainfall variability and sustainable agricultural development in the Ho Municipality,Ghana

El Niño‐Southern Oscillation (ENSO), which occurs in the Equatorial Pacific Ocean, has been identified to have significant influence on rainfall variability throughout the world, especially in the tropics. Such variability in rainfall has implications for agrarian economies, such as that in Ghana. This study therefore sought to demonstrate the effect of ENSO‐induced variability in annual and seasonal rainfall on the development of sustainable agriculture in the Ho Municipality of Ghana. Using 61 years of monthly rainfall data (1955–2015) for the Ho Municipality and ENSO indices, this study showed that 15% of the variability in total annual rainfall is explained by the ENSO phenomena. Mean annual rainfall and rainfall in the major rainy season decreased for El Niño years, in addition to a more variable rainfall compared to that received in La Niña years. The major growing season was observed to be longer in La Niña years and shorter in El Niño years. This means that the potential for crop cultivation will be severely hampered in an El Niño year. Farmers within the municipality are therefore encouraged to harness other complementary water sources for farming activities and also employ water management strategies during El Niño years.     

Conservation strategies for effective land management of protected areas using an erosion prediction information system (EPIS)

This research demonstrates the predictive modeling capabilities of a geographic information system (GIS)-based soil erosion potential model to assess the effects of implementing land use change within a tropical watershed. The Revised Universal Soil Loss Equation (RUSLE) was integrated with a GIS to produce an Erosion Prediction Information System (EPIS) and modified to reflect conditions found in the mountainous tropics. Research was conducted in the Zenzontia subcatchment of the Río Ayuquíla, located within the Sierra de Manantlán Biosphere Reserve (SMBR), México. Expanding agricultural activities within this area will accentuate the already high rate of soil erosion and resultant sediment loading occurring in the Río Ayuquíla. Two land-use change scenarios are modeled with the EPIS: (1) implementation of soil conservation practices in erosion prone locations; and (2) selection of sites for agricultural expansion which minimize potential soil loss. Confronted with limited financial resources and the necessity for expedient action, managers of the SMBR can draw upon the predictive capacity of the EPIS to facilitate rapid and informed land-use planning decisions.     

Modes of Variability of Annual and Seasonal Rainfall in Mexico

The aim of this study is to identify temporal and spatial variability patterns of annual and seasonal rainfall in Mexico. A set of 769 weather stations located in Mexico was examined. The country was divided into 12 homogeneous rainfall regions via principal component analysis. A Pettitt test was conducted to perform a homogeneity analysis for detecting abrupt changes in mean rainfall levels, and a Mann‐Kendall test was conducted to examine the presence of monotonically increasing/decreasing patterns in the data. In total, 14.4% of the annual series was deemed nonstationary. Fourteen percent of the samples were nonstationary in the winter and summer, and 9% were nonstationary in the spring and autumn. According to the results, the nonstationarity of some seasonal rainfall series may be associated with the presence of atmospheric phenomena (e.g., El Niño/Southern Oscillation, Pacific Decadal Oscillation, Atlantic Multidecadal Oscillation, and North Atlantic Oscillation). A rainfall frequency analysis was performed for the nonstationary annual series, and significant differences in the return levels can be expected for the scenarios analyzed. The identification of areas that are more susceptible to changes in rainfall levels will improve water resource management plans in the country, and it is expected that these plans will take into account nonstationary theory.     

TEMPORAL AND SPATIAL CHARACTERIZATION OF RAINFALL OVER CENTRAL AND SOUTH FLORIDA1

requency evaluation and spatial characterization of rainfall in Central and South Florida are presented. Point frequency analysis performed at all available sites has shown that the 2‐parameter Gamma probability density function is the best model for monthly rainfall frequency over Central and South Florida. The model's parameters estimated at 145 stations were used to provide monthly rainfall estimates for 10‐ and 100‐year dry and wet return periods. Experimental and theoretical variograms computed for these estimates, as well as the Kriging estimation variance maps, show that the existing rain gage network is less capable of resolving monthly rainfall variation in the wet season than the dry season. May is the dry‐to‐wet transition month, while October is the wet‐to‐dry transition month with average rainfall of 4.5 inches. Monthly average rainfall is above 7 inches during the wet season and below 3 inches during the dry season. Two‐thirds of the annual rainfall is accumulated in the wet season. Annual average rainfall is maximum (above 60 inches) in many areas along the east coast, and is minimum (below 45 inches) in many areas over Lake Okee‐chobee and Central Florida. Rainfall maps show a changing pattern between the wet and the dry seasons. Frontal rainfall occurs in the dry season, while convective rainfall, tropical depression, and hurricanes occur in the wet season. Average rainfall is higher along the east coast area in the dry season and it is higher along the west coast area in the wet season.     

Estimating Lag to Peak between Rainfall and Peak Streamflow with a Mixed‐Effects Model

We test the use of a mixed‐effects model for estimating lag to peak for small basins in Maine (drainage areas from 0.8 to 78 km²). Lag to peak is defined as the time between the center of volume of the excess rainfall during a storm event and the resulting peak streamflow. A mixed‐effects model allows for multiple observations at sites without violating model assumptions inherent in traditional ordinary least squares models, which assume each observation is independent. The mixed model includes basin drainage area and maximum 15‐min rainfall depth for individual storms as explanatory features. Based on a remove‐one‐site cross‐validation analysis, the prediction errors of this model ranged from ?42% to +73%. The mixed model substantially outperformed three published models for lag to peak and one published model for centroid lag for estimating lag to peak for small basins in Maine. Lag to peak estimates are a key input to rainfall–runoff models used to design hydraulic infrastructure. The improved accuracy and consistency with model assumptions indicates that mixed models may provide increased data utilization that could enhance models and estimates of lag to peak in other regions.     

Change in Fish Fauna as Indication of Aquatic Ecosystem Condition in Río Grande de Morelia–Lago de Cuitzeo Basin, Mexico

Chirostoma charari and C. compressum, and they are presumed extinct. Twelve (63%) of the remaining species had declines in distribution. Sixteen (80%) of the 20 localities sampled had lost species. The greatest declines occurred in Lago de Cuitzeo proper and in the lower portion of the Río Grande de Morelia watershed. Species losses from the lake were attributable to drying and hypereutrophication of the lake because of substantial reductions in the amount and quality of tributary inputs, whereas losses from the Río Grande de Morelia watershed were the result of pollution from agricultural, municipal, and industrial sources, especially in the region around the city of Morelia. Three localities in the upper portion of the Río Grande de Morelia watershed—Cointzio reservoir, La Mintzita spring, and Insurgente Morelos stream—contained most of the remaining fish species diversity in the basin and deserve additional protection. Fish faunal changes indicated major declines in the health of aquatic ecosystems in the Morelia–Cuitzeo basin.