Search NASA⌕ Search

SEARCH · Search NASA

Results for “correlation coefficient”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Investigating the ecological fallacy through sampling distributions constructed from finite populations

Correlation coefficients and linear regression values computed from group averages can differ from correlation coefficients and linear regression values computed using individual scores. This observation known as the ecological fallacy often assumes that all the individual scores are available from a population. In many situations, one must use a sample from the larger population. In such cases, the computed correlation coefficient and linear regression values will depend on the sample that is chosen and the underlying sampling distribution. The sampling distribution of correlation coefficients and linear regression values for group averages will be identical to the sampling distribution for individuals for normally distributed variables for random samples drawn from infinitely large continuous distributions. However, data that is acquired in practice is often acquired when sampling without replacement from a finite population. Our objective is to demonstrate through Monte Carlo simulations that the sampling distributions for correlation and linear regression will also be similar for individuals and group averages when sampling without replacement from normally distributed variables. These simulations suggest that when a random sample from a population is selected, the correlation coefficients and linear regression values computed from individual scores will not be more accurate in estimating the entire population values compared to samples when group averages are used as long as the sample size is the same.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

The pNab experiment and the quest for ever better neutron beam polarization

The Nab and pNAB collaborations are conducting a program of studies of free neutron beta decay, with the primary goal of testing the unitarity of the Cabbibo-Kobayashi-Maskawa matrix that describes quark mixing due to the weak interaction. For this purpose, a large, novel electromagnetic spectrometer, the Nab spectrometer, has been designed, built, and placed in use to determine the correlation coefficients in unpolarized neutron beta decay: "a", the neutrino-electron correlation coefficient; and "b", the Fierz term. The subject of this paper is pNAB, the second phase of the program, that will deploy the same spectrometer with a polarized neutron beam to determine "A", the beta asymmetry; and "B", the neutrino asymmetry coefficients. A focus of this paper will be on the strategies to provide a high and precisely known neutron beam polarization.

Baessler, Stefan↗
DOE OSTI2023

Advanced Fuel Cycle Cost Basis Report: Supporting Documents 3 Cost Correlations

This chapter provides an overview of the drivers for developing correlation coefficients, status on efforts to develop partial correlation coefficients, and recommended next steps. Mathematical methods for developing and using correlation coefficients were previously documented by the FCO EWG in Chapters 6 and 7 of [Ganda 2014].

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Statistical Relations among Solid Precipitation, Atmospheric Moisture and Cloud Parameters in the Arctic

Observations collected during cold-season precipitation periods at Utquagvik, Alaska and at the multidisciplinary drifting observatory for the study of Arctic climate (MOSAiC) are used to statistically analyze the relations among the atmospheric water cycle parameters including the columnar supercooled liquid and ice amounts (expressed as liquid-water and ice-water paths, i.e., LWP and IWP), the integrated water vapor (IWV) and the near-surface snowfall rate. Data come from radar and radiometer-based retrievals and from optical precipitation sensors. While the correlation between snowfall rate and LWP is rather weak, correlation coefficients between radar-derived snowfall rate and IWP are high (~0.8), which is explained, in part, by the generally low LWP/IWP ratios during significant precipitation. Correlation coefficients between snowfall rate and IWV are moderate (~0.45). Correlations are generally weaker if snowfall is estimated by optical sensors, which is, in part, due to blowing snow. Correlation coefficients between near-surface temperature and snowfall rates are low (r < 0.3). The results from the Alaska and MOSAiC sites are generally similar. These results are not very sensitive to the amount of time averaging (e.g., 15 min averaging versus daily averages). Observationally based relations among the water cycle parameters are informative about atmospheric moisture conversion processes and can be used for model evaluations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Uncertain of uncertainties? A comparison of uncertainty quantification metrics for chemical data sets

