Search

Demand for food and water are increasing while the extent of arable land and accessible fresh water are decreasing. This poses global challenges as economies continue to develop and the population grows. With agriculture as the leading consumer of water, better understanding how water is used to produce food may help support the increase of Crop Water Productivity (CWP; kg/m3), the ratio of crop output per unit of water input (or crop per drop). Previous large-scale CWP studies have been useful for broad water use modeling at coarser resolutions. However, obtaining more precise CWP, especially for specific crop types in a particular area and growing season as outlined here are important for informing farm-scale water management decision making. Therefore, this study focused on California's Central Valley utilizing high-spatial resolution satellite imagery of 30 m (0.09 hectares per pixel) to generate more precise CWP for commonly grown and water-intensive irrigated crops. First, two products were modeled and mapped. 1. Landsat based Actual Evapotranspiration (ETa; mm/d) to determine Crop Water Use (CWU; m3/m2), and 2. Crop Productivity (CP; kg/m2) to estimate crop yield per growing season. Then, CWP was calculated by dividing CP by CWU and mapped. The amount of water that can be saved by increasing CWP of each crop was further calculated. For example, in the 434 million m2 study area, a 10% increase in CWP across the 9 crops analyzed had a potential water savings of 31.5 million m3 of water. An increase in CWP is widely considered the best approach for saving maximum quantities of water. This paper proposed, developed, and implemented a workflow of combined methods utilizing cloud computing based remote sensing data. The environmental implications of this work in assessing water savings for food and water security in the 21st century are expected to be significant. [ABSTRACT FROM AUTHOR]

Simple Summary: Excessive use of groundwater resources in the Central Valley in California has led to major land sinking over the years. In this study, we rely on satellite imagery to monitor and assess the extent of this sinking. Specifically, we use images from a satellite that emits microwaves which allows us to directly obtain the deformation at a specific time. Then, we apply a machine learning algorithm to the resulting data in an attempt to effectively predict short-term future deformation. We find that the algorithm applied has a low error when compared to the actual data. This shows that machine learning algorithms could be incorporated into models that assess potential hazards associated with land sinking. Land subsidence caused by excessive groundwater pumping in Central Valley, California, is a major issue that has several negative impacts such as reduced aquifer storage and damaged infrastructures which, in turn, produce an economic loss due to the high reliance on crop production. This is why it is of utmost importance to routinely monitor and assess the surface deformation occurring. Two main goals that this paper attempts to accomplish are deformation characterization and deformation prediction. The first goal is realized through the use of Principal Component Analysis (PCA) applied to a series of Interferomtric Synthetic Aperture Radar (InSAR) images that produces eigenimages displaying the key characteristics of the subsidence. Water storage changes are also directly analyzed by the use of data from the Gravity Recovery and Climate Experiment (GRACE) twin satellites and the Global Land Data Assimilation System (GLDAS). The second goal is accomplished by building a Long Short-Term Memory (LSTM) model to predict short-term deformation after developing an InSAR time series using LiCSBAS, an open-source InSAR time series package. The model is applied to the city of Madera and produces better results than a baseline averaging model and a one dimensional convolutional neural network (CNN) based on a mean squared error metric showing the effectiveness of machine learning in deformation prediction as well as the potential for incorporation in hazard mitigation models. The model results can directly aid policy makers in determining the appropriate rate of groundwater withdrawal while maintaining the safety and well-being of the population as well as the aquifers' integrity. [ABSTRACT FROM AUTHOR]

Parameter identification for wildfire forecasting models often relies on case-by-case tuning or posterior diagnosis/analysis, which can be computationally expensive due to the complexity of the forward prediction model. In this paper, we introduce an efficient parameter flexible fire prediction algorithm based on machine learning and reduced order modelling techniques. Using a training dataset generated by physics-based fire simulations, the method forecasts burned area at different time steps with a low computational cost. We then address the bottleneck of efficient parameter estimation by developing a novel inverse approach relying on data assimilation techniques (latent assimilation) in the reduced order space. The forward and the inverse modellings are tested on two recent large wildfire events in California. Satellite observations are used to validate the forward prediction approach and identify the model parameters. By combining these forward and inverse approaches, the system manages to integrate real-time observations for parameter adjustment, leading to more accurate future predictions. [ABSTRACT FROM AUTHOR]

Evapotranspiration (ET) represents crop water use and is a key indicator of crop health. Accurate estimation of ET is critical for agricultural irrigation and water resource management. ET retrieval using energy balance methods with remotely sensed thermal infrared data as the key input has been widely applied for irrigation scheduling, yield prediction, drought monitoring and so on. However, limitations on the spatial and temporal resolution of available thermal satellite data combined with the effects of cloud contamination constrain the amount of detail that a single satellite can provide. Fusing satellite data from different satellites with varying spatial and temporal resolutions can provide a more continuous estimation of daily ET at field scale. In this study, we applied an ET fusion modeling system, which uses a surface energy balance model to retrieve ET using both Landsat and Moderate Resolution Imaging Spectroradiometer (MODIS) data and then fuses the Landsat and MODIS ET retrieval timeseries using the Spatial-Temporal Adaptive Reflectance Fusion Model (STARFM). In this paper, we compared different STARFM ET fusion implementation strategies over various crop lands in the central California. In particular, the use of single versus two Landsat-MODIS pair images to constrain the fusion is explored in cases of rapidly changing crop conditions, as in frequently harvested alfalfa fields, as well as an improved dual-pair method. The daily 30 m ET retrievals are evaluated with flux tower observations and analyzed based on land cover type. This study demonstrates improvement using the new dual-pair STARFM method compared with the standard one-pair STARFM method in estimating daily field scale ET for all the major crop types in the study area. [ABSTRACT FROM AUTHOR]

