Search NASA⌕ Search

SEARCH · Search NASA

Results for “Variable selection”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

IM3/HyperFACETS Thermodynamic Global Warming (TGW) Simulation Datasets

Publication For a thorough description of the methods, see the peer-reviewed paper: Jones, A.D., Rastogi, D., Vahmani, P. et al. Continental United States climate projections based on thermodynamic modification of historical weather. Sci Data 10, 664 (2023). https://doi.org/10.1038/s41597-023-02485-5 Overview The IM3 / HyperFACETS climate simulations provide 40-year historical (1980-2019) as well as four 80-year future simulations (2020-2099) over the U.S. The future simulations are split into near (2020-2059) and far future (2060-2099) segments. The future scenarios span a range of plausible changes in future climate (both Global Circulation Model (GCM) and Representative Concentration Pathways/Shared Socioeconomic Pathway (RCP/SSP) dimensions). The simulations provide climate variables with high spatiotemporal resolution (25 hourly variables and 207 3-hourly variables at 12 km2). The datasets are generated using dynamical downscaling with the WRF (Weather Research and Forecasting) model (version 4.2.1) and therefore preserve physical consistency across variables. WRF is a state-of-the-art, fully compressible, non-hydrostatic, mesoscale numerical weather prediction model. WRF is coupled with an urban canopy model (UCM), which resolves urban surfaces. The future scenarios were developed using a thermodynamic global warming approach where past events are replayed under a range of future warming conditions. These scenarios therefore provide a perspective on potential increases in extreme event intensity, geographic scope, and duration, with previously non-extreme conditions potentially crossing new thresholds to be considered extreme by today's standards. This approach is not intended to estimate future changes in extreme event frequency that might result from changes in large-scale atmospheric dynamics. This dataset has NOT been bias corrected. A bias corrected version of selected variables is under development and will be released here when available. Scenarios Files Data for each scenario is provided in weekly NetCDF files. 25 variables are available at hourly resolution, and 207 variables are available at three-hourly resolution. Spatial resolution is 12km and spans the conterminous United States (CONUS), including some areas of Canada and Mexico, resulting in a grid of 424 by 299 cells. The spatial projection is a Lambert Conformal Conic with the following proj-string: "+proj=lcc +lat_0=40.0000076293945 +lon_0=-97 +lat_1=30 +lat_2=45 +x_0=0 +y_0=0 +R=6370000 +units=m +no_defs". The available scenarios and simulation periods are listed below: historical | 1980 - 2019 rcp45cooler | 2020 - 2059 rcp45cooler | 2060 - 2099 rcp45hotter | 2020 - 2059 rcp45hotter | 2060 - 2099 rcp85cooler | 2020 - 2059 rcp85cooler | 2060 - 2099 rcp85hotter | 2020 - 2059 rcp85hotter | 2060 - 2099 * The first year (1979, 2019, and 2059) of data within each scenario represents a model warmup period and should not be used. These are located in the `spinup_files` directory. Historical year 2020 is considered an extra year of data beyond the simulation period and can be found in the `additional_files` directory. For information on specific variables and a more in-depth discussion of methodology, please refer to the data landing page at https://tgw-data.msdlive.org. Delta Warming Files The global and CONUS warming deltas for each scenario are provided in degrees Celsius annually and monthly. Restart Files Yearly restart files are provided for each scenario which can be used to restart the WRF model at a particular point in time. Spinup Files The first year of data within each simulation period represents a model warmup period and should not be used. The files are provided here for the sake of reproducibility. Additional Files Additional years of data are provided as an extension of the historic simulation.

Jones, Andrew D.↗
DOE OSTI2023

Predicting Postoperative Injury and Military Discharge Status After Knee Surgery in the US Army

