Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 577 records · Page 32

Comparison of Pre- and Post-Flight Estimations of Pitch Damping Coefficients for Unguided Entry Vehicles

Introduction: Pre-flight predictions of pitch damping coefficient (𝐶!!+𝐶!"̇)curves are used in entry simulations to determine flight readiness and if any dynamic remediation strategies are required. Accurate characterization of these pitch damping coefficient curves is crucial for mission success. Historically, pre-flight estimations of these curves are derived experimentally either via ballistic range testing or with 1-degree-of-freedom (DOF) free-or forced-oscillation wind tunnel testing. Methodology:With post-flight best-estimated trajectories, a fitting procedure as described in Karlgaard, et. al. [1] is used to compare pre-and post-flight pitch damping coefficients. This technique relies on an assumed analytical form of the trajectory, a limited 3-DOF representation described in Schoenenbergerand Queen [2], which results in anEuler-Cauchy solution form.This solution assumes constant density,thus limited amplitude multi-peak windows of the trajectory are considered to populate the pitch damping-amplitudes pace. Fig. 1 shows preliminary estimations of the pitch damping coefficients estimated using the post-flight best estimated trajectoryfor MER B, using a3-amplitude-peak window. Important to note is that the Euler-Cauchy solution also assumes planar motion, thus coning or 𝛽contribution is not included. Further,the estimated post-flight pitch damping approximations and the pre-flight predictions shown are effective representations of pitch damping, or the integrated effect of the pitch damping coefficient over one pitch cycle at the corresponding cycle’s peak amplitude. Effective pitch damping coefficient curves provide more intuitive insights about where the vehicle is stable and permits readily processed comparisons between data-sets. It is important to note that most pitch damping curves, such as those reported in aerodatabases and utilized in flight mechanics trajectory simulations, are represented in instantaneous angle-of-attack space distinct from what is presented in Fig 1. The extracted post-flight pitch damping coefficient points show a wide range of values across a narrow band of amplitudes (Fig.1a), albeit centered on the published pre-flight approximation. This also applies to the prediction as a function of Mach in Fig1b, where the post-flight scatter points are centered around the pre-flight predictions, but the pre-flight predictions notably lack data aboveMach 5.Further tuning of the technique, such as changing the number of peaks considered in each stencil, resolving density in higher resolution, and employing more accurate fitting procedures (beyond scipy.curve_fit()), will be performed. PosterFocus:This study will apply this reconstruction method to the as-flown trajectories for the following unguided entry vehicles: MER A, MER B, IRVE3, Phoenix, Insight, Pathfinder, and SIAD. The pitch damping coefficient will be derived from these trajectories and compared with the pre-flight predictions. Comparing across different missions, entry environments, and vehicle geometries can inform where our pre-flight testing methodologies accurately capture the flight dynamics. Results will provide valuable insight into the accuracy of pre-flight predicted pitch damping curves. Beyond that, this method can shed light on how accurate these curves need to be, which can inform uncertainty ranges applied to pre-flight curves for other missions.

pitch damping coefficient↗

Automated Machine Learning to Evaluate the Information Content of Tropospheric Trace Gas Columns for Fine Particle Estimates Over India: A Modeling Testbed

India is largely devoid of high-quality and reliable on-the-ground measurements of fine particulate matter (PM 2.5 ). Ground-level PM 2.5 concentrations are estimated from publicly available satellite Aerosol Optical Depth (AOD) products combined with other information. Prior research has largely overlooked the possibility of gaining additional accuracy and insights into the sources of PM using satellite retrievals of tropospheric trace gas columns. We evaluate the information content of tropospheric trace gas columns for PM 2.5 estimates over India within a modeling testbed using an Automated Machine Learning (AutoML) approach, which selects from a menu of different machine learning tools based on the data set. We then quantify the relative information content of tropospheric trace gas columns, AOD, meteorological fields, and emissions for estimating PM 2.5 over four Indian sub-regions on daily and monthly time scales. Our findings suggest that, regardless of the specific machine learning model assumptions, incorporating trace gas modeled columns improves PM 2.5 estimates. We use the ranking scores produced from the AutoML algorithm and Spearman’s rank correlation to infer or link the possible relative importance of primary versus secondary sources of PM 2.5 as a first step toward estimating particle composition. Our comparison of AutoML-derived models to selected baseline machine learning models demonstrates that AutoML is at least as good as user-chosen models. The idealized pseudo-observations (chemical-transport model simulations) used in this work lay the groundwork for applying satellite retrievals of tropospheric trace gases to estimate fine particle concentrations in India and serve to illustrate the promise of AutoML applications in atmospheric and environmental research.

Machine learning↗

Probabilistic Hydrological Estimation of LandSlides (PHELS): Global Ensemble Landslide Hazard Modelling

In this study we present a model for the global Probabilistic Hydrological Estimation of LandSlides (PHELS). PHELS estimates the daily hazard of hydrologically triggered landslides at a coarse spatial resolution of 36 km, by combining landslide susceptibility (LSS) and (percentiles of) hydrological variable(s). The latter include daily rainfall, a 7-day antecedent rainfall index (ARI7) or root-zone soil moisture content (rzmc) as hydrological predictor variables, or the combination of rainfall and rzmc. The hazard estimates with any of these predictor variables have areas under the receiver operating characteristic curve (AUC) above 0.68. The best performance was found with combined rainfall and rzmc predictors (AUC = 0.79), which resulted in the lowest number of missed alarms (especially during spring) and false alarms. Furthermore, PHELS provides hazard uncertainty estimates by generating ensemble simulations based on repeated sampling of LSS and the hydrological predictor variables. The estimated hazard uncertainty follows the behaviour of the input variable uncertainties, is about 13.6 % of the estimated hazard value on average across the globe and in time and is smallest for very low and very high hazard values.

Anne Felsberg↗

Hybrid Flush and Synthetic Air Data Filter for Entry Vehicle Atmospheric State Estimation

A hybrid flush/synthetic air data sensing filter utilizing Kalman-Schmidt and Rach-Tung-Striebel smoothers is developed to obtain entry vehicle atmosphere estimates. The filter/smoother blends information from pressure sensors distributed on the heatshield with measurements of the vehicle aerodynamic forces and moments computed from mass properties and inertial measurement unit data, and prior estimates of the atmosphere. The filter produces estimates of the atmospheric conditions along the entry trajectory, and systematic error estimates to reconcile differences between the pressure and aerodynamic data sources. The filter is applied to data acquired during the Mars Science Laboratory and Mars 2020 entry, descent, and landing at Gale crater and at Jezero crater, respectively. The results show that the hybrid filter produces estimates of the freestream flight condition with lower uncertainty than either the flush or synthetic air data algorithms. The filter accomplishes this result by incorporating additional data and computing estimates of systematic error parameters in the pressure data and the aerodynamic model to further reduce the uncertainties.

Christopher D. Karlgaard↗

Synergy of Satellite Radiation, Precipitation, and Other Meteorological Variable Observations for Global Mean Sea Surface Turbulent Heat Flux Estimation

Sea surface turbulent heat flux is one of the key components in global freshwater and energy balances and plays an important role in atmospheric dynamics, thermodynamics and general circulation. This turbulent heat flux is dominantly decided by sea surface latent heat release with some contribution from sensible heat exchange. The flux and its anomaly could significantly affect ocean heat storage and ocean circulation. They are part of and have great impacts on climate variability. Though the turbulent heat flux is extremely important for the climate and weather systems, only very limited ship and buoy observations of the flux are available over open oceans. There are no global, operational, direct measurements of this crucial meteorological variable. Global estimates are, basically, indirectly calculated from bulk turbulent flux parameterization with a combination of satellite column water vapor, sea surface water temperature, and wind speed observations. Critical parameters such as sea surface air temperature and humidity are estimated from empirical relations of theses variables with the water vapor and sea surface water temperature, respectively. Lacking accurate knowledge on surface air temperature and humidity, uncertainties in the parameterization and potential changes in the non-linear turbulent processes with long-term climate variations could cause large errors in estimated long-term turbulent fluxes from the indirect method as shown in current global sea surface turbulent flux datasets. This study uses synergized data of satellite global radiation, precipitation, and other meteorological variable observations to estimate sea surface turbulent heat fluxes. The radiation observations are made by the satellite Clouds and the Earth’s Radiant Energy System (CERES) sensors, while the global precipitation data is from NASA’s satellite Global Precipitation Climatology Project (GPCP). Other data includes satellite sea surface water temperature and wind speed observations. These datasets are obtained from a wide range of space sensors from passive to active instruments and from visible and near infrared to thermal infrared and microwave spectral sounders. One significant feature of these datasets is that they have multi-decades long climate records. Top-of-atmosphere (TOA) radiation and its anomaly represents the net heat energy input to the climate system, and the oceanic precipitation and ocean-land moisture transport can be used to quantify sea surface latent heat release. Based on the synergized datasets and the principle of global water and energy balances, global mean turbulent heat fluxes are estimated. The turbulent heat anomalies for the first two decades of the 21st century are, then, obtained mainly from global CERES radiation and GPCP precipitation anomalies, along with CERES derived ocean-land heat transports. Analysis indicates that the uncertainties in the estimated turbulent flux anomalies may be reduced considerably. The results suggest a strong needs in synchronized and synergized observations of atmospheric radiation, precipitation, and oceanic meteorological variables for long-term climate studies.

Bing Lin↗

Uncertainties in an Observation-Based Estimate of the Global Aerosol Direct Radiative Effect

Aerosols impact Earth’s radiation budget directly through interactions with radiation and indirectly via interactions with clouds. Aerosols have significant radiative impacts on Earth’s energy budget and represent the largest uncertainty in the radiative forcing of the current climate. Although the magnitude of the aerosol direct radiative effect (DRE) is estimated to be less than that of the indirect effect, uncertainties are large. While model-based estimates of aerosol radiative forcing are subject to uncertainties in modeled aerosol properties and uncertainties in simulated cloud cover and albedo, most observation-based studies of global aerosol DRE have their own difficulties with aerosol optical properties and have been limited to cloud-free skies in most cases. Estimates of the uncertainty of the aerosol DRE are themselves uncertain and have varied widely among different published studies. In this study we use the CERES-CALIPSO-CloudSat-MODIS (C3M) product to estimate the global clear-sky and all-sky aerosol DRE. C3M merges collocated data from CERES, CALIPSO, CloudSat, and MODIS with information from reanalysis and an aerosol transport model. With CALIOP observations of aerosol below optically thin clouds and above clouds, co-located with cloud albedo from MODIS, the C3M dataset allows a detailed exploration of observational uncertainties. By perturbing the aerosol properties used and comparing radiative effects in perturbed cases with the base case we characterize uncertainties in estimated aerosol DRE due to uncertainties in the aerosol properties involved. The results also provide a basis for determining measurement requirements to improve the accuracy of DRE estimates. We will describe the approach and present results.

Dave Winker↗

Estimating the concentration of total suspended solids in inland and coastal waters from Sentinel-2 MSI: A semi-analytical approach

Inland and coastal waters provide key ecosystem services and are closely linked to human well-being. In this study, we propose a semi-analytical method, which can be applied to Sentinel-2 MultiSpectral Instrument (MSI) images to retrieve high spatial-resolution total suspended solids (TSS) concentration in a broad spectrum of aquatic ecosystems ranging from clear to extremely turbid waters. The presented approach has four main steps. First, the remote sensing reflectance ( R rs ) at a band lacking in MSI (620 nm) is estimated through an empirical relationship from R rs at 665 nm. Second, waters are classified into four types (clear, moderately turbid, highly turbid, and extremely turbid). Third, semi-analytical algorithms are used to estimate the particulate backscattering coefficient ( b bp ) at a reference band depending on the water types. Last, TSS is estimated from bbp at the reference band. Validation and comparison of the proposed method with three existing methods are performed using a simulated dataset ( N = 1000), an in situ dataset collected from global inland and coastal waters ( N = 1265) and satellite matchups ( N = 40). Results indicate that the proposed method can improve TSS estimation and provide accurate retrievals of TSS from all three datasets, with a median absolute percentage error (MAPE) of 14.88 %, 31.50 % and 41.69 % respectively. We also present comparisons of TSS mapping between the Sentinel-3 Ocean and Land Colour Instrument (OLCI) and MSI in Lake Kasumigaura, Japan and the Tagus Estuary, Portugal. Results clearly demonstrate the advantages of using MSI for TSS monitoring in small water bodies such as rivers, river mouths and other nearshore waters. MSI can provide more detailed and realistic TSS estimates than OLCI in these water bodies. The proposed TSS estimation method was applied to MSI images to produce TSS time-series in Lake Kasumigaura, which showed good agreements with in situ and OLCI-derived TSS time-series.

Sentinel-2↗

Accounting for Point Estimate Uncertainty in Space Systems Reliability and Risk Analysis

Understanding and accounting for uncertainty in risk analysis is a critical step in the management and communication of risk in engineered systems. The component and system-level analysis to determine the probability of a negative outcome and its consequence is often quantified by a point estimate. Many Program and Enterprise decisions involving technical concerns and issues rely on reliability engineering activities to produce quantified risk analysis to inform the decision making process. At NASA, it is common to use a Probabilistic Risk Analysis (PRA) to inform the overall risk to Loss of Mission or Loss of Crew that involves integration across all spacecraft subsystem fault trees to produce an overall probability of mission failure. The point estimate is an estimate of this overall probability and is an immediate result of a fault tree model. It is the result of a model where the probability of each event is taken to be equal to its mean. The value provides an approximation of the overall mean without running any uncertainty calculations (e.g., no sampling). Using only the point estimate can lead to a false sense of precision and the point estimate may not match the resulting mean when uncertainty is taken into consideration. This paper will explore five conditions that can cause the PRA model mean to diverge from the point estimate and will provide engineers and managers insight into the importance of understanding uncertainty in the elements of PRA models.

Paul J Collier↗

Introducing the GeoRePORT Resource Size Tool: Reporting on Geothermal Resource Size Estimations Using the Geothermal Resource Portfolio Optimization and Reporting Technique (GeoRePORT): Preprint

The Geothermal Resource Portfolio Optimization and Reporting Technique (GeoRePORT) was developed with funding from the U.S. Department of Energy Geothermal Technologies Office to assist in identifying and pursuing long-term investment strategies through the development of a resource reporting protocol. The assessment protocols used in GeoRePORT allow for comparison of project attributes across locations and geological settings to understand the feasibility of geothermal development. This work introduces the Resource Size Tool, a new feature within the GeoRePORT package that compiles two independent methods for estimating geothermal resource size in terms of energy capacity in MW. Energy production potential for twenty-three case studies was estimated with the Resource Size Tool in order to 1) generate a reasonable range of resource size estimates for a particular geothermal field; 2) illustrate the advantages and limitations of each methodology (such as data input requirements, estimate accuracy and precision, and the appropriate circumstances of use); and 3) test the ability of the resource size tool to provide useful and accurate information for geothermal stakeholders. The tool employs two methods widely used in the geothermal industry: (1) USGS Volumetric and (2) Power Density. Results from our case studies show general overlap between these two methods in terms of resource size estimates; however, they also reveal key differences between the two approaches that should be considered when using such estimates to drive development. First, the two methods rely on different input parameters and therefore one method may be more appropriate and/or accurate for a given project than the other. Second, the Power Density method was found to generate wider ranges of resource size predictions, more consistently aligning with actual power production of the field but with larger scales of error; whereas the USGS Volumetric method predicts narrower ranges but tends to overestimate when compared to current MW production. Future work will refine variables used in the methods with input data from other sections of GeoRePORT and modify uncertainty levels based on the particular datasets used for a given project.

