Search NASA⌕ Search

SEARCH · Search NASA

Results for “linear estimation”

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

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

At least 1,729 records · Page 96
All sourcesNASA NTRSDOE OSTI
NASA NTRS2014

Analyzing Source Apportioned Methane in Northern California During DISCOVER-AQ-CA Using Airborne Measurements and Model Simulations

This study analyzes source apportioned methane (CH4) emissions and atmospheric concentrations in northern California during the Discover-AQ-CA field campaign using airborne measurement data and model simulations. Source apportioned CH4 emissions from the Emissions Database for Global Atmospheric Research (EDGAR) version 4.2 were applied in the 3-D chemical transport model GEOS-Chem and analyzed using airborne measurements taken as part of the Alpha Jet Atmospheric eXperiment over the San Francisco Bay Area (SFBA) and northern San Joaquin Valley (SJV). During the time period of the Discover-AQ-CA field campaign EDGAR inventory CH4 emissions were 5.30 Gg/day (Gg 1.0 109 grams) (equating to 1.9 103 Gg/yr) for all of California. According to EDGAR, the SFBA and northern SJV region contributes 30 of total emissions from California. Source apportionment analysis during this study shows that CH4 concentrations over this area of northern California are largely influenced by global emissions from wetlands and local/global emissions from gas and oil production and distribution, waste treatment processes, and livestock management. Model simulations, using EDGAR emissions, suggest that the model under-estimates CH4 concentrations in northern California (average normalized mean bias (NMB) -5 and linear regression slope 0.25). The largest negative biases in the model were calculated on days when hot spots of local emission sources were measured and atmospheric CH4 concentrations reached values 3.0 parts per million (model NMB -10). Sensitivity emission studies conducted during this research suggest that local emissions of CH4 from livestock management processes are likely the primary source of the negative model bias. These results indicate that a variety, and larger quantity, of measurement data needs to be obtained and additional research is necessary to better quantify source apportioned CH4 emissions in California and further the understanding of the physical processes controlling them.

Methane↗
NASA NTRS2016

Influence of Population Variation of Physiological Parameters in Computational Models of Space Physiology

The earliest manifestations of Visual Impairment and Intracranial Pressure (VIIP) syndrome become evident after months of spaceflight and include a variety of ophthalmic changes, including posterior globe flattening and distension of the optic nerve sheath. Prevailing evidence links the occurrence of VIIP to the cephalic fluid shift induced by microgravity and the subsequent pressure changes around the optic nerve and eye. Deducing the etiology of VIIP is challenging due to the wide range of physiological parameters that may be influenced by spaceflight and are required to address a realistic spectrum of physiological responses. Here, we report on the application of an efficient approach to interrogating physiological parameter space through computational modeling. Specifically, we assess the influence of uncertainty in input parameters for two models of VIIP syndrome: a lumped-parameter model (LPM) of the cardiovascular and central nervous systems, and a finite-element model (FEM) of the posterior eye, optic nerve head (ONH) and optic nerve sheath. Methods: To investigate the parameter space in each model, we employed Latin hypercube sampling partial rank correlation coefficient (LHSPRCC) strategies. LHS techniques outperform Monte Carlo approaches by enforcing efficient sampling across the entire range of all parameters. The PRCC method estimates the sensitivity of model outputs to these parameters while adjusting for the linear effects of all other inputs. The LPM analysis addressed uncertainties in 42 physiological parameters, such as initial compartmental volume and nominal compartment percentage of total cardiac output in the supine state, while the FEM evaluated the effects on biomechanical strain from uncertainties in 23 material and pressure parameters for the ocular anatomy. Results and Conclusion: The LPM analysis identified several key factors including high sensitivity to the initial fluid distribution. The FEM study found that intraocular pressure and intracranial pressure had dominant impact on the peak strains in the ONH and retro-laminar optic nerve, respectively; optic nerve and lamina cribrosa stiffness were also important. This investigation illustrates the ability of LHSPRCC to identify the most influential physiological parameters, which must therefore be well-characterized to produce the most accurate numerical results.

biomechanics↗
DOE OSTI2022

Correlation Calculations for the Russian Pu Metal Fast Experiments

