Search

Small Island Developing States (SIDS) are acknowledged as particularly vulnerable to extreme climate events; however, the realities for transport infrastructure and bridges are still poorly studied. Assessing bridges in this context can be challenging due to data scarcity, a lack of local standards, and uncertainty due to climate change. While bridges are designed to connect transport networks, they also carry energy, water, and communication networks, making them critical cascading failure points worthy of special attention in terms of risk assessment and resilience measures. We explore what resilience actually means for the design and management of SIDS bridge infrastructure by applying a post disaster forensics and systems approach that is not reliant on complex methods or large amounts of data. To demonstrate the practicality of our approach, we apply it to the island of Dominica, which is regularly impacted by both tropical storms and hurricanes. Our results document the extreme conditions for infrastructure and nearby settlements and the complex interrelated physical processes that occur during these events. We reflect on the implications for design approaches for bridges under these conditions and detail specific recommendations on how the resilience of existing and new bridges can be enhanced through practical measures that are achievable, even within the constraints experienced by those managing bridge infrastructure in SIDS contexts. This work adds to the growing number of studies exploring forensic disaster investigation and systems thinking, but is the first to explore bridge resilience in SIDS.

Understanding land use/land cover (LULC) effects on tropical soil infiltration is crucial for maximizing watershed scale hydro‐ecosystem services and informing land managers. This paper reports results from a multiyear investigation of LULC effects on soil bulk infiltration in steep, humid tropical, and lowland catchments. A rainfall simulator applied water at measured rates on 2 × 6 m plots producing infiltration through structured, granulated, and macroporous Ferralsols in Panama's central lowlands. Time‐lapse electrical resistivity tomography (ERT) helped to visualize infiltration depth and bulk velocity. A space‐for‐time substitution methodology allowed a land‐use history investigation by considering the following: (a) a continuously heavy‐grazed cattle pasture, (b) a rotationally grazed traditional cattle pasture, (c) a 4‐year‐old (y.o.) silvopastoral system with nonnative improved pasture grasses and managed intensive rotational grazing, (d) a 7 y.o. teak (Tectona grandis) plantation, (e) an approximately 10 y.o. secondary succession forest, (f) a 12 y.o. coffee plantation (Coffea canephora), (g) an approximately 30 y.o. secondary succession forest, and (h) a >100 y.o. secondary succession forest. Within a land cover, unique plot sites totalled two at (a), (c), (d), (e), and (g); three at (b); and one at (f) and (h). Our observations confirmed measured infiltration scale dependency by comparing our 12 m² plot‐scale measurements against 8.9 cm diameter core‐scale measurements collected by others from nearby sites. Preferential flow pathways (PFPs) significantly increased soil infiltration capacity, particularly in forests greater than or equal to 10 y.o. Time‐lapse ERT observations revealed shallower rapid bulk infiltration and increased rapid lateral subsurface flow in pasture land covers when compared with forest land covers and highlighted how much subsurface flow pathways can vary within the Ferralsol soil class. Results suggest that LULC effects on PFPs are the dominant mechanism by which LULC affects throughfall partitioning, runoff generation, and flow pathways.

Conventional markets can underprovide ecosystem services. Deliberate creation of a market for ecosystem services [e.g., a payments for ecosystem services (PES) scheme] can close the gap. The new ecosystem service market alters behaviors and quantities of ecosystem service provided and reveals prices for the ecosystems service: a market-clearing equilibrium. Assessing the potential for PES programs, which often act as ecological infrastructure investment mechanisms, requires forecasting the market-clearing equilibrium. Forecasting the equilibrium is complicated, especially at relevant social and ecological scales. It requires greater disciplinary integration than valuing ecosystem services or computing the marginal cost of making a land-use change to produce a service. We conduct anex antebenefit–cost assessment and forecast market-clearing prices and quantities for ecological infrastructure investment contracts in the Panama Canal Watershed. The Panama Canal Authority could offer contracts to private farmers to change land use to increase dry-season water flow and reduce sedimentation. A feasible voluntary contracting system yields a small program of about 1,840 ha of land conversion in a 279,000-ha watershed and generates a 4.9 benefit–cost ratio. Physical and social constraints limit market supply and scalability. Service delays, caused by lags between the time payments must be made and the time services stemming from ecosystem change are realized, hinder program feasibility. Targeting opportunities raise the benefit–cost ratio but reduce the hectares likely to be converted. We compare and contrast our results with prior state-of-the-art assessments on this system.

Cities in the humid- and seasonal-humid tropics depend on small watersheds for their water supplies. Under normal conditions with ample rainfall, water supplies are reliable. However, water shortages can occur during extended dry periods. The literature contains contradictory findings regarding the effectiveness of different land management strategies aimed at enhanced delivery of hydrologic ecosystem services during periods of significant rainfall deficit, so-called “green infrastructure”. Recent research results from field and modeling studies in the Panama Canal Watershed indicated that land-cover dependent flow paths play an important role in partitioning throughfall into subsurface stormflow and groundwater recharge. Land management practices considered included continuous and rotational grazing, silvipastoral treatments, and different ages of secondary succession including old regrowth forest. The effectiveness of land management was found to depend on both land use practice and annual rainfall as determined by orography in steep regions. These dependencies at least partially explain some of the discrepancies in the literature.

