Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Journal of Time Series Analysis”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5Linked to original sources
USGS2019

Insights into pāhoehoe lava emplacement using visible and thermal structure-from-motion photogrammetry

We present the evolution over 3 months of a 2016–2017 pāhoehoe flow at Kīlauea as it changed from a narrow sheet flow into a compound lava field fed by a stable system of tubes. The portion of the flow located on Kīlauea's coastal plain was characterized using helicopter‐based visible and thermal structure‐from‐motion photogrammetry to construct a series of georeferenced digital surface models and thermal maps on eight different days. Results reveal key influences on the emplacement and evolution of such long‐lived pāhoehoe flows. This region of the flow grew by ~12 × 10 6 m 3 with a near‐constant time‐average discharge rate of 1.2–2.7 m 3 /s. The development of two tube systems is captured and shows an initial nascent tube enhanced by a narrow topographic confinement, which later inflated and created a topographic inversion that modulated the emplacement of a second flow lobe with its own tube system. The analysis of breakouts at various stages of the field's life suggests that the evolution of the thermal and morphological properties of the flow surface reflect its maturity. Thermal properties of breakouts were used to expand the empirical relationship of breakout cooling to longer timescales. This study contributes to the long‐term development and validation of more accurate predictive models for pāhoehoe, required during the management of long‐lasting lava flow crises in Hawai'i and elsewhere.

Hawaii↗
USGS2010

Land cover characterization and mapping of continental southeast Asia using multi-resolution satellite sensor data

Land use/land cover change, particularly that of tropical deforestation and forest degradation, has been occurring at an unprecedented rate and scale in Southeast Asia. The rapid rate of economic development, demographics and poverty are believed to be the underlying forces responsible for the change. Accurate and up-to-date information to support the above statement is, however, not available. The available data, if any, are outdated and are not comparable for various technical reasons. Time series analysis of land cover change and the identification of the driving forces responsible for these changes are needed for the sustainable management of natural resources and also for projecting future land cover trajectories. We analysed the multi-temporal and multi-seasonal NOAA Advanced Very High Resolution Radiometer (AVHRR) satellite data of 1985/86 and 1992 to (1) prepare historical land cover maps and (2) to identify areas undergoing major land cover transformations (called ‘hot spots’). The identified ‘hot spot’ areas were investigated in detail using high-resolution satellite sensor data such as Landsat and SPOT supplemented by intensive field surveys. Shifting cultivation, intensification of agricultural activities and change of cropping patterns, and conversion of forest to agricultural land were found to be the principal reasons for land use/land cover change in the Oudomxay province of Lao PDR, the Mekong Delta of Vietnam and the Loei province of Thailand, respectively. Moreover, typical land use/land cover change patterns of the ‘hot spot’ areas were also examined. In addition, we developed an operational methodology for land use/land cover change analysis at the national level with the help of national remote sensing institutions.

International Journal of Remote Sensing↗
USGS2023

Modeling groundwater-level responses to multiple stresses using transfer-function models and wavelet analysis in a coastal aquifer system

In coastal aquifers, dynamic stresses such as climate forcings, groundwater withdrawals, and ocean tidal fluctuations cause nonlinear responses to groundwater levels. Such responses to the stresses impact groundwater resources and related flooding and infrastructure risks at multiple scales. We used time-series models such as transfer-function models and wavelet analysis to quantify the relative contribution of these stresses to groundwater-level fluctuation in wells from the unconfined and confined aquifers in an Atlantic coastal aquifer. Climate forcings, such as precipitation and temperature, explained most of the groundwater-level variation for wells in the unconfined aquifer, whereas groundwater withdrawals were the dominant driver of groundwater levels for wells in the confined aquifer. The impact of groundwater withdrawals also was detected in several wells in the unconfined aquifer. Although the influence of ocean tides on groundwater levels commonly is observed in coastal aquifers, we found that daily groundwater withdrawals can obscure the semi-diurnal coherence signal of the two series. The magnitude of groundwater-level fluctuation that could be explained solely by tides was minor compared to that explained by climate or withdrawal stresses. Transfer-function modeling showed seasonal withdrawals from wells in confined aquifers had a significant, yet heterogeneous influence on groundwater levels in coastal aquifers, which highlights climate and withdrawals as key compounding stresses in coastal hydrology. This study demonstrates the value of time-series approaches to advance characterization of groundwater systems in areas with limited hydrogeologic parameter information.