Abstract With the increasingly more important role of machine learning (ML) models in chemical research, the need for putting a level of confidence to the model predictions naturally arises. Several methods for obtaining uncertainty estimates have been proposed in recent years but consensus on the evaluation of these have yet to be established and different studies on uncertainties generally uses different metrics to evaluate them. We compare three of the most popular validation metrics (Spearman’s rank correlation coefficient, the negative log likelihood (NLL) and the miscalibration area) to the error-based calibration introduced by Levi et al. ( Sensors 2022 , 22 , 5540). Importantly, metrics such as the negative log likelihood (NLL) and Spearman’s rank correlation coefficient bear little information in themselves. We therefore introduce reference values obtained through errors simulated directly from the uncertainty distribution. The different metrics target different properties and we show how to interpret them, but we generally find the best overall validation to be done based on the error-based calibration plot introduced by Levi et al. Finally, we illustrate the sensitivity of ranking-based methods (e.g. Spearman’s rank correlation coefficient) towards test set design by using the same toy model ferent test sets and obtaining vastly different metrics (0.05 vs. 0.65).

Rasmussen, Maria H.↗
DOE OSTI2024

In Search of Optimum Fresh-Cut Raw Material: Using Computer Vision Systems as a Sensory Screening Tool for Browning-Resistant Romaine Lettuce Accessions

The popularity of ready-to-eat (RTE) salads has prompted novel technology to prolong the shelf life of their ingredients. Fresh-cut romaine lettuce is widely used in RTE salads; however, its tendency to quickly discolor continues to be a challenge for the industry. Selecting the ideal lettuce accessions for use in RTE salads is essential to ensure maximum shelf life, and it is critical to have a practical way to assess and compare the quality of multiple lettuce accessions that are being considered for use in fresh-cut applications. Thus, in this work we aimed to determine whether a computer vision system (CVS) composed of image acquisition, processing, and analysis could be effective to detect visual quality differences among 16 accessions of fresh-cut romaine lettuce during postharvest storage. The CVS involved a post-capturing color correction, effective image segmentation, and calculation of a browning index, which was tested as a predictor of quality and shelf life of fresh-cut romaine lettuce. The results demonstrated that machine vision software can be implemented to replace or supplement the scoring of a trained panel and instrumental quality measurements. Overall visual quality, a key sensory parameter that determines food preferences and consumer behavior, was highly correlated with the browning index, with a Pearson correlation coefficient of −0.85. Other important sensory decision parameters were also strongly or moderately correlated with the browning index, with Pearson correlation coefficients of −0.84 for freshness, 0.79 for off odor, and 0.57 for browning. The ranking of the accessions according to quality acceptability from the sensory evaluation produced a similar pattern to those obtained with the CVS. This study revealed that multiple lettuce accessions can be effectively benchmarked for their performance as fresh-cut sources via a CVS-based method. Future opportunities and challenges in using machine vision image processing to predict consumer preferences for RTE salad greens is also discussed.

Agriculture↗
DOE OSTI2025

Ensemble‐Based, Large‐Eddy Reconstruction of Wind Turbine Inflow in a Near‐Stationary Atmospheric Boundary Layer Through Generative Artificial Intelligence

ABSTRACT To validate the second‐by‐second dynamics of turbines in field experiments, it is necessary to accurately reconstruct the winds going into the turbine. Current time‐resolved inflow reconstruction techniques estimate wind behavior in unobserved regions using relatively simple spectral‐based models of the atmosphere. Here, we develop a technique for time‐resolved inflow reconstruction that is rooted in a large‐eddy simulation model of the atmosphere. Our “large‐eddy reconstruction” technique blends observations and atmospheric model information through a diffusion model machine learning algorithm, allowing us to generate probabilistic ensembles of reconstructions for a single 10‐min observational period. Our generated inflows can be used directly by aeroelastic codes or as inflow boundary conditions in a large‐eddy simulation. We verify the second‐by‐second reconstruction capability of our technique in three synthetic field campaigns, finding positive Pearson correlation coefficient values () between ground‐truth and reconstructed streamwise velocity, as well as smaller positive correlation coefficient values for unobserved fields (spanwise velocity, vertical velocity, and temperature). We validate our technique in three real‐world case studies by driving large‐eddy simulations with reconstructed inflows and comparing to independent inflow measurements. The reconstructions are visually similar to measurements, follow desired power spectra properties, and track second‐by‐second behavior ().

17 WIND ENERGY↗
DOE OSTI2024

Approximation of ice phenology of Maine lakes using Aqua MODIS surface temperature data

