Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Data Series”

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 37 records · Page 2Linked to original sources

Mapping rice-fallow cropland areas for short-season grain legumes intensification in South Asia using MODIS 250 m time-series data

The goal of this study was to map rainfed and irrigated rice-fallow cropland areas across South Asia, using MODIS 250 m time-series data and identify where the farming system may be intensified by the inclusion of a short-season crop during the fallow period. Rice-fallow cropland areas are those areas where rice is grown during the kharif growing season (June–October), followed by a fallow during the rabi season (November–February). These cropland areas are not suitable for growing rabi -season rice due to their high water needs, but are suitable for a short -season (≤3 months), low water-consuming grain legumes such as chickpea ( Cicer arietinum L.), black gram, green gram, and lentils. Intensification (double-cropping) in this manner can improve smallholder farmer’s incomes and soil health via rich nitrogen-fixation legume crops as well as address food security challenges of ballooning populations without having to expand croplands. Several grain legumes, primarily chickpea, are increasingly grown across Asia as a source of income for smallholder farmers and at the same time providing rich and cheap source of protein that can improve the nutritional quality of diets in the region. The suitability of rainfed and irrigated rice-fallow croplands for grain legume cultivation across South Asia were defined by these identifiers: (a) rice crop is grown during the primary ( kharif ) crop growing season or during the north-west monsoon season (June–October); (b) same croplands are left fallow during the second ( rabi ) season or during the south-east monsoon season (November–February); and (c) ability to support low water-consuming, short-growing season (≤3 months) grain legumes (chickpea, black gram, green gram, and lentils) during rabi season. Existing irrigated or rainfed crops such as rice or wheat that were grown during kharif were not considered suitable for growing during the rabi season, because the moisture/water demand of these crops is too high. The study established cropland classes based on the every 16-day 250 m normalized difference vegetation index (NDVI) time series for one year (June 2010–May 2011) of Moderate Resolution Imaging Spectroradiometer (MODIS) data, using spectral matching techniques (SMTs), and extensive field knowledge. Map accuracy was evaluated based on independent ground survey data as well as compared with available sub-national level statistics. The producers’ and users’ accuracies of the cropland fallow classes were between 75% and 82%. The overall accuracy and the kappa coefficient estimated for rice classes were 82% and 0.79, respectively. The analysis estimated approximately 22.3 Mha of suitable rice-fallow areas in South Asia, with 88.3% in India, 0.5% in Pakistan, 1.1% in Sri Lanka, 8.7% in Bangladesh, 1.4% in Nepal, and 0.02% in Bhutan. Decision-makers can target these areas for sustainable intensification of short-duration grain legumes.

International Journal of Digital Earth↗

Phenological classification of the United States: A geographic framework for extending multi-sensor time-series data

This study introduces a new geographic framework, phenological classification, for the conterminous United States based on Moderate Resolution Imaging Spectroradiometer (MODIS) Normalized Difference Vegetation Index (NDVI) time-series data and a digital elevation model. The resulting pheno-class map is comprised of 40 pheno-classes, each having unique phenological and topographic characteristics. Cross-comparison of the pheno-classes with the 2001 National Land Cover Database indicates that the new map contains additional phenological and climate information. The pheno-class framework may be a suitable basis for the development of an Advanced Very High Resolution Radiometer (AVHRR)-MODIS NDVI translation algorithm and for various biogeographic studies.

Remote Sensing↗

Prediction of forest canopy and surface fuels from Lidar and satellite time series data in a bark beetle-affected forest