Virginia↗
USGS2021

Incorporating water quality analysis into navigation assessments as demonstrated in the Mississippi River Basin

A description of historical and ambient water quality conditions is often required as part of navigational studies. This paper describes a series of tools developed by the USGS that can aid navigation managers in developing water quality assessments. The tools use R, a statistical software program, and provide methods to retrieve historical streamflow and water quality data, summarize observations, model concentrations and fluxes, and estimate seasonal, annual, and decadal trends. The utility of these tools is demonstrated by providing an analysis of the seasonal variability and long-term trends of nitrate plus nitrite, orthophosphate, and suspended sediment concentrations and fluxes at nine sites in the Mississippi River Basin. Trends in annual mean concentration and flux showed fairly stable nitrate plus nitrite at most of the nine sites, with increases in the Upper Mississippi and Missouri Rivers and decreases on the Illinois River over a 40-year period beginning in 1980. Orthophosphate concentration or flux increased at almost all sites over a similar time period. Conversely, a concurrent steady decline in suspended sediment concentrations and fluxes was noted at sites throughout the basin.

Mississippi River Basin↗
USGS2017

New techniques to measure cliff change from historical oblique aerial photographs and structure-from-motion photogrammetry

Oblique aerial photograph surveys are commonly used to document coastal landscapes. Here it is shown that adequate overlap may exist in these photographic records to develop topographic models with Structure-from-Motion (SfM) photogrammetric techniques. Using photographs of Fort Funston, California, from the California Coastal Records Project, imagery were combined with ground control points in a four-dimensional analysis that produced topographic point clouds of the study area’s cliffs for 5 years spanning 2002 to 2010. Uncertainty was assessed by comparing point clouds with airborne LIDAR data, and these uncertainties were related to the number and spatial distribution of ground control points used in the SfM analyses. With six or more ground control points, the root mean squared errors between the SfM and LIDAR data were less than 0.30 m (minimum 1⁄4 0.18 m), and the mean systematic error was less than 0.10 m. The SfM results had several benefits over traditional airborne LIDAR in that they included point coverage on vertical- to-overhanging sections of the cliff and resulted in 10–100 times greater point densities. Time series of the SfM results revealed topographic changes, including landslides, rock falls, and the erosion of landslide talus along the Fort Funston beach. Thus, it was concluded that SfM photogrammetric techniques with historical oblique photographs allow for the extraction of useful quantitative information for mapping coastal topography and measuring coastal change. The new techniques presented here are likely applicable to many photograph collections and problems in the earth sciences.

Journal of Coastal Research↗
USGS2008

Joint variability of global runoff and global sea surface temperatures

Global land surface runoff and sea surface temperatures (SST) are analyzed to identify the primary modes of variability of these hydroclimatic data for the period 1905-2002. A monthly water-balance model first is used with global monthly temperature and precipitation data to compute time series of annual gridded runoff for the analysis period. The annual runoff time series data are combined with gridded annual sea surface temperature data, and the combined dataset is subjected to a principal components analysis (PCA) to identify the primary modes of variability. The first three components from the PCA explain 29% of the total variability in the combined runoff/SST dataset. The first component explains 15% of the total variance and primarily represents long-term trends in the data. The long-term trends in SSTs are evident as warming in all of the oceans. The associated long-term trends in runoff suggest increasing flows for parts of North America, South America, Eurasia, and Australia; decreasing runoff is most notable in western Africa. The second principal component explains 9% of the total variance and reflects variability of the El Ni??o-Southern Oscillation (ENSO) and its associated influence on global annual runoff patterns. The third component explains 5% of the total variance and indicates a response of global annual runoff to variability in North Aflantic SSTs. The association between runoff and North Atlantic SSTs may explain an apparent steplike change in runoff that occurred around 1970 for a number of continental regions.

Journal of Hydrometeorology↗
USGS2010

Tidal calibration of Plate Boundary Observatory borehole strainmeters: Roles of vertical and shear coupling