Studies of lake ice phenology have historically relied on limited in situ data. Relatively few observations exist for ice out and fewer still for ice in, both of which are necessary to determine the temporal extent of ice cover. Satellite data provide an opportunity to better document patterns of ice phenology across landscapes and relate them to the climatological drivers behind changing ice phenology. We developed a model, the Cumulative Sum Method (CSM), that uses daytime and nighttime surface temperature observations from the Moderate Resolution Imaging Spectroradiometer (MODIS) sensor on board the Earth-observing Aqua satellite to approximate ice in (the onset of ice cover) and ice out from training datasets of 13 and 58 Maine lakes, respectively, during the 2002/2003 through 2017/2018 ice seasons. Ice in was signaled by reaching a threshold of cumulative negative degrees following the first day of the season below 0°C. Ice out was signaled by reaching a threshold of cumulative positive degrees following the first day of the year above 0°C. The comparison of observed and remotely sensed ice-in dates showed relative agreement with a correlation coefficient of 0.71 and a mean absolute error (MAE) of 9.8 days. Ice-out approximations had a correlation coefficient of 0.67 and an MAE of 8.8 days. Lakes smaller in surface area and nearer the Atlantic coast had the greatest error in approximation. Application of the CSM to 20 additional lakes in Maine produced a comparable ice-out MAE of 8.9 days. Ice-out model performance was weaker for the warmest years; there was a larger MAE of 12.0 days when the model was applied to the years 2019–2023 for the original 58 lakes. The development of this model, which utilizes daily satellite data, demonstrates the promise of remote sensing for quantifying ice phenology over short, temporal scales, and wider geographic regions than can be observed in situ, and allows exploration of the influence of surface temperature patterns on the process and timing of ice in and ice out.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Velocity reconstruction in the era of DESI and Rubin/LSST. II. Realistic samples on the light cone

Reconstructing the galaxy peculiar velocity field from the distribution of large-scale structure plays an important role in cosmology. On one hand, it gives us an insight into structure formation and gravity; on the other, it allows us to selectively extract the kinetic Sunyaev-Zel’dovich (kSZ) effect from cosmic microwave background maps. In this work, we employ high-accuracy synthetic galaxy catalogs on the light cone to investigate how well we can recover the velocity field when utilizing the three-dimensional spatial distribution of the galaxies in a modern large-scale structure experiment such as the Dark Energy Spectroscopic Instrument (DESI) and the Rubin Observatory Legacy Survey of Space and Time. In particular, we adopt the standard technique used in baryon acoustic oscillation analysis for reconstructing the Zel’dovich displacements of galaxies through the continuity equation, which yields a first-order approximation to their large-scale velocities. We investigate variations in the number density, bias, mask, area, redshift noise, and survey depth, as well as modifications to the settings of the standard reconstruction algorithm. Since our main goal is to provide guidance for planned kSZ analysis between DESI and the Atacama Cosmology Telescope, we apply velocity reconstruction to a faithful representation of DESI spectroscopic and photometric targets. We report the cross-correlation coefficient between the reconstructed and the true velocities along the line of sight. For the DESI Y1 spectroscopic survey, we expect the correlation coefficient to be r ≈ 0.64, while for a photometric survey with δ z /(1+z) = 0.02, as is approximately the case for the Legacy Survey used in the target selection of DESI galaxies, r shrinks by half to r ≈ 0.31. Here, we hope the results in this paper can be used to inform future kSZ stacking studies and other velocity reconstruction analyses planned with the next generation of cosmology experiments.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Evaluation of WSR-88D Level III and MRMS Rainfall Estimates Against Rain Gauge Observations at the Savannah River Site