As a result of the advocacy of Indigenous communities and increasing evidence of the ecological importance of fire, California has invested in the restoration of intentional burning (the practice of deliberately lighting low-severity fires) in an effort to reduce the occurrence and severity of wildfires. Recognizing the growing need to monitor the impacts of these smaller, low-severity fires, we leveraged Sentinel-2 imagery to reveal important inter- and intra-annual variation in grasslands before and after fires. Specifically, we explored three methodological approaches: (1) the complete time series of the normalized burn ratio (NBR), (2) annual summary metrics (mean, fifth percentile, and amplitude of NBR), and (3) maps depicting spatial patterns in these annual NBR metrics before and after fire. We also used a classification of pre-fire vegetation to stratify these analyses by three dominant vegetation cover types (grasses, shrubs, and trees). We applied these methods to a unique study area in which three adjacent grasslands had diverging fire histories and showed how grassland recovery from a low-severity intentional burn and a high-severity wildfire differed both from each other and from a reference site with no recent fire. On the low-severity intentional burn site, our results showed that the annual NBR metrics recovered to pre-fire values within one year, and that regular intentional burning on the site was promoting greater annual growth of both grass and shrub species, even in the third growing season following a burn. In the case of the high-severity wildfire, our metrics indicated that this grassland had not returned to its pre-fire phenological signals in at least three years after the fire, indicating that it may be undergoing a longer recovery or an ecological shift. These proposed methods address a growing need to study the effects of small, intentional burns in low-biomass ecosystems such as grasslands, which are an essential part of mitigating wildfires. [ABSTRACT FROM AUTHOR]

In 2020, wildfires scarred over 4,000,000 hectares in the western United States, devastating urban populations and ecosystems alike. The significant impact that wildfires have on plants, animals, and human environments makes wildfire adaptation, management, and mitigation strategies a critical task. This study uses satellite imagery from Landsat to calculate burn severity and map the fire progression for the Dolan Fire of central Coastal California which occurred in August 2020. Several environmental factors, such as temperature, humidity, fuel type, topography, surface conditions, and wind velocity, are known to affect wildfire spread and burn severity. The aim of this study is the investigation of the relationship between these environmental factors, estimates of burn severity, and fire spread patterns. Burn severity is calculated and classified using the Difference in Normalized Burn Ratio (dNBR) before being displayed as a time series of maps. The Dolan Fire had a moderate severity burn with an average dNBR of 0.292. The ignition site location, when paired with the patterns of fire spread, is consistent with wind speed and direction data, suggesting fire movement to the southeast of the fire ignition site. Patterns of increased burn severity are compared with both topography (slope and aspect) and fuel type. Locations that were found to be more susceptible to high burn severity featured Long Needle Timber Litter and Mature Timber fuels, intermediate slope angles between 15 and 35°, and north- and east-facing slopes. This study has implications for the future predictive modeling of wildfires that may serve to develop wildfire mitigation strategies, manage climate change impacts, and protect human lives. [ABSTRACT FROM AUTHOR]

There are an estimated 800,000 small reservoirs globally with a range of uses. Given the collective importance of these reservoirs to water resource management and wider society, it is essential that we can monitor and understand the hydrological dynamics of ungauged reservoirs, particularly in a changing climate. However, unlike large reservoirs, continuous and systematic hydrological observations of small reservoirs are often unavailable. In response, this study has developed a retrieval framework for water levels of small reservoirs using a deep learning algorithm and remotely sensed satellite data. Demonstrated at four reservoirs in California, satellite imagery from both Sentinel-1 and Sentinel-2 along with corresponding water level field measurements was collected. Post-processed images were fed into a water level inversion convolutional neural network model for water level inversion, while different combinations of these satellite images, sampling approaches for training/testing data, and attention modules were used to train the model and evaluated for accuracy. The results show that random sampling of training data coupled with Sentinel-2 satellite imagery was generally the most accurate initially. Performance is improved by incorporating a channel attention mechanism, with the average R2 increasing by 8.6% and the average RMSE and MAE decreasing by 15.5% and 36.4%, respectively. The proposed framework was further validated on three additional reservoirs in different regions. In conclusion, the retrieval framework proposed in this study provides a stable and accurate methodology for water level estimation of small reservoirs and can be a powerful tool for small reservoir monitoring over large spatial scales. [ABSTRACT FROM AUTHOR]

The photovoltaic (PV) industry boom has accelerated the need for accurately understanding the spatial distribution of PV energy systems. The synergy of remote sensing and artificial intelligence presents significant prospects for PV energy monitoring. Currently, numerous studies have focused on extracting rooftop PV systems from airborne or satellite imagery, but their small-scale and size-varying characteristics make the segmentation results suffer from PV internal incompleteness and small PV omission. To address these issues, this study proposed a size-aware deep learning network called Rooftop PV Segmenter (RPS) for segmenting small-scale rooftop PV systems from high-resolution imagery. In detail, the RPS network introduced a Semantic Refinement Module (SRM) to sense size variations of PV panels and reconstruct high-resolution deep semantic features. Moreover, a Feature Aggregation Module (FAM) enhanced the representation of robust features by continuously aggregating deeper features into shallower ones. In the output stage, a Deep Supervised Fusion Module (DSFM) was employed to constrain and fuse the outputs at different scales to achieve more refined segmentation. The proposed RPS network was tested and shown to outperform other models in producing segmentation results closer to the ground truth, with the F1 score and IoU reaching 0.9186 and 0.8495 on the publicly available California Distributed Solar PV Array Dataset (C-DSPV Dataset), and 0.9608 and 0.9246 on the self-annotated Heilbronn Rooftop PV System Dataset (H-RPVS Dataset). This study has provided an effective solution for obtaining a refined small-scale energy distribution database. [ABSTRACT FROM AUTHOR]

