SEARCH · Engineering Papers

Results for “inverse modeling”

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.

SourceResearch connectionFrom yearThrough year

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

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

NASA NTRS2019

TPSAS-NF1676L-34010-DND

We invert global black carbon (BC), organic carbon (OC) and desert dust (DD) aerosol emissions from POLDER/PARASOL spectral aerosol optical depth (AOD) and aerosol absorption optical depth (AAOD) using the GEOSChem inverse modelling framework. Our inverse modeling framework uses standard a priori emissions to provide a posteriori emissions that are constrained by POLDER/PARASOL AODs and AAODs. Our satellite-derived global aerosol emission database indicates an increase of BC and OC, while a decrease of DD with respect to the a priori emission inventories used in standard GEOS-Chem runs. The model simulations using a posteriori emissions (i.e. retrieved emissions) provide new estimate of global aerosol extinction and absorption constrained by observations. Additionally, the a posteriori model simulation of AOD, AAOD, single scattering albedo, Angstrom exponent, and absorption Angstrom exponent show better agreement with independent AERONET, MODIS, and OMI measurements than the a priori simulation. Thus, this study suggests that using satellite-constrained global aerosol emissions in aerosol transport models can improve the accuracy of simulated global aerosol properties and have potential to better understanding global aerosol effects on climate and human health.

Cheng Chen↗

NASA NTRS2019

Constraining global aerosol emissions using POLDER/PARASOL satellite remote sensing observations

We invert global black carbon (BC), organic carbon (OC) and desert dust (DD) aerosol emissions from POLDER/PARASOL spectral aerosol optical depth (AOD) and aerosol absorption optical depth (AAOD) using the GEOS-Chem inverse modeling framework. Our inverse modeling framework uses standard a priori emissions to provide a posteriori emissions that are constrained by POLDER/PARASOL AODs and AAODs. The following global emission values were retrieved for the three aerosol components: 18.4 Tg/yr for BC, 109.9 Tg/yr for OC and 731.6 Tg/yr for DD for the year 2010. These values show a difference of +166.7 %, +184.0 % and −42.4 %, respectively, with respect to the a priori values of emission inventories used in “standard” GEOS-Chem runs. The model simulations using a posteriori emissions (i.e., retrieved emissions) provide values of 0.119 for global mean AOD and 0.0071 for AAOD at 550 nm, which are +13.3 % and +82.1 %, respectively, higher than the AOD and AAOD obtained using the a priori values of emissions. Additionally, the a posteriori model simulation of AOD, AAOD, single scattering albedo, Ångström exponent and absorption Ångström exponent show better agreement with independent AERONET, MODIS and OMI measurements than the a priori simulation. Thus, this study suggests that using satellite-constrained global aerosol emissions in aerosol transport models can improve the accuracy of simulated global aerosol properties.

chemical transport model (CTM) simulations↗

DOE OSTI2023

Massively parallel modeling and inversion of electrical resistivity tomography data using PFLOTRAN

Abstract. Electrical resistivity tomography (ERT) is a broadly accepted geophysical method for subsurface investigations. Interpretation of field ERT data usually requires the application of computationally intensive forward modeling and inversion algorithms. For large-scale ERT data, the efficiency of these algorithms depends on the robustness, accuracy, and scalability on high-performance computing resources. In this regard, we present a robust and highly scalable implementation of forward modeling and inversion algorithms for ERT data. The implementation is publicly available and developed within the framework of PFLOTRAN, an open-source, state-of-the-art massively parallel subsurface flow and transport simulation code. The forward modeling is based on a finite-volume discretization of the governing differential equations, and the inversion uses a Gauss–Newton optimization scheme. To evaluate the accuracy of the forward modeling, two examples are first presented by considering layered (1D) and 3D earth conductivity models. The computed numerical results show good agreement with the analytical solutions for the layered earth model and results from a well-established code for the 3D model. Inversion of ERT data, simulated for a 3D model, is then performed to demonstrate the inversion capability by recovering the conductivity of the model. To demonstrate the parallel performance of PFLOTRAN's ERT process model and inversion capabilities, large-scale scalability tests are performed by using up to 131 072 processes on a leadership class supercomputer. These tests are performed for the two most computationally intensive steps of the ERT inversion: forward modeling and Jacobian computation. For the forward modeling, we consider models with up to 122 ×106 degrees of freedom (DOFs) in the resulting system of linear equations and demonstrate that the code exhibits almost linear scalability on up to 10 000 DOFs per process. On the other hand, the code shows superlinear scalability for the Jacobian computation, mainly because all computations are fairly evenly distributed over each process with no parallel communication.