The flow of water in partially saturated porous media is of importance in fields such as hydrology, agriculture, environment and waste management. It is also one of the most complex flows in nature. The Richards’ equation describes the flow of water in an unsaturated porous medium due to the actions of gravity and capillarity neglecting the flow of the non-wetting phase, usually air. Analytical solutions of Richards’ equation exist only for simplified cases, so most practical situations require a numerical solution in one- two- or three-dimensions, depending on the problem and complexity of the flow situation. Despite the fact that the first reasonably complete conservative numerical solution method was published in the early 1990s, the numerical solution of the Richards’ equation remains computationally expensive and in certain circumstances, unreliable. A universally robust and accurate solution methodology has not yet been identified that is applicable across the range of soils, initial and boundary conditions found in practice. Existing solution codes have been modified over years to attempt to increase robustness. Despite theoretical results on the existence of solutions given sufficiently regular data and constitutive relations, our numerical methods often fail to demonstrate reliable convergence behavior in practice, especially for higher-order methods. Because of robustness, the lack of higher-order accuracy and computational expense, alternative solution approaches or methods are needed. There is also a need for better documentation of improved solution methodologies and benchmark test problems to facilitate consistent advances and avoid re-inventing of the wheel.

Bartlett et al. [2016] performed a re-interpretation and modification of the space-time lumped USDA NRCS (formerly SCS) Curve Number (CN) method to extend its applicability to forested watersheds. We believe that the well documented limitations of the CN method severely constrains the applicability of the modifications proposed by Bartlett et al. [2016]. This forward-looking comment urges the research communities in hydrologic science and engineering to consider the CN method as a stepping stone that has outlived its usefulness in research. The CN method fills a narrow niche in certain settings as a parsimonious method having utility as an empirical equation to estimate runoff from a given amount of rainfall, which originated as a static functional form that fits rainfall-runoff data sets. Sixty five years of use and multiple reinterpretations have not resulted in improved hydrological predictability using the method. We suggest that the research community should move forward by (1) identifying appropriate dynamic hydrological model formulations for different hydro-geographic settings, (2) specifying needed model capabilities for solving different classes of problems (e.g. flooding, erosion/sedimentation, nutrient transport, water management, etc.) in different hydro-geographic settings, and (3) expanding data collection and research programs to help ameliorate the so-called “over-parameterization” problem in contemporary modeling. Many decades of advances in geo-spatial data and processing, computation, and understanding are being squandered on continued focus on the static CN regression method. It is time to truly “move beyond” the Curve Number method.

Preferential flow paths (PFPs) affect the hydrological response of humid tropical catchments but have not received sufficient attention. We consider PFPs created by tree roots and earthworms in a near-surface soil layer in steep, humid, tropical lowland catchments and hypothesize that observed hydrological behaviors can be better captured by reasonably considering PFPs in this layer. We test this hypothesis by evaluating the performance of four different physically-based distributed model structures without and with PFPs in different configurations. Model structures are tested both quantitatively and qualitatively using hydrological, geophysical, and geochemical data both from the Smithsonian Tropical Research Institute, Agua Salud Project experimental catchment(s) in central Panama and other sources in the literature. The performance of different model structures is evaluated using runoff volume error and three Nash-Sutcliffe efficiency measures against observed total runoff, storm- and baseflows along with visual comparison of simulated and observed hydrographs. Two of the four proposed model structures which include both lateral and vertical PFPs are plausible, but the one with explicit simulation of PFPs performs the best. A small number of vertical PFPs that fully extend below the root zone allow the model to reasonably simulate deep groundwater recharge, which plays a crucial role in baseflow generation. Results also show that the shallow lateral PFPs are the main contributor to the observed high flow characteristics. Their number and size distribution are found to be more important than the depth distribution. Our model results are corroborated by geochemical and geophysical observations.

Numerical solution of the one-dimensional Richards' equation is the recommended method for coupling groundwater to the atmosphere through the vadose zone in hyper-resolution Earth System Models, but requires fine spatial discretization, is computationally expensive, and may not converge due to mathematical degeneracy or when sharp wetting fronts occur. We transformed the one-dimensional Richards' equation into a new equation that describes the velocity of moisture content values in an unsaturated soil under the actions of capillarity and gravity. We call this new equation the Soil Moisture Velocity Equation (SMVE). The SMVE consists of two terms: an advection-like term that accounts for gravity and the integrated capillary drive of the wetting front, and a diffusion-like term that describes the flux due to the shape of the wetting front capillarity profile divided by the vertical gradient of the capillary pressure head. The SMVE advection-like term can be converted to a relatively easy to solve ordinary differential equation (ODE) using the method of lines and solved using a finite moisture-content discretization. Comparing against analytical solutions of Richards' equation shows that the SMVE advection-like term is >99% accurate for calculating infiltration fluxes neglecting the diffusion-like term. The ODE solution of the SMVE advection-like term is accurate, computationally efficient and reliable for calculating one-dimensional vadose zone fluxes in Earth System and large-scale coupled models of land-atmosphere interaction. It is also well suited for use in inverse problems such as when repeat remote sensing observations are used to infer soil hydraulic properties or soil moisture.