Wildfire behavior depends on the type, quantity, and condition of fuels, and the effect that bark beetle outbreaks have on fuels is a topic of current research and debate. Remote sensing can provide estimates of fuels across landscapes, although few studies have estimated surface fuels from remote sensing data. Here we predicted and mapped field-measured canopy and surface fuels from light detection and ranging (lidar) and Landsat time series explanatory variables via random forest (RF) modeling across a coniferous montane forest in Colorado, USA, which was affected by mountain pine beetles ( Dendroctonus ponderosae Hopkins) approximately six years prior. We examined relationships between mapped fuels and the severity of tree mortality with correlation tests. RF models explained 59%, 48%, 35%, and 70% of the variation in available canopy fuel, canopy bulk density, canopy base height, and canopy height, respectively (percent root-mean-square error (%RMSE) = 12–54%). Surface fuels were predicted less accurately, with models explaining 24%, 28%, 32%, and 30% of the variation in litter and duff, 1 to 100-h, 1000-h, and total surface fuels, respectively (%RMSE = 37–98%). Fuel metrics were negatively correlated with the severity of tree mortality, except canopy base height, which increased with greater tree mortality. Our results showed how bark beetle-caused tree mortality significantly reduced canopy fuels in our study area. We demonstrated that lidar and Landsat time series data contain substantial information about canopy and surface fuels and can be used for large-scale efforts to monitor and map fuel loads for fire behavior modeling at a landscape scale.

Colorado↗

Characteristics of vegetation phenology over the Alaskan landscape using AVHRR time-series data

Advanced Very High Resolution Radiometer (AVHRR) satellite data were acquired and composited into twice-a-month periods from 1 May 1991 to 15 October 1991 in order to map vegetation characteristics of the Alaskan landscape. Unique spatial and temporal qualities of the AVHRR data provide information that leads to a better understanding of regional biophysical characteristics of vegetation communities and patterns. These data provided synoptic views of the landscape and depicted phenological diversity, temporal vegetation phenology (green-up, peak of green, and senescence), photosynthetic activity, and regional landscape patterns. Products generated from the data included a phenological class map, phenological composite maps (onset, peak, and duration), and photosynthetic activity maps (mean and maximum greenness). The time-series data provide opportunities to study phenological processes at small landscape scales over time periods of weeks, months, and years. Regional patterns identified on some of the maps are unique to specific areas; others correspond to biophysical or ecoregional boundaries. The data provide new insights to landscape processes, ecology, and landscape physiognomy that allow scientists to look at landscapes in ways that were previously difficult to achieve.

Alaska↗

Dynamic rating method for computing discharge and stage from time-series data

Ratings are used for several reasons in water-resources investigations. The simplest rating relates discharge to the stage of a river (the stage-discharge relation). From a pure hydrodynamics perspective, all rivers and streams have some form of hysteresis in the relation between stage and discharge because flow becomes unsteady as a flood wave passes. The stage-discharge relation is unable to represent hysteresis. However, a dynamic rating method can capture hysteresis, which is driven by the variable energy slope of a flood wave. A dynamic rating method called DYNPOUND, which accommodates compact and compound channel geometry, was developed by simplifying the one-dimensional Saint-Venant equations. The DYNPOUND method was developed in the Python programming language and computes discharge from stage and stage from discharge. Stage and discharge time series computed with this dynamic rating method were compared to the U.S. Geological Survey (USGS) published stage and discharge time series. The results from the DYNPOUND method were also compared to in-person field measurements of stage and discharge made at 10 USGS streamgages. DYNPOUND was calibrated for 10 USGS streamgages using published discharge time-series data computed with a simple rating method. The calibration objective was to minimize the mean squared logarithmic error (MSLE) of the DYNPOUND-computed discharge with respect to the discharge time series computed by a simple rating method. For each site, the calibration process also included comparing all field measurements within a selected water year to the corresponding DYNPOUND-computed discharge data points. The MSLE of the DYNPOUND-computed discharge time series for the 10 sites ranged from 8.51×10 −4 to 1.36×10 −1 . For each site, an event-based period was selected to compare the discharge time series computed with the dynamic rating method to discharge field measurements made at the streamgages; the range of MSLE for the 10 DYNPOUND-computed discharge sites was from 4.79×10 −4 to 2.30×10 −2 .

continental United States↗

Investigating lake-area dynamics across a permafrost-thaw spectrum using airborne electromagnetic surveys and remote sensing time-series data in Yukon Flats, Alaska