58 GEOSCIENCES↗

DOE OSTI2022

A deep learning-based direct forecasting of CO 2 plume migration

Accurate and timely forecasts of CO 2 plume evolution in geological reservoirs are crucial for CO 2 migration detection, leakage risk assessment, and operation decision support. Conventional forecasting usually adopts a two-step strategy, first calibrating reservoir model parameters against observations using iterative inverse modeling (or history matching) and then applying the calibrated model for predictions. This method impedes real-time forecasts due to the heavy computational demand in inverse modeling and may suffer from poor prediction accuracy because of the limited observation data. In this work, we propose a deep learning-based latent space mapping framework to forecast CO 2 plume migration directly by avoiding the inverse modeling. We first use the convolutional autoencoder to map the high-dimensional complex plume extents onto low-dimensional latent space. Next, we use neural networks to learn the relationship between the observation variables and the prediction latent variables. And then for given observation data, we infer the prediction values directly. This one-step direct forecasting is computationally efficient which requires a few number of parallelizable reservoir simulations and it can provide accurate predictions with limited observations by learning the observation-prediction relationship in the reduced dimension. Therefore, our proposed method enables an in-time forecast of dynamic CO 2 plume distributions. In this work, we demonstrate the effectiveness and accuracy of our method in predicting the CO 2 plume migration using four metrics such as plume area, centroid movement distance, and plume spreading in the primary and secondary directions. And the spatio-temporal evolution patterns of plume migration under diverse geological complexities are also accurately quantified.

15 GEOTHERMAL ENERGY↗

NASA NTRS2019

Prior biosphere model impact on global terrestrial CO2 fluxes estimated from OCO-2 retrievals

This study assesses the impact of different state of the art global biospheric CO2 flux models, when applied as prior information, on inverse model “top-down” estimates of terrestrial CO2 fluxes obtained when assimilating Orbiting Carbon Observatory 2 (OCO-2) observations. This is done with a series of observing system simulation experiments (OSSEs) using synthetic CO2 column-average dry air mole fraction (XCO2) retrievals sampled at the OCO-2 satellite spatiotemporal frequency. The OSSEs utilized a 4-D variational (4D-Var) assimilation system with the GEOS-Chem global chemical transport model (CTM) to estimate CO2 net ecosystem exchange (NEE) fluxes using synthetic OCO-2 observations. The impact of biosphere models in inverse model estimates of NEE is quantified by conducting OSSEs using the NASA-CASA, CASA-GFED, SiB-4, and LPJ models as prior estimates and using NEE from the multi-model ensemble mean of the Multiscale Synthesis and Terrestrial Model Intercomparison Project as the “truth”. Results show that the assimilation of simulated XCO2 retrievals at OCO-2 observing modes over land results in posterior NEE estimates which generally reproduce “true” NEE globally and over terrestrial TransCom-3 regions that are well-sampled. However, we find larger spread among posterior NEE estimates, when using different prior NEE fluxes, in regions and seasons that have limited OCO-2 observational coverage and a large range in “bottom-up” NEE fluxes. Seasonally averaged posterior NEE estimates had standard deviations (SD) of ∼10 % to ∼50 % of the multi-model-mean NEE for different TransCom-3 land regions with significant NEE fluxes (regions/seasons with a NEE flux ≥0.5 PgC/yr). On a global average, the seasonally averaged residual impact of the prior model NEE assumption on the posterior NEE spread is ∼10 %–20 % of the posterior NEE mean. Additional OCO-2 OSSE simulations demonstrate that posterior NEE estimates are also sensitive to the assumed prior NEE flux uncertainty statistics, with spread in posterior NEE estimates similar to those when using variable prior model NEE fluxes. In fact, the sensitivity of posterior NEE estimates to prior error statistics was larger than prior flux values in some regions/times in the tropics and Southern Hemisphere where sufficient OCO-2 data were available and large differences between the prior and truth were evident. Overall, even with the availability of spatiotemporally dense OCO-2 data, noticeable residual differences (up to ∼20 %–30 % globally and 50 % regionally) in posterior NEE flux estimates remain that were caused by the choice of prior model flux values and the specification of prior flux uncertainties.

Sajeev Philip↗