The escalating frequency and severity of global wildfires necessitate an in-depth understanding and monitoring of wildfire smoke impacts, specifically its contribution to fine particulate matter (PM 2.5 ). We propose a data-fusion method to study wildfire contribution to PM 2.5 using satellite-derived smoke plume indicators and PM 2.5 monitoring data. Our study incorporates two types of monitoring data, the high-quality but sparse Air Quality System (AQS) stations and the abundant but less accurate PurpleAir (PA) sensors that are gaining popularity among citizen scientists. We propose a multi-resolution spatiotemporal model specified in the spectral domain to calibrate the PA sensors against accurate AQS measurements, and leverage the two networks to estimate wildfire contribution to PM 2.5 in California in 2020 and 2021. A Bayesian approach is taken to incorporate all uncertainties and our prior intuition that the dependence between networks, as well as the accuracy of PA network, vary by frequency. We find that 1% to 3% increase in PM 2.5 concentration due to wildfire smoke, and that leveraging PA sensors improves accuracy. [ABSTRACT FROM AUTHOR]

Earthquakes are known to generate disturbances in the ionosphere. Such disturbances, referred to as co-seismic ionospheric disturbances, or ionoquakes, were previously reported for large earthquakes with magnitudes Mw ≥ 6.6. This paper reports ionoquakes associated with the Ridgecrest earthquakes of magnitude (Mw = 6.4) , that occurred on 4 July 2019 in California, USA. The ionoquakes manifested in total electron content (TEC) in the form of traveling ionospheric disturbances (TIDs) within 1 h from the mainshock onset. These seismic-origin TIDs have unique wave characteristics that distinguish them from TIDs of non-seismic origin arising from a moderate geomagnetic activity on the same day. Moreover, in the space-time domain of the detection of seismic-origin TIDs, TIDs are absent on the day before and day after the earthquake day. Their spectral characteristics relate them to the Earth's normal modes and atmospheric resonance modes. We found the ground velocity associated with the mainshock, rather than the ground displacement, satisfies the threshold criteria for detectable ionoquakes in TEC measurements. Numerical simulation suggested that the coupled seismo–atmosphere–ionosphere (SAI) dynamics energized by the atmospheric waves are responsible for the generation of ionoquakes. This study's findings demonstrate the potential of using TEC measurement to detect the ionospheric counterparts of moderate earthquakes. [ABSTRACT FROM AUTHOR]

In the American West, wildfires and earthquakes are increasingly threatening the archaeological, historical, and tribal resources that define the collective identity and connection with the past for millions of Americans. The loss of said resources diminishes societal understanding of the role cultural heritage plays in shaping our present and future. This paper examines the viability of employing stationary and SLAM-based terrestrial laser scanning, close-range photogrammetry, automated surface change detection, GIS, and WebGL visualization techniques to enhance the preservation of cultural resources in California. Our datafication approach combines multi-temporal remote sensing monitoring of historic features with legacy data and collaborative visualization to document and evaluate how environmental threats affect built heritage. We tested our methodology in response to recent environmental threats from wildfire and earthquakes at Bodie, an iconic Gold Rush-era boom town located on the California and Nevada border. Our multi-scale results show that the proposed approach effectively integrates highly accurate 3D snapshots of Bodie's historic buildings before/after disturbance, or post-restoration, with surface change detection and online collaborative visualization of 3D geospatial data to monitor and preserve important cultural resources at the site. This study concludes that the proposed workflow enhances the monitoring of at-risk California's cultural heritage and makes a call to action to employ remote sensing as a pathway to advanced planning. [ABSTRACT FROM AUTHOR]

We investigate the utility of towed underwater camera systems with tightly coupled Global Navigation Satellite System (GNSS) positions to provide reef-scale bathymetric models with millimeter to centimeter resolutions and accuracies with Structure-from-Motion (SfM) photogrammetry. Successful development of these techniques would allow for detailed assessments of benthic conditions, including the accretion and erosion of reefs and adjacent sediment deposits, without the need for ground control points. We use a multi-camera system towed by a small vessel to map over 70,000 m2 of complex shallow (2–8 m water depth) bedrock reef, boulder fields, and fine (sand and gravel) sediments of Lake Tahoe, California. We find that multiple synchronized cameras increase overall mapping coverage and allow for wider survey line spacing. The accuracy of the techniques was sub-millimeter for local length measurements less than a meter, and the bathymetric reproducibility was found to scale with the accuracy of GNSS (3–5 cm), although this could be improved to sub-centimeter with the inclusion of one or more co-registered, but unsurveyed, control points. For future applications, we provide guidance on conducting field operations, correcting underwater image color, and optimizing the SfM workflows. We conclude that a GNSS-coupled underwater camera array is a promising technique to map shallow reefs at high accuracy and resolution without ground control. [ABSTRACT FROM AUTHOR]

Temporary roads are often placed in mountainous regions for logging purposes but then never decommissioned and removed. These abandoned forest roads often have unwanted environmental consequences. They can lead to altered hydrological regimes, excess erosion, and mass wasting events. These events can affect sediment budgets in streams, with negative consequences for anadromous fish populations. Maps of these roads are frequently non-existent; therefore, methods need to be created to identify and locate these roads for decommissioning. Abandoned logging roads in the Point Reyes National Seashore in California, an area partially under heavy forest canopy, were mapped using object-based image processing in concert with machine learning. High-resolution Q1 LiDAR point clouds from 2019 were used to create a bare earth model of the region, from which a slope model was derived. This slope model was then subjected to segmentation algorithms to identify and isolate regions of differing slopes. Regions of differing slopes were then used in a convolutional neural network (CNN), and a maximum likelihood classifier was used to delineate the historic road network. The accuracy assessment was conducted using historic aerial photos of the state of the region post-logging, along with ground surveys to verify the presence of logging roads in areas of question. This method was successfully able to identify road networks with a precision of 0.991 and an accuracy of 0.992. It was also found that the CNN was able to identify areas of highest disturbance to the slope gradient. This methodology is a valuable tool for decision makers who need to identify areas of high disturbance in order to mitigate adverse effects. [ABSTRACT FROM AUTHOR]

