Automated Calibration of the METRIC‐Landsat Evapotranspiration Process

A remaining challenge to applying satellite‐based energy‐balance algorithms for operational estimation of evapotranspiration (ET) is the calibration of the energy‐balance model. Customized calibration for each image date is generally required to overcome biases associated with radiometric accuracy of the image, uncertainties in aerodynamic features of the landscape, background thermal conditions, and model assumptions. The CIMEC process (calibration using inverse modeling at extreme conditions) is an endpoint calibration procedure where near extreme conditions in the image are identified where the ET can be estimated and assigned. In the Mapping EvapoTranspiration at high Resolution with Internalized Calibration (METRIC?) energy‐balance model, two endpoints represent the dry and wet ends of the ET spectrum. Generally, user‐intervention is required to select locations in the image to produce best accuracy. To bring the METRIC and similar processes into the domain of less experienced operators, a consistent, reproducible, and dependable statistics‐based procedure is introduced where relationships between vegetation amount and surface temperature are used to identify a subpopulation of locations (pixels) in an image that may best represent the calibration endpoints. This article describes the background and logic for the statistical approach, how the statistics were developed, area of interest requirements and assumptions, adjustment for dry conditions in desert climates, and implementation in a common image processing environment (ERDAS Imagine).     

Assessing Calibration Uncertainty and Automation for Estimating Evapotranspiration from Agricultural Areas Using METRIC

Agricultural irrigation accounts for a large fraction of the total water use in the western United States. The Mapping Evapotranspiration at high Resolution with Internalized Calibration (METRIC) remote sensing energy balance model is being used to estimate historical agricultural water use in western Nevada to evaluate basin‐wide water budgets. Each METRIC evapotranspiration (ET) estimate must be calibrated by a trained user, which requires some iterative time investment and results in variation in ET estimates between users. An automated calibration algorithm for the METRIC model was designed to generate ET estimates comparable to those from trained users by mimicking the manual calibration process. Automated calibration allows for rapid generation of METRIC ET estimates with minimal manual intervention, as well as uncertainty and sensitivity analysis of the model. The variation in ET estimates generated by the automated calibration algorithm was found to be similar to the variation in manual ET estimates. Results indicate that uncertainty was highest for fields with low ET levels and lowest for fields with high ET levels, with a seasonal mean uncertainty of approximately 5% for all fields. In addition, in a blind comparison, automated daily and seasonal ET estimates compared well with flux tower measurement ET data at multiple sites. Automated methods can generate first‐order ET estimates that are similar to time intensive manual efforts with less time investment.     

Experimental investigation on convection heat transfer characteristics of flat plate under environmental wind condition

By establishing wind tunnel and employing electrical heating method, the heat transfer characteristics of flat plate were investigated under environmental wind condition. Both the uniform and linear heat flux boundary conditions were adopted for comparison purpose. Besides, the impacts of heat flux q_wm, tilt angle α and wind incidence angle θ on heat transfer were explored in cases of windward and leeward facing positions. The local convection heat transfer coefficient h_cx and average convection Nusselt number Nu_cm were obtained. The results show that, when α is small, there are two maximum values of h_cx under linear heat flux boundary condition. As for Nu_cm, differences between the two boundary conditions seem indistinctive. At windward orientation, Nu_cm is not sensitive to α. While for leeward orientation, Nu_cm has a rapid decline progress with the increasing α, which indicates the tilt angle α is an important parameter to affect the heat transfer of plate. Finally, new correlations of Nu_cm have been developed, which were proven effective in engineering applications.     

A wind resource assessment around large mountain masses: The speed-up effect

The wind potential around an intensely mountainous area has been studied and an experimental analysis is presented the output of which could be used in the wind farm planning procedure aiming at maximization of the wind power production output of an area. The wind speed of a chosen site of Central Greece was studied based on field measurements around a large mountainous area of Central Greece. Understanding flow in the foothills and the wider area of the mountains is of great importance for estimating wind resource in rough terrain. In this article, special focus was given to the speed-up effect and forced air flow around mountainous masses.     