Tropical Storm Erika was a weakly organized tropical storm when its center of circulation passed more than 150 km north of the island of Dominica on August 27, 2015. Hurricane hunter flights had difficulty finding the center of circulation as the storm encountered a high shear environment. Satellite and radar observations showed gyres imbedded within the broader circulation. Radar observations from Guadeloupe show that one of these gyres formed in convergent mid-level flow triggered by orographic convection over the island of Dominica. Gauge-adjusted radar rainfall data indicated between 300 and 750 mm of rainfall on Dominica, most of it over a four hour period. The result was widespread flooding, destruction of property, and loss of life. The extremity of the rainfall on steep watersheds covered with shallow soils was hypothesized to result in near-equilibrium runoff conditions where peak runoff rates equal the watershed-average peak rainfall rate minus a small constant loss rate. Rain gauge adjusted radar rainfall estimates and indirect peak discharge (IPD) measurements from 16 rivers at watershed areas ranging from 0.9 to 31.4 km2 using the USGS Slope-Area method allowed testing of this hypothesis. IPD measurements were compared against the global envelope of maximum observed flood peaks versus drainage area and against simulations using the U.S. Army Corps of Engineers Gridded Surface/Subsurface Hydrologic Analysis (GSSHA) model to detect landslide-affected peak flows. Model parameter values were estimated from the literature. Reasonable agreement was found between GSSHA simulated peak flows and IPD measurements in some watersheds. Results showed that landslide dam failure affected peak flows in 5 of the 16 rivers, with peak flows significantly greater than the envelope curve values for the flood of record for like-sized watersheds on the planet. GSSHA simulated peak discharges showed that the remaining 11 peak flow values were plausible. Simulations of an additional 24 watersheds ranging in size from 2.2 to 75.4 km2 provided confirmation that the IPD measurements varied from 40 to nearly 100% of the envelope curve value depending on storm-total rainfall. Results presented in this paper support the hypothesis that on average, the peak discharges scaled linearly with drainage area, and the constant of proportionality was equivalent to 134 mm h−1, or a unit discharge of 37.22 m3 s−1 km−2. The results also indicate that after the available watershed storage was filled after approximately 450–500 mm of rain fell, runoff efficiencies exceeded 50–60%, and peak runoff rates were more than 80% of the peak rainfall rate minus a small constant loss rate of 20 mm h−1. These findings have important implications for design of resilient infrastructure, and means that rainfall rate was the primary determinant of peak flows once the available storage was filled in the absences of landslide dam failure.

Chemical analyses from 71 watersheds across an ∼450 km transect in west-central Panama provide insight into controls on weathering and rates of chemical denudation and CO2 consumption across an igneous arc terrain in the tropics. Stream and river compositions across this region of Panama are generally dilute, having a total dissolved solute value = 118 ± 91 mg/L, with bicarbonate and silica being the predominant dissolved species. Solute, stable isotope, and radiogenic isotope compositions are consistent with dissolution of igneous rocks present in Panama by meteoric precipitation, with geochemical signatures of rivers largely acquired in their upstream regions. Comparison of a headwater basin with its entire watershed observed considerably more runoff production from the high-elevation upstream portion of the catchment than in its much more spatially extensive downstream region. Rock alteration profiles document that weathering proceeds primarily by dissolution of feldspar and pyroxene, with base cations effectively leached in the following sequence: Na > Ca > Mg > K. Control on water chemistry by bedrock lithology is indicated through a linking of elevated ([Na + K]/[Ca + Mg]) ratios in waters to a high proportion of catchment area silicic bedrock and low ratios to mafic bedrock. Sr-isotope ratios are dominated by basement-derived Sr, with only very minor, if any, contribution from other sources. Cation weathering of Casil + Mgsil + Na + K spans about an order in magnitude, from 3 to 32 tons/km2/yr. Strong positive correlations of chemical denudation and CO2 consumption are observed with precipitation, mean watershed elevation, extent of land surface forest cover, and physical erosion rate.

Summary Process-based hydrological models have a long history dating back to the 1960s. Criticized by some as over-parameterized, overly complex, and difficult to use, a more nuanced view is that these tools are necessary in many situations and, in a certain class of problems, they are the most appropriate type of hydrological model. This is especially the case in situations where knowledge of flow paths or distributed state variables and/or preservation of physical constraints is important. Examples of this include: spatiotemporal variability of soil moisture, groundwater flow and runoff generation, sediment and contaminant transport, or when feedbacks among various Earth’s system processes or understanding the impacts of climate non-stationarity are of primary concern. These are situations where process-based models excel and other models are unverifiable. This article presents this pragmatic view in the context of existing literature to justify the approach where applicable and necessary. We review how improvements in data availability, computational resources and algorithms have made detailed hydrological simulations a reality. Avenues for the future of process-based hydrological models are presented suggesting their use as virtual laboratories, for design purposes, and with a powerful treatment of uncertainty.