Shallow estuarine habitats are globally undergoing rapid changes due to climate change and anthropogenic influences, resulting in spatiotemporal shifts in distribution and habitat extent. Yet, scientists and managers do not always have rapidly available data to track habitat changes in real-time. In this study, we apply a novel and a state-of-the-art image segmentation machine learning technique (DeepLab) to two years of high-resolution drone-based imagery of a marine flowering plant species (eelgrass, a temperate seagrass). We apply the model to eelgrass (Zostera marina) meadows in the Morro Bay estuary, California, an estuary that has undergone large eelgrass declines and the subsequent recovery of seagrass meadows in the last decade. The model accurately classified eelgrass across a range of conditions and sizes from meadow-scale to small-scale patches that are less than a meter in size. The model recall, precision, and F1 scores were 0.954, 0.723, and 0.809, respectively, when using human-annotated training data and random assessment points. All our accuracy values were comparable to or demonstrated greater accuracy than other models for similar seagrass systems. This study demonstrates the potential for advanced image segmentation machine learning methods to accurately support the active monitoring and analysis of seagrass dynamics from drone-based images, a framework likely applicable to similar marine ecosystems globally, and one that can provide quantitative and accurate data for long-term management strategies that seek to protect these vital ecosystems. [ABSTRACT FROM AUTHOR]

Sustainable groundwater management requires high-quality and low-latency estimates of changes in groundwater storage (∆Sgw). However, estimates of ∆Sgw produced using traditional methods, including groundwater models and well-based measurements, typically lag years behind the present because collecting the required on-the-ground data is a time consuming, expensive, and labor-intensive task. Satellite remote sensing measurements provide potential pathways to overcome these limitations by quantifying ∆Sgw through closing the water balance. However, the range of spatial scales over which ∆Sgw can be accurately estimated using remote sensing products remains unclear. To bridge this knowledge gap, this study quantified ∆Sgw for the period of 2002 through to 2021 using the water balance method and multiple remote sensing products in two subbasins (~2700 km2–3500 km2) within California's Central Valley: (1) the Kaweah–Tule Subbasin, a region where the pumping of groundwater to support agriculture has resulted in decades of decline in head levels, resulting in land subsidence, damage to infrastructure, and contamination of drinking water and (2) the Butte Subbasin, which receives considerably more rainfall and surface water and has not experienced precipitous drops in groundwater. The remote sensing datasets which we utilized included multiple sources for key hydrologic components in the study area: precipitation, evapotranspiration, and soil moisture. To assess the fidelity of the remote sensing-based model, we compared estimates of ∆Sgw to alternative estimates of ∆Sgw derived from independent sources of data: groundwater wells as well as a widely used groundwater flow model. The results showed strong agreement in the Kaweah–Tule Subbasin in long-term ∆Sgw trends and shorter-term trends during droughts, and modest agreement in the Butte Subbasin with remote sensing datasets suggesting more seasonal variability than validation datasets. Importantly, our analysis shows that the timely availability of remote sensing data can potentially enable ∆Sgw estimates at sub-annual latencies, which is timelier than estimates derived through alternate methods, and thus can support adaptive management and decision making. The models developed herein can aid in assessing aquifer dynamics, and can guide the development of sustainable groundwater management practices at spatial scales relevant for management and decision making. [ABSTRACT FROM AUTHOR]

Evapotranspiration (ET) is a crucial part of commercial grapevine production in California, and the partitioning of this quantity allows the separate assessment of soil and vine water and energy fluxes. This partitioning has an important role in agriculture since it is related to grapevine stress, yield quality, irrigation efficiency, and growth. Satellite remote sensing-based methods provide an opportunity for ET partitioning at a subfield scale. However, medium-resolution satellite imagery from platforms such as Landsat is often insufficient for precision agricultural management at the plant scale. Small, unmanned aerial systems (sUAS) such as the AggieAir platform from Utah State University enable ET estimation and its partitioning over vineyards via the two-source energy balance (TSEB) model. This study explores the assessment of ET and ET partitioning (i.e., soil water evaporation and plant transpiration), considering three different resistance models using ground-based information and aerial high-resolution imagery from the Grape Remote sensing Atmospheric Profile and Evapotranspiration eXperiment (GRAPEX). We developed a new method for temperature partitioning that incorporated a quantile technique separation (QTS) and high-resolution sUAS information. This new method, coupled with the TSEB model (called TSEB-2TQ), improved sensible heat flux (H) estimation, regarding the bias, with around 61% and 35% compared with the H from the TSEB-PT and TSEB-2T, respectively. Comparisons among ET partitioning estimates from three different methods (Modified Relaxed Eddy Accumulation—MREA; Flux Variance Similarity—FVS; and Conditional Eddy Covariance—CEC) based on EC flux tower data show that the transpiration estimates obtained from the FVS method are statistically different from the estimates from the MREA and the CEC methods, but the transpiration from the MREA and CEC methods are statistically the same. By using the transpiration from the CEC method to compare with the transpiration modeled by different TSEB models, the TSEB-2TQ shows better agreement with the transpiration obtained via the CEC method. Additionally, the transpiration estimation from TSEB-2TQ coupled with different resistance models resulted in insignificant differences. This comparison is one of the first for evaluating ET partitioning estimation from sUAS imagery based on eddy covariance-based partitioning methods. [ABSTRACT FROM AUTHOR]