Nuclear criticality experiments are often conducted in campaigns with multiple variations. These experiments reuse the same basic components, like the fuel, moderator, or positioning machines. The components have uncertainties in their geometry and composition that propagate to models of the experiments. Shared components create shared uncertainty between the $k_{eff}$ of benchmarks. The shared uncertainty is commonly quantified with a covariance, or correlation coefficient. These covariances can impact criticality safety and nuclear data validation applications. While benchmark evaluations tabulate an experiment’s uncertainty, they often lack a detailed calculation of correlations between experiments. Even some very commonly used benchmarks, like the Russian Pu Metal Fast (PMF) experiments, have missing correlations. This paper presents our approach to calculate the correlations for five of the Russian PMF experiments. The experiments share hemispherical Pu shells that induce a correlation between modeled $k_{eff}$ values. We estimated the correlations with simplified and detailed models of the experiments through linear-perturbation theory. The correlations between the experiments vary significantly between the detailed vs. simplified models. We also investigate how the correlations affect validation metrics of the experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2015

Coming of Age: Polarization as a Probe of Plant Canopy Water Status

We tested the hypothesis that the relative water content (RWC) of the sunlit leaves in a plant canopy may be estimated from polarized canopy imagery. Recently (IGARSS, July 27-31, 2015, Milan, Italy), we reported the results of laboratory polarization measurements of single detached leaves during dry down. We found that RWC was linearly related to the ratio of the reflectance of the interior of the leaf and the leaf transmittance. Here we report application of the laboratory results to estimate RWC for sunlit leaves in a plant canopy. Using a commercial-off-the-shelf (COTS) Nikon 810 camera with Nikkor 300 mm lens and Polaroid type HN-22 linear polarizer, we photographed in the principle plane a plant canopy displaying a gradient of water stress and collected, at each of multiple points along the gradient, two images, one with the polarization filter oriented for maximum scene response and a second with the filter oriented for minimum scene response. We converted the digital values in the two images to reflectance factor with reference to images of a white, flat, horizontal Spectralon surface. We classified the polarization imagery, identifying reflecting leaves, transmitting leaves, other sunlit vegetation and shadows. For each image pair we normalized the leaf internal reflectance by dividing by the cosine of the angle of incidence of the sunlight on the leaf, selected the leaf maximum transmittance in the scene and divided to obtain the ratio reflectance/transmittance, which we compared with leaf RWC. We determined the leaf relative water content by harvesting a section of leaf and immediately placing it in a sealed container in an ice chest. Later in the laboratory the leaf sample was weighed, rehydrated, weighed, dried and again weighed. RWC was determined using the standard formula. Our experimental results support our hypothesis, suggesting that the RWC of sunlit leaves in a plant canopy may be estimated from analysis of polarization imagery collected by a COTS camera system. Unlike remotely sensed estimates of canopy equivalent water thickness, our estimates of the RWC of sunlit canopy leaves provide leaf physiological information. We propose RWC estimates based upon sunlit leaves are more relevant to assessing the water status of a plant canopy than would be RWC estimates based upon large FOV canopy measurements.

estimates↗
DOE OSTI2023

Evaluation of dual-weighted residual and machine learning error estimation for projection-based reduced-order models of steady partial differential equations

Projection-based reduced-order models (pROMs) show great promise as a means to accelerate many-query applications such as forward error propagation, solving inverse problems, and design optimization. In order to deploy pROMs in the context of high-consequence decision making, accurate error estimates are required to determine the region(s) of applicability in the parameter space. The following paper considers the dual-weighted residual (DWR) error estimate for pROMs and compares it to another promising pROM error estimate, machine learned error models (MLEM). Here, we show how DWR can be applied to ROMs and then evaluate DWR on two partial differential equations (PDEs): a two-dimensional linear convection–reaction–diffusion equation, and a three-dimensional static hyper-elastic beam. It is shown that DWR is able to estimate errors for pROMs extrapolating outside of their training set while MLEM is best suited for pROMs used to interpolate within the pROM training set.

42 ENGINEERING↗
NASA NTRS2003

Deviation of Long-Period Tides from Equilibrium: Kinematics and Geostrophy