IMPACT OF VARIED DATA RESOLUTION ON HYDRAULIC MODELING AND FLOODPLAIN DELINEATION1

ABSTRACT: Current data collection technologies such as light detection and ranging (LIDAR) produce dense digital terrain data that result in more accurate digital terrain models (DTMs) for engineering applications. However, such data are redundant and often cumbersome for hydrologic and hydraulic modeling purposes. Data filtering provides a means of eliminating redundant points and facilitates model preparation. This paper demonstrates the impact of varied data resolution on a case study completed for a 2.3 mi² area with mild slopes (about 001 ft/ft) along Leith Creek near Laurinburg, North Carolina. For the original data set and seven filtered data sets, filtering induced changes in elevation, area, and hydraulic radius were determined for 10 water depths at 23 cross sections. Water surface elevations resulting from HEC‐RAS (Hydrologic Engineering Center‐River Analysis System) models for each data set were then compared. A hydraulic model sensitivity analysis was also conducted to compare filtering error to error introduced by variation in flow rates and roughness values. Finally, automated floodplain delineation was performed for each filter level based on the computed hydraulic model results and the filtered LIDAR elevations. Data filtering results indicate that significant time savings are achieved throughout the modeling process and that filtering to four degrees can be performed without compromising cross‐sectional geometry, hydraulic model results, or floodplain delineation results.     

The effects of sinusoidal leading edge of turbine blades on the power coefficient of horizontal-axis wind turbine (HAWT)

In order to improve the aerodynamic performance of horizontal-axis wind turbine (HAWT), a sinusoidal shape is applied to turbine blade. In this study, four types of modified blades were chosen based on variations in amplitude and wavelength of protuberance along the leading edge. Compared with the baseline model, the power coefficients (C_p) of HAWT with modified blades were improved, especially at low tip speed ratios. At low wind speed (V = 6 m/s), blades with short wavelength obtain significant improvement in C_p compared with the baseline model. As wind speed increases, this improvement decreases. In addition, turbine blade with large amplitude and long wavelength obtains better C_p values at higher wind speeds than lower ones, which have a great potential to be more superior at relatively higher wind speeds.     

Blade thickness effect on the aerodynamic performance of an asymmetric NACA six series blade vertical axis wind turbine in low wind speed

ABSTRACT

Vertical axis wind turbine (VAWT) is an economic and widely used energy converter for converting wind energy into useful form of energy, like mechanical and electrical energy. For efficient energy conversion in low wind speed and to have improved power coefficient of asymmetric blade VAWT, selection of optimum blade thickness is needed thus entailing its detailed investigation with respect to different operating wind speed conditions. Present study methodically explores the impact of thickness to chord (t/c) ratio on aerodynamic performance of a three bladed asymmetrical blade H-Darrieus VAWT at different low wind speed conditions by using 2D unsteady CFD simulations. The optimal t/c is obtained on the basis of maximum power coefficient and average moment coefficient of the turbine. The aerodynamic performance curves are obtained at different operating and t/c conditions and the performance insights are corroborated with the findings from the flow physics study to come to some concrete conclusions on the effects of the thickness to chord ratio. The present study identifies large blade curvature to create a large diverging passage on the blade suction surface as the prominent reason for aerodynamic performance drop at a high t/c ratio.     

Global Daily Reference Evapotranspiration Modeling and Evaluation1