Accurate characterization of evapotranspiration (ET) is imperative in water-limited cropping systems such as California vineyards and almond orchards. Satellite-based ET modeling techniques, including the atmosphere–land exchange inverse model (ALEXI) and associated flux disaggregation technique (DisALEXI), have proven reliable in determining field scale ET. However, validation efforts typically focus on ET and omit an evaluation of partitioned evaporation (E) and transpiration (T). ALEXI/DisALEXI is based on the two-source energy balance (TSEB) model, making it uniquely qualified to derive E and T individually. The current study evaluated E and T estimates derived using two formulations of DisALEXI; one based on Priestley-Taylor (DisALEXI-PT) and the other on Penman-Monteith (DisALEXI-PM). The modeled values were validated against partitioned fluxes derived from the conditional eddy covariance (CEC) approach using EC flux towers in three wine grape vineyards and three almond orchards for the year 2021. Modeled estimates were derived using Landsat 8 Collection 2 thermal infrared and surface reflectance imagery as well as Harmonized Landsat and Sentinel-2 surface reflectance datasets as input into DisALEXI. The results indicated that the modeled total ET fluxes were similar between the two methods, but the partitioned values diverged, with DisALEXI-PT overestimating E and slightly underestimating T when compared to CEC estimates. Conversely, DisALEXI-PM agreed better with CEC-derived E and overestimated T estimates under non-advective conditions. Compared to one another, DisALEXI-PM estimated canopy temperatures ~5 °C cooler and soil temperatures ~5 °C warmer than DisALEXI-PT, causing differences in E and T of −2.6 mm day−1 and +2.6 mm day−1, respectively. The evaluation of the iterative process required for DisALEXI indicates DisALEXI-PM ET values converge on ALEXI ET with proportionate adjustments to E and T, while DisALEXI-PT convergence is driven by adjustments to E. The analysis presented here can potentially drive improvements in the modeling framework to provide specific soil and canopy consumptive water use information in unique canopy structures, allowing for improved irrigation and water use efficiencies in these water-limited systems. [ABSTRACT FROM AUTHOR]

The objectives of this study were to evaluate spatial sampling and statistical aspects of landscape-level wildfire rate of spread (ROS) estimates derived from airborne thermal infrared imagery (ATIR). Wildfire progression maps and ROS estimates were derived from repetitive ATIR image sequences collected during the 2017 Thomas and Detwiler wildfire events in California. Three separate landscape sampling unit (LSU) sizes were used to extract remotely sensed environmental covariates known to influence fire behavior. Statistical relationships between fire spread rates and landscape covariates were analyzed using (1) bivariate regression, (2) multiple stepwise regression, (3) geographically weighted regression (GWR), (4) eigenvector spatial filtering (ESF) regression, (5) regression trees (RT), and (6) and random forest (RF) regression. GWR and ESF regressions reveal that relationships between covariates and ROS estimates are substantially non-stationary and suggest that the global association of fire spread controls are locally differentiated on landscape scales. Directional slope is by far the most strongly associated covariate of ROS for the imaging sequences analyzed and the size of LSUs has little influence on any of the covariate relationships. [ABSTRACT FROM AUTHOR]

The fusion of global navigation satellite system (GNSS) and interferometric synthetic aperture radar (InSAR) deformation data can leverage the advantages of GNSS high temporal resolution and InSAR high spatial resolution, and obtain more abundant deformation data for constraints on geophysical structural and mechanical parameters. Existing studies seldom consider the spatial heterogeneity of largescale deformation data, which easily leads to obvious spatial aggregation of errors in the results of fusion. Here, we propose a novel multiresolution segmentation fusion (MRSF) method that uses a multiresolution segmentation algorithm to automatically classify the spatial heterogeneity of InSAR deformation data with similar deformation characteristics. We applied the MRSF method to the fusion of GNSS and InSAR deformation data covering the central valley aquifer system (CVAS) in southern California to verify its precision and robustness. Results show that the MRSF method can accurately reflect spatiotemporal evolution characteristics of displacement data and reliably estimate deformation for the times and locations of missing data. We then tested this method for geophysical parameter estimation by constructing three different sets of data, including dense GNSS sites, sparse GNSS sites, and sparse GNSS sites fused with InSAR data using MRSF, to invert the slip distribution of the Cascadia subduction zone. Results show that the inverted slip of the fused InSAR and GNSS data is comparable to that of the dense GNSS sites. Therefore, the MRSF method can obtain deformation results with high precision and high spatiotemporal resolution and effectively compensate for the lack of data caused by sparse GNSS sites during the geophysical inversion process. [ABSTRACT FROM AUTHOR]

Remote sensing techniques offer significant potential for generating accurate thick oil slick maps critical for marine oil spill response. However, field validation and methodology assessment challenges remain. Here, we report on an approach to leveraging oil emissions from the Coal Oil Point (COP) natural marine hydrocarbon seepage offshore of southern California, where prolific oil seepage produces thick oil slicks stretching many kilometers. Specifically, we demonstrate and validate a remote sensing approach as part of the Seep Assessment Study (SAS). Thick oil is sufficient for effective mitigation strategies and is set at 0.15 mm. The brightness temperature of thick oil, TBO, is warmer than oil-free seawater, TBW, allowing segregation of oil from seawater. High spatial-resolution airborne thermal and visible slick imagery were acquired as part of the SAS; including along-slick "streamer" surveys and cross-slick calibration surveys. Several cross-slick survey-imaged short oil slick segments that were collected by a customized harbor oil skimmer; termed "collects". The brightness temperature contrast, ΔTB (TBO − TBW), for oil pixels (based on a semi-supervised classification of oil pixels) and oil thickness, h, from collected oil for each collect provided the empirical calibration of ΔTB(h). The TB probability distributions provided TBO and TBW, whereas a spatial model of TBW provided ΔTB for the streamer analysis. Complicating TBW was the fact that streamers were located at current shears where two water masses intersect, leading to a TB discontinuity at the slick. This current shear arose from a persistent eddy down current of the COP that provides critical steering of oil slicks from the Coal Oil Point. The total floating thick oil in a streamer observed on 23 May and a streamer observed on 25 May 2016 was estimated at 311 (2.3 bbl) and 2671 kg (20 bbl) with mean linear floating oil 0.14 and 2.4 kg m−1 with uncertainties by Monte Carlo simulations of 25% and 7%, respectively. Based on typical currents, the average of these two streamers corresponds to 265 g s−1 (~200 bbl day−1) in a range of 60–340 bbl day−1, with significant short-term temporal variability that suggests slug flow for the seep oil emissions. Given that there are typically four or five streamers, these data are consistent with field emissions that are higher than the literature estimates. [ABSTRACT FROM AUTHOR]