A multicomponent borehole strainmeter directly measures changes in the diameter of its cylindrical housing at several azimuths. To transform these measurements to formation strains requires a calibration matrix, which must be estimated by analyzing the installed strainmeter's response to known strains. Typically, theoretical calculations of Earth tidal strains serve as the known strains. This paper carries out such an analysis for 12 Plate Boundary Observatory (PBO) borehole strainmeters, postulating that each of the strainmeters' four gauges responds ("couples") to all three horizontal components of the formation strain tensor, as well as to vertical strain. Orientation corrections are also estimated. The fourth extensometer in each PBO strainmeter provides redundant information used to reduce the chance that coupling coefficients could be misleadingly fit to inappropriate theoretical tides. Satisfactory fits between observed and theoretically calculated tides were obtained for three PBO strainmeters in California, where the calculated tides are corroborated by other instrumentation, as well as for six strainmeters in Oregon and Washington, where no other instruments have ever recorded Earth tidal strain. Several strainmeters have unexpectedly large coupling coefficients for vertical strain, which increases the strainmeter's response to atmospheric pressure. Vertical coupling diminishes, or even changes the sign of, the apparent response to areal strain caused by Earth tides or deep Earth processes because near the free surface, vertical strains are opposite in sign to areal strain. Vertical coupling does not impair the shear strain response, however. PBO borehole strainmeters can provide calibrated shear strain time series of transient strain associated with tectonic or magmatic processes.

Journal of Geophysical Research↗
USGS2018

Attribution analysis of the Ethiopian drought of 2015

In northern and central Ethiopia, 2015 was a very dry year. Rainfall was only from one-half to three-quarters of the usual amount, with both the “belg” (February–May) and “kiremt” rains (June–September) affected. The timing of the rains that did fall was also erratic. Many crops failed, causing food shortages for many millions of people. The role of climate change in the probability of a drought like this is investigated, focusing on the large-scale precipitation deficit in February–September 2015 in northern and central Ethiopia. Using a gridded analysis that combines station data with satellite observations, it is estimated that the return period of this drought was more than 60 years (lower bound 95% confidence interval), with a most likely value of several hundred years. No trend is detected in the observations, but the large natural variability and short time series means large trends could go undetected in the observations. Two out of three large climate model ensembles that simulated rainfall reasonably well show no trend while the third shows an increased probability of drought. Taking the model spread into account the drought still cannot be clearly attributed to anthropogenic climate change, with the 95% confidence interval ranging from a probability decrease between preindustrial and today of a factor of 0.3 and an increase of a factor of 5 for a drought like this one or worse. A soil moisture dataset also shows a nonsignificant drying trend. According to ENSO correlations in the observations, the strong 2015 El Niño did increase the severity of the drought.

Journal of Climate↗
USGS2017

Evaluation of remote-sensing techniques to measure decadal-scale changes of Hofsjokull ice cap, Iceland

Dynamic surficial changes and changes in the position of the firn line and the areal extent of Hofsjökull ice cap, Iceland, were studied through analysis of a time series (1973–98) of synthetic-aperture radar (SAR) and Landsat data. A digital elevation model of Hofsjökull, which was constructed using SAR interferometry, was used to plot the SAR backscatter coefficient (σ°) vs elevation and air temperature along transects across the ice cap. Seasonal and daily σ° patterns are caused by freezing or thawing of the ice-cap surface, and abrupt changes in σ° are noted when the air temperature ranges from ∼−5° to 0°C. Late-summer 1997 σ° (SAR) and reflectance (Landsat) boundaries agree and appear to be coincident with the firn line and a SAR σ° boundary that can be seen in the January 1998 SAR image. In January 1994 through 1998, the elevation of this σ° boundary on the ice capwas quite stable, ranging from 1000 to 1300 m, while the equilibrium-line altitude, as measured on the ground, varied considerably. Thus the equilibrium line may be obscured by firn from previous years. Techniques are established to measure long-term changes in the elevation of the firn line and changes in the position of the ice margin.

Hofsjokull ice cap↗
USGS2012

Analysis of time series of glacier speed: Columbia Glacier, Alaska

During the summers of 1984 and 1985, laser measurements were made of the distance from a reference location to markers on the surface of the lower reach of Columbia Glacier, Alaska. The data have numerous gaps, mostly because of inclement weather. The laser measurements were corrected for variations in atmospheric temperature and pressure and for the trajectory of the marker. The marker speed was calculated by using a combination of cubic splines and digital filters to arrive at a data set with 1-hour time intervals. The speed varies from 7 to 15 m/d and has three noteworthy components: (1) a low-frequency perturbation in speed with a time scale of days related to increased precipitation, (2) semidiurnal and diurnal variations related to sea tides, and (3) diurnal variations related to glacier surface melt. The low frequency and tidal period variations are separated by low-pass filtering and subtracting the filtered data from the original data. The variations from melt runoff and tides have the same frequency range and are examined using harmonic analysis and multiple regression with the various forcing functions. The high-frequency portion of the ice speed signal is dominated by tidal influences and to a lesser extent by meltwater influences. The low-frequency portion of the signal is dominated by the effects of precipitation.