New empirical estimates of the long-period fortnightly (Mf) tide obtained from TOPEX/Poseidon (T/P) altimeter data confirm significant basin-scale deviations from equilibrium. Elevations in the low-latitude Pacific have reduced amplitude and lag those in the Atlantic by 30 deg or more. These interbasin amplitude and phase variations are robust features that are reproduced by numerical solutions of the shallow-water equations, even for a constant-depth ocean with schematic interconnected rectangular basins. A simplified analytical model for cooscillating connected basins also reproduces the principal features observed in the empirical solutions. This simple model is largely kinematic. Zonally averaged elevations within a simple closed basin would be nearly in equilibrium with the gravitational potential, except for a constant offset required to conserve mass. With connected basins these offsets are mostly eliminated by interbasin mass flux. Because of rotation, this flux occurs mostly in a narrow boundary layer across the mouth and at the western edge of each basin, and geostrophic balance in this zone supports small residual offsets (and phase shifts) between basins. The simple model predicts that this effect should decrease roughly linearly with frequency, a result that is confirmed by numerical modeling and empirical T/P estimates of the monthly (Mm) tidal constituent. This model also explains some aspects of the anomalous nonisostatic response of the ocean to atmospheric pressure forcing at periods of around 5 days.

Egbert, Gary D.↗
DOE OSTI2017

Distributed Estimation of All the Eigenvalues and Eigenvectors of Matrices Associated With Strongly Connected Digraphs

This letter considers the problem of estimating all the eigenvalues and eigenvectors of an irreducible matrix, corresponding to a strongly connected digraph, in the absence of knowledge on the global network topology. To this end, in this work we propose a unified distributed strategy performed by each node in the network and relies only on the local information. The key idea is to transform the nonlinear problem of computing both the eigenvalues and eigenvectors of an irreducible matrix into a linear one. Specifically, we first transform distributively the irreducible matrix into a nonsingular irreducible matrix. Each node in the network then estimates in a distributed fashion the inverse of the nonsingular matrix by solving a set of linear equations based on a consensus-type algorithm. The eigenvalues and the corresponding eigenvectors are finally computed by exploiting the relations between the eigenvalues and eigenvectors of both the inverse and the original irreducible matrices. A numerical example is provided to demonstrate the effectiveness of the proposed distributed strategy.

14 SOLAR ENERGY↗
DOE OSTI2016

A perturbative analysis of synchrotron spectral index variation over the microwave sky

In this paper, we implement a perturbative approach, first proposed by Bouchet and Gispert, to estimate the variation of the spectral index of galactic polarized synchrotron emission, using a linear combination of simulated Stokes Q polarization maps of selected frequency bands from WMAP and Planck observations on a region of sky dominated by the synchrotron Stokes Q signal. We find that a first order perturbative analysis recovers the input spectral index map well. Along with the spectral index variation map, our method provides a fixed reference index, β-hat {sub 0s}, over the sky portion being analyzed. Using Monte-Carlo simulations, we find that 〈 β-hat {sub 0s}〉=−2.84±0.01, which matches very closely with the position of a peak at β{sub s}(p)=−2.85 of the empirical probability density function of input synchrotron indices obtained from the same sky region. For thermal dust, the mean recovered spectral index 〈 β-hat {sub d}〉=2.00±0.004 from simulations, matches very well with the spatially fixed input thermal dust spectral index β{sub d}=2.00. As accompanying results of the method, we reconstruct cosmic microwave background, thermal dust, and a synchrotron template components with fixed spectral indices over the entire sky region. We use, in our analysis, full pixel–pixel noise covariance matrices of all frequency bands, estimated from the sky region being analyzed. The perturbative technique of this work (1) can build a model with an arbitrary but sufficient degree of accuracy (and precession) as allowed by the data and (2) can produce maximum likelihood estimators for reference indices and templates asymptotically.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2014

Benchmarking a Soil Moisture Data Assimilation System for Agricultural Drought Monitoring

Despite considerable interest in the application of land surface data assimilation systems (LDAS) for agricultural drought applications, relatively little is known about the large-scale performance of such systems and, thus, the optimal methodological approach for implementing them. To address this need, this paper evaluates an LDAS for agricultural drought monitoring by benchmarking individual components of the system (i.e., a satellite soil moisture retrieval algorithm, a soil water balance model and a sequential data assimilation filter) against a series of linear models which perform the same function (i.e., have the same basic inputoutput structure) as the full system component. Benchmarking is based on the calculation of the lagged rank cross-correlation between the normalized difference vegetation index (NDVI) and soil moisture estimates acquired for various components of the system. Lagged soil moistureNDVI correlations obtained using individual LDAS components versus their linear analogs reveal the degree to which non-linearities andor complexities contained within each component actually contribute to the performance of the LDAS system as a whole. Here, a particular system based on surface soil moisture retrievals from the Land Parameter Retrieval Model (LPRM), a two-layer Palmer soil water balance model and an Ensemble Kalman filter (EnKF) is benchmarked. Results suggest significant room for improvement in each component of the system.

soil moisture↗
DOE OSTI2018