Background: Researchers have assessed postoperative injury or disability predictors in the military setting but typically focused on 1 type of surgical procedure at a time, used relatively small sample sizes, or investigated mixed cohorts with civilian populations. Purpose: To identify the relationship between baseline variables and injury incidence or military discharge status in US Army soldiers after knee surgery. Study Design: Case-control study; Level of evidence, 3. Methods: Data were obtained from a repository containing personnel, performance, and medical records for all active-duty US Army soldiers. Multivariate logistic regressions were used to estimate the effects of numerous variables on postoperative injury or on medical discharge. Variable selection and model validation were conducted using the k-fold method. Results: A total of 7567 soldiers underwent knee surgery between 2017 and 2019. Meniscal procedures were the most common type of surgery (39%), and approximately 71% of the cohort had a postoperative injury. Significant predictors for sustaining a postoperative injury included having a previous nonknee injury (odds ratio [OR], 1.5), female sex (OR, 1.3), and Black race (OR, 1.2). Within 4 years after surgery, 17% of soldiers were discharged from the military because of knee-related disability. Significant predictors for discharge from duty included enlisted rank (OR, 2.3), recent fitness test failure (OR, 1.9), number of previous knee surgeries (OR, 1.7), and having a previous nonknee injury (OR, 1.6). Conclusion: After knee surgery, nearly three-fourths of the soldiers in this cohort sustained a postoperative injury and almost one-fifth of soldiers were medically discharged from the military within 4 years. This study identified variables that indicate statistically increased risk for these postoperative outcomes and highlighted potentially modifiable factors.

Orthopedics↗
DOE OSTI2020

An updated estimate of the Mu2e experiment sensitivity

The Mu2e experiment at Fermilab will search for the conversion of a negative muon into an electron inside the field of a nucleus. This process does not conserve charged-lepton flavour and is heavily suppressed in the Standard Model (SM), with a branching ratio < 10-50. Any evidence of it would be an undeniable evidence of new physics beyond the SM. The project sets out to achieve a single event sensitivity of $\sim 3 × 10^{-17}$ on the ratio between the probability for a conversion of a negative muon into an electron and the one for a muon capture by the nucleus. Such a sensitivity would represent a 4 orders of magnitude improvement on the previous upper limit for the process, making possible to test predictions of different extensions of the SM. Mu2e uses three superconducting solenoids to produce and measure the muon conversions. In the first solenoid, the Production Solenoid, pions and kaons are produced, together with other secondary products, by 8 GeV kinetic energy proton interactions in a tungsten target. A gradient magnetic field is specifically designed to direct low momentum particles into the Transport Solenoid, an S-shaped magnet that filters out particles with unwanted charge and momentum. Muons, produced by pion and kaon decays, finally reach the aluminum Stopping Target in the Detector Solenoid, where they eventually stop and convert. The result of the conversion process is a monochromatic electron of ~105 MeV/c momentum. The Detector Solenoid also hosts the two main detectors: a straw tube tracker and two CsI calorimeter disks, both providing measurement of event kinematics. A germanium detector and a LaBr crystal are located downstream of the Detector Solenoid to measure the X and gamma rays produced by the muon captures in the Stopping Target. A veto system of scintillators covering the Detector Solenoid and half of the Transport Solenoid is used to identify and reject cosmic rays interactions. With respect to the initial project, the Mu2e running plan has evolved to a staged configuration with 2 years at reduced intensity before the 2025 accelerator shutdown for the neutrino beam upgrade and 2 or 3 years at full intensity after that. This, together with geometry changes and a better knowledge of detector performances obtained by the first slice tests, has required a full revision of the signal over background selection that is the subject of this thesis. In order to achieve a new estimate of Mu2e sensitivity, the simulation of the data corresponding to the first 2 years of data acquisition has been performed. This includes both conversion electrons (CE) and the main sources of background: cosmics, decay-in orbits (DIO), radiative pion captures (RPC) and antiprotons. The characteristics of the signal and of the main backgrounds have been studied to define the best selection variables for CE. A special effort has been devoted to the evaluation of the antiproton background. The lack of experimental data for antiproton production cross section makes the systematic uncertainty on this background significant. A new parameterization of the cross section has been developed to fit the existing data and to provide a more reliable estimate of the systematic uncertainty by comparing the results of the old and the new model. A special effort has been devoted to the optimization of the antiprotons Monte Carlo generator. This study has also revealed that the dominant component within this background is represented by antiprotons produced in the opposite direction with respect to the Transport Solenoid entrance and then redirected to it by back-scattering processes in the Production Target. This ultimately highlights the sensitivity of the background estimate to the G4 handling of antiproton interactions in the 1 to 3 GeV/c momentum range. Finally, the final experimental sensitivity has been studied. The momentum and time selection have been optimized to obtain the 5-sigma discovery reach or, in case of no signal, the upper limit on conversion probability. The results confirm that, in the firs t two years of data taking, Mu2e will be able to improve the current experimental sensitivity for muon-to-electron conversion in an atomic field by more than 3 order of magnitudes.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Optimization of direct air capture processes using reactive transport models of adsorption-desorption cycles