Summary Vadose zone flow problems are usually solved from the Richards equation. Solution to the Richards equation is generally challenging because the hydraulic conductivity and diffusivity in the equation are strongly non-linear functions of water content. The finite water-content method was proposed as an alternative general solution method of the vadose zone flow problem for infiltration, falling slugs, and vadose zone response to water table dynamics based on discretizing the water content domain into numerous bins instead of the traditional spatial discretization. In this study, we develop an improved approach to the original finite water-content method (referred to as TO method hereinafter) that better simulates diffusive effects but retains the robustness of the TO method. The approach treats advection and diffusion separately and considers diffusion on a bin by bin basis. After discretizing into water content bins, we treat the conductivity and diffusivity in individual bins as water content dependent constant evaluated at given water content corresponding to each bin. For each bin, we can solve the flow equations analytically since the hydraulic conductivity and diffusivity can be treated as a constant. We then develop solutions for each bin to determine the diffusive water amounts at each time step. The water amount ahead of the convective front for each bin is redistributed among water content bins to account for diffusive effects. The application of developed solution is straightforward only involving algebraic manipulations at each time step. The method can mainly improve water content profiles, but has no significant difference for the total infiltration rate and cumulative infiltration compared to the TO method. Although the method separately deals with advection and diffusion, it can account for the coupling effects of advection and diffusion reasonably well.

Generally, one expects evapotranspiration (ET) maps derived from optical/thermal Landsat and MODIS satellite imagery to improve decision support tools and lead to superior decisions regarding water resources management. However, there is lack of supportive evidence to accept or reject this expectation. We “benchmark” three existing hydrologic decision support tools with the following benchmarks: annual ET for the ET Toolbox developed by the United States Bureau of Reclamation, predicted rainfall-runoff hydrographs for the Gridded Surface/Subsurface Hydrologic Analysis model developed by the U.S. Army Corps of Engineers, and the average annual groundwater recharge for the Distributed Parameter Watershed Model used by Daniel B. Stephens & Associates. The conclusion of this benchmark study is that the use of NASA/USGS optical/thermal satellite imagery can considerably improve hydrologic decision support tools compared to their traditional implementations. The benefits of improved decision making, resulting from more accurate results of hydrologic support systems using optical/thermal satellite imagery, should substantially exceed the costs for acquiring such imagery and implementing the remote sensing algorithms. In fact, the value of reduced error in estimating average annual groundwater recharge in the San Gabriel Mountains, California alone, in terms of value of water, may be as large as $1 billion, more than sufficient to pay for one new Landsat satellite.

Stream hydrograph separation using naturally occurring geochemical tracers holds great potential for elucidating mineral weathering and solute transport. This study addresses a critical need to characterize catchment runoff generation in the humid tropics using multiple natural tracers for hydrograph separation and concentration/discharge (C/Q) hysteresis analysis. We use hydrometric and geochemical data collected at the start of the wet season from three small, steep catchments located in the humid seasonal tropics of central Panama that differ primarily in land cover. We apply a dual source hydrograph separation model between two end-members: new event water precipitation and pre-event water stored in the catchment. We compare the effectiveness of electrical conductivity (EC) and stable water isotopes (δD and δ18O) tracers for identifying precipitation event water in stream runoff using across forested (1.43 km2), mixed land use ‘mosaic’ (1.82 km2) and pasture (0.42 km2) catchments. Hysteretic C/Q loops are analyzed for flowpath interpretation using δD, Ca2+, Mg2+, Na+, K+, Cl−, and SO42−. During a medium-large magnitude event on May 23, 2013, forest and mosaic stream δD, Ca2+, Mg2+, and Na+ exhibited clockwise hysteresis, SO42− exhibited anticlockwise hysteresis, and K+ and Cl− each showed no hysteresis. EC as a surrogate for total dissolved solids agrees acceptably with stable water isotope hydrograph separations during small peak runoff events ( 10 mm/h) in the mosaic catchment. Early wet-season events indicate lower event water fractions than events farther into the wet season. Despite previous work showing land cover strongly controls storm runoff efficiencies, hydrograph separation and hysteresis analyses only indicate weak event water delivery differences between the paired forest and mosaic catchments.

The Isthmus of Panama comprises a lithologically diverse andesitic oceanic arc of Late Cretaceous to Holocene age; it has large spatial variation in rainfall, displays a large range of physical erosion rates, and, therefore, is an ideal location to examine silicate weathering in the tropics. We use a multiyear data set of river chemistry for a 450 km transect across the Cordillera Central of west-central Panama to investigate controls on chemical weathering in tropical small mountainous rivers. Sea-salt corrected cation weathering yields (Casil + Mgsil + Na + K) range over more than an order in magnitude from 3.1 to 31.7 t/km2/yr, while silicate weathering yields (Casil + Mgsil + Na + K + Si) range from 6.9 to 69.5 t/km2/yr. Watershed lithology is the primary control on riverine chemistry, but landscape topographic character and land cover and/or land use also influence solute delivery potential. Strong statistical links of small mountainous river chemical weathering fluxes with rainfall and physical weathering rates attest to the importance of runoff and erosion in maintaining elevated bedrock weathering rates. CO2 consumption ranges from 155 × 103 mol/km2/yr to 1566 × 103 mol/km2/yr, in the upper range of global rates, leading us to suggest that andesite terrains should be considered separately when calculating removal of CO2 from the atmosphere via silicate weathering.