Abstract: Accurate and reliable evapotranspiration (ET) datasets are crucial in regional water and energy balance studies. Due to the complex instrumentation requirements, actual ET values are generally estimated from reference ET values by adjustment factors using coefficients for water stress and vegetation conditions, commonly referred to as crop coefficients. Until recently, the modeling of reference ET has been solely based on important weather variables collected from weather stations that are generally located in selected agro‐climatic locations. Since 2001, the National Oceanic and Atmospheric Administration’s Global Data Assimilation System (GDAS) has been producing six‐hourly climate parameter datasets that are used to calculate daily reference ET for the whole globe at 1‐degree spatial resolution. The U.S. Geological Survey Center for Earth Resources Observation and Science has been producing daily reference ET (ETo) since 2001, and it has been used on a variety of operational hydrological models for drought and streamflow monitoring all over the world. With the increasing availability of local station‐based reference ET estimates, we evaluated the GDAS‐based reference ET estimates using data from the California Irrigation Management Information System (CIMIS). Daily CIMIS reference ET estimates from 85 stations were compared with GDAS‐based reference ET at different spatial and temporal scales using five‐year daily data from 2002 through 2006. Despite the large difference in spatial scale (point vs. ～100 km grid cell) between the two datasets, the correlations between station‐based ET and GDAS‐ET were very high, exceeding 0.97 on a daily basis to more than 0.99 on time scales of more than 10 days. Both the temporal and spatial correspondences in trend/pattern and magnitudes between the two datasets were satisfactory, suggesting the reliability of using GDAS parameter‐based reference ET for regional water and energy balance studies in many parts of the world. While the study revealed the potential of GDAS ETo for large‐scale hydrological applications, site‐specific use of GDAS ETo in complex hydro‐climatic regions such as coastal areas and rugged terrain may require the application of bias correction and/or disaggregation of the GDAS ETo using downscaling techniques.     

A comparative analysis of wind resource parameters using WAsP and windPRO

This paper investigates the accuracy of the wind resource estimation for a site in a central India region using a latest licensed version of WAsP 11 and windPRO 3.1. Whole one year measured met mast wind data has been taken using anemometer and wind vane at 10 m and 25 m height, respectively above ground level. The digitized elevation and roughness model of the corresponding site shows the roughness class 4 (roughness length 1.2525 m). The wind data has been extrapolated up to 80 m height by using power and log law models which provide the power density near about 120 W/m². As per the micro sitting guidelines for the virtual wind farm installation 5D X 7D mapping has been selected which Indicates the total power output by installing 8 Vestas V-90 1.8 MW wind turbine from WAsP is 31.561 GWh and from windPRO is 28.083 GWh.     

Best location of the small wind turbine around a single high-rise building

This research is a three-dimensional investigation about the aerodynamic interaction between the wind flow and a single high-rise building. In order to find location(s) with high potential of velocity around the building, a wide variety of wind speeds ranging from 2 to 10 m/s is studied. On the other hand, a high-rise building with the ratio of height to width of H/W = 3 is considered. Computations are performed numerically by means of the finite volume approach. Several results are obtained in the present numerical study. For example, it is found that due to wind-structure vertical interaction, locations with enhanced velocities are developed on the building roof in which the rate of this enhancement increases with increasing the wind speed. In addition, over the building, “lines C and D” are realized as the best locations having high power potentials and low turbulence intensities. In addition, lateral wind-structure interaction revealed that for all wind speeds, location of L/W = 0.5 is the best for the small wind turbine installation.     

EVAPOTRANSPIRATION MEASUREMENT AND ESTIMATION OF THREE WETLAND ENVIRONMENTS IN THE UPPER ST. JOHNS RIVER BASIN,FLORIDA1