geological↗

Spatial Transferability of Machine Learning Based Volume Estimation Models

High-quality traffic volume data is essential for efficient transportation planning and operations. However, such high-quality data is expensive to collect, owing primarily to the high capital cost of installing and maintaining continuous counting stations (CCSs). Recent availability of probe-based vehicle data offers a cost-effective solution for increasing the observability of traffic volumes. However, having ample ground truth traffic data is a prerequisite for developing robust volume estimation models. Though this might not be a big issue in many states, states with scarce CCS data might be able to benefit from robust volume estimation models developed in (adjacent) data-rich states. While there is a reasonable amount of spatial transferability research in the transportation domain, there is a dearth of knowledge on the spatial transferability of probe-based volume estimation models. To address this gap, this paper explores spatial transferability of volume estimation models developed from data in three states (Colorado, North Carolina, and Pennsylvania). Results indicate that it is extremely important to maintain temporal consistency when attempting spatial transferability of volume estimation models. It was also found that models trained on regions with lower peak traffic volumes will limit the performance of models transferred to states with higher peak hourly traffic volumes. Corroborating findings from existing spatial transferability research on other topics, it was found that a meta-model (developed using data from multiple states) performs better than volume estimation models developed within any one of the states.

ADVANCED PROPULSION SYSTEMS↗

