Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “time series analysis”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 271 records · Page 15

Southwest Water Resources: Monitoring Surface Water Extents of Remote Stock Ponds in the Southwestern United States Using Earth Observing Systems for Enhanced Water Resources Management

Due to increasingly frequent and severe drought conditions in the southwestern US, land managers and livestock producers need to monitor stock ponds with increasing regularity. The ability to assess stock pond water levels with Earth observing satellite systems would enhance monitoring efforts of partners at the US Forest Service, Arizona Department of Game and Fish, and the Diablo Trust. This study employed Landsat 8 Operational Land Imager (OLI), Sentinel-1 C-band Synthetic Aperture Radar (C-SAR), and Sentinel-2 Multispectral Instrument (MSI) to monitor surface water extent for hundreds of critical stock ponds in Arizona. Using methods adapted from previously developed image processing workflows, this project conducted a time-series analysis to capture seasonal and interannual variations in surface water area between 2013 to 2021. In addition, end users can monitor the surface water extent of stock ponds through the developed Google Earth Engine software tool called Surface Water Identification and Forecasting Tool (SWIFT). SWIFT incorporates the Automated Water Extraction Index, Modified Normalized Difference Water Index, and Tasseled Cap-Wetness Index for optical imagery and the incidence angle, VV and VH polarization bands for Sentinel-1 imagery to detect small water bodies in the study area with an overall accuracy range of 88-93%. These tools will empower our partners to monitor the extents of water in their stock ponds remotely, enabling them to develop data-informed and sustainable management solutions for decades to come.

Rainey Aberle↗

Powder River Basin Water Resources: Mapping Russian Olive in the Powder River Basin to Inform Invasive Species Management

Since its introduction in the late 1800s, Elaeagnus augustifolia (Russian olive) has become a widespread invasive shrub that poses a threat to native riparian species in the United States by competing with native riparian plants for space and resources. To date, limited information on the distribution of Russian olive in the Powder River Basin of Montana and Wyoming have hampered management efforts and decision making. Here, we detect and model the distribution of Russian Olive using field surveys, ocular sampling, and variables from Landsat 8 Operational Land Imager (OLI), Sentinel-2 MultiSpectral Instrument (MSI), and Shuttle Radar Topography Mission (SRTM) using the Random Forest algorithm. We derived topographic, spectral, and hydrological variables from Landsat 8 OLI, Sentinel-2 MSI, and SRTM to utilize as model inputs. The team was able to successfully create a spectral Russian olive detection map for the Powder River Basin (RMSE =15.44%, R2 = 0.6482). The team also examined change in stream channel geomorphology from 1984-2020 in a time-series analysis using Landsat visible imagery and the RivMap MATLAB package and found little change. Our results will help our partners at the Powder River County Weed Board, Gay Ranch, United States Geological Survey, and University of Northern Colorado to locate and prioritize areas for riparian habitat restoration and to understand the region’s hydrology and geomorphology.

Catherine Buczek↗

Mesoamerica Ecological Forecasting: Assessing Land Cover Change to Inform Management Planning for the Mesoamerican Biological Corridor

In 1992, Central America and Mexico drew up an agreement to establish the Mesoamerican Biological Corridor (MBC) which defines natural corridors to connect nearly 600 protected areas. The MBC is home to 9% of the world's terrestrial species on 0.7% of the world's landmass, yet this biodiverse area has been impacted by great levels of deforestation. The MBC supports protected areas and the important conservation efforts that are tied into the area’s economic and sustainable development. The NASA DEVELOP team partnered with NASA SERVIR, Sistema de la Integración Centroamericana (SICA), Tropical Agriculture Research and High Education Center (CATIE), and Ministries of the Environment for Costa Rica, El Salvador, and Guatemala to assess forest cover change in the MBC. While the southern states of Mexico are included in the MBC, the team excluded Mexico in this study. The team acquired data from Landsat 5 Thematic Mapper (TM), Landsat 7 Enhanced Thematic Mapper Plus (ETM+), Landsat 8 Operational Land Imager (OLI), and Landsat 9 OLI-2 to develop a forest versus non-forest classification. This classification was used to create a Land Use Land Cover Change (LULC) trend map and Deforestation Detection Time Series analysis between 1992 and 2022. The team found that minimum distance classification was the most effective classifier for the project scope. Analysis showed that 6.18% of the study area experienced forest loss and 10.99% experienced forest growth. These observations will help partners visualize the evolution and severity of deforestation and allow decision making for future land management and transboundary conservation efforts.