ABSTRACT: Accurate estimates of evapotranspiration from areas dominated by wetland vegetation are needed in the water budget of the Upper St. Johns River Basin. However, local data on evapotranspiration rates, especially in wetland environments, were lacking in the project area. In response to this need, the St. Johns River Water Management District collected evapotranspiration field data in Fort Drum Marsh Conservation Area over the period 1996 through 1999. Three large lysimeters were installed to measure the evapotranspiration from different wetland environments: sawgrass (Cladium jamaicense), cattail (Typha domingensis), and open water. In addition, pan evaporation was measured with a standard class “A” pan. Concurrently, meteorological data including rainfall, solar radiation, wind speed, relative humidity, air temperature, and atmospheric pressure were collected. By comparing computed evapotranspiration rates with those measured in the lysimeters, parameters in the Penman‐Monteith, the Priestley‐Taylor, and Reference‐ET methods, and evaporation pan coefficients were estimated for monthly and seasonal cycles. The results from the data collected in this study show that mean monthly evapotranspiration rates, computed by the different methods, are relatively close. From a practical point of view, results indicate that the evaporation pan can be used equally well as the more complex and data‐intensive methods. This paper presents the measured evapotranspiration rates, evaporation pan coefficients, and the estimated parameter values for three different methods to compute evapotranspiration in the project area. Since local data on evaporation are often scarce or lacking, this information may be useful to watershed hydrologists for practical application in other project regions.     

MODELING RUNOFF FROM VARIABLE SOURCE AREAS IN HUMID,SHALLOW WATER TABLE ENVIRONMENTS1

ABSTRACT: Variable Source Areas (VSAs) are zones with water saturated soils in forested wetlands fringing streams and creeks. Runoff from these areas is generated by saturation excess after a shallow water table rises and inundates the ground surface. In humid regions, like Florida and the Southeast, VSAs are believed to produce most of the runoff in shallow water table environments. Modeling the spatial extent and temporal fluctuation of a VSA is difficult because the formation of a VSA depends on a number of hydrological and morphological factors like rainfall intensity, soil texture, water table depth, and topographic attributes of the terrain. In this paper, we couple a digital elevation model with a two‐dimensional variable saturation model to illustrate the formation of a VSA at the hillside scale. The topography derived from the digital elevation model forms the upper domain geometry for the two‐dimensional finite element simulations of variable saturated flow. The objectives are: (1) to model the spatial and dynamic fluctuation of a VSA, and (2) to understand the roles of rainfall variability and terrain attributes on the formation of a VSA. Results show that hillsides with shallow water table depths, low saturated hydraulic conductivity, mild slopes, and concave slope curvature were more susceptible to runoff from a variable source. Runoff from a variable source showed little sensitivity to rainfall intensity. In general, landscapes with steep slopes generated a small VSA and a seepage face that vanished rapidly with time. In contrast, flat terrains are more amenable to VSA and retain ground surface inundation for longer periods of time.     

Watershed Evapotranspiration Increased due to Changes in Vegetation Composition and Structure Under a Subtropical Climate1

Abstract: Natural forests in southern China have been severely logged due to high human demand for timber, food, and fuels during the past century, but are recovering in the past decade. The objective of this study was to investigate how vegetation cover changes in composition and structure affected the water budgets of a 9.6‐km² Dakeng watershed located in a humid subtropical mountainous region in southern China. We analyzed 27 years (i.e., 1967‐1993) of streamflow and climate data and associated vegetation cover change in the watershed. Land use/land cover census and Normalized Difference of Vegetation Index (NDVI) data derived from remote sensing were used to construct historic land cover change patterns. We found that over the period of record, annual streamflow (Q) and runoff/precipitation ratio did not change significantly, nor did the climatic variables, including air temperature, Hamon’s potential evapotranspiration (ET), pan evaporation, sunshine hours, and radiation. However, annual ET estimated as the differences between P and Q showed a statistically significant increasing trend. Overall, the NDVI of the watershed had a significant increasing trend in the peak spring growing season. This study concluded that watershed ecosystem ET increased as the vegetation cover shifted from low stock forests to shrub and grasslands that had higher ET rates. A conceptual model was developed for the study watershed to describe the vegetation cover‐streamflow relationships during a 50‐year time frame. This paper highlighted the importance of eco‐physiologically based studies in understanding transitory, nonstationary effects of deforestation or forestation on watershed water balances.     