Updating Nuclear Energy Cost Estimates for Net Zero World Initiative

Energy modeling of decarbonized scenarios in integrated energy systems requires nuclear energy parameters that are critical for forecasting, modeling and cost structure analysis. Using updated real-world data has always been a challenge to estimate current nuclear reactors costs and deployment scenarios. Given this, an updated set of parameters for overnight capital costs and operation and maintenance costs are estimated for the Net Zero World initiative using recent reports that provided a vast set of open sources data inputs. This paper follows the methodology developed in the Net Zero World report and applies the new ranges estimated in the Gateway for Accelerated Innovation in Nuclear report that address many of the current challenges in obtaining accurate cost data for advanced nuclear concepts. The final goal is to provide new estimates of the overnight capital costs and operational costs for different countries. The present paper improves the earlier capital cost estimations, building on recent literature that aims to obtain accurate data for modeling and simulation to enhance energy system evaluations and support decision-making in areas like de-carbonization and capacity expansion. Finally, the paper compares the new cost estimates with the old cost results.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Ensuring Solution Uniqueness in Three-Phase Power System State Estimation

This paper is concerned with the issue of potential non-unique solutions in three-phase state estimation. Theory of observability analysis for positive sequence power system state estimation is based on certain assumptions that avoid possibility of multiple solutions. Also, it is shown that observability of a positive sequence network remains independent of the network parameters or the operating state. When extending single-phase observability analysis directly to the three-phase case, this paper considers the possibility of converging to multiple solutions, i.e. solution non-uniqueness, even for cases where state estimator successfully converges. The study illustrates via numerical examples the likelihood of converging to entirely different solutions for certain network parameters. It also examines how the operating state, particularly under unbalanced loading, leads to solution non-uniqueness. The paper then describes an alternative approach to ensure a unique solution in three-phase state estimation. This method aims to accurately and uniquely estimate the state of any unbalanced three-phase system, irrespective of load imbalance, network configuration, existence of synchronous generators or transformers.

