Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “covariates”

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 1,243 records · Page 69
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Abiotic influences on continuous conifer forest structure across a subalpine watershed

Understanding the abiotic drivers of high-elevation forest physiognomy is essential for forecasting how mountain ecosystems will respond to emerging environmental pressures. Most prior studies of these relationships have relied on small samples of the full landscape, resulting in limited power to detect dominant covariates and their interactions. Here we report the first evaluation of abiotic influences on a complement of accurate, wall-to-wall estimates of conifer forest structure and composition at the watershed scale. In a subalpine conifer domain in the Colorado Rocky Mountains (USA), we developed a novel method for deriving stand structure metrics from waveform LiDAR data, which showed high fidelity with field inventory. We quantified the relationships between structural and compositional metrics and climate, topographic, edaphic, and geologic factors. Our results showed that peak snow water equivalent (SWE), snow disappearance rate, and elevation explained most of the variation in forest structure. The highest stand density, basal area, maximum canopy height, and quadratic mean diameter occurred in sites with SWE around one standard deviation below mean, but with long snow residence times. Stand density decreased linearly with elevation, while other metrics peaked between 3000 m and 3200 m. Substrate properties had weaker influence. Continuous mapping of through-canopy forest structure enabled our novel findings of the dominant role of snowpack in explaining structural and compositional variation, and of elevation thresholds. Our reproducible approach facilitates assessment of forest-topoclimate relationships in other conifer-dominated landscapes and improves understanding of the baseline patterns controlling forest structure, which is needed for predicting long-term ecological change.

forest ecology↗
DOE OSTI2025

Contributions of vegetation heterogeneity within tower footprint to CO 2 flux estimations through graph neural network modeling

Net ecosystem exchange of CO 2 (Fc) measured directly by eddy covariance towers is based on various assumptions, including large, flat and homogenous land cover type. In reality, often a tower site is not large enough for flux measurements, and landscapes consist of patches of different land cover types within the flux footprint. In addition, some portions of fluxes are contributed by different cover types when a footprint exceeds the size of the target ecosystem. The contributions of non-dominant patches to Fc are often ignored. Here, in this study, we propose a novel integrated modeling framework that combines random forest (RF) and XGBoost with a residual correction module based on a deep graph convolutional network (DeeperGCN) to simulate Fc for seven flux measurement sites in southwest Michigan. High-resolution remote sensing vegetation indices, soil properties, meteorological variables, and footprint-weighted spatial features were used as model inputs at three spatial resolutions (10 m, 20 m, 30 m), and their importance in predicting Fc with DeeperGCN was assessed. We found that residual correction using DeeperGCN significantly improved prediction accuracy, with the R 2 increasing from 0.9098 to 0.9479 for RF and from 0.9235 to 0.9433 for XGBoost. At site level, the maximum improvement in R 2 reached 0.1617. Paired t-tests confirmed that these improvements were statistically significant (p < 0.05). Among all predictors, leaf area index and incoming shortwave radiation emerged as the dominant drivers of spatial residual variation, followed by precipitation, relative humidity, and selected vegetation indices. The 20 m resolution yielded the best balance between model performance and computational efficiency. In conclusion, our modeling framework effectively captures both spatial heterogeneity and nonlinear interactions, offering a robust solution for spatially explicit flux modeling in structurally diverse ecosystems beyond the study sites.

footprint model↗
DOE OSTI2024

Changes in inundation drive carbon dioxide and methane fluxes in a temperate wetland