We have developed an alternative to the one-dimensional partial differential equation (PDE) attributed to Richards (1931) that describes unsaturated porous media flow in homogeneous soil layers. Our solution is a set of three ordinary differential equations (ODEs) derived from unsaturated flux and mass conservation principles. We used a hodograph transformation, the Method of Lines, and a finite water-content discretization to produce ODEs that accurately simulate infiltration, falling slugs, and groundwater table dynamic effects on vadose zone fluxes. This formulation, which we refer to as “finite water-content”, simulates sharp fronts and is guaranteed to conserve mass using a finite-volume solution. Our ODE solution method is explicitly integrable, does not require iterations and therefore has no convergence limits and is computationally efficient. The method accepts boundary fluxes including arbitrary precipitation, bare soil evaporation, and evapotranspiration. The method can simulate heterogeneous soils using layers. Results are presented in terms of fluxes and water content profiles. Comparing our method against analytical solutions, laboratory data, and the Hydrus-1D solver, we find that predictive performance of our finite water-content ODE method is comparable to or in some cases exceeds that of the solution of Richards' equation, with or without a shallow water table. The presented ODE method is transformative in that it offers accuracy comparable to the Richards (1931) PDE numerical solution, without the numerical complexity, in a form that is robust, continuous, and suitable for use in large watershed and land-atmosphere simulation models, including regional-scale models of coupled climate and hydrology.

Data from laboratory experiments on a 143 cm tall and 14.5 cm diameter column, packed with Wedron sand with varied constant upper boundary fluxes and water table velocities for both falling and rising water tables are used to validate a finite water-content vadose zone simulation methodology. The one-dimensional finite water-content Talbot and Ogden (2008) (T-O) infiltration and redistribution method was improved to simulate groundwater table dynamic effects and compared against the numerical solution of the Richards equation using Hydrus-1D. Both numerical solutions agreed satisfactorily with time series measurements of water content. Results showed similar performance for both methods, with the T-O method on average having higher Nash-Sutcliffe efficiencies and smaller absolute biases. Hydrus-1D was more accurate in predicting deponding times in the case of a falling water table, while Hydrus-1D and the T-O method had similar errors in predicted ponding times in the case of a rising water table in six of nine tests. The improved T-O method was able to predict general features of vadose zone moisture dynamics with moving water table and surface infiltration using an explicit, mass-conservative formulation. The advantage of an explicit formulation is that it is numerically simple, using forward Euler solution methodology, and is guaranteed to converge and to conserve mass. These properties make the improved T-O method presented in this paper a robust and computationally efficient alternative to the numerical solution of the Richards equation in hydrological modeling applications involving groundwater table dynamic effects on vadose zone soil moistures.

Summary Numerical solution of the Richards’ equation (RE) in variably saturated soils continues to be a challenge due to its highly non-linear behavior. This is particularly true as soils approach saturation and the behavior of the fundamental partial differential equation changes from elliptic to parabolic. In this paper, a finite volume predictor–corrector method with adaptive time-stepping was developed to solve the 1D vertical RE. The numerical method was mass-conservative and non-iterative. In the predictor step, the pressure head-based form of the RE was solved using the cell-centered finite volume method and the pressure head was updated. In the corrector step, the soil water content was calculated by solving the mixed form RE. Five different schemes to evaluate the inter-cell hydraulic conductivity were investigated. The robustness and accuracy of the numerical model were demonstrated through simulation of experimental tests, including free drainage, field infiltration into wet and dry soils, and laboratory infiltration with falling water table. Numerical results were compared against laboratory measurements, simulation results from the Hydrus-1D program, or analytical solution when available. Results showed that the developed scheme is robust and accurate in simulating variably saturated flows with various boundary conditions. The arithmetic mean and Szymkiewicz’s mean of inter-cell hydraulic conductivity performed better than other methods especially in the case of infiltration into very dry soil.

We have developed a one-dimensional numerical method to simulate infiltration and redistribution in the presence of a shallow dynamic water table. This method builds upon the Green-Ampt infiltration with Redistribution (GAR) model and incorporates features from the Talbot-Ogden (T-O) infiltration and redistribution method in a discretized moisture content domain. The redistribution scheme is more physically meaningful than the capillary weighted redistribution scheme in the T-O method. Groundwater dynamics are considered in this new method instead of hydrostatic groundwater front. It is also computationally more efficient than the T-O method. Motion of water in the vadose zone due to infiltration, redistribution, and interactions with capillary groundwater are described by ordinary differential equations. Numerical solutions to these equations are computationally less expensive than solutions of the highly nonlinear Richards' (1931) partial differential equation. We present results from numerical tests on 11 soil types using multiple rain pulses with different boundary conditions, with and without a shallow water table and compare against the numerical solution of Richards' equation (RE). Results from the new method are in satisfactory agreement with RE solutions in term of ponding time, deponding time, infiltration rate, and cumulative infiltrated depth. The new method, which we call “GARTO” can be used as an alternative to the RE for 1-D coupled surface and groundwater models in general situations with homogeneous soils with dynamic water table. The GARTO method represents a significant advance in simulating groundwater surface water interactions because it very closely matches the RE solution while being computationally efficient, with guaranteed mass conservation, and no stability limitations that can affect RE solvers in the case of a near-surface water table.