Power System State Estimation, Three-Phase, Distri↗

Resilient Distributed Field Estimation

We study resilient distributed field estimation under measurement attacks. A network of agents or devices measures a large, spatially distributed physical field parameter. An adversary arbitrarily manipulates the measurements of some of the agents. Each agent's goal is to process its measurements and information received from its neighbors to estimate only a few specific components of the field. Here, we present SAFE, the saturating adaptive field estimator, a consensus + innovations distributed field estimator that is resilient to measurement attacks. Under sufficient conditions on the compromised measurement streams, the physical coupling between the field and the agents' measurements, and the connectivity of the cyber communication network, SAFE guarantees that each agent's estimate converges almost surely to the true value of the components of the parameter in which the agent is interested. Finally, we illustrate the performance of SAFE through numerical examples.

97 MATHEMATICS AND COMPUTING↗

First-Principles Cost Estimation of a Sodium Fast Reactor Nuclear Plant

A multi-tiered cost analysis is performed to estimate full costs of a nuclear power plant (NPP) based on sodium-cooled fast reactor (SFR) technology. To address the lack of fully transparent cost estimations from past undertakings for NPPs, we have developed a detailed and first-principles-based cost estimate for a generalized SFR NPP. Our intent is to achieve a high degree of transparency with our cost assumptions and develop a cost model that is flexible and easily extendable to variations in NPP design and other nuclear reactor types. Furthermore, we strive to achieve a clear organization of costs and complete identification of key cost drivers based on first principles. To this end, the cost results of our analysis as given in Table 24 and Table 25 are organized and categorized into a code of accounts (COA) under development at Idaho National Laboratory (INL). Varying degrees of first-principles methods are employed, such as design for manufacture and assembly® (DFMA® ), to elucidate costs in all process levels of the plant equipment, buildings and site structures, personnel, and other miscellaneous but significant cost elements. These approaches have been successfully applied in past cost analysis projects and are designed to enable rapid and flexible cost estimation. Application of these techniques for evaluating NPP costs is similar in concept to the full, detailed estimation of construction and fabrication costs determined in a later stage of NPP development. Note that our approach tries to avoid use of other past analysis results and data such as those from the legacy Energy Economic Data Base (EEDB) Program, as these resources are based on historical NPP costs and thus may not be indicative of new reactor technologies or construction and fabrication/manufacturing techniques. However, we provide a comparison of our SFR NPP cost results in Table 90 against those included in the EEDB for a representative pressurized water reactor (PWR).

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Estimating the State of Charge in Lithium Primary Batteries: Recent Advances and Critical Insights