DOE OSTI2024

Permafrost Carbon: Progress on Understanding Stocks and Fluxes Across Northern Terrestrial Ecosystems

Significant progress in permafrost carbon science made over the past decades include the identification of vast permafrost carbon stocks, the development of new pan-Arctic permafrost maps, an increase in terrestrial measurement sites for CO 2 and methane fluxes, and important factors affecting carbon cycling, including vegetation changes, periods of soil freezing and thawing, wildfire, and other disturbance events. Process-based modeling studies now include key elements of permafrost carbon cycling and advances in statistical modeling and inverse modeling enhance understanding of permafrost region C budgets. By combining existing data syntheses and model outputs, the permafrost region is likely a wetland methane source and small terrestrial ecosystem CO 2 sink with lower net CO 2 uptake toward higher latitudes, excluding wildfire emissions. For 2002–2014, the strongest CO2 sink was located in western Canada (median: -52 g C m -2 y -1 ) and smallest sinks in Alaska, Canadian tundra, and Siberian tundra (medians: -5 to -9 g C m -2 y -1 ). Eurasian regions had the largest median wetland methane fluxes (16–18 g CH4 m -2 y -1 ). Quantifying the regional scale carbon balance remains challenging because of high spatial and temporal variability and relatively low density of observations. More accurate permafrost region carbon fluxes require: (a) the development of better maps characterizing wetlands and dynamics of vegetation and disturbances, including abrupt permafrost thaw; (b) the establishment of new year-round CO 2 and methane flux sites in underrepresented areas; and (c) improved models that better represent important permafrost carbon cycle dynamics, including non-growing season emissions and disturbance effects.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2023

Using Orbiting Carbon Observatory-2 (OCO-2) column CO2 retrievals to rapidly detect and estimate biospheric surface carbon flux anomalies

The global carbon cycle is experiencing continued perturbations via increases in atmospheric carbon concentrations, which are partly reduced by terrestrial biosphere and ocean carbon uptake. Greenhouse gas satellites have been shown to be useful in retrieving atmospheric carbon concentrations and observing surface and atmospheric CO2 seasonal-to-interannual variations. However, limited attention has been placed on using satellite column CO2 retrievals to evaluate surface CO2 fluxes from the terrestrial biosphere without advanced inversion models at low latency. Such applications could be useful to monitor, in near real time, biosphere carbon fluxes during climatic anomalies like drought, heatwaves, and floods, before more complex terrestrial biosphere model outputs and/or advanced inversion modelling estimates become available. Here, we explore the ability of Orbiting Carbon Observatory-2 (OCO-2) column-averaged dry air CO2 (XCO2) retrievals to directly detect and estimate terrestrial biosphere CO2 flux anomalies using a simple mass-balance approach. An initial global analysis of surface–atmospheric CO2 coupling and transport conditions reveals that the western US, among a handful of other regions, is a feasible candidate for using XCO2 for detecting terrestrial biosphere CO2 flux anomalies. Using the CarbonTracker model reanalysis as a test bed, we first demonstrate that a well-established mass-balance approach can estimate monthly surface CO2 flux anomalies from XCO2 enhancements in the western United States. The method is optimal when the study domain is spatially extensive enough to account for atmospheric mixing and has favorable advection conditions with contributions primarily from one background region. We find that errors in individual soundings reduce the ability of OCO-2 XCO2 to estimate more frequent, smaller surface CO2 flux anomalies. However, we find that OCO-2 XCO2 can often detect and estimate large surface flux anomalies that leave an imprint on the atmospheric CO2 concentration anomalies beyond the retrieval error/uncertainty associated with the observations. OCO-2 can thus be useful for low-latency monitoring of the monthly timing and magnitude of extreme regional terrestrial biosphere carbon anomalies.

Andrew F. Feldman↗

DOE OSTI2022

Multimodal Freight Energy Model for Emerging Freight Technology Analysis