Wetlands cycle carbon by being net sinks for carbon dioxide (CO 2 ) and net sources of methane (CH 4 ). Daily and seasonal temporal patterns, dissolved oxygen (DO) availability, inundation status (flooded or dry/partially flooded), water depth, and vegetation can affect the magnitude of carbon uptake or emissions, but the extent and interactive effects of these variables on carbon gas fluxes are poorly understood. We characterized the linkages between carbon fluxes and these environmental and temporal drivers at the Old Woman Creek National Estuarine Research Reserve (OWC), OH. We measured diurnal gas flux patterns in an upstream side channel (called the cove) using chamber measurements at six sites (three vegetated and three non-vegetated). We sampled hourly from 7 AM to 7 PM and monthly from July to October 2022. DO concentrations and water levels were measured monthly. Water inundation status had the most influential effect on carbon fluxes with flooded conditions supporting higher CH 4 fluxes (0.39 μmol CH 4 m –2 s –1 ; –1.23 μmol CO 2 m –2 s –1 ) and drier conditions supporting higher CO 2 fluxes (0.03 μmol CH 4 m –2 s –1 ; 0.86 μmol CO 2 m –2 s –1 ). When flooded, the wetland was a net CO 2 sink; however, it became a source for both CH 4 and CO 2 when water levels were low. We compared chamber-based gas fluxes from the cove in flooded (July) and dry (August) months to fluxes measured with an eddy covariance tower whose footprint covers flooded portions of the wetland. The diurnal pattern of carbon fluxes at the tower did not vary with changing water levels but remained a CO 2 sink and a CH 4 source even when the cove where we performed the chamber measurements dried out. Furthermore, these results emphasize the role of inundation status on wetland carbon cycling and highlight the importance of fluctuating hydrologic patterns, especially hydrologic drawdowns, under changing climatic conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Carbon acquisition ecological strategies to connect soil microbial biodiversity and carbon cycling

Soil carbon feedbacks to global change are uncertain, and the biological processes that govern soil organic matter decomposition are not resolved in current ecosystem models. Though it is recognized that microbial biodiversity influences decomposition rates, incorporating this relationship into ecosystem models is challenging because microbial communities are prohibitively diverse. It is likely necessary to distill microbial biodiversity by focusing on functional groups or ecological strategies. The ecological strategies that currently dominate the microbial ecology literature derive from macroecological theory, have clear weaknesses, and have had limited success when applied to predict soil carbon dynamics. Here, we present a new framework for soil microorganisms: Carbon Acquisition Ecological Strategies (CAES), and we outline a path toward incorporating microbial biodiversity into ecosystem models using this framework to enhance predictions of soil carbon feedbacks to global change. Because a microorganism's diet is central to its ecological niche and likely to covary with other ecologically significant traits, we posit that carbon acquisition may serve as a tractable foundation for developing ecological strategies. Further, we describe four candidate ecological strategies for soil microorganisms: 1° decomposers that assimilate complex plant polymers, 2° decomposers that assimilate microbial necromass, passive consumers that assimilate dissolved organic carbon, and predatory microbes that assimilate live microbial biomass. These strategies are directly linked to soil carbon pools currently represented in ecosystem models and may provide a foundation for greater integration of microbial community dynamics into ecosystem models.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Ion Clusters Reveal the Sources, Impacts, and Drivers of Freshwater Salinization

Population growth, land use change, climate change, and natural resource extraction are driving the salinization of freshwater resources worldwide. Reversing these trends will require data-centric approaches that identify salt sources, environmental drivers, and ecosystem responses. In this study, we applied principal component analysis and hierarchical clustering to identify ion covariance patterns, or “ion clusters,” in Broad Run, an urban stream in the Mid-Atlantic United States. These clusters correspond to distinct hydrologic regimes and reveal specific salinization risks: (1) phosphorus pollution mobilized during summer storms (Cluster 1); (2) elevated concentrations of sulfate and bicarbonate during baseflow (Cluster 2), likely reflecting groundwater discharge; and (3) elevated specific conductance and sodium, chloride, and potassium ion concentrations during snowmelt and rain-on-snow events (Cluster 3), driven by deicer and anti-icer wash-off. These ion fingerprints offer a transferable framework for diagnosing salt sources, assessing ecological risk, and identifying management targets. Our findings underscore the need for next-generation stormwater infrastructure and smart growth policies to protect aquatic life in rapidly urbanizing watersheds.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Investigating ultrafine aerosol turbulent fluxes during atmospheric new particle formation events