Lithium primary batteries (LPBs) remain essential in critical applications such as military, aerospace, medical and emergency devices, and portable electronics. Their superior energy density over lithium-ion batteries offers a significant advantage for long-duration use. Therefore, accurate estimation of the state of charge (SoC) is essential for ensuring the reliable and safe operation of these batteries. While extensive research has been conducted on SoC estimation techniques for lithium-ion secondary batteries, LPBs present unique challenges that complicate accurate SoC estimation. Moreover, research on nondestructive testing techniques for SoC estimation in LPBs is significantly lacking. In this review article, it is aimed to provide a comprehensive overview of recent advancements in SoC estimation for LPBs and generates new insights and directions for future research. Herein, existing methods are discussed and their effectiveness and mechanisms are identified, and areas for further optimization are outlined. More theoretical/experimental efforts to advance SoC detection in LPBs is recommended due to challenges identified with existing techniques.

25 ENERGY STORAGE↗

Parameter estimation with maximal updated densities

A recently developed measure-theoretic framework solves a stochastic inverse problem (SIP) for models where uncertainties in model output data are predominantly due to aleatoric (i.e., irreducible) uncertainties in model inputs (i.e., parameters). The subsequent inferential target is a distribution on parameters. Another type of inverse problem is to quantify uncertainties in estimates of “true” parameter values under the assumption that such uncertainties should be reduced as more data are incorporated into the problem, i.e., the uncertainty is considered epistemic. A major contribution of this work is the formulation and solution of such a parameter identification problem (PIP) within the measure-theoretic framework developed for the SIP. The approach is novel in that it utilizes a solution to a stochastic forward problem (SFP) to update an initial density only in the parameter directions informed by the model output data. In other words, this method performs “selective regularization” only in the parameter directions not informed by data. Further, the solution is defined by a maximal updated density (MUD) point where the updated density defines the measure-theoretic solution to the PIP. Another significant contribution of this work is the full theory of existence and uniqueness of MUD points for linear maps with Gaussian distributions. Data-constructed Quantity of Interest (QoI) maps are also presented and analyzed for solving the PIP within this measure-theoretic framework as a means of reducing uncertainties in the MUD estimate. We conclude with a demonstration of the general applicability of the method on two problems involving either spatial or temporal data for estimating uncertain model parameters. The first problem utilizes spatial data from a stationary partial differential equation to produce a MUD estimate of an uncertain boundary condition. The second problem utilizes temporal data obtained from the state-of-the-art ADvanced CIRCulation (ADCIRC) model to obtain a MUD estimate of uncertain wind drag coefficients for a simulated extreme weather event near the Shinnecock Inlet located in the Outer Barrier of Long Island, NY, USA.