Rainfall data for the Savannah River Site (SRS) has been historically measured by rain gauges. These instruments serve as ground truth for most climatological and weather applications; however, gauge measurements are prone to errors or biases under certain weather conditions. Rainfall estimates from radar reflectivity values have been developed and improved over the years and serve as an alternative or supplement for gauge measurements. This study compares measurements from tipping bucket rain gauges with rainfall estimates from the NOAA NWS WSR-88D Level III hourly rainfall product and the Multi-Radar Multi-Sensor (MRMS) gauge-corrected precipitation estimate. Results show good agreement between radar derived amounts and ground measurements, with MRMS values showing correlation coefficients of 0.60-0.95 and RMSE values of less than 1.0 cm (0.4 in). The WSR-88D Level III estimates result in correlation coefficients of 0.46-0.86 and RMSE values of less than 1.3 cm (0.5 in). Few outliers are observed for each data pair and are evaluated against precipitation classification products (WSR-88D Hybrid Hydrometeor Classification product and MRMS Precipitation Flag product). The rain gauges used in this study are not part of the Hydrometeorological Automated Data System (HADS) network used to correct the MRMS estimates and therefore, results of this work provide an independent validation of the MRMS gauge correction scheme.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Correlations between azimuthal anisotropy and mean transverse momentum in pp, pPb, and peripheral PbPb collisions

Correlations between azimuthal anisotropy and mean transverse momentum of charged particles in proton-proton (pp), proton-lead (pPb), and peripheral lead-lead (PbPb) collisions are presented as a function of charged particle multiplicity. The pp, pPb and PbPb collision data were collected using the CMS detector at the LHC with a center-of-mass energy per nucleon pair of 13, 8.16 and 5.02 TeV, respectively. The two- and four-particle cumulants for the second- and third-order Fourier anisotropy harmonics are correlated with the mean transverse momentum of charged particles on an event-by-event basis. In pp and pPb systems, the observed correlation coefficients based on two-particle cumulants are found to change from negative to positive values as the charged particle multiplicity decreases. The sign changes disappear when the correlated particles are required to be further apart in pseudorapidity. Additionally, no sign changes in correlation coefficients are observed when employing four-particle cumulants. Models incorporating initial-state gluon saturation and final-state hydrodynamic evolutions are compared to pPb data and the predicted sign changes are not observed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Comparison of Machine Learning-Based Predictive Models of the Nutrient Loads Delivered from the Mississippi/Atchafalaya River Basin to the Gulf of Mexico

Predicting nutrient loads is essential to understanding and managing one of the environmental issues faced by the northern Gulf of Mexico hypoxic zone, which poses a severe threat to the Gulf’s healthy ecosystem and economy. The development of hypoxia in the Gulf of Mexico is strongly associated with the eutrophication process initiated by excessive nutrient loads. Due to the complexities in the excessive nutrient loads to the Gulf of Mexico, it is challenging to understand and predict the underlying temporal variation of nutrient loads. The study was aimed at identifying an optimal predictive machine learning model to capture and predict nonlinear behavior of the nutrient loads delivered from the Mississippi/Atchafalaya River Basin (MARB) to the Gulf of Mexico. For this purpose, monthly nutrient loads (N and P) in tons were collected from US Geological Survey (USGS) monitoring station 07373420 from 1980 to 2020. Machine learning models—including autoregressive integrated moving average (ARIMA), gaussian process regression (GPR), single-layer multilayer perceptron (MLP), and a long short-term memory (LSTM) with the single hidden layer—were developed to predict the monthly nutrient loads, and model performances were evaluated by standard assessment metrics—Root Mean Square Error (RMSE) and Correlation Coefficient (R). The residuals of predictive models were examined by the Durbin–Watson statistic. The results showed that MLP and LSTM persistently achieved better accuracy in predicting monthly TN and TP loads compared to GPR and ARIMA. In addition, GPR models achieved slightly better test RMSE score than ARIMA models while their correlation coefficients are much lower than ARIMA models. Moreover, MLP performed slightly better than LSTM in predicting monthly TP loads while LSTM slightly outperformed for TN loads. Furthermore, it was found that the optimizer and number of inputs didn’t show effects on the LSTM performance while they exhibited impacts on MLP outcomes. This study explores the capability of machine learning models to accurately predict nonlinearly fluctuating nutrient loads delivered to the Gulf of Mexico. Further efforts focus on improving the accuracy of forecasting using hybrid models which combine several machine learning models with superior predictive performance for nutrient fluxes throughout the MARB.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Characterization of the quantum state of top quark pairs produced in proton-proton collisions at $\sqrt{s}$ = 13 TeV using the beam and helicity bases