Lakes in boreal lowlands cycle carbon and supply an important source of freshwater for wildlife and migratory waterfowl. The abundance and distribution of these lakes are supported, in part, by permafrost distribution, which is subject to change. Relationships between permafrost thaw and lake dynamics remain poorly known in most boreal regions. Here, new airborne electromagnetic (AEM) data collected during June 2010 and February 2016 were used to constrain deep permafrost distribution. AEM data were coupled with Landsat-derived lake surface-area data from 1979 through 2011 to inform temporal lake behavior changes in the 35 500- km 2 Yukon Flats ecoregion of Alaska. Together, over 1500 km of AEM data, and roughly 30 years of Landsat data were used to explore processes that drive lake dynamics across a variety of permafrost thaw states not possible in studies conducted with satellite imagery or field measurements alone. Clustered time-series data identified lakes with similar temporal dynamics. Clusters possessed similarities in lake permanence (i.e. ephemeral versus perennial), subsurface permafrost distribution, and proximity to rivers and streams. Of the clustered lakes, ~66% are inferred to have at least intermittent connectivity with other surface-water features, ~19% are inferred to have shallow subsurface connectivity to other surface water features that served as a low-pass filter for hydroclimatic fluctuations, and ~15% appear to be isolated by surrounding permafrost (i.e. no connectivity). Integrated analysis of AEM and Landsat data reveals a progression from relatively synchronous lake dynamics among disconnected lakes in the most spatially continuous, thick permafrost to quite high spatiotemporal heterogeneity in lake behavior among variably-connected lakes in regions with notably less continuous permafrost. Variability can be explained by the preferential development of thawed permeable gravel pathways for lateral water redistribution in this area. The general spatial progression in permafrost thaw state and lake area behavior may be extended to the temporal dimension. However, extensive permafrost thaw, beyond what is currently observed, is expected to promote ubiquitous subsurface connectivity, eventually evolving to a state of increased lake synchronicity.

Alaska↗

Monitoring landscape dynamics in central U.S. grasslands with harmonized Landsat-8 and Sentinel-2 time series data

Remotely monitoring changes in central U.S. grasslands is challenging because these landscapes tend to respond quickly to disturbances and changes in weather. Such dynamic responses influence nutrient cycling, greenhouse gas contributions, habitat availability for wildlife, and other ecosystem processes and services. Traditionally, coarse-resolution satellite data acquired at daily intervals have been used for monitoring. Recently, the harmonized Landsat-8 and Sentinel-2 (HLS) data increased the temporal frequency of the data. Here we investigated if the increased data frequency provided adequate observations to characterize highly dynamic grassland processes. We evaluated HLS data available for 2016 to (1) determine if data from Sentinel-2 contributed to an improvement in characterizing landscape processes over Landsat-8 data alone, and (2) quantify how observation frequency impacted results. Specifically, we investigated into estimating annual vegetation phenology, detecting burn scars from fire, and modeling within-season wetland hydroperiod and growth of aquatic vegetation. We observed increased sensitivity to the start of the growing season (SOST) with the HLS data. Our estimates of the grassland SOST compared well with ground estimates collected at a phenological camera site. We used the Continuous Change Detection and Classification (CCDC) algorithm to assess if the HLS data improved our detection of burn scars following grassland fires and found that detection was considerably influenced by the seasonal timing of the fires. The grassland burned in early spring recovered too quickly to be detected as change events by CCDC; instead, the spectral characteristics following these fires were incorporated as part of the ongoing time-series models. In contrast, the spectral effects from late-season fires were detected both by Landsat-8 data and HLS data. For wetland-rich areas, we used a modified version of the CCDC algorithm to track within-season dynamics of water and aquatic vegetation. The addition of Sentinel-2 data provided the potential to build full time series models to better distinguish different wetland types, suggesting that the temporal density of data was sufficient for within-season characterization of wetland dynamics. Although the different data frequency, in both the spatial and temporal dimensions, could cause inconsistent model estimation or sensitivity sometimes; overall, the temporal frequency of the HLS data improved our ability to track within-season grassland dynamics and improved results for areas prone to cloud contamination. The results suggest a greater frequency of observations, such as from harmonizing data across all comparable Landsat and Sentinel sensors, is still needed. For our study areas, at least a 3-day revisit interval during the early growing season (weeks 14–17) is required to provide a >50% probability of obtaining weekly clear observations.