New particle formation (NPF) is an important source of atmospheric aerosols, yet quantitatively predicting its occurrence remains challenging, partly because ground-based measurements cannot capture NPF accurately if they occur in the upper atmosphere. While observations have confirmed the presence of new particles near the top of the planetary boundary layer (PBL), their origins and the direction of their vertical transport have remained ambiguous. Here, we propose and validate an analytical framework using airborne eddy-covariance (EC) and continuous wavelet transform (CWT) to directly quantify the vertical flux of newly formed particles and determine its dominant vertical transportation direction. Analyzing data from a dedicated airborne campaign over the Southern Great Plains, we observed a persistent and strong downward particle flux during NPF events, with a mean value of 133.8 cm?3 m s?¹ downward flux in the entrainment zone, whereas fluxes on non-NPF days were negligible. Our spectral analysis further confirms that these directional fluxes can be reliably captured using standard 1 Hz aerosol instrumentation. These findings suggest that new particle formation, driven by the entrainment of air from the overlying residual/stable layer during PBL growth, is a significant and potentially understudied source of boundary layer aerosols. The framework presented here provides a methodology to correctly attribute NPF events to specific altitudes, thereby improving the mechanistic understanding required for accurate atmospheric models.

Zhang, Ruoyu↗
DOE OSTI2020

Tracking Seasonal and Interannual Variability in Photosynthetic Downregulation in Response to Water Stress at a Temperate Deciduous Forest

Abstract The understanding and modeling of photosynthetic dynamics affected by climate variability can be highly uncertain. In this paper, we examined a well‐characterized eddy covariance site in a drought‐prone temperate deciduous broadleaf forest combining tower measurements and satellite observations. We find that an increase in spring temperature usually leads to enhanced spring gross primary production (GPP), but a GPP reduction in late growing season due to water limitation. We evaluated how well a coupled fluorescence‐photosynthesis model (SCOPE) and satellite data sets track the interannual and seasonal variations of tower GPP from 2007 to 2016. In SCOPE, a simple stress factor scaling of Vcmax as a linear function of observed predawn leaf water potential ( ψ pd ) shows a good agreement between modeled and measured interannual variations in both GPP and solar‐induced chlorophyll fluorescence (SIF) from the Global Ozone Monitoring Experiment‐2 (GOME‐2). The modeled and satellite‐observed changes in SIF yield are ~30% smaller than corresponding changes in light use efficiency (LUE) under severe stress, for which a common linear SIF to GPP scaling would underestimate the stress reduction in GPP. Overall, GOME‐2 SIF tracks interannual tower GPP variations better than satellite vegetations indices (VIs) representing canopy “greenness.” However, it is still challenging to attribute observed SIF variations unequivocally to greenness or physiological changes due to large GOME‐2 footprint. Higher‐resolution SIF data sets (e.g., TROPOMI) already show the potential to well capture the downregulation of late‐season GPP and could pave the way to better disentangle canopy structural and physiological changes in the future.

He, Liyin↗
DOE OSTI2020

A Nonequilibrium Thermodynamic Approach for Surface Energy Balance Closure

Abstract The surface energy imbalance, viz. that the turbulent dissipation does not fully account for available energy, has for long been an outstanding challenge in geophysical studies. In this study, we developed a novel approach based on nonequilibrium thermodynamics by representing the atmospheric boundary layer as a heat engine. In addition, an analytically tractable approach was used to estimate the ground heat flux based on Green's function approach, which in turn determines the available energy that drives the atmospheric heat engine. The proposed model was evaluated using heat fluxes measured by eight AmeriFlux eddy covariance towers with atmospheric temperature profiles recorded at adjacent radiosonde sites. The surface energy balance closure can be improved by ~11% over various landscapes, by including the estimated power production from the atmospheric heat engine.

Li, Peiyuan↗
DOE OSTI2019

The Correlation of Mesoscale Humidity Anomalies With Mesoscale Organization of Marine Stratocumulus From Observations Over the ARM Eastern North Atlantic Site

The structure of mesoscale cellular organization in marine cloud-topped boundary layers is found to be well characterized by using a novel compositing approach based on mesoscale variations in the column water vapor path (WVP). The approach is applied to ground-based observations from the Atmospheric Radiation Measurement (ARM) Eastern North Atlantic site. Based on a set of satellite and ground-based observational criteria, 381 hr of closed-cell and 227 hr of open-cell cases were selected from late 2015 to early 2018. Strong correlations are found between mesoscale-filtered cloud properties and WVP for both closed-cell and open-cell regimes. In the moist atmospheric columns, the clouds are thicker with higher tops and lower bases, and stronger precipitation compared to that in the dry columns. Overall, cloud properties of open and closed cells covary similarly with mesoscale moisture perturbations, except that the correlation of liquid water path with WVP is much weaker for open cells. The relations of subcloud properties with WVP are also examined. In the moist columns, surface equivalent potential temperature and relative humidity are higher than in the dry columns. A marginally negative correlation between mesoscale-filtered subcloud turbulence eddy dissipation rate and WVP is found for closed cells. These results are consistent with the conceptual model of closed cellular mesoscale circulations of Zhou and Bretherton (2019).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Wildfire‐Smoke Aerosols Lead to Increased Light Use Efficiency Among Agricultural and Restored Wetland Land Uses in California's Central Valley