Enhancing statistical moment calculations for stochastic Galerkin solutions with Monte Carlo techniques

Here, we provide a method for enhancing stochastic Galerkin moment calculations to the linear elliptic equation with random diffusivity using an ensemble of Monte Carlo solutions. This hybrid approach combines the accuracy of low-order stochastic Galerkin and the computational efficiency of Monte Carlo methods to provide statistical moment estimates which are significantly more accurate than performing each method individually. The hybrid approach involves computing a low-order stochastic Galerkin solution, after which Monte Carlo techniques are used to estimate the residual. We show that the combined stochastic Galerkin solution and residual is superior in both time and accuracy for a one-dimensional test problem and a more computational intensive two-dimensional linear elliptic problem for both the mean and variance quantities.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

A General Method for Estimating Zonal Transmission Interface Limits from Nodal Network Data: Preprint

Capacity expansion models for the electric power system often employ zonal (rather than nodal) resolution, necessitating estimates of aggregate power transfer limits across the interfaces between model zones. Interface limits between planning areas are sometimes published, but they are not generalizable to arbitrary zone shapes. There is thus a need for a reproducible method for estimating interface transfer limits (ITLs) between user-defined zones directly from nodal transmission system data. Here, we present a simple method for estimating ITLs using a DC power flow approximation via the power transfer distribution factor (PTDF) matrix. Linear optimization is performed to identify the distribution of power flows that maximizes the total flow on interface-crossing lines, subject to individual line ratings, limits on bus injection/withdrawal, and the relationships among flows, injections, and withdrawals imposed by the PTDF matrix. We demonstrate the application of the method on a 134-zone ~65000-bus system, and we explore the influence of flow direction, contingency level, and zone size on the estimated ITLs. There is significant heterogeneity in the ratio of the ITL to the sum of interface-crossing line ratings, which highlights the importance of accounting for the physical constraints on power flows imposed by Kirchhoff's laws when estimating zonal ITLs.

17 WIND ENERGY↗
NASA NTRS1974

On the design of optimal input signals in system identification

The problem of designing optimal inputs in the identification of multi-input multi-output linear systems with unknown time-varying parameters is considered using a Bayesian approach. A sensitivity index gives a measure of performance for the closed-loop system inputs. The computation of the optimal closed-loop mappings is shown to be a nontrivial exercise in stochastic control with no analytic solution, but optimal open-loop and affine laws yield much more tractable problems. For time-invariant systems, the sensitivity index considered is shown to be equivalent to the trace of the (strictly positive definite) information matrix associated with the system. Numerical examples are given. A Kalman filter is used to estimate the parameters. A necessary condition for the Kalman filter not to diverge when applying linear feedback is also given.

Lopez-Toledo, A. A.↗
DOE OSTI2015

Decentralizing Network Inference Problems With Multiple-Description Fusion Estimation (MDFE)

Network inference (or tomography) problems, such as traffic matrix estimation or completion and link loss inference, have been studied rigorously in different networking applications. These problems are often posed as under-determined linear inverse (UDLI) problems and solved in a centralized manner, where all the measurements are collected at a central node, which then applies a variety of inference techniques to estimate the attributes of interest. This paper proposes a novel framework for decentralizing these large-scale under-determined network inference problems by intelligently partitioning it into smaller subproblems and solving them independently and in parallel. The resulting estimates, referred to as multiple descriptions, can then be fused together to compute the global estimate. We apply this Multiple Description and Fusion Estimation (MDFE) framework to three classical problems: traffic matrix estimation, traffic matrix completion, and loss inference. Using real topologies and traces, we demonstrate how MDFE can speed up computation while maintaining (even improving) the estimation accuracy and how it enhances robustness against noise and failures. In conclusion, we also show that our MDFE framework is compatible with a variety of existing inference techniques used to solve the UDLI problems.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2012

Constructing an Efficient Self-Tuning Aircraft Engine Model for Control and Health Management Applications