While freight movement represents a small portion of the volume in the transportation sector, they are a critical contributor for energy consumed in that sector. To reduce logistics costs, energy consumption and negative environmental impact, emerging technologies, such as digitalization of logistics and alternative powertrain, have been developed and extended applications in freight. These trends are expected to grow and provide opportunities for greater efficiencies in freight movement and corresponding energy use. Still, the complexity of freight systems present challenges in evaluating these benefits, especially in the multimodal inter-city freight. Addressing the research need, this paper develops a multimodal freight energy modeling framework for the analysis of emerging freight technology scenarios. The framework is a bi-level optimization problem: network cost minimization problem (lower-level) and energy minimization problem (upper-level). The lower-level problem is a mode-path assignment problem in multimodal inter-city freight networks, where commodity-specific congestion effects on trans-shipment links are considered. For this model, an inverse modeling approach is applied to infer parameters of the lower-level model. The upper-level problem is designed to search for an optimal scenario that has the lowest energy consumption among different levels of technology applications. The proposed model is empirically tested to analyze truck load-pooling and multimodal load-pooling scenarios through stand-alone and mixed applications using freight shipments originating from or destined to the Chicago region. This framework can be used to explore the impact of emerging freight technologies on mode-path freight flow and energy consumption in the national multimodal freight network.

ADVANCED PROPULSION SYSTEMS↗

NASA NTRS1984

Inversion of vegetation canopy reflectance models for estimating agronomic variables. III - Estimation using only canopy reflectance data as illustrated by the suits model. IV - Total inversion of the SAIL model

The possibility of estimating agronomic and spectral parameters for a vegetation canopy from the canopy reflectance (CR) data in the infrared region is investigated for a set of solar/view directions. It is shown that such an estimation is possible, in principle, for the Suits (1972) model for a homogeneous canopy. The technique is then applied to a more complex model, the SAIL (scattering by arbitrarily inclined leaves) model which explicitly includes the leaf angle distribution in the computation of the canopy reflectance. It is concluded that, given the expected accuracy of CR measurements and the accuracy of the SAIL model in representing CR in the infrared region, the agronomic parameters, leaf area index, and leaf angle distribution can be estimated fairly accurately using ancillary data on spectral parameters.

Goel, N. S.↗

NASA NTRS2021

Improved representation of the global dust cycle using observational constraints on dust properties and abundance

Even though desert dust is the most abundant aerosol by mass in Earth's atmosphere, atmospheric models struggle to accurately represent its spatial and temporal distribution. These model errors are partially caused by fundamental difficulties in simulating dust emission in coarse-resolution models and in accurately representing dust microphysical properties. Here we mitigate these problems by developing a new methodology that yields an improved representation of the global dust cycle. We present an analytical framework that uses inverse modeling to integrate an ensemble of global model simulations with observational constraints on the dust size distribution, extinction efficiency, and regional dust aerosol optical depth. We then compare the inverse model results against independent measurements of dust surface concentration and deposition flux and find that errors are reduced by approximately a factor of 2 relative to current model simulations of the Northern Hemisphere dust cycle. The inverse model results show smaller improvements in the less dusty Southern Hemisphere, most likely because both the model simulations and the observational constraints used in the inverse model are less accurate. On a global basis, we find that the emission flux of dust with a geometric diameter up to 20 µm (PM20) is approximately 5000 Tg yr−1, which is greater than most models account for. This larger PM20 dust flux is needed to match observational constraints showing a large atmospheric loading of coarse dust. We obtain gridded datasets of dust emission, vertically integrated loading, dust aerosol optical depth, (surface) concentration, and wet and dry deposition fluxes that are resolved by season and particle size. As our results indicate that this dataset is more accurate than current model simulations and the MERRA-2 dust reanalysis product, it can be used to improve quantifications of dust impacts on the Earth system.

Dust cycle↗

NASA NTRS2021

Improved Representation of the Global Dust Cycle Using Observational Constraints on Dust Properties and Abundance

Even though desert dust is the most abundant aerosol by mass in Earth’s atmosphere, atmospheric models 25 struggle to accurately represent its spatial and temporal distribution. These model errors are partially caused byfundamental difficulties in simulating dust emission in coarse-resolution models and in accurately representing dust microphysical properties. Here we mitigate these problems by developing a new methodology that yields an improved representation of the global dust cycle. We present an analytical framework that uses inverse modeling to integrate an ensemble of global model simulations with observational constraints on the dust size distribution, 30 extinction efficiency, and regional dust aerosol optical depth. We then compare the inverse model results against independent measurements of dust surface concentration and deposition flux and find that errors are reduced by approximately a factor of two relative to current model simulations of the Northern Hemisphere dust cycle. The inverse model results show smaller improvements in the less dusty Southern Hemisphere, most likely because both the model simulations and the observational constraints used in the inverse model are less accurate. On a global 35 basis, we find that the emission flux of dust with geometric diameter up to 20 μm (PM20) is approximately 5,000 Tg/year, which is greater than most models account for. This larger PM20 dust flux is needed to match observational constraints showing a large atmospheric loading of coarse dust. We obtain gridded data sets of dust emission, vertically integrated loading, dust aerosol optical depth, (surface) concentration, and wet and dry deposition fluxes that are resolved by season and particle size. As our results indicate that this data set is more accurate than current 40 model simulations and the MERRA-2 dust reanalysis product, it can be used to improve quantifications of dust impacts on the Earth system.