Minnesota↗

Removing tidal-period variations from time-series data using low-pass digital filters

Several low-pass, digital filters are examined for their ability to remove tidal Period Variations from a time-series of water surface elevation for San Francisco Bay. The most efficient filter is the one which is applied to the Fourier coefficients of the transformed data, and the filtered data recovered through an inverse transform. The ability of the filters to remove the tidal components increased in the following order: 1) cosine-Lanczos filter, 2) cosine-Lanczos squared filter; 3) Godin filter; and 4) a transform fitter. The Godin fitter is not sufficiently sharp to prevent severe attenuation of 2–3 day variations in surface elevation resulting from weather events.

Journal of Physical Oceanography↗

Analysis of hydrologic and geochemical time-series data at James Cave, Virginia: Implications for epikarst influence on recharge in Appalachian karst aquifers

The epikarst, which consists of highly weathered rock in the upper vadose zone of exposed karst systems, plays a critical role in determining the hydrologic and geochemical characteristics of recharge to an underlying karst aquifer. This study utilized time series (2007–2014) of hydrologic and geochemical data of drip water collected within James Cave, Virginia, to examine the influence of epikarst on the quantity and quality of recharge in a mature, doline-dominated karst terrain. Results show a strong seasonality of both hydrology and geochemistry of recharge, which has implications for management of karst aquifers in temperate climatic zones. First, recharge (discharge from the epikarst to the underlying aquifer) reaches a maximum between late winter and early spring, with the onset of the recharge season ranging from as early as December to as late as March during the study period. The timing and duration of the recharge season were found to be a function of precipitation in excess of evapotranspiration on a seasonal time scale. Secondly, seasonally variable residence times for water in the epikarst influence rock-water interaction and, hence, the geochemical characteristics of recharge. Overall, results highlight the strong and complex influence that the epikarst has on karst recharge, which requires long-term and high-resolution data sets to accurately understand and quantify.

Virginia↗

Global cropland-extent product at 30-m resolution (GCEP30) derived from Landsat satellite time-series data for the year 2015 using multiple machine-learning algorithms on Google Earth Engine cloud