Measurements of the spin correlation coefficients in the beam basis are presented for top quark-antiquark $(t\bar{t}⁢)$ systems produced in proton-proton collisions at $\sqrt{s}$ = 13 TeV collected by the CMS experiment in 2016–2018, and corresponding to an integrated luminosity of 138 fb −1 . The $t\bar{t}$ system is reconstructed from final states containing an electron or muon and jets. Together with the previously reported results in the helicity basis, these measurements are used to decompose the system into the Bell and spin eigenstates in various kinematic regions. The spin correlation coefficients are also used to evaluate properties of the $t\bar{t}$ quantum state, such as the purity, von Neumann entropy, and entanglement. All results are consistent with standard model predictions.

Hayrapetyan, Aram [Yerevan Physics Institute]↗
DOE OSTI2024

Unique impacts of strong and westward-extended western Pacific subtropical high on ozone pollution over eastern China

As a subtropical anticyclonic high-pressure system that typically forms over the northwestern Pacific Ocean in summer, the Western Pacific subtropical high (WPSH) affects meteorological conditions and ozone pollution in China. The relationship between maximum daily 8-h average ozone (MDA8 O 3 ) concentrations and the extremely strong and westward-extended WPSH occurred in 2022 is investigated using observations, reanalysis data and GEOS-Chem model simulations. The intensity of WPSH has a significant positive correlation with MDA8 O3 over southern China during July-August in 2022, with a correlation coefficient of +0.44, but the correlation is negative (–0.40) in northern China. During the strong WPSH days, MDA8 O 3 increased by 16.5µgm -3 (16.4% relative to July-August average) over southern China and decreased by 19.0µgm -3 (14.5%) in northern China compared to the weak WPSH days. The unique dipole pattern in the relationship between ozone levels and the WPSH in 2022 exhibited a contrast to that during 2015–2021. The difference is primarily due to the extremely strong WPSH intensity and its unusual westward expansion in 2022. In this case, an anomalous anticyclone at 500 hPa dominates over southern China, which creates conditions conducive for ozone formation and accumulation. The anticyclone weakened horizontal winds and reduced the dispersion of ozone, alongside a high temperature and low relative humidity, which favored the chemical production of ozone. In contrast, abnormal northerly winds enhanced ozone diffusion in northern China and the low temperature reduced ozone chemical production. Here, this study reveals the mechanism for the significant impact of strong and westward-extended WPSH on ozone concentrations over China, emphasizing the role of the WPSH location in modulating meteorology and ozone levels.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Empirical radius formulas for canonical neutron stars from bidirectionally selecting features of equations of state in extended Bayesian analyses of observational data

Significant advancement in Bayesian inference of nuclear equation of state (EOS) from gravitational wave and x-ray observations of neutron stars (NSs) has been made by the nuclear astrophysics community especially since GW170817. By extending the traditional Bayesian analysis which normally ends at presenting the marginalized posterior probability distribution functions (PDFs) of individual EOS parameters and their correlations (or sometimes only the Pearson correlation coefficients which are only reliably useful when the variables are linearly correlated while they are actually often not), we search for a data-driven and robust empirical formula for the radius 𝑅 1.4 of canonical NSs in terms of the characteristic EOS parameters (features). We also identify the single most important but currently poorly known EOS parameter for determining the 𝑅 1.4 . Using three regression-model-building methodologies: bidirectional stepwise feature selection, least absolute shrinkage selection operator (LASSO) regression, and neural network regression on a large set of posterior EOSs and the corresponding 𝑅 1.4 values inferred from earlier comprehensive Bayesian analyses of NS observational data, we systematically and rigorously develop the most probable 𝑅 1.4 formulas with varying statistical accuracy and technical complexity. Here, the most important EOS parameters for determining 𝑅 1.4 are found consistently in each of the feature selection processes to be (in order of decreasing importance): curvature 𝐾 sym , slope 𝐿, skewness 𝐽 sym of nuclear symmetry energy, skewness 𝐽 0 , incompressibility 𝐾 0 of symmetric nuclear matter, and the magnitude 𝐸 sym ⁡(𝜌 0 ) of symmetry energy at the saturation density 𝜌 0 of nuclear matter.

Bayesian methods↗
DOE OSTI2026

Fast baryonic field painting for Sunyaev-Zel’dovich analyses: Transfer function vs hybrid effective field theory