In this study, we develop and implement a reactive transport model in COMSOL Multiphysics® to address the challenges of direct air carbon capture. The model is validated against experimental data and used to simulate the cyclic steady state of the adsorption-desorption process. The optimization of this model is achieved through advanced trust-region methods integrated with Gaussian Processes. Key decision variables, including adsorption and desorption times, desorption temperature and pressure, input velocity, bed porosity, column length, and radius were optimized to minimize the capture cost. After optimization, a sensitivity analysis revealed the complex interplay between the decision variables and their effect on the specific energy and cost of removing the CO 2 . We optimized the capture cost while taking into account the trade-off between energy consumption and productivity. The resulting minimum capture cost was determined to be 265.2 $/t-CO 2 , which aligns with expected values reported in the literature. Numerical results suggest the effectiveness of the optimization strategies applied, and underscore the importance of simultaneous decision variable selection in improving the performance in direct air capture processes. We also extend the modeling approach to a 2D axisymmetric model to better visualize CO₂ uptake and temperature profiles, revealing significant radial gradients during the regeneration step. As a main drawback, this enhanced model comes with a computational cost approximately 40 times higher than that of the 1D model.

Adsorption-desorption process↗
DOE OSTI2022

Evaluating the performance of random forest and iterative random forest based methods when applied to gene expression data

Gene-to-gene networks, such as Gene Regulatory Networks (GRN) and Predictive Expression Networks (PEN) capture relationships between genes and are beneficial for use in downstream biological analyses. There exists multiple network inference tools to produce these gene-to-gene networks from matrices of gene expression data. Random Forest-Leave One Out Prediction (RF-LOOP) is a method that has been shown to be efficient at producing these gene-to-gene networks, frequently known as GEne Network Inference with Ensemble of trees (GENIE3). Random Forest can be replaced in this process by iterative Random Forest (iRF), which performs variable selection and boosting. Here we validate that iterative Random Forest-Leave One Out Prediction (iRF-LOOP) produces higher quality networks than GENIE3 (RF-LOOP). We use both synthetic and empirical networks from the Dialogue for Reverse Engineering Assessment and Methods (DREAM) Challenges by Sage Bionetworks, as well as two additional empirical networks created from Arabidopsis thaliana and Populus trichocarpa expression data.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Pavement condition and climatic data in southeast Texas: A dataset for evaluating flood impacts on pavement performance

Effective pavement maintenance is essential for economic stability, optimal network performance, and roadway safety. Achieving this requires thorough evaluation of pavement conditions, including structural integrity, surface roughness, and distress characteristics. Pavement performance indicators play a critical role in influencing vehicle safety and ride quality. Recent advances have emphasized the use of data-driven modeling to anticipate pavement behavior, with the goal of optimizing resource allocation and refining Maintenance and Rehabilitation (M&R) strategies through accurate condition assessment. A foundational requirement for these modeling efforts is the availability of standardized, high-quality datasets that can support robust and reproducible infrastructure analysis. This data article presents a comprehensive dataset assembled to facilitate pavement performance prediction, with a geographic focus on Southeast Texas, particularly the flood-vulnerable area of Beaumont. The dataset encompasses pavement and traffic attributes, meteorological records, flood simulation outputs, ground deformation measurements, and topographic indices, enabling detailed examination of both load-associated and non-load-associated degradation mechanisms. Data preprocessing was performed using ArcGIS Pro, Microsoft Excel, and Python to ensure consistency and usability in data-driven modeling applications, including machine learning workflows. Key contributions of this dataset include its utility in analyzing the climatic and environmental factors affecting pavement conditions, identifying critical predictive features, and enabling in-depth correlation analysis across diverse variables. By filling existing gaps in input variable selection resources, this dataset supports the development of predictive tools for estimating future maintenance demand and enhancing the resilience of pavement networks in flood-impacted areas. The resource highlights the importance of standardized datasets for advancing pavement management practices and provides a robust foundation for ongoing infrastructure performance modeling.