Reflectance-based vegetation indices can be valuable for assessing crop nitrogen (N) status and predicting grain yield. While proximal sensors have been widely studied in agriculture, there is increasing interest in utilizing aerial sensors. Given that few studies have compared aerial and proximal sensors, the objective of this study was to quantitatively compare the sensitivity of aerially sensed Normalized Difference Vegetation Index (NDVI) and Normalized Difference Red-Edge Index (NDRE) and proximally sensed NDVI for assessing total N uptake at panicle initiation (PI-NUP) and predicting grain yield in rice. Nitrogen response trials were established over a 3-year period (10 site-years) at various locations throughout the Sacramento Valley rice growing region of California. At PI, a multispectral unmanned aircraft system (UAS) was used to measure NDVIUAS and NDREUAS (average ground sampling distance: 3.7 cm pixel−1), and a proximal GreenSeeker (GS) sensor was used to record NDVIGS. To enable direct comparisons across the different indices on an equivalent numeric scale, each index was normalized by calculating the Sufficiency-Index (SI) relative to a non-N-limiting plot. Kernel density distributions indicated that NDVIUAS had a narrower range of values that were poorly differentiated compared to NDVIGS and NDREUAS. The critical PI-NUP where yields did not increase with higher PI-NUP averaged 109 kg N ha−1 (±4 kg N ha−1). The relationship between SI and PI-NUP for the NDVIUAS saturated lower than this critical PI-NUP (96 kg N ha−1), whereas NDVIGS and NDREUAS saturated at 111 and 130 kg N ha−1, respectively. This indicates that NDVIUAS was less suitable for making N management decisions at this crop stage than NDVIGS and NDREUAS. Linear mixed effects models were developed to evaluate how well each SI measured at PI was able to predict grain yield. The NDVIUAS was least sensitive to variation in yields as reflected by having the highest slope (2.4 Mg ha−1 per 0.1 SI). In contrast, the slopes for NDVIGS and NDREUAS were 0.9 and 1.1 Mg ha−1 per 0.1 SI, respectively, indicating greater sensitivity to yields. Altogether, these results indicate that the ability of vegetation indices to inform crop management decisions depends on the index and the measurement platform used. Both NDVIGS and NDREUAS produced measurements sensitive enough to inform N fertilizer management in this system, whereas NDVIUAS was more limited. [ABSTRACT FROM AUTHOR]

Off-channel areas are one of the most impacted aquatic habitats by humans globally, as extensive agricultural and urban development has limited them to roughly 10% of historical extent. This is also true for California's Sacramento River Valley, where historically frequent widespread inundation has been reduced to a few off-channel water bodies along the mid-Sacramento River. This remaining shallow-water habitat provides crucial ecological benefits to multiple avian and fish species, but especially to floodplain-adapted species such as Chinook salmon (Oncorhynchus tshawytscha). Characterizing spatiotemporal off-channel dynamics, including inundation extent and residence time, is fundamental to better understanding the intrinsic value of such habitats and their potential to support recovery actions. Remote sensing techniques have been increasingly used to map surface water at regional and local scales, with improved resolutions. As such, this study maps off-channel inundation areas and describes their temporal dynamics by analyzing pixel-based time- series of multiple water indices, modified Normalized Difference Water Index (mNDWI) and the Automated Water Extraction Index (AWEI), generated from LandSat-8 and Sentinel-2 data between 2013–2021. Quantified off-channel area was similar with each water index and method used, but improved performance was associated with Sentinel-2 products and AWEI index to identify wetted areas under lower mainstem discharges. Results indicate an uneven distribution of off-channel habitat in the study area, with limited inundated areas in upstream reaches (<16% of total off-channel area for greater flows). In addition, much less habitat exists for flows under 400 m3/s, an important migration cue for endangered winter-run Chinook salmon, limiting juvenile access to areas with enhanced rearing conditions. Off-channel habitat residence times averaged between 7 and 16 days, primarily defined by the rate of receding flows, with rapid flow recession providing marginal off-channel habitat. This study shows reasonable performance of moderate resolution LandSat-8 and Sentinel-2 remote sensing imagery to characterize shallow-water inundated habitat in higher-order rivers, and as a method to inform restoration and native fish recovery efforts. [ABSTRACT FROM AUTHOR]

This study proposes a method to detect ocean fronts from in situ temperature and density glider measurements. This method is applied to data collected along the CalCOFI Line 90, south of the California Current System (CCS), over the 2006–2013 period. It is based on image-processing techniques commonly applied to sea surface temperature (SST) satellite data. Front detection results using glider data are consistent with those obtained in other studies carried out in the CCS. SST images of the Multi-scale Ultra-high Resolution (MUR) dataset were also used to compare the probability of occurrence or front frequency (FF) obtained with the two datasets. Glider and MUR temperatures are highly correlated. Along Line 90, frontal frequency exhibited the same maxima near the transition zone (~130 km offshore) as derived from MUR and glider datasets. However, marked differences were found in the bimonthly FF probability with high (low) front frequency in spring-summer for glider (MUR) data. Methodological differences explaining these contrasting results are investigated. Thermohaline-compensated fronts are more abundant towards the oceanic zone, although most fronts are detected using both temperature and density criteria, indicating a significant contribution of temperature to density in this region. [ABSTRACT FROM AUTHOR]

Because so many ancestral populations, at various times, have lived within the constraints of earthen architectural technology, it has significant representation in building traditions across large temporal and geographic expanses. Adobe, known also as dagga, ferey, cob, and by other names, is a variant in which clays and other sediments are combined with organic materials and formulated into discrete construction components, often in communities of practice for which adobe recipes, preparation, and application are integral to daily intersections of home and community. For archaeologists, community partners, and interested publics who wish to learn more about it, a large portion of this architectural culture is no longer visible above the surface but is accessible through archaeology. Yet low impact sampling, such as probing and test excavation, rarely reveals adobe features. Ground-penetrating radar (GPR) allows a tantalizing and non-invasive perspective on subsurface adobe features; increasingly there is a buried structural landscape emerging at sites where adobe architecture prevailed. The case studies presented here from 18th and 19th century sites in central and southern California serve as a guide for further survey. [ABSTRACT FROM AUTHOR]