Self-tuning aircraft engine models can be applied for control and health management applications. The self-tuning feature of these models minimizes the mismatch between any given engine and the underlying engineering model describing an engine family. This paper provides details of the construction of a self-tuning engine model centered on a piecewise linear Kalman filter design. Starting from a nonlinear transient aerothermal model, a piecewise linear representation is first extracted. The linearization procedure creates a database of trim vectors and state-space matrices that are subsequently scheduled for interpolation based on engine operating point. A series of steady-state Kalman gains can next be constructed from a reduced-order form of the piecewise linear model. Reduction of the piecewise linear model to an observable dimension with respect to available sensed engine measurements can be achieved using either a subset or an optimal linear combination of "health" parameters, which describe engine performance. The resulting piecewise linear Kalman filter is then implemented for faster-than-real-time processing of sensed engine measurements, generating outputs appropriate for trending engine performance, estimating both measured and unmeasured parameters for control purposes, and performing on-board gas-path fault diagnostics. Computational efficiency is achieved by designing multidimensional interpolation algorithms that exploit the shared scheduling of multiple trim vectors and system matrices. An example application illustrates the accuracy of a self-tuning piecewise linear Kalman filter model when applied to a nonlinear turbofan engine simulation. Additional discussions focus on the issue of transient response accuracy and the advantages of a piecewise linear Kalman filter in the context of validation and verification. The techniques described provide a framework for constructing efficient self-tuning aircraft engine models from complex nonlinear simulations.Self-tuning aircraft engine models can be applied for control and health management applications. The self-tuning feature of these models minimizes the mismatch between any given engine and the underlying engineering model describing an engine family. This paper provides details of the construction of a self-tuning engine model centered on a piecewise linear Kalman filter design. Starting from a nonlinear transient aerothermal model, a piecewise linear representation is first extracted. The linearization procedure creates a database of trim vectors and state-space matrices that are subsequently scheduled for interpolation based on engine operating point. A series of steady-state Kalman gains can next be constructed from a reduced-order form of the piecewise linear model. Reduction of the piecewise linear model to an observable dimension with respect to available sensed engine measurements can be achieved using either a subset or an optimal linear combination of "health" parameters, which describe engine performance. The resulting piecewise linear Kalman filter is then implemented for faster-than-real-time processing of sensed engine measurements, generating outputs appropriate for trending engine performance, estimating both measured and unmeasured parameters for control purposes, and performing on-board gas-path fault diagnostics. Computational efficiency is achieved by designing multidimensional interpolation algorithms that exploit the shared scheduling of multiple trim vectors and system matrices. An example application illustrates the accuracy of a self-tuning piecewise linear Kalman filter model when applied to a nonlinear turbofan engine simulation. Additional discussions focus on the issue of transient response accuracy and the advantages of a piecewise linear Kalman filter in the context of validation and verification. The techniques described provide a framework for constructing efficient self-tuning aircraft engine models from complex nonlinear simulatns.

Armstrong, Jeffrey B.↗
NASA NTRS2004

CDGPS-Based Relative Navigation for Multiple Spacecraft

This thesis investigates the use of Carrier-phase Differential GPS (CDGPS) in relative navigation filters for formation flying spacecraft. This work analyzes the relationship between the Extended Kalman Filter (EKF) design parameters and the resulting estimation accuracies, and in particular, the effect of the process and measurement noises on the semimajor axis error. This analysis clearly demonstrates that CDGPS-based relative navigation Kalman filters yield good estimation performance without satisfying the strong correlation property that previous work had associated with "good" navigation filters. Several examples are presented to show that the Kalman filter can be forced to create solutions with stronger correlations, but these always result in larger semimajor axis errors. These linear and nonlinear simulations also demonstrated the crucial role of the process noise in determining the semimajor axis knowledge. More sophisticated nonlinear models were included to reduce the propagation error in the estimator, but for long time steps and large separations, the EKF, which only uses a linearized covariance propagation, yielded very poor performance. In contrast, the CDGPS-based Unscented Kalman relative navigation Filter (UKF) handled the dynamic and measurement nonlinearities much better and yielded far superior performance than the EKF. The UKF produced good estimates for scenarios with long baselines and time steps for which the EKF would diverge rapidly. A hardware-in-the-loop testbed that is compatible with the Spirent Simulator at NASA GSFC was developed to provide a very flexible and robust capability for demonstrating CDGPS technologies in closed-loop. This extended previous work to implement the decentralized relative navigation algorithms in real time.

Mitchell, Megan Leigh↗
NASA NTRS1987

Deconvolution of spectra for intimate mixtures