42 ENGINEERING↗
DOE OSTI2024

Identifying Disadvantaged Communities in the United States: An Energy-Oriented Mapping Tool that Aggregates Environmental and Socioeconomic Burdens

This paper defines a policy-relevant nationwide composite index to identify communities disproportionately impacted by environmental, energy, and climate injustices in the United States. We review existing vulnerability indicators and indices to assess the tradeoffs of different design parameters, including variable selection, geographic unit, dimensionality reduction, weighting, and aggregation methods. From this methodological review, we create the first nationwide, census tract-level index of cumulative burden that includes energy-relevant indicators alongside climate, social, environmental, and economic indicators, and is flexible to the inclusion of additional data sources. We provide a summary of the sources of inputs used to develop a definition for "disadvantaged communities" that can be used to prioritize energy investments. We discuss use-cases for this index including the implementation of the Justice40 Initiative, which calls for 40% of certain federal clean energy benefits to flow to disadvantaged communities in the United States. We use our results to examine historic allocations of federal energy investments and show that communities that we identify as disadvantaged received about 37% fewer funds per capita than non-disadvantaged communities.

cumulative burden↗
DOE OSTI2026

Revealing the complex chemistry of grain boundaries in K-doped BaFe 2 As 2 with atom probe tomography

Iron-based superconductors have attractive properties for high-field applications, but there is a lack of understanding of the effect of grain boundary chemistry on the in-field performance. The near atomic-scale resolution, ppm sensitivity and 3D analysis offered by atom probe tomography make it a powerful tool to investigate the nanoscale structure and chemistry of these defects in fine-grained K-doped BaFe 2 As 2 samples. A computational method to systematically extract and compare the Gibbsian interfacial excess of chemical species across grain boundaries has been explored in this work. The robustness of the method has been tested by evaluating the effects of selected variables on simulated APT datasets. The accuracy and precision of the calculated Gibbsian interfacial excess were found to be stable over a range of analysis conditions: varying grain boundary widths and detection efficiencies, spatial precisions below 1.5 nm, and bin widths between 1.2 and 1.6 nm. For the K-doped BaFe 2 As 2 samples studied, segregation of As, Ba, K and impurities of O, Na, and Sb were found at grain boundaries. The Gibbsian excess values were found to vary widely between different boundaries, showing the complexity of the grain boundary chemistry in this material. Possible links between the observed critical current density (Jc) of these samples and their nano- and micro-structure have also been investigated and discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Geology and elevation shape bacterial assembly in Antarctic endolithic communities

Ice free areas of continental Antarctica are among the coldest and driest environments on Earth, and yet, they support surprisingly diverse and highly adapted microbial communities. Endolithic growth is one of the key adaptations to such extreme environments and often represents the dominant life-form. Despite growing scientific interest, little is known of the mechanisms that influence the assembly of endolithic microbiomes across these harsh environments. Here, we used metagenomics to examine the diversity and assembly of endolithic bacterial communities across Antarctica within different rock types and over a large elevation range. While granite supported richer and more heterogeneous communities than sandstone, elevation had no apparent effect on taxonomic richness, regardless of rock type. Conversely, elevation was clearly associated with turnover in community composition, with the deterministic process of variable selection driving microbial assembly along the elevation gradient. The turnover associated with elevation was modulated by geology, whereby for a given elevation difference, turnover was consistently larger between communities inhabiting different rock types. Overall, selection imposed by elevation and geology appeared stronger than turnover related to other spatially-structured environmental drivers. Our findings indicate that at the cold-arid limit of life on Earth, geology and elevation are key determinants of endolithic bacterial heterogeneity. This also suggests that warming temperatures may threaten the persistence of such extreme-adapted organisms.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Microbial and Environmental Processes Shape the Link between Organic Matter Functional Traits and Composition