Flow-measurement weirs in natural channels experience sedimentation. Removing accumulated sediment is often impractical, prohibitively expensive, or both. This paper describes the results o...

Recent upgrades to operational radar-rainfall products in terms of quality and resolution call for reexamination of the factors that contribute to the uncertainty of radar-rainfall estimation. Remapping or regridding of radar observations onto Cartesian coordinates is implemented by practitioners when radar estimates are compared against rain gauge observations, in hydrologic applications, or for merging data from different radars. However, assuming perfect radar observations, many of the widely used remapping methodologies do not conserve mass for the rainfall rate field. The most popular remapping approaches used are those based on extracting information from radar bins whose centers fall within a certain distance from the center of the Cartesian grid. This paper develops a mass-conserving method for remapping, which is called “precise remapping,” which is compared against two other commonly used remapping methods. Results show that the choice of the remapping method can make a substantial difference in grid-averaged rainfall accumulations (up to more than 100%). Differences were quantified using observations from two radars, collected during a field experiment. The interpolation grid resolution was also found to affect interpolated rainfall estimates. Approximate remapping methods tend to be much more sensitive to the interpolation grid resolution than precise remapping. High-resolution radar data such as those from radars with short gate spacing or narrow beams, or the super-resolution Weather Surveillance Radar-1988 Doppler (WSR-88D) sampling format, are significantly more sensitive (by up to 100%) to the remapping method and the interpolation grid resolution than the legacy WSR-88D rainfall data resolution of 1° × 1 km.

Statistical analyses were performed using five years (1996–2000) of Weather Surveillance Radar–1988 Doppler (WSR-88D) weather-radar rainfall estimates that were produced for the Global-Continental Scale International Project over the entire Mississippi River basin. The project radar rainfall estimates were adjusted using a Z-R relation optimized to improve performance at the time scale of individual storms. The accuracy of radar-rainfall estimates were analyzed during the warm season considering a number of factors, including number of overlapping radars, distance from gauge to nearest radar, gauge elevation, and the geographic location of the radar. Results were used to identify optimal radar-rainfall estimation areas (ORREAs) within the Mississippi River basin with high correlation between storm-total radar and rain gauge rainfall. Additionally, estimation of the bias between storm-total radar-rainfall and rain gauge rainfall accumulations in areas without gauges was assessed to identify the app...

Summary Tropical hydrology is poorly understood for a number of reasons. Intense biological activity in the tropics introduces complexities to the hydrologic process. Bioturbation, rapid rates of decay, and intensive insect activity all tend to promote rapid flow paths in the upper soil. Aggressive weathering leads to clays depleted of light cations and deep soil profiles. Processes in the seasonal tropics are further complicated by seasonal transitions, and very large changes in catchment storage between seasons. Beginning in 2005, we installed a suite of hydrologic sensors in a 16.7 ha first-order catchment in the Panama Canal Watershed to observe hydrologic variables and identify the dominant streamflow generation processes. The site is located near the village of Gamboa, which is located on the east bank of the Panama Canal at the confluence of Lake Gatun and the Chagres River. The study catchment is located on the north side of a ridge off the eastern flank of a 230 m tall hill known as Cerro Pelado, and is covered by 70–120 year old re-growth triple-canopy forest. Measurements included: rainfall above the canopy, throughfall, stemflow, evapotranspiration, shallow groundwater levels and streamflow. Deep groundwater storage was not measured. This paper describes measurements made, data collected, and the worth of those data in estimating the mass balance closure of a first-order catchment during the wet season. We compare measurements of the different components of the water cycle with observations from other published studies from the tropics. Data analysis results indicate water balance closure errors of approximately 8%.

Hydrological processes in the humid tropics differ from other regions in having greater energy inputs and faster rates of change. In this Review it is argued that understanding of the key hydrological interactions there remains limited, and a vision of future research designed to address these shortcomings is outlined.

The TalbotOgden model is a mass conservative method to simulate flow of a wetting liquid in variably-saturated porous media. The principal feature of this model is the discretization of the moisture content domain into bins. This paper gives an analysis of the relationship between the number of bins and the computed flux. Under the circumstances of discrete bins and discontinuous wetting fronts, we show that fluxes increase with the number of bins. We then apply this analysis to the continuous case and get an upper bound of the difference of infiltration rates when the number of bins tends to infinity. We also extend this model by creating a two dimensional moisture content domain so that there exists a probability distribution of the moisture content for different soil systems. With these theoretical and experimental results and using a Dynamic Data Driven Application System (DDDAS), sensors can be put in soils to detect the infiltration fluxes, which are important to compute the proper number of bins for a specific soil system and predict fluxes. Using this feedback control loop, the extended TalbotOgden model can be made more efficient for estimating infiltration into soils.