Here, we present two approaches for “painting” baryonic properties relevant to the Sunyaev-Zel’dovich (SZ) effect—optical depth and Compton-y—onto three-dimensional N-body simulations, using the MillenniumTNG suite as a benchmark. The goal of these methods is to produce fast and accurate reconstruction methods to aid future analyses of baryonic feedback using the SZ effect. The first approach employs a Gaussian process emulator to model the SZ quantities via a transfer function, while the second utilizes hybrid effective field theory (HEFT) to reproduce these quantities within the simulation. Our analysis involves comparing both methods to the true MillenniumTNG optical depth and Compton-y fields using several metrics, including the cross-correlation coefficient, power spectrum, and power spectrum error. Additionally, we assess how well the reconstructed fields correlate with dark matter haloes across various mass thresholds. The results indicate that the transfer function method yields more accurate reconstructions for fields with initially high correlations (r ≈ 1), such as between the optical depth and dark matter fields. Conversely, the HEFT-based approach proves more effective in enhancing correlations for fields with weaker initial correlations (r ∼ 0.5), such as between the Compton-y and dark matter fields. Lastly, we discuss extensions of our methods to improve the reconstruction performance at the field level.

Liu, R. Henry [University of California, Berkeley,↗
DOE OSTI2024

Identification of Suitable Vacuum Gas Oils as Plasticizers Using HT-GC × GC-HRMS

Vacuum gas oils (VGOs) have long been heralded as effective plasticizers due to their high boiling-points and lubricating properties. One of these VGOs, HyVac Oil 93050, is a paraffinic plasticizer that is currently utilized in polymer–plasticizer explosive formulations, while this oil effectively increases the elasticity and decreases the sensitivity for current formulations. We seek to identify alternative VGOs with similar density, viscosity, molecular composition, and impurities for future formulations. In this study, 18 VGOs, including HyVac Oil 93050, were initially evaluated for their density and viscosity. The oils were then ranked based on physical characteristics and analyzed for molecular composition using high-temperature comprehensive two-dimensional gas chromatography with high-resolution time-of-flight mass spectrometry (HT-GC × GC-HRMS). Further, the HT-GC × GC-HRMS chromatograms of the top 10 most similar VGOs to HyVac Oil 93050, according to density, were compared utilizing the Pearson correlation coefficient. Three of the oils with densities similar to those of HyVac Oil 93050 had a Pearson correlation within the lot-to-lot variation of HyVac Oil 93050. For the top 3 candidates, Pearson correlation was then utilized as a feature selection technique to discover significant chemical differences. Positive chemical ionization (PCI) and negative chemical ionization (NCI) HT-GC × GC-HRMS chromatograms were also evaluated and aided in the discovery and identification of several key compounds including additives that acted as stabilizers for the VGOs. After this further chemical analysis, two of the original 17 potential VGOs were determined to be physically and chemically similar to HyVac Oil 93050.

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

Bayesian merging of numerical modeling and remote sensing for saltwater intrusion quantification in the Vietnamese Mekong Delta

Saltwater intrusion has become one of the most concerning issues in the Vietnamese Mekong Delta (VMD) due to its increasing impacts on agriculture and food security of Vietnam. Reliable estimation of salinity plays a crucial role to mitigate the impacts of saltwater intrusion. Here, this study developed a hybrid technique that merges satellite imagery with numerical simulations to improve the estimation of salinity in the VMD. The salinity derived from Landsat images and by numerical simulations was fused using the Bayesian inference technique. The results indicate that our technique significantly reduces the uncertainties and improves the accuracy of salinity estimates. The Nash–Sutcliffe coefficient is 0.74, which is much higher than that of numerical simulation (0.63) and Landsat estimation (0.6). The correlation coefficient between the ensemble and measured salinity is relatively high (0.88). The variance of the ensemble salinity errors (5.0 ppt 2 ) is lower than that of Landsat estimation (10.4 ppt 2 ) and numerical simulations (9.6 ppt 2 ). The proposed approach shows a great potential to combine multiple data sources of a variable of interest to improve its accuracy and reliability wherever these data are available.

54 ENVIRONMENTAL SCIENCES↗