Dissolved organic matter (DOM) is a large and complex mixture of molecules that fuels biogeochemical reaction in virtually all ecosystems on Earth. However, the relative importance of deterministic and stochastic processes in structuring DOM composition remains poorly characterized. Here we develop a framework for partitioning molecular composition based on key molecular traits, including lability vs. recalcitrance and activity vs. inactivity. Within this framework, we examine the ecological processes governing the assembly of DOM fractions by deploying aquatic microcosms on mountainsides that span gradients of temperature and nutrient loading in subtropical and subarctic ecosystems. Across study regions, deterministic and stochastic processes primarily structure active and inactive fractions, respectively. However, recalcitrant molecules are more deterministically assembled than labile molecules in the inactive fraction. Deterministic processes leading to variable selection generally exhibit more variation across the energy supply gradient for inactive fractions, and their importance increases with energy supply for recalcitrant molecules in both active and inactive fractions. Together, our results indicate that active and inactive fractions of DOM assemblages are structured by contrasting ecological processes, and their recalcitrant components are sensitive to global change. In conclusion, our framework opens new avenues to understand the assembly and turnover of DOM in a changing world, which can be used to predict carbon cycling at local to global scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Application of artificial intelligence methods in the international roughness index prediction of rigid and composite pavements: a systematic review

The International Roughness Index (IRI) is a widely adopted metric for quantifying pavement roughness, directly influencing vehicle safety, ride comfort, and overall roadway performance. In recent years, the use of Machine Learning (ML) models for IRI prediction has gained momentum, with the goal of improving the allocation of maintenance and rehabilitation resources by enabling accurate assessments of pavement conditions. Most prior reviews, however, have concentrated on flexible pavements, leaving a notable gap regarding rigid and composite pavements. To address this gap, the present study conducts a systematic review of Artificial Intelligence (AI) methods applied to IRI prediction for rigid and composite pavements. Literature published between 2004 and 2025 is synthesized to highlight prevailing trends, methodological contributions, and directions for future research. Particular attention is given to the types of models employed, the datasets used for training and validation, and the role of input variables and data-processing strategies. Across the included studies, ensemble learning methods (especially gradient boosting variants such as XGBoost), artificial neural networks, and hybrid architectures frequently achieved high predictive skill, with several models reporting test-set coefficients of determination approaching 0.9–0.96, indicating strong potential for capturing the influence of traffic, pavement structure, and climatic factors. Since these results are obtained from heterogeneous datasets and evaluation protocols, they are interpreted qualitatively rather than as strict cross-study rankings. Analysis of input variables revealed that pavement age and initial IRI were included in 91% (21 of 23) and 78% (18 of 23) of studies, respectively. Climatic variables such as the freezing index appeared in 57% (13 of 23), while traffic-related factors were considered in 65% (15 of 23). The findings underscore the importance of standardized, high-quality datasets, such as those from the Long-Term Pavement Performance (LTPP) program, along with data consistency, model interpretability, computational efficiency, and replicability in enhancing IRI prediction. Future research should focus on incorporating input variable selection techniques to identify the most influential predictors, thereby improving accuracy and robustness. Integrating these approaches with advanced non-linear data-driven models, coupled with robust hyperparameter optimization, holds considerable promise for strengthening the reliability of IRI prediction and supporting resilient pavement management strategies.

42 ENGINEERING↗
DOE OSTI2024

Data and scripts associated with the manuscript "Organic Molecules are Deterministically Assembled in River Sediments"

This data package is associated with the publication "Organic Molecules are Deterministically Assembled in River Sediments" submitted to Scientific Reports (Stegen et al., 2024). The study applies community ecology methods to dissolved organic matter (DOM) chemistry from variably inundated riverbed sediments to uncover principles governing DOM composition at a reach-scale. This data package documents the workflow used to process and generate the main findings in the manuscript. The R scripts reference the raw, unprocessed Fourier transform ion cyclotron resonance mass spectrometry (FTICR-MS) data from another data package, available on ESS-DIVE at https://data.ess-dive.lbl.gov/view/doi:10.15485/1834208. The scripts then process the raw FTICR-MS data and generate the findings and figures presented in the associated manuscript. In brief, this study demonstrates that DOM assemblages in variably inundated sediments are primarily governed by deterministic variable selection, including sediment moisture effecting the degree of deterministic assembly. See the manuscript for more details pertaining to interpretation and implications of the findings. This data package is associated with the GitHub repository found at https://github.com/WHONDRS-Hub/ECA_2020_Sed.This data package is comprised of 6 scripts and 7 folders. The file-level metadata file (file ending in "flmd.csv") lists all files contained in this data package and descriptions for each. The data dictionary (file ending in "dd.csv) describes all tabular data columns and their respective definitions and units. The FTICR_Processing_Scripts produce the outputs found in the "Processed_Data" folder. The remaining scripts (located in the parent directory) produce the outputs found in the following four folders: (1) "MCD_Dendrograms", "MCD_Randomizations", "MCD_bNTI_Outcomes", and "OM_Null_Modeling". The fifth script additionally takes the three comma-separated values (CSV) files found in the parent directory as input ("VGC_texture.csv", "merged_weights.csv", and "ECA2_FTICR_BetaDisp.csv"). The outputs of each of the five scripts serve as the input to the following script, with the final outputs stored in the folder "OM_Null_Modeling".

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Leaf gas exchange data, Nov 2018 and May 2019, Daintree Rainforest, Australia