aerosols, dust↗

NASA NTRS1991

Computer modeling of inversion layer MOS solar cells and arrays

A two dimensional numerical model of the inversion layer metal insulator semiconductor (IL/MIS) solar cell is proposed by using the finite element method. The two-dimensional current flow in the device is taken into account in this model. The electrostatic potential distribution, the electron concentration distribution, and the hole concentration distribution for different terminal voltages are simulated. The results of simple calculation are presented. The existing problems for this model are addressed. Future work is proposed. The MIS structures are studied and some of the results are reported.

Ho, Fat Duen↗

DOE OSTI2025

Diffusion Model-Guided Inverse Design of Bimetallic Catalysts for Ammonia Decomposition

In the past decade, artificial intelligence and deep learning have played increasingly prominent roles in materials design and discovery. Among these, generative AI models, known for their ability to create unique and complex structures, have emerged as state-of-the-art tools for materials screening due to their high efficiency and low computational cost. In catalysis, one of the major challenges is identifying promising material candidates within an immense chemical space. This challenge can be addressed using generative approaches, such as diffusion-based inverse design models. In this study, we present a machine learning-guided workflow that employed a diffusion model for the inverse design of bimetallic alloy catalysts for low-carbon ammonia decomposition, a key reaction for ammonia emission control and sustainable hydrogen production. Catalyst candidates were evaluated using nitrogen adsorption energy as the key descriptor, inspired by multiscale modeling. The proposed workflow identified low-cost, environmentally friendly catalysts with excellent catalytic performance, which have been validated theoretically and experimentally. Our framework decoupled the generative and property-prediction components, enhancing both flexibility and accuracy in the catalytic material design process.

Adsorption↗

NASA NTRS1984

Continental and oceanic crustal magnetization modelling

Inversion of magnetic data from the MAGSAT satellite, to arrive at intensities of magnetization of the Earth's crust, was performed by two different methods. The first method uses a spherical harmonic model of the magnetic field. The coefficients believed to represent sources in the Earth's crust can then be inverted to arrive at vertical dipole moments per unit area at the Earth's surface. The spherical harmonic models contain coefficients of degrees of harmonics up to 23. The dipole moment per unit area for a surface element can then be determined by summing the contribution for each individual degree of harmonic. The magnetic moments were calculated for continental and oceanic areas separately as well as over certain latitudinal segments. Of primary concern was to determine whether there are any differences between continental and oceanic areas. The second analysis with magnetization intensities was made using narrower ranges of degrees of harmonics, assuming that higher degrees are present in the core field signal.

Harrison, C. G. A.↗

DOE OSTI2021

Utilizing Earth Observations of Soil Freeze/Thaw Data and Atmospheric Concentrations to Estimate Cold Season Methane Emissions in the Northern High Latitudes

The northern wetland methane emission estimates have large uncertainties. Inversion models are a qualified method to estimate the methane fluxes and emissions in northern latitudes but when atmospheric observations are sparse, the models are only as good as their a priori estimates. Thus, improving a priori estimates is a competent way to reduce uncertainties and enhance emission estimates in the sparsely sampled regions. Here, we use a novel way to integrate remote sensing soil freeze/thaw (F/T) status from SMOS satellite to better capture the seasonality of methane emissions in the northern high latitude. The SMOS F/T data provide daily information of soil freezing state in the northern latitudes, and in this study, the data is used to define the cold season in the high latitudes and, thus, improve our knowledge of the seasonal cycle of biospheric methane fluxes. The SMOS F/T data is implemented to LPX-Bern DYPTOP model estimates and the modified fluxes are used as a biospheric a priori in the inversion model CarbonTracker Europe-CH 4 . The implementation of the SMOS F/T soil state is shown to be beneficial in improving the inversion model’s cold season biospheric flux estimates. Our results show that cold season biospheric CH4 emissions in northern high latitudes are approximately 0.60 Tg lower than previously estimated, which corresponds to 17% reduction in the cold season biospheric emissions. This reduction is partly compensated by increased anthropogenic emissions in the same area (0.23 Tg), and the results also indicates that the anthropogenic emissions could have even larger contribution in cold season than estimated here.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Modeling of transient conduction in building envelope assemblies: A review

