Engineering topics
Wood, Eric F.
Publications and source records attributed to Wood, Eric F..
The Future of Evapotranspiration: Global Requirements for Ecosystem Functioning, Carbon and Climate Feedbacks, Agricultural Management, and Water Resources
The fate of the terrestrial biosphere is highly uncertain given recent and projected changes in climate. This is especially acute for impacts associated with changes in drought frequency and intensity on the distribution and timing of water availability. The development of effective adaptation strategies for these emerging threats to food and water security are compromised by limitations in our understanding of how natural and managed ecosystems are responding to changing hydrological and climatological regimes. This information gap is exacerbated by insufficient monitoring capabilities from local to global scales. Here, we describe how evapotranspiration (ET) represents the key variable in linking ecosystem functioning, carbon and climate feedbacks, agricultural management, and water resources, and highlight both the outstanding science and applications questions and the actions, especially from a space-based perspective, necessary to advance them.
Technical Report Series on Global Modeling and Data Assimilation, Volume 41 : GDIS Workshop Report
The workshop "An International Global Drought Information System Workshop: Next Steps" was held on 10-13 December 2014 in Pasadena, California. The more than 60 participants from 15 countries spanned the drought research community and included select representatives from applications communities as well as providers of regional and global drought information products. The workshop was sponsored and supported by the US National Integrated Drought Information System (NIDIS) program, the World Climate Research Program (WCRP: GEWEX, CLIVAR), the World Meteorological Organization (WMO), the Group on Earth Observations (GEO), the European Commission Joint Research Centre (JRC), the US Climate Variability and Predictability (CLIVAR) program, and the US National Oceanic and Atmospheric Administration (NOAA) programs on Modeling, Analysis, Predictions and Projections (MAPP) and Climate Variability & Predictability (CVP). NASA/JPL hosted the workshop with logistical support provided by the GEWEX program office. The goal of the workshop was to build on past Global Drought Information System (GDIS) progress toward developing an experimental global drought information system. Specific goals were threefold: (i) to review recent research results focused on understanding drought mechanisms and their predictability on a wide range of time scales and to identify gaps in understanding that could be addressed by coordinated research; (ii) to help ensure that WRCP research priorities mesh with efforts to build capacity to address drought at the regional level; and (iii) to produce an implementation plan for a short duration pilot project to demonstrate current GDIS capabilities. See http://www.wcrp-climate.org/gdis-wkshp-2014-objectives for more information.
The North American Multi-Model Ensemble (NMME): Phase-1 Seasonal to Interannual Prediction, Phase-2 Toward Developing Intra-Seasonal Prediction
The recent US National Academies report "Assessment of Intraseasonal to Interannual Climate Prediction and Predictability" was unequivocal in recommending the need for the development of a North American Multi-Model Ensemble (NMME) operational predictive capability. Indeed, this effort is required to meet the specific tailored regional prediction and decision support needs of a large community of climate information users. The multi-model ensemble approach has proven extremely effective at quantifying prediction uncertainty due to uncertainty in model formulation, and has proven to produce better prediction quality (on average) then any single model ensemble. This multi-model approach is the basis for several international collaborative prediction research efforts, an operational European system and there are numerous examples of how this multi-model ensemble approach yields superior forecasts compared to any single model. Based on two NOAA Climate Test Bed (CTB) NMME workshops (February 18, and April 8, 2011) a collaborative and coordinated implementation strategy for a NMME prediction system has been developed and is currently delivering real-time seasonal-to-interannual predictions on the NOAA Climate Prediction Center (CPC) operational schedule. The hindcast and real-time prediction data is readily available (e.g., http://iridl.ldeo.columbia.edu/SOURCES/.Models/.NMME/) and in graphical format from CPC (http://origin.cpc.ncep.noaa.gov/products/people/wd51yf/NMME/index.html). Moreover, the NMME forecast are already currently being used as guidance for operational forecasters. This paper describes the new NMME effort, presents an overview of the multi-model forecast quality, and the complementary skill associated with individual models.
Hyperresolution Global Land Surface Modeling: Meeting a Grand Challenge for Monitoring Earth's Terrestrial Water
Monitoring Earth's terrestrial water conditions is critically important to many hydrological applications such as global food production; assessing water resources sustainability; and flood, drought, and climate change prediction. These needs have motivated the development of pilot monitoring and prediction systems for terrestrial hydrologic and vegetative states, but to date only at the rather coarse spatial resolutions (approx.10-100 km) over continental to global domains. Adequately addressing critical water cycle science questions and applications requires systems that are implemented globally at much higher resolutions, on the order of 1 km, resolutions referred to as hyperresolution in the context of global land surface models. This opinion paper sets forth the needs and benefits for a system that would monitor and predict the Earth's terrestrial water, energy, and biogeochemical cycles. We discuss six major challenges in developing a system: improved representation of surface-subsurface interactions due to fine-scale topography and vegetation; improved representation of land-atmospheric interactions and resulting spatial information on soil moisture and evapotranspiration; inclusion of water quality as part of the biogeochemical cycle; representation of human impacts from water management; utilizing massively parallel computer systems and recent computational advances in solving hyperresolution models that will have up to 10(exp 9) unknowns; and developing the required in situ and remote sensing global data sets. We deem the development of a global hyperresolution model for monitoring the terrestrial water, energy, and biogeochemical cycles a grand challenge to the community, and we call upon the international hydrologic community and the hydrological science support infrastructure to endorse the effort.
Comparison of Two Methods for Estimating the Sampling-Related Uncertainty of Satellite Rainfall Averages Based on a Large Radar Data Set
The uncertainty of rainfall estimated from averages of discrete samples collected by a satellite is assessed using a multi-year radar data set covering a large portion of the United States. The sampling-related uncertainty of rainfall estimates is evaluated for all combinations of 100 km, 200 km, and 500 km space domains, 1 day, 5 day, and 30 day rainfall accumulations, and regular sampling time intervals of 1 h, 3 h, 6 h, 8 h, and 12 h. These extensive analyses are combined to characterize the sampling uncertainty as a function of space and time domain, sampling frequency, and rainfall characteristics by means of a simple scaling law. Moreover, it is shown that both parametric and non-parametric statistical techniques of estimating the sampling uncertainty produce comparable results. Sampling uncertainty estimates, however, do depend on the choice of technique for obtaining them. They can also vary considerably from case to case, reflecting the great variability of natural rainfall, and should therefore be expressed in probabilistic terms. Rainfall calibration errors are shown to affect comparison of results obtained by studies based on data from different climate regions and/or observation platforms.
Estimation of Rainfall Sampling Uncertainty: A Comparison of Two Diverse Approaches
The spatial and temporal intermittence of rainfall causes the averages of satellite observations of rain rate to differ from the "true" average rain rate over any given area and time period, even if the satellite observations are perfectly accurate. The difference of satellite averages based on occasional observation by satellite systems and the continuous-time average of rain rate is referred to as sampling error. In this study, rms sampling error estimates are obtained for average rain rates over boxes 100 km, 200 km, and 500 km on a side, for averaging periods of 1 day, 5 days, and 30 days. The study uses a multi-year, merged radar data product provided by Weather Services International Corp. at a resolution of 2 km in space and 15 min in time, over an area of the central U.S. extending from 35N to 45N in latitude and 100W to 80W in longitude. The intervals between satellite observations are assumed to be equal, and similar In size to what present and future satellite systems are able to provide (from 1 h to 12 h). The sampling error estimates are obtained using a resampling method called "resampling by shifts," and are compared to sampling error estimates proposed by Bell based on earlier work by Laughlin. The resampling estimates are found to scale with areal size and time period as the theory predicts. The dependence on average rain rate and time interval between observations is also similar to what the simple theory suggests.
An Agenda for Land-Surface Hydrology Research and a Call for the Second International Hydrological Decade
An agenda for land-surface hydrology research is proposed to open the debate for more comprehensive prioritization of science and application activities in the hydrologic sciences. A set of science questions are posed and the observational requirements to achieve substantial progress are identified. In this context, the proposal to initiate the 2nd International Hydrologic Decade (IHD) is put forth. The benefits of this initiative for enhanced scientific understanding and improved capability in meeting societal needs are also identified.
Simulation of Boreal Ecosystem Carbon and Water Budgets: Scaling from Local to Regional Extents
A coupled water and energy balance model is developed. This model can predict the partitioning of water and energy between major source, sink and storage elements within the Boreal-Ecosystem-Atmospheric Study (BOREAS) areas. The results of testing the model against data collected at BOREAS tower sites during Intensive Field Campaigns and remotely sensed data collected across the BOREAS region are presented.
A simple hydrologically based model of land surface water and energy fluxes for general circulation models
A generalization of the single soil layer variable infiltration capacity (VIC) land surface hydrological model previously implemented in the Geophysical Fluid Dynamics Laboratory (GFDL) general circulation model (GCM) is described. The new model is comprised of a two-layer characterization of the soil column, and uses an aerodynamic representation of the latent and sensible heat fluxes at the land surface. The infiltration algorithm for the upper layer is essentially the same as for the single layer VIC model, while the lower layer drainage formulation is of the form previously implemented in the Max-Planck-Institut GCM. The model partitions the area of interest (e.g., grid cell) into multiple land surface cover types; for each land cover type the fraction of roots in the upper and lower zone is specified. Evapotranspiration consists of three components: canopy evaporation, evaporation from bare soils, and transpiration, which is represented using a canopy and architectural resistance formulation. Once the latent heat flux has been computed, the surface energy balance is iterated to solve for the land surface temperature at each time step. The model was tested using long-term hydrologic and climatological data for Kings Creek, Kansas to estimate and validate the hydrological parameters, and surface flux data from three First International Satellite Land Surface Climatology Project Field Experiment (FIFE) intensive field campaigns in the summer-fall of 1987 to validate the surface energy fluxes.
Scaling water and energy fluxes in climate systems - Three land-atmospheric modeling experiments
Three numerical experiments that investigate the scaling of land-surface processes - either of the inputs or parameters - are reported, and the aggregated processes are compared to the spatially variable case. The first is the aggregation of the hydrologic response in a catchment due to rainfall during a storm event and due to evaporative demands during interstorm periods. The second is the spatial and temporal aggregation of latent heat fluxes, as calculated from SiB. The third is the aggregation of remotely sensed land vegetation and latent and sensible heat fluxes using TM data from the FIFE experiment of 1987 in Kansas. In all three experiments it was found that the surface fluxes and land characteristics can be scaled, and that macroscale models based on effective parameters are sufficient to account for the small-scale heterogeneities investigated.
Research into the influence of spatial variability and scale on the parameterization of hydrological processes
The objectives of the research were as follows: (1) Extend the Representative Elementary Area (RE) concept, first proposed and developed in Wood et al, (1988), to the water balance fluxes of the interstorm period (redistribution, evapotranspiration and baseflow) necessary for the analysis of long-term water balance processes. (2) Derive spatially averaged water balance model equations for spatially variable soil, topography and vegetation, over A RANGE OF CLIMATES. This is a necessary step in our goal to derive consistent hydrologic results up to GCM grid scales necessary for global climate modeling. (3) Apply the above macroscale water balance equations with remotely sensed data and begin to explore the feasibility of parameterizing the water balance constitutive equations at GCM grid scale.
On hydrologic similarity: A dimensionless flood frequency model using a generalized geomorphologic unit hydrograph and partial area runoff generation
One of the shortcomings of the original theory of the geomorphologic unit hydrograph (GUH) is that it assumes that runoff is generated uniformly from the entire catchment area. It is now recognized that in many catchments much of the runoff during storm events is produced on partial areas which usually form on narrow bands along the stream network. A storm response model that includes runoff generation on partial areas by both Hortonian and Dunne mechanisms was recently developed by the authors. In this paper a methodology for integrating this partial area runoff generation model with the GUH-based runoff routing model is presented; this leads to a generalized GUH. The generalized GUH and the storm response model are then used to estimate physically based flood frequency distributions. In most previous work the initial moisture state of the catchment had been assumed to be constant for all the storms. In this paper we relax this assumption and allow the initial moisture conditions to vary between storms. The resulting flood frequency distributions are cast in a scaled dimensionless framework where issues such as catchment scale and similarity can be conveniently addressed. A number of experiments are performed to study the sensitivity of the flood frequency response to some of the 'similarity' parameters identified in this formulation. The results indicate that one of the most important components of the derived flood frequency model relates to the specification of processes within the runoff generation model; specifically the inclusion of both saturation excess and Horton infiltration excess runoff production mechanisms. The dominance of these mechanisms over different return periods of the flood frequency distribution can significantly affect the distributional shape and confidence limits about the distribution. Comparisons with observed flood distributions seem to indicate that such mixed runoff production mechanisms influence flood distribution shape. The sensitivity analysis also indicated that the incorporation of basin and rainfall storm scale also greatly influences the distributional shape of the flood frequency curve.
Similarity and scale in catchment storm response
Until recently, very little progress had been made in understanding the relationship between small-scale variability of topography, soil, and rainfalls and the storm response seen at the catchment scale. The work reviewed here represents the first attempt at a systematic theoretical framework for such understanding in the context of surface runoff generation by different processes. The parameterization of hydrological processes over a range of scales is examined, and the concept of the 'representative elementary area' (REA) is introduced. The REA is a fundamental scale for catchment modeling at which continuum assumptions can be applied for the spatially variable controls and parameters, and spatial patterns no longer have to be considered explicitly. The investigation of scale leads into the concept of hydrologic similarity in which the effects of the environmental controls on runoff generation and flood frequency response be investigated independently of catchment scale. The paper reviews the authors' initial results and hopefully will motivate others to also investigate the issues of hydrologic scale and similarity.
Evapotranspiration and runoff from large land areas: Land surface hydrology for atmospheric general circulation models
A land surface hydrology parameterization for use in atmospheric GCM's is presented. The parameterization incorporates subgrid scale variability in topography, soils, soil moisture and precipitation. The framework of the model is the statistical distribution of a topography-soils index, which controls the local water balance fluxes, and is therefore taken to represent the large land area. Spatially variable water balance fluxes are integrated with respect to the topography-soils index to yield our large topography-soils distribution, and interval responses are weighted by the probability of occurrence of the interval. Grid square averaged land surface fluxes result. The model functions independently as a macroscale water balance model. Runoff ratio and evapotranspiration efficiency parameterizations are derived and are shown to depend on the spatial variability of the above mentioned properties and processes, as well as the dynamics of land surface-atmosphere interactions.
Heterogeneity and scaling land-atmospheric water and energy fluxes in climate systems
The effects of small-scale heterogeneity in land surface characteristics on the large-scale fluxes of water and energy in land-atmosphere system has become a central focus of many of the climatology research experiments. The acquisition of high resolution land surface data through remote sensing and intensive land-climatology field experiments (like HAPEX and FIFE) has provided data to investigate the interactions between microscale land-atmosphere interactions and macroscale models. One essential research question is how to account for the small scale heterogeneities and whether 'effective' parameters can be used in the macroscale models. To address this question of scaling, three modeling experiments were performed and are reviewed in the paper. The first is concerned with the aggregation of parameters and inputs for a terrestrial water and energy balance model. The second experiment analyzed the scaling behavior of hydrologic responses during rain events and between rain events. The third experiment compared the hydrologic responses from distributed models with a lumped model that uses spatially constant inputs and parameters. The results show that the patterns of small scale variations can be represented statistically if the scale is larger than a representative elementary area scale, which appears to be about 2 - 3 times the correlation length of the process. For natural catchments this appears to be about 1 - 2 sq km. The results concerning distributed versus lumped representations are more complicated. For conditions when the processes are nonlinear, then lumping results in biases; otherwise a one-dimensional model based on 'equivalent' parameters provides quite good results. Further research is needed to fully understand these conditions.
Evaluation of a distributed catchment scale water balance model
The validity of some of the simplifying assumptions in a conceptual water balance model is investigated by comparing simulation results from the conceptual model with simulation results from a three-dimensional physically based numerical model and with field observations. We examine, in particular, assumptions and simplifications related to water table dynamics, vertical soil moisture and pressure head distributions, and subsurface flow contributions to stream discharge. The conceptual model relies on a topographic index to predict saturation excess runoff and on Philip's infiltration equation to predict infiltration excess runoff. The numerical model solves the three-dimensional Richards equation describing flow in variably saturated porous media, and handles seepage face boundaries, infiltration excess and saturation excess runoff production, and soil driven and atmosphere driven surface fluxes. The study catchments (a 7.2 sq km catchment and a 0.64 sq km subcatchment) are located in the North Appalachian ridge and valley region of eastern Pennsylvania. Hydrologic data collected during the MACHYDRO 90 field experiment are used to calibrate the models and to evaluate simulation results. It is found that water table dynamics as predicted by the conceptual model are close to the observations in a shallow water well and therefore, that a linear relationship between a topographic index and the local water table depth is found to be a reasonable assumption for catchment scale modeling. However, the hydraulic equilibrium assumption is not valid for the upper 100 cm layer of the unsaturated zone and a conceptual model that incorporates a root zone is suggested. Furthermore, theoretical subsurface flow characteristics from the conceptual model are found to be different from field observations, numerical simulation results, and theoretical baseflow recession characteristics based on Boussinesq's groundwater equation.
A detailed model for simulation of catchment scale subsurface hydrologic processes
A catchment scale numerical model is developed based on the three-dimensional transient Richards equation describing fluid flow in variably saturated porous media. The model is designed to take advantage of digital elevation data bases and of information extracted from these data bases by topographic analysis. The practical application of the model is demonstrated in simulations of a small subcatchment of the Konza Prairie reserve near Manhattan, Kansas. In a preliminary investigation of computational issues related to model resolution, we obtain satisfactory numerical results using large aspect ratios, suggesting that horizontal grid dimensions may not be unreasonably constrained by the typically much smaller vertical length scale of a catchment and by vertical discretization requirements. Additional tests are needed to examine the effects of numerical constraints and parameter heterogeneity in determining acceptable grid aspect ratios. In other simulations we attempt to match the observed streamflow response of the catchment, and we point out the small contribution of the streamflow component to the overall water balance of the catchment.