In this study, we present a nationwide machine learning model for hourly PM 2.5 estimation for the continental United States (US) using high temporal resolution Geostationary Operational Environmental Satellites (GOES-16) Aerosol Optical Depth (AOD) data, meteorological variables from the European Center for Medium Range Weather Forecasting (ECMWF) and ancillary data collected between May 2017 and December 2020. A model sensitivity analysis was conducted on predictor variables to determine the optimal model. It turns out that GOES16 AOD, variables from ECMWF, and ancillary data are effective variables in PM 2.5 estimation and historical reconstruction, which achieves an average mean absolute error (MAE) of 3.0 μ g/m 3 , and a root mean square error (RMSE) of 5.8 μ g/m 3 . This study also found that the model performance as well as the site measured PM 2.5 concentrations demonstrate strong spatial and temporal patterns. Specifically, in the temporal scale, the model performed best between 8:00 p.m. and 11:00 p.m. (UTC TIME) and had the highest coefficient of determination (R 2 ) in Autumn and the lowest MAE and RMSE in Spring. In the spatial scale, the analysis results based on ancillary data show that the R 2 scores correlate positively with the mean measured PM 2.5 concentration at monitoring sites. Mean measured PM 2.5 concentrations are positively correlated with population density and negatively correlated with elevation. Water, forests, and wetlands are associated with low PM 2.5 concentrations, whereas developed, cultivated crops, shrubs, and grass are associated with high PM 2.5 concentrations. In addition, the reconstructed PM 2.5 surfaces serve as an important data source for pollution event tracking and PM 2.5 analysis. For this purpose, from May 2017 to December 2020, hourly PM 2.5 estimates were made for 10 km by 10 km and the PM 2.5 estimates from August through November 2020 during the period of California Santa Clara Unite (SCU) Lightning Complex fires are presented. Based on the quantitative and visualization results, this study reveals that a number of large wildfires in California had a profound impact on the value and spatial-temporal distributions of PM 2.5 concentrations. [ABSTRACT FROM AUTHOR]

Fire severity represents fire-induced environmental changes and is an important variable for modeling fire emissions and planning post-fire rehabilitation. Remotely sensed fire severity is traditionally evaluated using the differenced normalized burn ratio (dNBR) derived from multispectral imagery. This spectral index is based on bi-temporal differenced reflectance changes caused by fires in the near-infrared (NIR) and short-wave infrared (SWIR) spectral regions. Our study aims to evaluate the spectral sensitivity of the dNBR using hyperspectral imagery by identifying the optimal bi-spectral NIR SWIR combination. This assessment made use of a rare opportunity arising from the pre- and post-fire airborne image acquisitions over the 2013 Rim and 2014 King fires in California with the Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) sensor. The 224 contiguous bands of this sensor allow for 5760 unique combinations of the dNBR at a high spatial resolution of approximately 15 m. The performance of the hyperspectral dNBR was assessed by comparison against field data and the spectral optimality statistic. The field data is composed of 83 in situ measurements of fire severity using the Geometrically structured Composite Burn Index (GeoCBI) protocol. The optimality statistic ranges between zero and one, with one denoting an optimal measurement of the fire-induced spectral change. We also combined the field and optimality assessments into a combined score. The hyperspectral dNBR combinations demonstrated strong relationships with GeoCBI field data. The best performance of the dNBR combination was derived from bands 63, centered at 0.962 µm, and 218, centered at 2.382 µm. This bi-spectral combination yielded a strong relationship with GeoCBI field data of R2 = 0.70 based on a saturated growth model and a median spectral index optimality statistic of 0.31. Our hyperspectral sensitivity analysis revealed optimal NIR and SWIR bands for the composition of the dNBR that are outside the ranges of the NIR and SWIR bands of the Landsat 8 and Sentinel-2 sensors. With the launch of the Precursore Iperspettrale Della Missione Applicativa (PRISMA) in 2019 and several planned spaceborne hyperspectral missions, such as the Environmental Mapping and Analysis Program (EnMAP) and Surface Biology and Geology (SBG), our study provides a timely assessment of the potential and sensitivity of hyperspectral data for assessing fire severity. [ABSTRACT FROM AUTHOR]

A new approach is proposed to derive evapotranspiration (E) and irrigation requirements by implementing the combination equation models of Penman–Monteith and Shuttleworth and Wallace with surface parameters and resistances derived from Sentinel-2 data. Surface parameters are derived from Sentinel-2 and used as an input in these models; namely: the hemispherical shortwave albedo, leaf area index and water status of the soil and canopy ensemble evaluated by using a shortwave infrared-based index. The proposed approach has been validated with data acquired during the GRAPEX (Grape Remote-sensing Atmospheric Profile and Evapotranspiration eXperiment) in California irrigated vineyards. The E products obtained with the combination equation models are evaluated by using eddy covariance flux tower measurements and are additionally compared with surface energy balance models with Landsat-7 and -8 thermal infrared data. The Shuttleworth and Wallace (S-W S-2) model provides an accuracy comparable to thermal-based methods when using local meteorological data, with daily E errors < 1 mm/day, which increased from 1 to 1.5 mm/day using meteorological forcing data from atmospheric models. The advantage of using the S-W S-2 modeling approach for monitoring ET is the high temporal revisit time of the Sentinel-2 satellites and the finer pixel resolution. These results suggest that, by integrating the thermal-based data fusion approach with the S-W S-2 modeling scheme, there is the potential to increase the frequency and reliability of satellite-based daily evapotranspiration products. [ABSTRACT FROM AUTHOR]