AN EVALUATION OF PRECIPITATION GAGE DENSITY IN A MOUNTAINOUS TERRAIN1

ABSTRACT: Inter-station analysis was employed to evaluate the adequacy of the precipitation network in topographically complex West Virginia. A 25-year period was determined as the minimum lingth of record needed for relatively stable and fairly accurate estimates of long-term (50-year) precipitation and in frequency analysis. Data from the 83 National Weather Service stations with 25-year records were adjusted for consistency and evaluated separately by zones east (31 stations) and west (52 stations) of the Appalachian divide. Correlation coefficients (r) and average standard errors of estimate were computed for all station pairs within 50 miles distance and 1000 feet elevation difference of each other. The third polynomial equation of inter-station distance eliminated using elevation and land slope as the criteria in network design in this mountainous terrain. A network with (r) = 0.9 estimates annual precipitation with accuracy as great as 5 percent, but requires about 250 additional gages (i.e., about 200 percent of the present density).     

A COMPARISON OF HAND- AND SPLINE-DRAWN PRECIPITATION MAPS FOR MOUNTAINOUS MONTANA1

ABSTRACT: Average annual precipitation for the period 1961–1990 was estimated for a mountainous region in Montana with a Laplacian thin-plate spline (ANUSPLIN) and compared to a hand-drawn map. Input data included latitude, longitude, and elevation from a three-arc-second U.S. Geological Survey Digital Elevation Model of the Bozeman and Billings 1 × 2 topographic quadrangles and precipitation data at 96 stations. The two maps are similar in appearance. Digital comparison of the two maps with ARC/INFO's Grid tools shows that mean annual precipitation for the hand-contoured map is 22.9 inches and for the ANUSPLIN map is 23.7 inches. Of the 5,760,000 cells, 53 percent showed no difference between ANUSPLIN and hand-drawn maps; 19 percent showed a two-inch difference, and 28 percent showed more than 2 inches difference. Input data and model output at the same location are not different (standard deviation 1.77, p-value 0.76). Hand-drawn maps show two inches more precipitation during the 1961–90 period than during the 1941–1970 period. Similarly, measured data at 73 sites for the period 1961–1990 are on average 2.4 inches higher than the same stations during the 1941–1970 period. The difference is significant (p-value > 0.0001).     

A System Dynamics Model for Conjunctive Management of Water Resources in the Snake River Basin

The Pacific Northwest is expected to witness changes in temperature and precipitation due to climate change. In this study, we enhance the Snake River Planning Model (SRPM) by modeling the feedback loop between incidental recharge and surface water supply resulting from surface water and groundwater extraction for irrigation and provide a case study involving climate change impacts and management scenarios. The new System Dynamics‐Snake River Planning Model (SD‐SRPM) is calibrated to flow at Box Canyon Springs located along a major outlet of the East Snake Plain Aquifer. A calibration of the model to flow at Box Canyon Springs, based on historic diversions (1950‐1995) resulted in an r² value of 0.74 and a validation (1996‐2005) r² value of 0.60. After adding irrigation entities to the model an r² value of 0.91, 0.88, and 0.87 were maintained for modeled vs. observed (1991‐2005) end‐of‐month reservoir content in Jackson Lake, Palisades, and American Falls, the three largest irrigation reservoirs in the system. The scenarios that compared the impacts of climate change were based on ensemble mean precipitation change scenarios and estimated changes to crop evapotranspiration (ET). Increased ET, despite increased precipitation, generally increased surface water shortages and discharge of springs. This study highlights the need to develop and implement models that integrate the human‐natural system to understand the impacts of climate change.     

Estimating turf pesticide volatilization from simple evapotranspiration models