42 ENGINEERING↗

How to estimate soil organic carbon stocks of agricultural fields? perspectives using ex-ante evaluation

Estimating soil organic carbon (SOC) stocks of agricultural fields has a range of important applications from development of sustainable management practices to monitoring carbon stocks. There are many estimation strategies with the potential for more reliable estimates of SOC stock and more efficient use of soil sampling and analysis resources, especially by leveraging readily available auxiliary information such as remote sensing. However, concrete guidance for strategy selection is lacking. This study narrows this gap with a comparison of strategies for estimating deep SOC stock (0–60 cm) in a prototypical field. Using high density SOC stock measurements and simulation, we built on past studies by 1) ex-ante evaluating a large number of strategy options, 2) using a Bayesian approach to quantify the uncertainty of the comparison, and 3) considering multiple Bayesian models to assess sensitivity to this modeling choice. We found that, using readily available auxiliary information, both balanced and stratified sampling offer substantial improvements over simple random sampling. The auxiliary information most important for this improvement is a Sentinel-2 SOC index = blue / (green × red), followed by the topographic wetness index. We found that these results are robust to the choice of mapping method, but that there is uncertainty in the magnitude of improvement. Here, we recommend future studies implement this Bayesian approach for simulated ex-ante evaluation of SOC stock estimation strategies across more fields to investigate the generalizability of these findings.

54 ENVIRONMENTAL SCIENCES↗

Data-Driven Distributed Algorithms for Estimating Eigenvalues and Eigenvectors of Interconnected Dynamical Systems

Here, the paper presents data-driven algorithms to estimate in a distributed manner the eigenvalues, right and left eigenvectors of an unknown linear (or linearized) interconnected dynamic system. In particular, the proposed algorithms do not require the identification of the system model in advance before performing the estimation. As a first step, we consider interconnected dynamical system with distinct eigenvalues. The proposed strategy first estimates the eigenvalues using the well-known Prony method. The right and left eigenvectors are then estimated by solving distributively a set of linear equations. One important feature of the proposed algorithms is that the topology of communication network used to perform the distributed estimation can be chosen arbitrarily, given that it is connected, and is also independent of the structure or sparsity of the system (state) matrix. The proposed distributed algorithms are demonstrated via a numerical example.

97 MATHEMATICS AND COMPUTING↗