Abstract There are few observational studies measuring the ecosystem‐scale productivity effects of changes in incident diffuse photosynthetically active radiation (PAR diffuse ), especially related to wildfire smoke. Climate change‐induced increases to the duration and intensity of fire conditions have made smoke a common occurrence across western North America, with largely unquantified ecosystem feedbacks. Under equivalent amounts of radiation, increased atmospheric particulate matter could lead to a boost in productivity as scattering redistributes photons throughout multilayer canopies. In this work, we leverage a meso‐network of eddy covariance measurement sites across a unique array of managed and restored C 3 and C 4 canopy types to understand how recent wildfire smoke affected ecosystem productivity during the summer of 2018, an especially smoky year in the agriculturally productive Central Valley. We find that diffuse PAR diffuse increased by more than a third compared to the previous growing season, while total PAR was only slightly diminished. These conditions caused nearly a doubling of light use efficiency over the range of diffuse fraction observed, with the highest sensitivity to diffuse fraction exhibited by corn and alfalfa crops. We utilized an empirical model to assess the trade‐off between enhanced diffuse fraction and reduced total PAR. Under mean radiation conditions, daily integrated gross ecosystem productivity increased by 1.2–4.2% compared to the previous growing season. Finally, we explore the potential negative effect of heightened ozone, a copollutant often associated with wildfire. In addition to the effects of wildfire smoke, the results of this natural experiment can help validate future predictions of aerosol‐productivity feedbacks.

Hemes, Kyle S.↗
DOE OSTI2020

Synoptic Meteorology Explains Temperate Forest Carbon Uptake

Abstract While substantial attention has been paid to the effects of both global climate oscillations and local meteorological conditions on the interannual variability of ecosystem carbon exchange, the relationship between the interannual variability of synoptic meteorology and ecosystem carbon exchange has not been well studied. Here we use a clustering algorithm to identify a summertime cyclonic precipitation system northwest of the Great Lakes to determine (a) the association at a daily scale between the occurrence of this system and the local meteorology and net ecosystem exchange at three Great Lakes region forested eddy covariance sites and (b) the association between the seasonal prevalence of this system and the summertime net ecosystem exchange of these sites. We find that temperature, in addition to precipitation and cloud cover, is an important explanatory factor for the suppression of net ecosystem productivity that occurs during these cyclonic events in this region. In addition, the prevalence of this cyclonic system can explain a significant proportion of the interannual variability in summertime forest ecosystem exchange in this region. This explanatory power is not due to a simple accumulation of low‐productivity days that cooccur with this meteorological event, but rather a broader association between the frequency of these events and several aspects of prevailing seasonal conditions. This work demonstrates the usefulness of conceptualizing meteorology in terms of synoptic systems for explaining the interannual variability of regional carbon fluxes.

Randazzo, Nina A.↗
DOE OSTI2020

Intra-Annual and Interannual Dynamics of Evaporation Over Western Lake Erie

Evaporation (E) is a critical component of the water and energy budget in lake systems yet is challenging to quantify directly and continuously. We examined the magnitude and changes of E and its drivers over Lake Erie - the shallowest and most southern lake of the Laurentian Great Lakes. We deployed two eddy-covariance tower sites in the western Lake Erie Basin - one located nearshore (CB) and one offshore (LI) - from September 2011 through May 2016. Monthly E varied from 5 to 120 mm, with maximum E occurring in August–October. The annual E was 635 ± 42 (±SD) mm at CB and 604 ± 32 mm at LI. Mean winter (October–March) E was 189 ± 61 mm at CB and 178 ± 25 mm at LI, accounting for 29.8% and 29.4% of annual E. Mean daily E was 1.8 mm during the coldest month (January) and 7.4 mm in the warmest month (July). Monthly E exhibited a strong positive linear relationship to the product of wind speed and vapor pressure deficit. Pronounced seasonal patterns in surface energy fluxes were observed with a 2-month lag in E from R n , due to the lake's heat storage. This lag was shorter than reports regarding other Great Lakes. Difference in E between the offshore and nearshore sites reflected within-lake spatial heterogeneity, likely attributable to climatic and bathymetric differences between them. These findings suggest that predictive models need to consider lake-specific heat storage and spatial heterogeneity in order to accurately simulate lake E and its seasonal dynamics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Joint CO 2 Mole Fraction and Flux Analysis Confirms Missing Processes in CASA Terrestrial Carbon Uptake Over North America