Executive Summary Global food and water security analysis and management require precise and accurate global cropland-extent maps. Existing maps have limitations, in that they are (1) mapped using coarse-resolution remote-sensing data, resulting in the lack of precise mapping location of croplands and their accuracies; (2) derived by collecting and collating national statistical data that are often subjective, leading to substantial uncertainties in cropland-area estimates, as well as their locations; and (3) extracted from one or more classes of a land use–land cover product in which cropland classes are not the focus of mapping, leading to their mixing with other classes and creating significant errors of omission and commission. These limitations can be overcome by producing high-resolution cropland-extent maps using satellite-sensor data, such as Landsat 30-m resolution or higher. The most fundamental cropland product is the high-resolution cropland-extent map because all higher level cropland products, such as crop-watering method (that is, whether crops are irrigated or rainfed), crop types, cropping intensities, cropland fallows, crop productivity, and crop-water productivity, are dependent on a precise and accurate cropland-extent product. Given these realities, the overarching goal of this study was to produce a Landsat satellite-derived global cropland-extent product at 30-m resolution. The work, which involved a paradigm shift in how global cropland-extent maps are produced, involved the following five key steps: (1) petabyte-scale computing that involved multiyear, 8- to 16-day, time-series Landsat 30-m resolution data for the global land surface; (2) composition of analysis-ready data (ARD) cubes; (3) creation of a large global-reference data hub for machine learning; (4) use of multiple machine-learning algorithms (MLAs) by writing software and computing in the cloud; and (5) Google Earth Engine (GEE) cloud computing. The five key steps involved nine distinct phases. First, the world was segmented into 74 agroecological zones (AEZs). Second, Landsat 8- to 16-day data were used to time-composite 10-band (blue, green, red, near-infrared, short-wave infrared band 1, short-wave infrared band 2, thermal infrared, enhanced vegetation index, normalized difference water index, and normalized difference vegetation index) Landsat 30-m resolution data cubes for every 2- to 4-month time period during 3- to 4-year periods (stated as nominal-year 2015 or, simply, 2015), along with two additional 30-m resolution bands (Shuttle Radar Topography Mission elevation, and slope) in each of the 74 AEZs. Third, more than 100,000 reference-training data samples were collected using ground data (some of which were collected using a mobile application), as well as submeter- to 5-m-resolution, very high-resolution imagery sourced from other reliable sources. Fourth, reference-training data were used to create a knowledge base for separating cropland from noncropland. Fifth, MLAs such as the pixel-based supervised random forest and support-vector machines were written on the GEE using Python and JavaScript. Sixth, object-based recursive hierarchical segmentation algorithm was used, in addition to MLAs, to overcome uncertainties. Seventh, MLAs used the knowledge base to classify and separate cropland from noncropland. Eighth, accuracy assessment was conducted by generating error matrices for each of the 74 AEZs using 19,171 independent validation-data samples. Ninth, cropland areas were computed for all countries of the world and compared with United Nation’s (UN’s) Food and Agricultural Organization (FAO) and other national statistics. The outcome was a Landsat-derived global cropland-extent product at 30-m resolution (GCEP30), which has an overall accuracy of 91.7 percent. For the cropland class, producer’s accuracy was 83.4 percent, and user’s accuracy was 78.3 percent. GCEP30 calculated (using direct pixel count) the global net-cropland area (GNCA) for the year 2015 as 1.873 billion hectares (~12.6 percent of the Earth’s terrestrial area). The continental cropland distribution as a percentage of GNCA was Asia, 33 percent; Europe, 25.5 percent; Africa, 16.7 percent; North America, 14.4 percent; South America, 8.1 percent; and Australia and Oceania, 2.4 percent. The worldwide cropland areas in GCEP30 for 2015 were higher by 236 to 299 million hectares (Mha) compared to national statistics reported elsewhere for the same year (for example, in Food and Agriculture Organization’s corporate statistical database [FAOSTAT] and in the monthly irrigated and rainfed crop areas [MIRCA] database). The global cropland area reported for 2015 increased by 344 Mha (22.5 percent), compared to the year 2000. During the same period (2000–2015), the world’s population increased by 20 percent. Whereas some of these areal increases are real increases in cropland areas, others are due to the types of data, methods, and approaches used. Using the highest known resolution (compared to previous coarse-resolution global products) enabled this study to capture fragmented croplands. Coarse-resolution data compute areas on the basis of subpixels, which, for a large proportion of certain land use–land cover classes, will show only a certain percentage of the total pixel area as actual area. Subpixel areas can lead to substantial uncertainties in area computation, as determining the exact fraction of cropland areas within a coarse-resolution pixel is resource intensive and subject to errors. Other innovations in GCEP30 include reference-data hubs, machine learning, and cloud computing. Cropland areas in 214 countries, territories, departments, and regions were calculated for the year 2015 using GCEP30, on the basis of UN’s global administrative unit layers (GAUL) boundaries. The 10 leading countries in terms of cropland area (as a percentage of the GNCA) were India (9.6 percent), United States (8.95 percent), China (8.82 percent), Russia (8.32 percent), Brazil (3.42 percent), Ukraine (2.32 percent), Canada (2.29 percent), Argentina (2.05 percent), Indonesia (2 percent), and Nigeria (1.91 percent). Together, these 10 countries occupy 50 percent of the global cropland, and they have 52 percent of the global population. Their combined cropland area increased by 2 percent between 2000 and 2015, compared to the substantial increase in population of 517 million (15.5 percent). Together, India, United States, China, and Russia encompass 36 percent of the total area. In the United States and Canada, from 2000 to 2015, cropland decreased by about 2 percent, whereas their populations increased by 14 and 13 percent, respectively. The additional food requirements in these 10 countries, which are caused by increased populations, as well as increasing nutritional demands, are met by production increases in existing cropland or through virtual food trade, or both. More than 18 countries, territories, departments, or regions had 60 percent or more of their geographic area as cropland: Republic of Moldova, San Marino, and Hungary had more than 80 percent of the country’s area as cropland; Denmark, Ukraine, Ireland, and Bangladesh, 70 to 80 percent; and Uruguay, Netherlands, United Kingdom, Spain, Lithuania, Poland, Gaza Strip, Czechia, Italy, India, and Azerbaijan, 60 to 70 percent. Europe and South Asia can be considered agricultural capitals of the world, on the basis of their percentages of geographic area as cropland. United States, China, and Russia, which all have high cropland areas, are ranked second, third, and fourth in the world; India is ranked first. However, the amount of cropland as a percentage of the country’s geographic area is relatively very low for United States (18.3 percent), China (17.7 percent), and Russia (9.5 percent), whereas it is 60.5 percent for India. Most African and South American countries, territories, departments, or regions have less than 15 percent of their geographic area as cropland. China and India together house 36 percent of the world’s population; however, between 2000 and 2015, the amount of China’s cropland area fell by 18.9 percent, owing to urban expansion and the abandonment of farmlands caused by demographic changes (that is, the movement of population from villages to cities). In contrast, China’s population grew by 10 percent. The amount of India’s cropland increased by 8.5 percent, whereas its population grew by 20 percent. This study showed that, out of the 10 leading cropland countries, Ukraine, Nigeria, Russia, and Indonesia showed an 18 to 31 percent increase in cropland areas, on the basis of GCEP30 by the year 2015, compared to 2000. Nigeria’s cropland area increased by 25 percent, and its population increased by 31 percent in the same period. In these countries, food security is maintained by cropland expansion, productivity increases, and virtual food trade. Nevertheless, this trend of increasing net-cropland area and productivity will likely become difficult to maintain, owing to diminishing arable lands and plateauing of 50 years of continual yield increases, requiring policymakers to explore novel and data-supported approaches to solving future food security issues. The GCEP30 product, which can be browsed at full resolution at www.croplands.org , has been released for public download and use through U.S. Geological Survey (USGS)–National Aeronautics and Space Administration (NASA) Land Processes Distributed Active Archive Center (see https://lpdaac.usgs.gov/news/release-of-gfsad-30-meter-cropland-extent-products/ ).

Professional Paper↗

Perspectives on monitoring gradual change across the continuity of Landsat sensors using time-series data

There are many types of changes occurring over the Earth's landscapes that can be detected and monitored using Landsat data. Here we focus on monitoring “within-state,” gradual changes in vegetation in contrast with traditional monitoring of “abrupt” land-cover conversions. Gradual changes result from a variety of processes, such as vegetation growth and succession, damage from insects and disease, responses to shifts in climate, and other factors. Despite the prevalence of gradual changes across the landscape, they are largely ignored by the remote sensing community. Gradual changes are best characterized and monitored using time-series analysis, and with the successful launch of Landsat 8 we now have appreciable data continuity that extends the Landsat legacy across the previous 43 years. In this study, we conducted three related analyses: (1) comparison of spectral values acquired by Landsats 7 and 8, separated by eight days, to ensure compatibility for time-series evaluation; (2) tracking of multitemporal signatures for different change processes across Landsat 5, 7, and 8 sensors using anniversary-date imagery; and (3) tracking the same type of processes using all available acquisitions. In this investigation, we found that data representing natural vegetation from Landsats 5, 7, and 8 were comparable and did not indicate a need for major modification prior to use for long-term monitoring. Analyses using anniversary-date imagery can be very effective for assessing long term patterns and trends occurring across the landscape, and are especially good for providing insights regarding trends related to long-term and continuous trends of growth or decline. We found that use of all available data provided a much more comprehensive level of understanding of the trends occurring, providing information about rate, duration, and intra- and inter-annual variability that could not be readily gleaned from the anniversary date analyses. We observed that using all available clear Landsat 5–8 observations with the new Continuous Change Detection and Classification (CCDC) algorithm was very effective for illuminating vegetation trends. There are a number of potential challenges for assessing gradual changes, including atmospheric impacts, algorithm development and visualization of the changes. One of the biggest challenges for studying gradual change will be the lack of appropriate data for validating results and products.

Remote Sensing of Environment↗

Modelling population change from time series data

Information on change in population size over time is among the most basic inputs for population management. Unfortunately, population changes are generally difficult to identify, and once identified difficult to explain. Sources of variald (patterns) in population data include: changes in environment that affect carrying capaciyy and produce trend, autocorrelative processes, irregular environmentally induced perturbations, and stochasticity arising from population processes. In addition. populations are almost never censused and many surveys (e.g., the North American Breeding Bird Survey) produce multiple, incomplete time series of population indices, providing further sampling complications. We suggest that each source of pattern should be used to address specific hypotheses regarding population change, but that failure to correctly model each source can lead to false conclusions about the dynamics of populations. We consider hypothesis tests based on each source of pattern, and the effects of autocorrelated observations and sampling error. We identify important constraints on analyses of time series that limit their use in identifying underlying relationships.

Book chapter↗

The use of continuous water-quality time-series data to compute nutrient loadings for selected Iowa streams, 2008–17

In support of nutrient reduction efforts, nitrate (as nitrate plus nitrite) and phosphorus loads and yields were computed for selected streams in Iowa based on continuously monitored sensor data for 2008–17 and 2014–17, respectively. Sample data were used to assess nitrate sensor bias and to create phosphorus-turbidity surrogate models. Where needed, nitrate loads were corrected for site-specific sensor bias, which was determined to be as high as 9.25 percent. Nitrate loads presented in this report using continuous (generally 15-minute interval) data were on average 4 percent less, but as much as 38 percent less, than annual loads computed from daily mean nitrate concentrations not corrected for sensor bias. Streamflow-based phosphorus models had poorer fit (adjusted coefficient of determination values less than 0.75) than turbidity-based models (adjusted coefficient of determination approximately 0.9). However, alternate models based on streamflow were used to obtain a more complete annual phosphorus load despite seasonal and fragmentary sensor data. Mean annual nitrate yields for 18 selected sites (96 site-years) ranged from 1.68 to 164 pounds per square mile per day (lb/mi 2 /d), compared to 19.4 lb/mi 2 /d average statewide yield needed to achieve the nitrate-reduction goal. Mean annual phosphorus yields for selected sites on the Maquoketa River, South Raccoon River, and West Nishnabotna River range from 1.57 to 7.19 lb/mi 2 /d, compared to 1.06 lb/mi 2 /d average statewide yield needed to achieve the phosphorus-reduction goal.

Iowa↗

The use of continuous water-quality time-series data to compute total phosphorus loadings for the Turkey River at Garber, Iowa, 2018–20

In support of nutrient reduction efforts, total phosphorus loads and yields were computed for the Turkey River at Garber, Iowa (U.S. Geological Survey station 05412500), for January 1, 2018, to April 30, 2020, based on continuously monitored turbidity sensor data. Sample data were used to create a total phosphorus turbidity-surrogate model. Streamflow-based total phosphorus models were used during periods of missing sensor data to obtain a more complete annual total phosphorus load. This report presents methods needed to accurately compute site-specific loads and track annual progress toward nutrient reduction goals within the State. Annual total phosphorus loads for the Turkey River at Garber, Iowa, were 1,740 and 1,490 U.S. short tons for 2018 and 2019, respectively, with annual yields ranging from 3.01 to 3.53 pounds per acre per year, compared to a mean statewide yield of 0.73 pound per acre per year needed to achieve the total phosphorus-reduction goal.

Iowa↗

Use of continuous water-quality time-series data to compute total phosphorus concentrations and loads for the Missouri River at St. Joseph and Hermann, Missouri, 2007–22

In support of Missouri’s Nutrient Loss Reduction Strategy, which was created to reduce the nutrient contamination of Missouri’s waterways from point and nonpoint sources, total phosphorus concentrations and loads were computed for the Missouri River at St. Joseph, Missouri, streamgage (U.S. Geological Survey station 06818000) and the Missouri River at Hermann, Mo., streamgage (U.S. Geological Survey station 06934500) for October 2007 to September 2022 using surrogate models and continuous turbidity sensor data. To obtain a more complete total phosphorus record for the study period, LOAD ESTimator (LOADEST) regression models using flow were used when turbidity sensor data were unavailable to estimate daily total phosphorus loads. This report presents the methods and results for the computed total phosphorus concentrations, loads, and yields for the two study sites on the Missouri River. With continued data collection and ongoing model evaluation and maintenance, the surrogate models may be useful into the future for computing total phosphorus concentrations and loads. Daily mean total phosphorus concentrations calculated using a surrogate model at the Missouri River at St. Joseph, Mo., streamgage during the 15-year study period (water years 2008 through 2022) ranged from 0.104 to 4.56 milligrams per liter (mg/L; median of 0.272 mg/L), and computed total phosphorus daily loads (with gaps in the daily record filled using the LOADEST regression model) ranged from 5.19 to 1,760 tons per day (tons/d; median of 36.5 tons/d). Annual loads ranged from 9,570 tons in water year 2022 to 50,500 tons in water year 2019. The total load for the study period was 437,000 tons. For the Missouri River at Hermann, Mo., streamgage during the same 15-year study period, daily mean total phosphorus concentrations, calculated using surrogate models applied to low and high turbidity values, ranged from 0.183 to 1.97 mg/L (median of 0.319 mg/L), and computed total phosphorus daily loads (with gaps in the daily record filled using the LOADEST regression model) ranged from 12.7 to 1,970 tons/d (median of 76.8 tons/d). Annual loads ranged from 22,600 tons in water year 2022 to 101,000 tons in water year 2019. The total load for the study period was 833,000 tons, which is nearly twice that at the Missouri River at St. Joseph, Mo., streamgage.

Missouri↗

VIMTS: Variational-based Imputation for Multi-modal Time Series

Multi-modal time series data in real applications often contain data of different dimensionalities, e.g., high-dimensional modality such as image data series, and low-dimensional univariate time series. Multi-modal time series data with missing high-dimensional modal values are ubiquitous in real-world classification and regression applications. To accurately predict the target labels, it is important to appropriately impute the high-dimensional modal missing values. However, most existing imputation methods focus on multivariate time series, fail to simultaneously consider temporal dependencies within each series and the correlations across the series, and also lack a probabilistic interpretation. In this paper, we propose a novel method, which uses a new structured variational approximation technique for the imputation of missing values in multi-modal time series. Instead of directly imputing high-dimensional modal missing values, we use the variational approximation technique to impute intermediate lower-dimensional feature representations of high-dimensional modal missing values from simple modalities related to high-dimensional modality and then feed them into a dynamical model. The dynamical model captures the temporal dependencies of the feature representations and finally predicts the target labels. In order to address the optimization difficulties caused by the lack of ground truth values of lower-dimensional feature representations, we also propose a two-stage isolated optimization strategy for better convergence. We evaluate our method on a real-world stream monitoring dataset. Our extensive experiments demonstrate that the proposed method outperforms several state-of-the-art methods in both data imputation and prediction performance.

Conference Paper↗