Journal of Geophysical Research Solid Earth↗
USGS1994

Regional Aquifer-System Analysis Program of the U.S. Geological Survey, 1978-1992

The major ground-water systems of the United States have been investigated by the U.S. Geological Survey (USGS) through its Regional Aquifer-System Analysis (RASA) Program. During the first 15 years of the program (1978-92), 25 regional aquifer systems, including the most heavily pumped aquifers in the Nation, were intensively studied. As of mid-1992, 18 of the regional aquifer studies are completed or nearly so; 7 of the regional aquifer studies are ongoing, and compilation of a national ground-water atlas is in progress. This report summarizes the status of each RASA study and briefly describes the hydrology of the 25 regional aquifer systems. Important study results and examples of applications of study results are presented for some of the completed RASA investigations. The major contributions of the RASA Program are (1) assembly of data from numerous local studies and long-time data collections into systematic regional data bases; (2) comprehensive descriptions of the geologic, hydrologic, and geochemical characteristics of the regional aquifer systems; and (3) an understanding of how the regional ground-water-flow systems function under natural (predevelopment) and current (developed) conditions. To provide the comprehensive system descriptions, many of the RASA studies present, for the first time, maps depicting the hydrogeologic frameworks, water chemistry, potentiometric surfaces, and other aspects of entire regional aquifer systems. To provide an understanding of how the flow systems function, several of the completed RASA studies provide, for the first time, hydrologic budgets for both predevelopment and developed conditions. The results of the RASA Program are contained in nearly 900 reports published by the USGS, as well as various State and local agencies, and in articles published in scientific journals. A series of U.S. Geological Survey Professional Papers are being published to summarize and synthesize the results of each RASA study.

Circular↗
USGS2012

Age estimates and uplift rates for late Pleistocene marine terraces: Southern Oregon portion of the Cascadia forearc

Interest in the Cascadia subduction zone has increased because recent investigations have suggested that slip along plates at certain types of convergent margins is characteristically accompanied by large earthquakes. In addition, other investigations have suggested that convergent margins can be broadly classified by the magnitude of their uplift rates. The authors generated new uranium series, amino acid, and stable isotope data for southern Oregon marine terrace fossils. These data, along with terrace elevations and two alternative estimates of sea level at the time of terrace formation, allow one to determine terrace ages and uplift rates. Uranium series analysis of fossil coral yields an age of 83 {plus minus} 5 ka for the Whisky Run terrace at Coquille Point in Bandon, Oregon. A combination of amino acid and oxygen isotope data suggest ages of about 80 and 105 ka for the lowest two terraces at Cape Blanco. These ages indicate uplift rates of 0.45-1.05 and 0.81-1.49 m/kyr for Coquille Point and Cape Blanco, respectively. In order to assess the utility of the southern Oregon uplift rates for predicting the behavior of the Cascadia subduction zone, the authors compared late Quaternary uplift rates derived from terrace data from subduction zones around the world. On the basis of this comparison the southern Oregon rates of vertical deformation are not usually high or low. Furthermore, late Quaternary uplift rates show little relationship to the type of convergent margin. In the case of the southern Oregon coast, variability in uplift rate probably reflects local structures in the overriding plate, and the rate of uplift cannot be used as a simple index of the potential for great earthquakes along the southern Cascadia subduction zone.

California, Oregon↗
USGS2021

Multi-year hydroclimatic droughts and pluvials across the conterminous United States

Time series of water‐year runoff for 2,109 hydrologic units (HUs) across the conterminous United States (CONUS) for the 1900 through 2014 period were used to identify drought and pluvial (i.e., wet) periods. Characteristics of the drought and pluvial events including frequency, duration, and severity were examined and compared. Additionally, a similar analysis was performed using gridded tree‐ring reconstructions of the Palmer Drought Severity Index (PDSI) for the period 1475 through 2005 to place the drought and pluvial characteristics determined using water‐year runoff for 1900 through 2014 in the context of multi‐century climate variability. The temporal and spatial variability of droughts and pluvials determined using runoff for the 1900 through 2014 period indicated that most drought events in the CONUS occurred before about 1970, whereas most pluvial periods occurred after about 1970. This change in the frequencies of drought and pluvial events around 1970 was largely related to an increase in fall (October through December) precipitation across much of the central United States. Also, the duration and severity of droughts and pluvials identified using runoff for the 1900 through 2014 period generally were not significantly different from the drought and pluvial characteristics identified using the PDSI for the 1475 through 2005 period.