A previously developed model by Haith et al. (2002) related pesticide volatilization from turf to evapotranspiration (ET) by scaling factors determined from vapor pressures and heats of vaporization. Although the model provided volatilization estimates that compared well with field measurements, it relied on the Penman ET equation, requiring hourly temperature, wind speed, and solar radiation data, none of which are routinely available at field sites. The current study determined that the volatilization model works equally well with a simpler ET equation requiring only daily temperatures. Three daily temperature-based ET models were evaluated as vehicles for estimating pesticide volatilization from turf: Hamon, Hargreaves-Samani, and a modified Priestley-Taylor. When compared with field volatilization measurements for eight pesticides, volatilization estimates produced from the Hargreaves-Samani model most closely approximated both the field observations and the previous estimates based on the more data-intensive Penman model. Mean estimated volatilization exceeded mean observations by 15% and the coefficient of variation (R2) between estimates and observations was 0.65. The comparable values based on Penman ET were 17% and 0.63, respectively.     

Predictors of Ips confusus Outbreaks During a Record Drought in Southwestern USA: Implications for Monitoring and Management

In many ecosystems the effects of disturbance can be cryptic and disturbance may vary in subtle spatiotemporal ways. For instance, we know that bark beetle outbreaks are more frequent in temperate forests during droughts; however, we have little idea about why they occur in some locations and not others. Understanding biotic and abiotic factors promoting bark beetle outbreaks can be critical to predicting and responding to pest outbreaks. Here we address the environmental factors which are associated with Ips confusus outbreaks during the 2002 widespread drought within the distribution range of pinyon pine woodlands in Arizona. We used univariate statistics to test if whether tree characteristics, other herbivores, stand properties, soil type, wind, and topography were associated with I. confusus outbreak, and logistic regression to create a predictive model for the outbreaks. We found that I. confusus attacks occur in low elevation stands on steeper slopes, where favorable winds for I. confusus dispersion occur. I. confusus select larger trees, in high density stands with understory shrubs that exhibit phenotypic traits characteristic of resistance to stem-boring moths. The model was highly accurate, and explained 95% of the variability in occurrence (98% of the absences and 95% of the presences). Accurate prediction of the impacts of disturbance allow us to anticipate, minimize or mitigate for and eventually counteract its effects, especially those affecting diversity and ecosystem function. Identification of outbreak risk areas can guide regional and national management towards the reduction of infestation risk and enhancing conservation of pinyon-juniper woodlands.     

Integrating limnological characteristics of high mountain lakes into the landscape of a natural area

A general conceptual watershed-lake model of the complex interactions among climatic conditions, watershed location and characteristics, lake morphology, and fish predation was used to evaluate limnological characteristics of high mountain lakes. Our main hypothesis was that decreasing elevation in mountainous terrain corresponds to an increase in diversity of watershed size and lake area, depth, temperature, nutrient concentrations, and productivity. A second hypothesis was that watershed location and aspect relative to climatic gradients within mountainous terrain influences the limnological characteristics of the lakes. We evaluated these hypotheses by examining watershed location, aspect and size; lake morphology; water quality; and phytoplankton and zooplankton community characteristics among high mountain forest and subalpine lakes in Mount Rainier National Park. Although many of the comparisons between all forest and subalpine lakes were statistically insignificant, the results revealed trends that were consistent with our hypotheses. The forest lake group included more lakes with larger watersheds, larger surface areas, greater depths, higher concentrations of nutrients, and higher algal biovolumes than did the group of subalpine lakes. Deep lakes, which were mostly of the forest lake type, exhibited thermal stratification and relatively high values of some of the water-quality variables near the lake bottoms. However, the highest near-surface water temperatures and phytoplankton densities and the taxonomic structures of the phytoplankton and zooplankton assemblages were more closely related to geographical location, which corresponded to a west-east climate gradient in the park, than to lake type. Some crustacean and rotifer taxa, however, were limited in distribution by lake type. Fish predation did not appear to play an important role in the structure of the crustacean zooplankton communities at the genus level with the exception of Mowich Lake, where crustacean taxa were absent from the zooplankton community. This was the only lake inhabited by a true zooplanktivourous species of fish.