Hanna Jung↗

Evaluating SAR Radiometric Terrain Correction Solutions: Optimal products for applied users

Operational applications of Synthetic Aperture Radar (SAR) are under development around the world, driven by the regularly-acquired, free-and-open source C-band SAR observations provided by ESA’s Sentinel-1 sensor constellation since 2014. Groups like SERVIR, a joint NASA and USAID initiative, are at the forefront of remote sensing applications for societal benefit. A takeaway from SERVIR’s experience is the need for appropriately geocoded and fully calibrated SAR data that is ready to use for a range of ecosystems-related applications. Radiometric Terrain Corrected (RTC) data are key entry-level products for multiple applications that range from ecosystems to hazards. This work fills a gap in current research by evaluating several RTCs produced by open-source software solutions (SNAP-7 and ISCE-2), the gold standard commercial software (GAMMA), a Google Earth Engine (GEE) based workflow, and the uncorrected GRD products currently available in GEE. RTCs were analyzed for geolocation quality, absolute radiometric calibration, and fidelity of the radiometric terrain flattening over ten sites representing varied terrains. In addition, a time series analysis was conducted over two locations. Overall, no significant differences for radiometric calibration were found across RTC products. However, all RTCs performed better than uncorrected GRD products. The main differences between products were found in geolocation quality. These results not only demonstrate the need for the uptake and distribution of RTC products for ecosystems applications, but demonstrate the ability to do so with open source methods, adding value to developing affordable operational applications.

Helen Blue Parache↗

Mangroves Cover Change Trajectories 1984-2020: The Gradual Decrease of Mangroves in Colombia

Awareness of the significant benefits of mangroves to human lives and their role in regulating environmental processes has increased during the recent decades. Yet there remains significant uncertainty about the mangrove change trajectories and the drivers of change at national scales. In Colombia, the absence of historical satellite imagery and persistent cloud cover have impeded the accurate mapping of mangrove extent and change over time. We create a temporally consistent Landsat-derived dataset using the LandTrendr algorithm to track the historical land cover and mangrove conversion from 1984-2020 across Colombia. Over this period, mangrove extent decreased by ~48.000ha (14% of total mangrove area). We find a gradual reduction of mangrove extent along the Pacific coast since 2004, whereas, in the Caribbean, mangrove cover declined around during 1984-1988 and also after 2012. Our time-series analysis matches with drivers of mangrove change at three local sites. For instance, hydroclimatic events, dredging activities, and high sediment loads transported by the rivers have collectively improved mangrove recovery in some sites. In contrast, human activities pressure linked to agricultural expansion and road construction have degraded mangroves. The transition from dense mangrove to other vegetation types is the most significant conversion affecting mangrove cover in Colombia, impacting an area of 38,469 ± 2,829 ha. We anticipate increased mangrove loss, especially along the Pacific coast, resulting from intensified human activity. Prioritization of conservation areas is needed to support local institutions, maintain currently protected areas, and develop strategies (e.g. payment for ecosystem services) to preserve one of the most pristine mangrove regions in the Western Hemisphere.

Landsat↗

Shoshone River Water Resources II​: Quantifying Sediment Input in the Shoshone River in Wyoming using the Soil and Water Assessment Tool for Enhanced Water Quality Monitoring

The Willwood Dam, an irrigation diversion dam located on the Shoshone River (Wyoming, USA), has faced ongoing issues with sediment accumulation and needs frequent sediment flushing to remain operable. However, high suspended sediment levels during flushing events have negatively impacted downstream aquatic ecology and recreational opportunities. To address these problems, DEVELOP partnered with the Wyoming Department of Environmental Quality, Shoshone River Partners, and United States Geologic Survey. During term one, the team developed a workflow to map turbidity using PlanetScope imagery, analyze time series precipitation data, and create landcover maps. For this term, we focused on three methods to gain a better understanding of sediment sources: 1) improving remote sensing of turbidity, 2) modeling sediment transport within the watershed using the Soil and Water Assessment Tool (SWAT), and 3) conducting a snow cover time series analysis using Suomi NPP VIIRS imagery. Through remote sensing, we found Dry Creek/Homesteader Creek and Penney Gulch/Rough Gulch had the highest concentration sediment plumes. The SWAT+ model created a high-resolution grid model of the watershed, identifying high sedimentation in the western and southern subbasins but displayed low-correlated calibration and validation results due to limited observed data. We analyzed the snow cover extent alongside other hydrologic variables and found that snow melt events correlated with increases in suspended sediment concentration and turbidity values. Coupling remote sensing with hydrological modeling will give watershed managers a new perspective on high-priority regions for implementation of remediation to reduce the sediment build-up at the Willwood Dam.