Terrestrial biosphere models (TBMs) play a key role in the detection and attribution of carbon cycle processes at local to global scales and in projections of the coupled carbon-climate system. TBM evaluation commonly involves direct comparison to eddy-covariance flux measurements. This study uses atmospheric CO 2 mole fraction ([CO 2 ]) measured in situ from aircraft and tower, in addition to flux-measurements from summer 2016 to evaluate the CASA TBM. WRF-Chem is used to simulate [CO 2 ] using biogenic CO 2 fluxes from a CASA parameter-based ensemble and CarbonTracker version 2017 (CT2017) in addition to transport and CO 2 boundary condition ensembles. The resulting “super ensemble” of modeled [CO 2 ] demonstrates that the biosphere introduces the majority of uncertainty to the simulations. Both aircraft and tower [CO 2 ] data show that the CASA ensemble net ecosystem exchange (NEE) of CO 2 is biased high (NEE too positive) and identify the maximum light use efficiency E max a key parameter that drives the spread of the CASA ensemble in summer 2016. These findings are verified with flux-measurements. The direct comparison of the CASA flux ensemble with flux-measurements confirms missing sink processes in CASA. Separating the daytime and nighttime flux, we discover that the underestimated net uptake results from missing sink processes that result in overestimation of respiration. NEE biases are smaller in the CT2017 posterior biogenic fluxes, which assimilates observed [CO 2 ]. Flux tower analyses, however, reveal an unrealistic overestimation of nighttime respiration in CT2017 due to the limitation of inversion strategy.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Enlarged Nonclosure of Surface Energy Balance With Increasing Atmospheric Instabilities Linked to Changes in Coherent Structures

The nonclosure of surface energy balance remains an outstanding problem in eddy covariance (EC) measurements of land-surface fluxes of heat, water vapor, and CO 2 . Our data collected from an EC tower over a semiarid sagebrush ecosystem indicate that under unstable atmospheric conditions, the nonclosure becomes increased with increasing instability, consistent with many other studies. It is demonstrated that the increased nonclosure is not caused by the inadequate sampling of large-scale turbulent motions using a 30-min averaging interval, though the scales of turbulent motions dominating sensible and latent heat fluxes become enlarged with increasing instability. Quadrant analysis is then used to reveal that the flux contributions from ejections remain nearly constant with increasing instability, whereas the flux contributions from sweeps are reduced and their time fractions increase. Our results imply that the increased nonclosure of surface energy balance is associated with changes in turbulent structures including their dominant time scales and flux contributions of ejections and sweeps as the atmospheric instability increases, which require further studies using vertically distributed observations and/or large eddy simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Groundwater regulates interannual variations in evapotranspiration in a riparian semiarid ecosystem

Dryland ecosystems are an important component of Earth’s carbon, water, and energy budgets. Changes in hydroclimatic conditions can significantly alter land surface fluxes in these ecosystems. In regions with groundwater-river water exchange, land surface fluxes are strongly modulated by groundwater table fluctuations in response to river stage variations; however, interactions between fluxes and groundwater remain poorly understood in drylands. Here, we use eddy covariance data collected in paired upland and riparian semiarid ecosystems during two years with contrasting river flow conditions to examine interactions among river stage conditions, water table dynamics, meteorological drivers, and ecosystem fluxes. Elevated groundwater availability resulting from groundwater-surface water exchange at the riparian site enhanced the latent heat flux (LE, the energy form of evapotranspiration), which was sustained longer in the growing season during both years. At the upland site, mean LE during the dry season differed by less than 2% between the two years. In contrast, at the riparian site contrasting river flow conditions resulted in a 41% difference in peak monthly LE between the two years. The two-year mean LE in the summer at the riparian site was twice that at the upland site. Carbon uptake (NEP) and LE were well correlated at the upland site during both years, but at the riparian site enhanced summertime LE during the high river flow year did not correspond t enhanced NEP. Our results demonstrate that enhanced groundwater availability modulated by river stage variability is a critical driver of land surface fluxes in semiarid riparian ecosystems.