Summary HBV is a conceptual hydrological model extensively used in operational hydrological forecasting and water balance studies. In this paper, we apply the HBV model on the 1924 km 2 semi-humid Leaf River catchment near Collins, Mississippi. We analyze individual sensitivity of the parameters by calibrating the model using the Multi-Objective Shuffled Complex Evolution (MOSCEM) algorithm, perform Monte-Carlo based identifiability analysis and investigate the dynamic behavior of the parameters using the Dynamic Identifiability Analysis (DYNIA) approach in reference to the hydrological process in the catchment. The sensitivity analysis using two objective measures showed that there are distinct groups of parameters that control total runoff volume errors and errors from the high-flow series. The DYNIA analysis revealed that parameters have specific periods where they show higher identifiability and play a crucial role in representing the predicted stream flow. Temporal changes of parameter optima were observed due either to inadequacies in the model structure or possible time-varying catchment response subject to unsteady hydrodynamic and hydroclimatic conditions.

This research develops a one-parameter model of saturated source area dynamics and the spatial distribution of soil moisture. The single required parameter is the maximum soil moisture deficit within the catchment. The concept behind the development of the model comes from the fact that the complexity of topographically-driven runoff generation can be reduced through the use of geomorphological scaling relations. The scaling formulation allows the prediction of the dynamics of saturated source areas as a function of basin-wide soil moisture state. This model offers a number of potential advantages. Firstly, the model parameter is independent of topographic index distribution and its associated scale effects. Secondly, it may be possible to measure this single parameter using field measurements or perhaps remote sensing, which gives the model significant potential for application in ungauged basins. Finally, the fact that this parameter is a physical characteristic of the basin, estimation of this parameter avoids regionalization and parameter transferability problems. The model is tested using rainfall–runoff data from the 10.4 ha experimental catchment known as Tarrawara in Australia, the 37 km2 Town Creek catchment in U.S.A., and the 620 km2 Balaphi and the 850 km2 Likhu sub-catchments of the Koshi river in Nepal. In sub-catchments of Koshi river, the simulation results compare favorably against the calibrated TOPMODEL both in terms of direct runoff and the spatial distribution of soil moisture state. In the Tarrawara and Town Brook catchments, simulation results compare favorably against observed storm runoff using all observed data, without calibration.

This article presents a methodology for estimating total incoming solar radiation from Triangular Irregular Network (TIN) topographic meshes. The algorithm also computes terrain slope degree and aspect (slope orientation) and accounts for self shading and cast shadows, sky view fractions for diffuse radiation, remote albedo and atmospheric backscattering, by using a vectorial approach within a topocentric coordinate system establishing geometric relations between groups of TIN elements and the sun position. A normal vector to the surface of each TIN element describes its slope and aspect while spherical trigonometry allows computing a unit vector defining the position of the sun at each hour and day of the year. Sky view fraction, useful to determine diffuse and backscattered radiation, is computed for each TIN element at prescribed azimuth intervals targeting the steepest elevation gradient. A comparison between the sun zenith angle and the steepest gradient allows deciding whether or not the pivot element is shaded. Finally, remote albedo is computed from the sky view fraction complementary functions for observed albedo values of the surrounding terrain. The sensitivity of the different radiative components to seasonal changes in atmospheric transmissivitties and surrounding albedo is tested in a mountainous watershed in Wyoming. This methodology represents an improvement on the current algorithms to compute terrain and radiation values on unstructured-mesh terrain models. All terrain-related features (e.g., slope, aspect, sky view fraction) can be pre-computed and stored for easy access into a subsequent, progressive-in-time, numerical simulation.

The Upper Rio Chagres drains 414 km2 of steep, mountainous terrain in central Panama. A tropical air mass thunderstorm on 10 July 2007 produced a flood across the basin that peaked at 720 m3 s− 1 at a headwaters gage draining 17.5 km2 and 1710 m3 s− 1 at a downstream gage draining 414 km2. The storm also triggered numerous landslides in the upper basin, which facilitated the formation of large logjams along portions of the channel where transport capacity of wood was reduced by a change in channel geometry such as a bend or channel expansion. During field work in February 2008, we characterized three jams with surface areas of 400–2450 m2; two of these jams resulted in storage of substantial (1100–8200 m3) sediment wedges upstream. We returned to these sites in March 2009 to document changes in the logjams and sediment storage. Drawing on observations made in the basin since 2002, and site visits during 2008 and 2009, we suggest that jams such as these last two years or less. We propose that wood dynamics in the Upper Chagres alternate between brief periods of moderate wood load in the form of large logjams and much longer periods of essentially no wood load, a situation that contrasts with the more consistent wood loads in catchments of similar size in temperate environments and with limited studies of more consistent wood load in tropical catchments with no landslides.