International Journal of Climatology↗
USGS1988

Pumping tests in nonuniform aquifers - The radially symmetric case

Traditionally, pumping-test-analysis methodology has been limited to applications involving aquifers whose properties are assumed uniform in space. This work attempts to assess the applicability of analytical methodology to a broader class of units with spatially varying properties. An examination of flow behavior in a simple configuration consisting of pumping from the center of a circular disk embedded in a matrix of differing properties is the basis for this investigation. A solution describing flow in this configuration is obtained through Laplace-transform techniques using analytical and numerical inversion schemes. Approaches for the calculation of flow properties in conditions that can be roughly represented by this simple configuration are proposed. Possible applications include a wide variety of geologic structures, as well as the case of a well skin resulting from drilling or development. Of more importance than the specifics of these techniques for analysis of water-level responses is the insight into flow behavior during a pumping test that is provided by the large-time form of the derived solution. The solution reveals that drawdown during a pumping test can be considered to consist of two components that are dependent and independent of near-well properties, respectively. Such an interpretation of pumping-test drawdown allows some general conclusions to be drawn concerning the relationship between parameters calculated using analytical approaches based on curve-matching and those calculated using approaches based on the slope of a semilog straight line plot. The infinite-series truncation that underlies the semilog analytical approaches is shown to remove further contributions of near-well material to total drawdown. In addition, the semilog distance-drawdown approach is shown to yield an expression that is equivalent to the Thiem equation. These results allow some general recommendations to be made concerning observation-well placement for pumping tests in nonuniform aquifers. The relative diffusivity of material on either side of a discontinuity is shown to be the major factor in controlling flow behavior during the period in which the front of the cone of depression is moving across the discontinuity. Though resulting from an analysis of flow in an idealized configuration, the insights of this work into flow behavior during a pumping test are applicable to a wide class of nonuniform units.

Journal of Hydrology↗
USGS2020

Maximum entropy derived statistics of sound speed structure in a fine-grained sediment inferred from sparse broadband acoustic measurements on the New England continental shelf

Marginal probability distributions for parameters representing an effective sound-speed structure of a fine-grained sediment are inferred from a data ensemble maximum entropy method that utilizes a sparse spatially distributed set of received pressure time series resulting from multiple explosive sources in a shallow-water ocean environment possessing significant spatial variability of the seabed. A remote sensing seabed acoustics experiment undertaken in March 2017 off the New England Shelf was designed so that multiple independent analyses could infer the statistical properties of the seabed. The current analysis incorporates the measured horizontal variability from interpretations of a subbottom profiling survey of the experimental area. An idealized range- and azimuth-dependent parameterization of the seabed is derived from identification of horizons within the seabed that define multiple sediment layers. A sparse set of explosive charges were deployed on circular tracks with radii of about 2, 4, and 6.5 km with an acoustic array at the center to correlate a set of random measurements to physical acoustic processes that characterize the seabed. The mean values of a surface sound speed ratio and a linear sound speed gradient for the fine-grained sediment layer derived from 12 data samples processed in the 25–275-Hz band provide an estimate of the effective sound-speed structure in a 130-km $^2$ area. The inferred sediment sound speed values are evaluated by predicting measured time series data not used in the statistical inference, and are also compared to historical measurements. Finally, the low-frequency maximum entropy estimate of the sediment sound speed along with physical measurements derived from piston core measurements are utilized to estimate the sediment grain bulk modulus. predictions made by the viscous grain shearing model.

Maine, Vermont, New Hampshire, Massachusetts, Rhod↗
USGS2017

Improved efficiency of maximum likelihood analysis of time series with temporally correlated errors