Suspended Sediment↗

Aircraft Parameter Estimation Considering Process and Measurement Noise

A practical formulation is proposed for parameter estimation using the filter-error method, which is a maximum-likelihood estimator for dynamic systems having both process and measurement noise inputs. The novelty of the proposed formulation is that by accurately estimating the measurement noise covariance matrix using a time series analysis method, the remaining unknowns (which include the unknown parameters in the state-space matrices and the process noise covariance matrix) become decorrelated and can be estimated simultaneously in a straightforward manner. The approach is demonstrated using simulation data and flight test data from a subscale airplane. Results indicate that proposed algorithm can obtain accurate modeling results when both measurement noise and process noise are present in the data.

Kalman filter↗

Integration of Total-Sky Imager Data with a Physics-Based Smart Persistence Model for Intra-Hour Forecasting of Solar Radiation

Short-term solar forecasting models based solely on global horizontal irradiance (GHI) measurements are often unable to discriminate the forecasting of the factors affecting GHI from those that can be precisely computed by atmospheric models. Our previous study introduced a Physics-based Smart Persistence model for Intra-hour forecasting of solar radiation (PSPI) that decomposed the forecasting of GHI into the computation of extraterrestrial solar radiation and solar zenith angle and the forecasting of cloud albedo and cloud fraction. The extraterrestrial solar radiation and solar zenith angle were accurately computed by the Solar Position Algorithm (SPA) developed at the National Renewable Energy Laboratory (NREL). A cloud retrieval technique was used to estimate cloud albedo and cloud fraction from surface-based observations of GHI. With the assumption of persistent cloud structures, the cloud albedo and cloud fraction were predicted for future time steps using a two-stream approximation and a 5-minute exponential weighted moving average, respectively. The model evaluation indicated the estimation and forecast of cloud fraction mostly contributed to the uncertainty of the PSPI though it overcame the persistence and smart persistence models in all forecast time horizons between 5 and 60 minutes. This study aims to enhance the PSPI by ingesting surface-based observations of cloud fraction from a total sky imager (TSI). The estimation and forecast of cloud albedo is correspondingly improved by utilizing the cloud fraction observations and thus leads to more accurate GHI forecast. Various time-series analysis methods are also investigated on the forecasting of cloud fraction and cloud albedo for further improving the GHI forecast. These improvements are valuable for many applications, such as forecasting energy use for buildings, grid operations, and ultimately bringing down the cost of solar energy.

14 SOLAR ENERGY↗

Quantifying Error in Photovoltaic Installation Metadata: Preprint

In this research, we quantify the level of metadata error for a fleet of 2860 photovoltaic (PV) systems, using metadata values provided by fleet owners. Using satellite imagery and time series analysis techniques available in open-source Python packages Panel-Segmentation and PVAnalytics, respectively, we evaluate the accuracy of PV system metadata such as location, azimuth, tilt, and mounting configuration (fixed tilt vs. tracking). We find that approximately 75% of provided latitude-longitude coordinates are within 190 meters of the actual solar installation. We were unable to link 7.8% of latitude-longitude coordinates to any solar installation via satellite imagery analysis. We evaluate the level of error in owner-provided mounting configuration (fixed tilt vs. single-axis tracking), finding only 8 systems with an incorrect mounting configuration. When evaluating azimuth and tilt parameters, we find that approximately 64% of the data is correct, with data for 860 systems (approximately 30%) not provided by system owners. To illustrate the importance of having correct solar metadata, we evaluate how incorrect metadata affects solar performance estimates by modeling system AC energy output at ground-truth vs. incorrect latitude-longitude coordinates, mounting configurations, and azimuth-tilt configurations. Energy output estimates can vary significantly if incorrect metadata parameters are used, with incorrect mounting configuration leading to the largest discrepancy with over 20% variation in expected energy output.