This data package contains leaf gas exchange data for 11 canopy tree species in the control and drought treatments at the Daintree Rainforest Observatory in Queensland, Australia (AU-DRO) with the purpose of testing for differences in gas exchange between treatments. Data were collected at the end of the dry season (November 2018) and end of the wet season (May 2019). Using a canopy crane, large branches (>1 m) were collected from the top of the canopy and immediately placed in a bucket of water, upon which the LI-6800 leaf chamber was clamped onto a recently mature leaf with minimal herbivory. Measurements were recorded every 3 seconds for 3 minutes total. The species included in the dataset are: Castanospermum australe, Cleistanthus myrianthus, Cryptocarya mackinnoniana, Dysoxylum papuanum, Elaeocarpus angustifolius (Elaeocarpus grandis), Endiandra microneura, Mallotus paniculatus, Myristica globosa, Rockinghamia angustifolia, Synima cordierorum, Syzygium spp. The two data files included are: 1. LICOR_Data_AU_DRO_2019_GasEx.csv (a compilation of individual LICOR output files, with correction of user input errors such as incorrect species spelling) 2. ESS_DIVE_Data_AU_DRO_2019_GasEx.csv (selected variables and converted units from LICOR_Data_AU_DRO_2019_GasEx.csv to follow the standard gas exchange reporting format from Ely et al. 2021 Ecological Informatics). The metadata files included are: (3) ESS_DIVE_Metadata_AU_DRO_2019_GasEx.xlsx (4) File_Submission_Metadata_AU_DRP_2019_GasEx.xlsx (5) E-Field_Log_AU_DRO_2019_GasEx.xlsx (6) Data_Column_Description_AU_DRO_2019_GasEx.xlsx.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

sPHENIX heavy flavor jet tagging studies in p+p at $\sqrt{s_{NN}}=200~GeV$

Heavy-flavor jets, which are initiated from heavy quarks, are ideal probes for studying flavor dependent parton energy loss. We report on the performance of jet flavor tagging using two Neural Network Machine Learning (ML) models: the Long Short-Term Memory (LSTM) model and an Attention-based Neural Network, in simulations of 200 GeV p + p collisions. The tagging performance of bottom quark initiated jets with both ML models surpasses that of the traditional cut-based method. Technical details, including sample and kinematic variable selections, the machine learning training and testing setup with parameter tuning, and outcome comparisons, will be discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Inverse Modeling of Hydrologic Parameters in CLM4 via Generalized Polynomial Chaos in the Bayesian Framework

In this work, generalized polynomial chaos (gPC) expansion for land surface model parameter estimation is evaluated. We perform inverse modeling and compute the posterior distribution of the critical hydrological parameters that are subject to great uncertainty in the Community Land Model (CLM) for a given value of the output LH. The unknown parameters include those that have been identified as the most influential factors on the simulations of surface and subsurface runoff, latent and sensible heat fluxes, and soil moisture in CLM4.0. We set up the inversion problem in the Bayesian framework in two steps: (i) building a surrogate model expressing the input–output mapping, and (ii) performing inverse modeling and computing the posterior distributions of the input parameters using observation data for a given value of the output LH. The development of the surrogate model is carried out with a Bayesian procedure based on the variable selection methods that use gPC expansions. Our approach accounts for bases selection uncertainty and quantifies the importance of the gPC terms, and, hence, all of the input parameters, via the associated posterior probabilities.

97 MATHEMATICS AND COMPUTING↗