Missik, Justine E.↗
DOE OSTI2021

Atmospheric Rivers Impacting Northern California Exhibit a Quasi-Decadal Frequency

Periods of water surplus and deficit in Northern California follow a pronounced quasi-decadal cycle. This cycle is largely driven by the frequency of atmospheric rivers (ARs), affecting the region’s wet and dry periods. Our analyses demonstrate that the quasi-decadal cycle of AR frequency relies on moisture transport associated with the position and intensity of the Aleutian Low. In observations, the Aleutian Low is shown to covary with tropical Pacific sea surface temperature anomalies. A modeling experiment, which incorporates ocean observations from the equatorial Pacific into the fully coupled climate model, provides support that the quasi-decadal cycle of the Aleutian Low is forced by the tropical Pacific. Subsequently, the tropical Pacific modulates the wet season moisture transport toward California on decadal time scales, affecting AR frequency. These results provide metrics for improving interannual-to-decadal prediction of AR activity, which drives hydrological cycles in Northern California.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Tundra Underlain By Thawing Permafrost Persistently Emits Carbon to the Atmosphere Over 15 Years of Measurements

Warming of the Arctic can stimulate microbial decomposition and release of permafrost soil carbon (C) as greenhouse gases, and thus has the potential to influence climate change. At the same time, plant growth can be stimulated and offset C release. This study presents a 15-year time series comprising chamber and eddy covariance measurements of net ecosystem C exchange in a tundra ecosystem in Alaska where permafrost has been degrading due to regional warming. The site was a carbon dioxide source to the atmosphere with a cumulative total loss of 781.6 g C m -2 over the study period. Both gross primary productivity (GPP) and ecosystem respiration (R eco ) were already likely higher than historical levels such that increases in R eco losses overwhelmed GPP gains in most years. This shift to a net C source to the atmosphere likely started in the early 1990s when permafrost was observed to warm and thaw at the site. Shifts in the plant community occur more slowly and are likely to constrain future GPP increases as compared to more rapid shifts in the microbial community that contribute to increased Reco. Observed rates suggest that cumulative net soil C loss of 4.18–10.00 kg C m -2 —8%–20% of the current active layer soil C pool—could occur from 2020 to the end of the century. This amount of permafrost C loss to the atmosphere represents a significant accelerating feedback to climate change if it were to occur at a similar magnitude across the permafrost region.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Scaling Analysis of Two–Phase Flow in Fractal Permeability Fields

Fluid mixing in permeable media is essential in many practical applications. The mixing process is a consequence of velocity fluctuations owing to geological heterogeneities and mobility contrast of fluids. Heterogeneities in natural rocks are often spatially correlated, and their properties, such as permeability, may be described using fractal distributions. This work models the fractal characteristics of such permeability fields in which the covariance function is expressed as a power-law function. A generalized scaling relation is derived relating various fractal permeability fields using the magnitude of their fluctuations. Here, this relation reveals the self-similar behavior of two-phase flow in such permeable media. To that end, a recently developed, high-resolution numerical simulator is employed to validate the analytically derived scaling relations. Two flow problems are considered in which flow is governed by 1) a linear, and 2) a nonlinear transport equation. Due to the probabilistic representation of the fractal permeability fields, a sensitivity study is conducted for each flow scenario to determine the number of realizations required for statistical convergence. Scaling analysis is performed using ensemble averages of simulated saturation profiles and their mixing lengths. Results support the validity of the developed scaling relation across the range of investigated flow conditions at intermediate times. The dynamics of linear flow in the asymptotic regime is affected by the correlation structure of heterogeneity. In nonlinear flow, scaling behavior appears to be dominated by the degree of nonlinearity.

58 GEOSCIENCES↗