azimuth↗

Nasa Space Flight Instruments: Cost Time Trends

Are NASA’s space flight instruments becoming cheaper or more expensive as time marches forward? After analyzing the costs of hundreds of instruments launched over the last 30 years, the short answer to this question is no… and yes. This paper gives a visual analysis of the cost time trends for various NASA space flight instrument types, such as optical, particles detectors, fields detectors and microwave instruments. In addition to the statistical approaches utilized, such as significance tests, cluster analysis and principle components analysis (PCA), we will also discuss the intangibles which are likely at play, including technological progress, NASA policy and the luck of the draw associated with mission manifests. This analysis was performed as the main driver for the NASA Instrument Cost Model (NICM) recent cost estimating model redesign. Started in 2004, the first version of NICM was based off of instruments launched from 1985-2005, or 20 years’ worth of data. As NICM hit its 10-year anniversary, we wanted to know: should NICM continue to only use the most recent 20 years’ worth of data (1995-2015)? Are instruments becoming cheaper or more expensive as time marches forward? There is evidence in favor of a drop in the median dollar-per-kg value across some instrument types, but little in others. Whereas further research is needed to substantiate, Particles and Optical-Planetary instrument types show moderate to strong evidence of a downward trend in dollar-per-kg. Further research is required to study the nature of this trend (shift, taper, cyclic, etc.). Little evidence for a similar downward trend was detected for Fields or Microwave instruments, or Optical instruments on Earth Orbiting spacecraft. We presented evidence in favor of a drop in the median dollar-per-kg value for Particles and Optical-Planetary instrument types. While similar evidence was weak at best for Fields and Microwave instruments. We can speculate as to the causes for this effect, but we are also equipped to begin to rule out, or at least prioritize, some of the suspected drivers. We observed, for Particles and Optical-Planetary instruments, that perhaps a launch manifest effect was playing part of the role in the observed decrease in dollar-per-kg over the years, noting that the more flagship class missions, which have more money to spend on their instruments, were seen in the earlier years in our data, versus the later years which were dominated by less expensive class missions. However, if this were a dominating driver, would we not have seen the downward trend in the Fields and Microwave instruments as well, which were drawn from that same launch manifest? The fact that we did not observe this helps us rule out the launch manifest effect, and other drivers, such as advances in technology, that seem to be more likely suspects. In that case, however, why would technology advances be helping the Particles and Optical-Planetary instruments only? Why would it not be impacting Optical-Earth Orbiting instruments? Further suspects were looked at as well and ruled out, such as the “Faster, Better, Cheaper”era of NASA development which did not seem to actually impact trends by instrument type on a dollar-perkg scale. VI. Future Work A. Time Series Detailed Statistical Assessment The analysis discussed above sets the foundation for a more rigorous time series analysis of the data. Time series analysis will further explore evidence to-date of time trends for the instrument types which showed the strongest indicators for a decrease in dollar-per-kg: Optical (Planetary) and Particles instruments. More than providing evidence and top-level significance tests, time series analysis would help elucidate what kind of trend that exists in the data, their significance and allow statistically based forecasting (see Figure 10)

Mrozinksi, Joseph↗

Solar Cycle Fine Structure and Surface Rotation from Ca II K-Line Time Series Data

Analysis of three and a half decades of data from the NSO/AFRL/Sac Peak K-line monitoring program yields evidence for four components to the variation: (a) the solar cycle, with considerable fine structure and a quasi-periodicity of 122.4 days; (b) a stochastic process, faster than (a) and largely independent of it, (c) a quasi-periodic signal due to rotational modulation, and of course (d) observational errors (shown to be quite small). Correlation and power spectrum analyses elucidate periodic and aperiodic variation of these chromospheric parameters. Time-frequency analysis is especially useful for extracting information about differential rotation, and in particular elucidates the connection between its behavior and fine structure of the solar cycle on approximately one-year time scales. These results further suggest that similar analyses will be useful at detecting and characterizing differential rotation in stars from stellar light-curves such as those being produced at NASA's Kepler observatory.

Scargle, Jeff↗

Time‐series multi‐omics analysis of micronutrient stress in Sorghum bicolor reveals iron and zinc crosstalk and regulatory network conservation