Most time series of geophysical phenomena have temporally correlated errors. From these measurements, various parameters are estimated. For instance, from geodetic measurements of positions, the rates and changes in rates are often estimated and are used to model tectonic processes. Along with the estimates of the size of the parameters, the error in these parameters needs to be assessed. If temporal correlations are not taken into account, or each observation is assumed to be independent, it is likely that any estimate of the error of these parameters will be too low and the estimated value of the parameter will be biased. Inclusion of better estimates of uncertainties is limited by several factors, including selection of the correct model for the background noise and the computational requirements to estimate the parameters of the selected noise model for cases where there are numerous observations. Here, I address the second problem of computational efficiency using maximum likelihood estimates (MLE). Most geophysical time series have background noise processes that can be represented as a combination of white and power-law noise, 1 / f α "> 1 / f α 1/fα with frequency, f . With missing data, standard spectral techniques involving FFTs are not appropriate. Instead, time domain techniques involving construction and inversion of large data covariance matrices are employed. Bos et al. (J Geod, 2013 . doi: 10.1007/s00190-012-0605-0 ) demonstrate one technique that substantially increases the efficiency of the MLE methods, yet is only an approximate solution for power-law indices >1.0 since they require the data covariance matrix to be Toeplitz. That restriction can be removed by simply forming a data filter that adds noise processes rather than combining them in quadrature. Consequently, the inversion of the data covariance matrix is simplified yet provides robust results for a wider range of power-law indices.

Journal of Geodesy↗
USGS2008

Long-term water quality and biological responses to multiple best management practices in Rock Creek, Idaho

Water quality and macroinvertebrate assemblage data from 1981 to 2005 were assessed to evaluate the water quality and biological responses of a western trout stream to the implementation of multiple best management practices (BMPs) on irrigated cropland. Data from Rock Creek near Twin Falls, Idaho, a long-term monitoring site, were assembled from state and federal sources to provide the evaluation. Seasonal loads of the nonpoint source pollutants suspended sediment (SS), total phosphorus (TP), and nitrate-nitrite (NN) were estimated using a regression model with time-series streamflow data and constituent concentrations. Trends in the macroinvertebrate assemblages were evaluated using a number of biological metrics and nonmetric multidimensional scaling ordination. Regression analysis found significant annual decreases in TP and SS flow-adjusted concentrations during the BMP implementation period from 1983 to 1990 of about 7 and 10%, respectively. These results are coincident with the implementation of multiple BMPs on about 75% of the irrigated cropland in the watershed. Macroinvertebrate assemblages during this time also responded with a change in taxa composition resulting in improved biotic index scores. Taxon specific TP and SS optima, empirically derived from a large national dataset, predicted a decrease in SS concentrations of about 37% (52 to 33 mg/l) and a decrease in TP concentrations of about 50% (0.20 to 0.10 mg/l) from 1981 to 1987. Decreasing trends in TP, SS, and NN pollutant loads were primarily the result of naturally low streamflow conditions during the BMP post-implementation period from 1993 to 2005. Trends in macroinvertebrate responses during 1993 to 2005 were confounded by the introduction of the New Zealand mudsnail (Potamopyrgus antipodarum), which approached densities of 100,000 per m 2 in riffle habitat. The occurrence of this invasive species appears to have caused a major shift in composition and function of the macroinvertebrate assemblages. ?? 2008 American Water Resources Association.

Journal of the American Water Resources Associatio↗
USGS2017

Impact of tile drainage on evapotranspiration in South Dakota, USA, based on high spatiotemporal resolution evapotranspiration time series from a multi-satellite data fusion system

Soil drainage is a widely used agricultural practice in the midwest USA to remove excess soil water to potentially improve the crop yield. Research shows an increasing trend in baseflow and streamflow in the midwest over the last 60 years, which may be related to artificial drainage. Subsurface drainage (i.e., tile) in particular may have strongly contributed to the increase in these flows, because of its extensive use and recent gain in the popularity as a yield-enhancement practice. However, how evapotranspiration (ET) is impacted by tile drainage on a regional level is not well-documented. To explore spatial and temporal ET patterns and their relationship to tile drainage, we applied an energy balance-based multisensor data fusion method to estimate daily 30-m ET over an intensively tile-drained area in South Dakota, USA, from 2005 to 2013. Results suggest that tile drainage slightly decreases the annual cumulative ET, particularly during the early growing season. However, higher mid-season crop water use suppresses the extent of the decrease of the annual cumulative ET that might be anticipated from widespread drainage. The regional water balance analysis during the growing season demonstrates good closure, with the average residual from 2005 to 2012 as low as -3 mm. As an independent check of the simulated ET at the regional scale, the water balance analysis lends additional confidence to the study. The results of this study improve our understanding of the influence of agricultural drainage practices on regional ET, and can affect future decision making regarding tile drainage systems.

IEEE Journal of Selected Topics in Applied Earth O↗