Visible to near infrared reflectance spectra of macroscopic mixtures have been shown to be linear combinations of the reflections of the pure mineral components in the mixture. However, for microscopic mixtures the mixing systematics are in general nonlinear. The systematics may be linearized by conversion of reflectance to single scattering albedo (w), where the equations which relate reflectance to w depend on the method of data collection. Several proposed mixing models may be used to estimate mineral abundances from the reflectance spectra of intimate mixtures. These models are summarized and a revised model is presented. A noniterative (linear) least squares approach was used for curve fitting and the data, measured as bi-directional reflectance with incidence and emergence angles of 30 and 0 deg were converted to w by a simplified version of Hapke's equation for bi-directional reflectance. This model was tested with two mixture series composed of 45 to 75 micron particles: an anorthite-enstatite series and an olivine-magnetite series. The data indicate that the simplified Hapke's equation may be used to convolve reflectance spectra into mineral abundances if appropriate endmembers are known or derived from other techniques. For surfaces that contain a significant component of very low albedo material, a somewhat modified version of this technique will need to be developed. Since the abundances are calculated using a noniterative approach, the application of this method is especially efficient for large spectral data sets, such as those produced by mapping spectrometers.

Mustard, John F.↗
NASA NTRS2022

Robust Controller Synthesis for Vision-based Spacecraft Guidance and Control

This work develops a method for Robust Controller for Vision-based Spacecraft (RCVS) guidance and control, integral to the robust autonomy framework for multi-spacecraft for- mation control and reconfiguration applications. The method is built around the use of a photo-realistic simulator, where a camera is deployed on a tracking spacecraft (ego) in order to observe an uncontrolled spacecraft (target) in a Low Earth Orbit (LEO). In this direction, the proposed approach performs the relative state (attitude and position) estimation of the target spacecraft using Convolutional Neural Network (CNN). The state estimation error is then modeled and the corresponding error-bounds are obtained around a nominal trajectory of the ego and target spacecraft. Next, this work proposes a linear matrix inequalities (LMIs) based approach to controller synthesis, guaranteed to be robust against both model uncertainties and measurement errors, resulting from vision-based estimation. This controller is comprised of two distinct components, one synthesized based on the nominal trajectory, while the “robust” component corrects for deviations from the nominal trajectory. Finally, a tracking scenario that directly utilize the image data for spacecraft guidance and control, is presented to showcase the performance of the proposed robust autonomy framework.

Rahmani, Amir↗
NASA NTRS2023

Assessing High Resolution Thermal Monitoring of Complex Intertidal Environments From Space: the Case of ECOSTRESS at Rias Baixas, NW Iberia

Currently, thermal sensors on satellites and other orbital platforms with high spatial resolution designed to monitor the oceans are often insufficient to assess the surface temperature of small bodies of water or within the coastal boundary layer. As the accuracy and sensitivity of remote sensing satellites improve, the demand for more accurate and up-to-date basic data sets for calibration increases. The quality of the thermal data collected by ECOSTRESS on the International Space Station allows the characterization of thermal stress levels in coastal ecosystems with a high spatial resolution of 70 m and a return time from hours to 5 days. This study focused on the calibration of ECOSTRESS estimates with in situ surface temperature data from sensors installed at 3 cm depth in the sediment on the intertidal muddy sandflats of three of the Rías Baixas in Galicia, NW Iberian Peninsula, from 2019 to 2021. A final number of 45 ECOSTRESS temperature images were analyzed. From these, 20% contained substantial georeferencing errors which had to be corrected manually with GIS software tools. We applied the Fourier's law of Heat Conduction to derive surface estimates from loggers sub-surface measurements that could be directly compared with ECOSTRESS data. Overall, a good calibration which explained more than 80% of ECOSTRESS temperature estimates for the whole dataset was obtained, but with an intrinsic cold bias around 1.39 °C. When temperatures at a depth of 3 cm were used, the linear fit became worse and the negative bias increased to 1.49 °C. Closer inspection revealed that night measurements were responsible for this larger deviation, as ECOSTRESS estimates became much colder compared to within sediment measurements because of the combined effect of the instrument intrinsic bias and nocturnal surface cooling. The best calibration was obtained when surface estimates were calculated just for the nighttime, as the cold bias decreased to 0.93 °C. More importantly, during hot daytime conditions in emersion above 20 °C, ECOSTRESS data matched surface temperature estimates, probably because of a better performance of ECOSTRESS algorithm at dry surfaces with lower emissivity. Thus, during the most ecologically relevant periods when high temperatures could drive thermal stress in many commercially-important bivalve species, ECOSTRESS provides accurate surface estimates that can be used to derive sub surface temperatures at the depths at which the different burrowing organisms live. We thus conclude that this instrument constitutes an important global tool to examine thermal stress at an unprecedented spatial scale for complex sea-land boundary systems.

ECOSTRESS↗