Micronutrient stress impacts growth, biomass production, and grain yield in crops. Multi-omics studies are valuable resources in identifying genes for functional studies and trait improvement, such as accumulation of Fe or Zn under deficient or excess conditions for bioenergy or grain agriculture. We conducted transcriptomics and ionomics analyses on Sorghum bicolor BTx623, grown under Fe and Zn limited and excess conditions over a 21-day period. To identify early and late transcriptional response in roots and leaves, 180 RNAseq libraries were sequenced for differential expression and co-expression network analyses. Fe and Zn accumulation was measured using ICP-MS at each time point, and a fluorometer was used to estimate chlorophyll content in leaves. Among the four treatments, Fe limitation and Zn excess resulted in the largest phenotypic effects and transcriptional response in roots and leaves. Several of the reduction (Strategy I) and chelation (Strategy II) strategy genes that improve bioavailability of Fe and Zn in plant roots often used by non-grass and grass species, respectively, were differentially expressed. Gene regulatory network (GRN) analysis of roots revealed enrichment of genes from Fe limiting and Zn excess which strongly connect to homologues of SbFIT, SbPYE, and SbBTS as hub genes. The GRN for leaf responses showed homologues of SbPYE and SbBTS as hubs connecting genes for chloroplast biosynthesis, Fe-S cluster assembly, photosynthesis, and ROS scavenging. Expression analyses suggest sorghum uses Strategy II genes for Fe and Zn uptake, as expected, but can also utilize Strategy I genes, which may be advantageous in variable moisture environments. We found strong overlap between Fe and Zn responsive GRNs, indicative of micronutrient crosstalk. We also found conservation of root and leaf GRNs, and known homologous genes suggest strong constraints on homeostasis networks in plants. These data will provide a resource for functional genetics to enhance micronutrient transport in sorghum, and opportunities to conduct further comparative GRN analysis across diverse crops species.

59 BASIC BIOLOGICAL SCIENCES↗

Clustering Analysis of Commercial Vehicles Using Automatically Extracted Features from Time Series Data

Standard of practice approaches to time series cluster analysis involve careful feature engineering, often utilizing expert input to tune and select features by hand. In many cases, expert input may not be readily available, or there may not yet exist a community consensus on the ideal features for a given application. This paper compares the results of several cluster analysis methods, using both hand selected features and those extracted automatically, when applied to large geospatial time series telematics data from commercial trucking fleets. The impacts of feature selection, dimensionality reduction, and choice of clustering algorithm on the quality of clustering results are explored. Results from this analysis confirm prior results that domain agnostic features are competitive with the hand engineered features with respect to clustering quality metrics. These results also provide new insight into the most successful strategies for identifying structure in large unstructured vehicle telematics data, and suggest that time series clustering using automatic feature extraction can be an effective approach to extract structure from large scale geospatial time series data in cases when hand selected features are not available.

33 ADVANCED PROPULSION SYSTEMS↗

Deep Cellular Recurrent Network for Efficient Analysis of Time-Series Data With Spatial Information

Efficient processing of large-scale time series data is an intricate problem in machine learning. Conventional sensor signal processing pipelines with hand engineered feature extraction often involve huge computational cost with high dimensional data. Deep recurrent neural networks have shown promise in automated feature learning for improved time-series processing. However, generic deep recurrent models grow in scale and depth with increased complexity of the data. This is particularly challenging in presence of high dimensional data with temporal and spatial characteristics. Consequently, this work proposes a novel deep cellular recurrent neural network (DCRNN) architecture to efficiently process complex multi-dimensional time series data with spatial information. Here, the cellular recurrent architecture in the proposed model allows for location-aware synchronous processing of time series data from spatially distributed sensor signal sources. Extensive trainable parameter sharing due to cellularity in the proposed architecture ensures efficiency in the use of recurrent processing units with high-dimensional inputs. This study also investigates the versatility of the proposed DCRNN model for classification of multi-class time series data from different application domains. Consequently, the proposed DCRNN architecture is evaluated using two time-series datasets: a multichannel scalp EEG dataset for seizure detection, and a machine fault detection dataset obtained in-house. The results suggest that the proposed architecture achieves state-of-the-art performance while utilizing substantially less trainable parameters when compared to comparable methods in the literature.

60 APPLIED LIFE SCIENCES↗

Variability of African Farming Systems from Phenological Analysis of NDVI Time Series