Since the initial experiments nearly 50 years ago, techniques for detecting caves using airborne and spacecraft acquired thermal imagery have improved markedly. These advances are largely due to a combination of higher instrument sensitivity, modern computing systems, and processor-intensive analytical techniques. Through applying these advancements, our goals were to: (1) Determine the efficacy of methods designed for terrain analysis and applied to thermal imagery; (2) evaluate the usefulness of predawn and midday imagery for detecting caves; and (3) ascertain which imagery type (predawn, midday, or the difference between those two times) was most informative. Using forward stepwise logistic (FSL) and Least Absolute Shrinkage and Selection Operator (LASSO) regression analyses for model selection, and a thermal imagery dataset acquired from the Mojave Desert, California, we examined the efficacy of three well-known terrain descriptors (i.e., slope, topographic position index (TPI), and curvature) on thermal imagery for cave detection. We also included the actual, untransformed thermal DN values (hereafter "unenhanced thermal") as a fourth dataset. Thereafter, we compared the thermal signatures of known cave entrances to all non-cave surface locations. We determined these terrain-based analytical methods, which described the "shape" of the thermal landscape, hold significant promise for cave detection. All imagery types produced similar results. Down-selected covariates per imagery type, based upon the FSL models, were: Predawn— slope, TPI, curvature at 0 m from cave entrance, as well as slope at 1 m from cave entrance; midday— slope, TPI, and unenhanced thermal at 0 m from cave entrance; and difference— TPI and slope at 0 m from cave entrance, as well as unenhanced thermal and TPI at 3.5 m from cave entrance. We provide recommendations for future research directions in terrestrial and planetary cave detection using thermal imagery. [ABSTRACT FROM AUTHOR]

Water scarcity during severe droughts has profound hydrological and ecological impacts on rivers. However, the drying dynamics of river surface extent during droughts remains largely understudied. Satellite remote sensing enables surveys and analyses of rivers at fine spatial resolution by providing an alternative to in-situ observations. This study investigates the seasonal drying dynamics of river extent in California where severe droughts have been occurring more frequently in recent decades. Our methods combine the use of Landsat-based Global Surface Water (GSW) and global river bankful width databases. As an indirect comparison, we examine the monthly fractional river extent (FrcSA) in 2071 river reaches and its correlation with streamflow at co-located USGS gauges. We place the extreme 2012–2015 drought into a broader context of multi-decadal river extent history and illustrate the extraordinary change between during- and post-drought periods. In addition to river extent dynamics, we perform statistical analyses to relate FrcSA with the hydroclimatic variables obtained from the National Land Data Assimilation System (NLDAS) model simulation. Results show that Landsat provides consistent observation over 90% of area in rivers from March to October and is suitable for monitoring seasonal river drying in California. FrcSA reaches fair (>0.5) correlation with streamflow except for dry and mountainous areas. During the 2012–2015 drought, 332 river reaches experienced their lowest annual mean FrcSA in the 34 years of Landsat history. At a monthly scale, FrcSA is better correlated with soil water in more humid areas. At a yearly scale, summer mean FrcSA is increasingly sensitive to winter precipitation in a drier climate; and the elasticity is also reduced with deeper ground water table. Overall, our study demonstrates the detectability of Landsat on the river surface extent in an arid region with complex terrain. River extent in catchments of deficient water storage is likely subject to higher percent drop in a future climate with longer, more frequent droughts. [ABSTRACT FROM AUTHOR]

Daily evapotranspiration (ETd) plays a key role in irrigation water management and is particularly important in drought-stricken areas, such as California and high-value crops. Remote sensing allows for the cost-effective estimation of spatial evapotranspiration (ET), and the advent of small unmanned aerial systems (sUAS) technology has made it possible to estimate instantaneous high-resolution ET at the plant, row, and subfield scales. sUAS estimates ET using "instantaneous" remote sensing measurements with half-hourly/hourly forcing micrometeorological data, yielding hourly fluxes in W/m2 that are then translated to a daily scale (mm/day) under two assumptions: (a) relative rates, such as the ratios of ET-to-net radiation (Rn) or ET-to-solar radiation (Rs), are assumed to be constant rather than absolute, and (b) nighttime evaporation (E) and transpiration (T) contributions are negligible. While assumption (a) may be reasonable for unstressed, full cover crops (no exposed soil), the E and T rates may significantly vary over the course of the day for partially vegetated cover conditions due to diurnal variations of soil and crop temperatures and interactions between soil and vegetation elements in agricultural environments, such as vineyards and orchards. In this study, five existing extrapolation approaches that compute the daily ET from the "instantaneous" remotely sensed sUAS ET estimates and the eddy covariance (EC) flux tower measurements were evaluated under different weather, grapevine variety, and trellis designs. Per assumption (b), the nighttime ET contribution was ignored. Each extrapolation technique (evaporative fraction (EF), solar radiation (Rs), net radiation-to-solar radiation (Rn/Rs) ratio, Gaussian (GA), and Sine) makes use of clear skies and quasi-sinusoidal diurnal variations of hourly ET and other meteorological parameters. The sUAS ET estimates and EC ET measurements were collected over multiple years and times from different vineyard sites in California as part of the USDA Agricultural Research Service Grape Remote Sensing Atmospheric Profile and Evapotranspiration eXperiment (GRAPEX). Optical and thermal sUAS imagery data at 10 cm and 60 cm, respectively, were collected by the Utah State University AggieAir sUAS Program and used in the Two-Source Energy Balance (TSEB) model to estimate the instantaneous or hourly sUAS ET at overpass time. The hourly ET from the EC measurements was also used to validate the extrapolation techniques. Overall, the analysis using EC measurements indicates that the Rs, EF, and GA approaches presented the best goodness-of-fit statistics for a window of time between 1030 and 1330 PST (Pacific Standard Time), with the Rs approach yielding better agreement with the EC measurements. Similar results were found using TSEB and sUAS data. The 1030–1330 time window also provided the greatest agreement between the actual daily EC ET and the extrapolated TSEB daily ET, with the Rs approach again yielding better agreement with the ground measurements. The expected accuracy of the upscaled TSEB daily ET estimates across all vineyard sites in California is below 0.5 mm/day, (EC extrapolation accuracy was found to be 0.34 mm/day), making the daily scale results from TSEB reliable and suitable for day-to-day water management applications. [ABSTRACT FROM AUTHOR]