This article presents the development of an approach to calibrate weather radar, by applying the method of the Marching–Jury Backward Beam Equation (J. Atmadja and A.C. Bagtzoglou, Pollution source identification in heterogeneous porous media, Water Resour. Res. 37, (2001), pp. 2113–2125; A.C. Bagtzoglou and J. Atmadja, The marching-jury backward beam equation and quasi-reversibility methods for hydrologic inversion: Application to contaminant plume spatial distribution recovery, Water Resour. Res. 39 (2003), p. 1038; A.C. Bagtzoglou and S. Baun, Near real-time atmospheric contamination source identification by an optimization-based inverse method, Inv. Prob. Sc. Eng. 13 (2005), pp. 261–278), coupled to physically-based discrete Fourier image processing techniques. We have enhanced the method by integrating an optimization scheme that takes as input parameters the stabilization parameter and the pseudo coefficient of diffusion. The objective function is set as a weighted sum of different mass and peak err...

At the watershed scale, soil moisture is the major control for rainfall–runoff response, especially where saturation excess runoff processes dominate. From the ecological point of view, the pools of soil moisture are fundamental ecosystem resources providing the transpirable water for plants. In drylands particularly, soil moisture is one of the major controls on the structure, function, and diversity in ecosystems. In terms of the global hydrological cycle, the overall quantity of soil moisture is small, ∼0.05%; however, its importance to the global energy balance and the distribution of precipitation far outweighs its physical amount. In soils it governs microbial activity that affects important biogeochemical processes such as nitrification and CO2 production via respiration. During the past 20 years, technology has advanced considerably, with the development of different electrical sensors for determining soil moisture at a point. However, modeling of watersheds requires areal averages. As a result, point measurements and modeling grid cell data requirements are generally incommensurate. We review advances in sensor technology, particularly emerging geophysical methods and distributed sensors, aimed at bridging this gap. We consider some of the data analysis methods for upscaling from a point to give an areal average. Finally, we conclude by offering a vision for future research, listing many of the current scientific and technical challenges.

Mansfield Hollow Lake (MHL) and Willimantic Reservoir (WR) are two reservoir lakes located in eastern Connecticut in the northeastern United States. MHL formed behind the Mansfield Hollow Dam constructed by the U.S. Army Corps of Engineers in 1952 and is primarily fed by the Fenton, Mount Hope and Natchaug Rivers. The WR lies approximately 1-km downstream from the Mansfield Hollow Dam. Total dissolved nitrogen, phosphorus and chlorophyll a measurements indicate the water bodies could be classified as borderline mesotrophic/eutrophic. A steady-state numerical software package (Bathtub) designed to facilitate application of empirical eutrophication models to morphometrically complex reservoirs was used to determine the trophic status in MHL and WR based on different phosphorus and nitrogen loading budgets. The short hydraulic residence times and rapid flushing rates in MHL and WR are directly related to the flow rates in the streams discharging into MHL. The low flow period could significantly incr...

The Mansfield Hollow Lake (MHL) and Willimantic Reservoir (WR) are two reservoirs located in eastern Connecticut in the north-eastern USA. The MHL was constructed by the US Army Corps of Engineers in 1952, being primarily fed by the Fenton, Mount Hope and Natchaug Rivers. The WR lies downstream from the Mansfield Hollow Dam. The physical and chemical characteristics of MHL, WR, and the three streams discharging into them were intensively evaluated during four sampling quarters in 2001 and 2002. This study focused on possible sources of eutrophication, and occurrence of trace amounts of selected priority pollutants, in the sediments and waters of MHL and WR. Analytical results for total dissolved nitrogen and phosphorus and chlorophyll-a indicated that, in terms of eutrophication status, the waterbodies can be classified as mesotrophic. Comparison of the analytical results in this study to historical limnological data reported for lakes in southern New England and eastern Connecticut indicates that the median levels of total nitrogen, total phosphorus, dissolved inorganic carbon, calcium, magnesium, and iron measured in this study were higher than previously reported median values.

Summary The worth of archival, long-term, high spatial and temporal resolution radar-rainfall estimates for hydrologic modeling has not yet been established. In particular, uncertainties remain regarding accuracy at the event time-scale, which is a modeling requirement for mass conservation. To explore this issue, an archival radar-rainfall precipitation dataset produced from the National Reflectivity Composite for the GEWEX Continental-Scale International Project (GCIP) and US National Weather Service, Weather Surveillance Radar 1988-D (WSR-88D) Stage III (SIII) data were compared with data from 982 National Climatic Data Center (NCDC) rain gages throughout the Mississippi River Basin. These two datasets are derived from the same radar-reflectivity observations but they are different. The GCIP dataset was developed from low resolution (5 dBZ bins) national radar mosaic imagery using a systematic procedure that employed one climatologically-derived reflectivity-rainfall relationship, while the SIII radar-rainfall estimates are produced using a processing algorithm developed by the US National Weather Service that employs near-real time adjustment and varied parameterizations in different River Forecast Centers around the Mississippi River Basin. Statistical analyses were performed to compare both radar-rainfall datasets with rain gage data on an event time-scale for the entire year of 1997. The investigation shows that the GCIP dataset, given a change in the equation used to calculate rainfall from reflectivity, has a lower root mean square difference for the majority of gages in four of the five river forecast centers that cover the Mississippi River Basin compared to SIII in 1997.