Food security exists when people have access to sufficient, safe and nutritious food at all times to meet their dietary needs. The natural resource base is one of the many factors affecting food security. Its variability and decline creates problems for local food production. In this study we characterize for sub-Saharan Africa vegetation phenology and assess variability and trends of phenological indicators based on NDVI time series from 1982 to 2006. We focus on cumulated NDVI over the season (cumNDVI) which is a proxy for net primary productivity. Results are aggregated at the level of major farming systems, while determining also spatial variability within farming systems. High temporal variability of cumNDVI occurs in semiarid and subhumid regions. The results show a large area of positive cumNDVI trends between Senegal and South Sudan. These correspond to positive CRU rainfall trends found and relate to recovery after the 1980's droughts. We find significant negative cumNDVI trends near the south-coast of West Africa (Guinea coast) and in Tanzania. For each farming system, causes of change and variability are discussed based on available literature (Appendix A). Although food security comprises more than the local natural resource base, our results can perform an input for food security analysis by identifying zones of high variability or downward trends. Farming systems are found to be a useful level of analysis. Diversity and trends found within farming system boundaries underline that farming systems are dynamic.

Vrieling, Anton↗

Optimal Frequency-domain Analysis for Spacecraft Time Series: Introducing the Missing-data Multitaper Power Spectrum Estimator

While the Lomb–Scargle periodogram is foundational to astronomy, it has a significant shortcoming: the variance in the estimated power spectrum does not decrease as more data are acquired. Statisticians have a 60 yr history of developing variance-suppressing power spectrum estimators, but most are not used in astronomy because they are formulated for time series with uniform observing cadence and without seasonal or daily gaps. Here we demonstrate how to apply the missing-data multitaper power spectrum estimator to spacecraft data with uniform time intervals between observations but missing data during thruster fires or momentum dumps. The F-test for harmonic components may be applied to multitaper power spectrum estimates to identify statistically significant oscillations that would not rise above a white noise–based false alarm probability. Multitapering improves the dynamic range of the power spectrum estimate and suppresses spectral window artifacts. We show that the multitaper–F-test combination applied to Kepler observations of KIC 6102338 detects differential rotation without requiring iterative sinusoid fitting and subtraction. Significant signals reside at harmonics of both fundamental rotation frequencies and suggest an antisolar rotation profile. Next we use the missing-data multitaper power spectrum estimator to identify the oscillation modes responsible for the complex "scallop-shell" shape of the K2 light curve of EPIC 203354381. We argue that multitaper power spectrum estimators should be used for all time series with regular observing cadence.

79 ASTRONOMY AND ASTROPHYSICS↗

Multiple-time scales analysis of physiological time series under neural control

We discuss multiple-time scale properties of neurophysiological control mechanisms, using heart rate and gait regulation as model systems. We find that scaling exponents can be used as prognostic indicators. Furthermore, detection of more subtle degradation of scaling properties may provide a novel early warning system in subjects with a variety of pathologies including those at high risk of sudden death.

Non-NASA Center↗

Characterizing the Effect of Capillary Heterogeneity on Multiphase Flow Pulsation in an Intermediate–Scale Beadpack Experiment using Time Series Clustering and Frequency Analysis

An intermediate-scale beadpack drainage experiment was conducted to investigate how simple layered lamination heterogeneity affects CO 2 flow. Two simple layers of capillary barriers are manually packed in the tank and slow drainage was carried out using analog fluids to mimic the capillary- and gravity-dominated CO 2 upward migration process in deep saline aquifers. Nonwetting phase saturation time series clustering analysis and frequency analysis have been conducted on the experimental data. Additionally, modified invasion percolation numerical simulations were done on a digital model of the beadpack to compare to experimental results. Results show that capillary barriers can lead to strong pulsation behavior, which in turn can cause unexpected early breaching through other barriers. The inlet pressure is found to be able to respond to saturation changes in far regions of the domain, indicating that the wetting phase can transmit pressure changes from the other phase. Although static simulations were not able to capture all the dynamic behavior observed in the experiment, Monte Carlo composite simulation results combining many different realizations can better illustrate how the nonwetting phase will behave in the heterogeneous domain. Furthermore, our results suggest the need for CO 2 storage site selection with preference given to aquifers with more capillary barriers with finer grain sizes to avoid flow pulsation and to retard plume upward migration.

58 GEOSCIENCES↗