As buildings age, retrofits are becoming an increasingly important topic for the ever-growing and aging existing building stock. To compare designs or evaluate in-service building envelopes, thermal modeling is utilized to evaluate the thermal performance of envelope assemblies; however, it can be difficult to model the thermal performance of as-built assemblies due to degradation or missing documentation. To address this issue, inverse modeling can be applied to infer the properties of as-built envelope assemblies. This paper presents a review of the literature and published methods to model and infer the transient conductive performance of building envelopes. This review serves as a survey of existing transient conduction algorithms to evaluate performance, computational speed, and relevance for inverse modeling applications. In addition to the literature review, this work also evaluates the computational performance of the most prevalent transient conduction algorithms against the ASHRAE 1052RP toolkit to assess inverse modeling potential. This methodology serves as a foundation for future research to characterize the transient thermal performance of as-built building envelope assemblies.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2023

Arctic Cloud‐Base Ice Precipitation Properties Retrieved Using Bayesian Inference

Abstract Cloud‐climate feedbacks are still the greatest source of uncertainty in current climate projections. Arctic clouds, which are predominantly stratiform and supercooled, often long‐lived, and nearly continuously precipitate ice particles, contribute roughly 10% of the uncertainty attributed to the global cloud feedback. This Arctic cloud uncertainty is driven by incomplete observational and theoretical knowledge required to estimate and explain the state and active processes occurring in those clouds. A focus on ice precipitation properties at Arctic cloud base rather than the surface deconfounds the product of cloud condensate sink processes from the influence of the atmospheric thermodynamic state below cloud base, rendering cloud‐base properties a more appealing target for inference and evaluation of model simulations. Here I describe an inverse model for the estimation of cloud base ice precipitation properties over Utqiagvik, North Slope of Alaska, using the synthesis of ground‐based radar and lidar measurements. By leveraging a Markov Chain Monte Carlo algorithm as the core of the inverse model, a wide range of particle size distributions are sampled, and different combinations of ice habit models are examined, both of which are typically fixed in other retrieval methods. Results show intriguing links between different cloud base thermodynamic and ice precipitation properties. Apparent ice number concentration enhancements at temperatures of −5 and −15°C suggest possible secondary ice production (SIP). The analysis alludes to an overestimation of SIP occurrence and intensity, especially in studies relying only on radar or lidar measurements. Finally, reflectivity‐dependent ice precipitation rate and ice water content parameterizations are presented.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Estimating Watershed Subsurface Permeability From Stream Discharge Data Using Deep Neural Networks

Subsurface permeability is a key parameter in watershed models that controls the contribution from the subsurface flow to stream flows. Since the permeability is difficult and expensive to measure directly at the spatial extent and resolution required by fully distributed watershed models, estimation through inverse modeling has had a long history in subsurface hydrology. The wide availability of stream surface flow data, compared to groundwater monitoring data, provides a new data source to infer soil and geologic properties using integrated surface and subsurface hydrologic models. As most of the existing methods have shown difficulty in dealing with highly nonlinear inverse problems, we explore the use of deep neural networks for inversion owing to their successes in mapping complex, highly nonlinear relationships. We train various deep neural network (DNN) models with different architectures to predict subsurface permeability from stream discharge hydrograph at the watershed outlet. The training data are obtained from ensemble simulations of hydrographs corresponding to an permeability ensemble using a fully-distributed, integrated surface-subsurface hydrologic model. The trained model is then applied to estimate the permeability of the real watershed using its observed hydrograph at the outlet. Our study demonstrates that the permeabilities of the soil and geologic facies that make significant contributions to the outlet discharge can be more accurately estimated from the discharge data. Their estimations are also more robust with observation errors. Compared to the traditional ensemble smoother method, DNNs show stronger performance in capturing the nonlinear relationship between permeability and stream hydrograph to accurately estimate permeability. Our study sheds new light on the value of the emerging deep learning methods in assisting integrated watershed modeling by improving parameter estimation, which will eventually reduce the uncertainty in predictive watershed models.

54 ENVIRONMENTAL SCIENCES↗