Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “network inference”

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,333 records · Page 74

The Lusi seismic experiment: An initial study to understand the effect of seismic activity to Lusi

The spectacular Lumpur Sidoarjo (Lusi) eruption started in northeast Java on the 29 of May 2006 following a M6.3 earthquake striking the island [1,2]. Initially, several gas and mud eruption sites appeared along the reactivated strike-slip Watukosek fault system [3] and within weeks several villages were submerged by boiling mud. The most prominent eruption site was named Lusi. The Lusi seismic experiment is a project aims to begin a detailed study of seismicity around the Lusi area. In this initial phase we deploy 30 seismometers strategically distributed in the area around Lusi and along the Watukosek fault zone that stretches between Lusi and the Arjuno Welirang (AW) complex. The purpose of the initial monitoring is to conduct a preliminary seismic campaign aiming to identify the occurrence and the location of local seismic events in east Java particularly beneath Lusi.This network will locate small event that may not be captured by the existing BMKG network. It will be crucial to design the second phase of the seismic experiment that will consist of a local earthquake tomography of the Lusi-AW region and spatial and temporal variations of vp/vs ratios. The goal of this study is to understand how the seismicity occurring along the Sunda subduction zone affects to the behavior of the Lusi eruption. Our study will also provide a large dataset for a qualitative analysis of earthquake triggering studies, earthquake-volcano and earthquake-earthquake interactions. In this study, we will extract Green’s functions from ambient seismic noise data in order to image the shallow subsurface structure beneath LUSI area. The waveform cross-correlation technique will be apply to all of recordings of ambient seismic noise at 30 seismographic stations around the LUSI area. We use the dispersive behaviour of the retrieved Rayleigh waves to infer velocity structures in the shallow subsurface.

58 GEOSCIENCES↗

Using GPS to Detect Imminent Tsunamis

A promising method of detecting imminent tsunamis and estimating their destructive potential involves the use of Global Positioning System (GPS) data in addition to seismic data. Application of the method is expected to increase the reliability of global tsunami-warning systems, making it possible to save lives while reducing the incidence of false alarms. Tsunamis kill people every year. The 2004 Indian Ocean tsunami killed about 230,000 people. The magnitude of an earthquake is not always a reliable indication of the destructive potential of a tsunami. The 2004 Indian Ocean quake generated a huge tsunami, while the 2005 Nias (Indonesia) quake did not, even though both were initially estimated to be of the similar magnitude. Between 2005 and 2007, five false tsunami alarms were issued worldwide. Such alarms result in negative societal and economic effects. GPS stations can detect ground motions of earthquakes in real time, as frequently as every few seconds. In the present method, the epicenter of an earthquake is located by use of data from seismometers, then data from coastal GPS stations near the epicenter are used to infer sea-floor displacements that precede a tsunami. The displacement data are used in conjunction with local topographical data and an advanced theory to quantify the destructive potential of a tsunami on a new tsunami scale, based on the GPS-derived tsunami energy, much like the Richter Scale used for earthquakes. An important element of the derivation of the advanced theory was recognition that horizontal sea-floor motions contribute much more to generation of tsunamis than previously believed. The method produces a reliable estimate of the destructive potential of a tsunami within minutes typically, well before the tsunami reaches coastal areas. The viability of the method was demonstrated in computational tests in which the method yielded accurate representations of three historical tsunamis for which well-documented ground-motion measurements were available. Development of a global tsunami-warning system utilizing an expanded network of coastal GPS stations was under consideration at the time of reporting the information for this article.

Song, Y. Tony↗

Predicting bulge to total luminosity ratio of galaxies using deep learning

ABSTRACT We present a deep learning model to predict the r-band bulge-to-total luminosity ratio (B/T) of nearby galaxies using their multiband JPEG images alone. Our Convolutional Neural Network (CNN) based regression model is trained on a large sample of galaxies with reliable decomposition into the bulge and disc components. The existing approaches to estimate the B/T ratio use galaxy light-profile modelling to find the best fit. This method is computationally expensive, prohibitively so for large samples of galaxies, and requires a significant amount of human intervention. Machine learning models have the potential to overcome these shortcomings. In our CNN model, for a test set of 20 000 galaxies, 85.7 per cent of the predicted B/T values have absolute error (AE) less than 0.1. We see further improvement to 87.5 per cent if, while testing, we only consider brighter galaxies (with r-band apparent magnitude <17) with no bright neighbours. Our model estimates the B/T ratio for the 20 000 test galaxies in less than a minute. This is a significant improvement in inference time from the conventional fitting pipelines, which manage around 2–3 estimates per minute. Thus, the proposed machine learning approach could potentially save a tremendous amount of time, effort, and computational resources while predicting B/T reliably, particularly in the era of next-generation sky surveys such as the Legacy Survey of Space and Time (LSST) and the Euclid sky survey which will produce extremely large samples of galaxies.

Grover, Harsh (ORCID:0000000321338142)↗

Properties of the Binary Black Hole Merger GW150914

On September 14, 2015, the Laser Interferometer Gravitational-Wave Observatory (LIGO) detected a gravitational-wave transient (GW150914); we characterize the properties of the source and its parameters. The data around the time of the event were analyzed coherently across the LIGO network using a suite of accurate waveform models that describe gravitational waves from a compact binary system in general relativity. GW150914 was produced by a nearly equal mass binary black hole of masses 36(+5/-4) solar mass and 29(+4/-4) solar mass; for each parameter we report the median value and the range of the 90% credible interval. The dimensionless spin magnitude of the more massive black hole is bound to be less than 0.7 (at 90% probability). The luminosity distance to the source is 410(+160/-180) Mpc, corresponding to a redshift 0.09(+0.03/−0.04) assuming standard cosmology. The source location is constrained to an annulus section of 610 sq deg, primarily in the southern hemisphere. The binary merges into a black hole of mass 62(+4/−4) solar mass and spin 0.67(+0.05/−0.07). This black hole is significantly more massive than any other inferred from electromagnetic observations in the stellar-mass regime.

Abbott, B. P.↗

Simple spatial scaling rules behind complex cities

Although most of wealth and innovation have been the result of human interaction and cooperation, we are not yet able to quantitatively predict the spatial distributions of three main elements of cities: population, roads, and socioeconomic interactions. By a simple model mainly based on spatial attraction and matching growth mechanisms, we reveal that the spatial scaling rules of these three elements are in a consistent framework, which allows us to use any single observation to infer the others. All numerical and theoretical results are consistent with empirical data from ten representative cities. In addition, our model can also provide a general explanation of the origins of the universal super- and sub-linear aggregate scaling laws and accurately predict kilometre-level socioeconomic activity. Our work opens a new avenue for uncovering the evolution of cities in terms of the interplay among urban elements, and it has a broad range of applications.

97 MATHEMATICS AND COMPUTING↗

Leveraging Hydropower Multi-Sensor Data for Inference and Age-Informed Modeling

Increased demand of operational flexibility such as faster ramp up/down in generation, and more frequent start/stops are putting hydropower plants and their associated components in unprecedented stress. Consequently, these plants are at the high risk of extended and more frequent outage to accommodate unscheduled, and unexpected maintenance. Therefore, hydropower plants are in critical need of data driven and age-informed analysis for their regular and unscheduled operation. Yet not all hydropower plants are exhaustively equipped with sensors and/or measurement streams for their respective components – demanding solutions on how to detect, identify, and locate the cause of any event from the unobservable. Idaho National Laboratory (INL) analyzed the anonymized measurements and event records from the Hydropower Research Institute (HRI) to address this issue, as part of the Water Power Technologies Office (WPTO) funded one year multi-lab project. First, we investigated how time series of multiple sensor measurements can be leveraged to identify an event “root cause” as well as to develop an inference (i.e., estimate the unobservable) problem. INL also investigated how individual hydropower components’ reaction or response times vary across the pre-event, during event, and post-event conditions – enabling the hydropower dynamic models to be age-informed. Finally, the impact of clustering multi-sensor time series on short-term vibration prediction is analyzed. INL will present key findings from these analyses and recommend next steps for stakeholder adoption.

13 HYDRO ENERGY↗

Statistical characterization of experimental magnetized liner inertial fusion stagnation images using deep-learning-based fuel–background segmentation

Significant variety is observed in spherical crystal x-ray imager (SCXI) data for the stagnated fuel–liner system created in Magnetized Liner Inertial Fusion (MagLIF) experiments conducted at the Sandia National Laboratories Z-facility. As a result, image analysis tasks involving, e.g., region-of-interest selection (i.e. segmentation), background subtraction and image registration have generally required tedious manual treatment leading to increased risk of irreproducibility, lack of uncertainty quantification and smaller-scale studies using only a fraction of available data. We present a convolutional neural network (CNN)-based pipeline to automate much of the image processing workflow. This tool enabled batch preprocessing of an ensemble of N scans = 139 SCXI images across N exp = 67 different experiments for subsequent study. The pipeline begins by segmenting images into the stagnated fuel and background using a CNN trained on synthetic images generated from a geometric model of a physical three-dimensional plasma. The resulting segmentation allows for a rules-based registration. Our approach flexibly handles rarely occurring artifacts through minimal user input and avoids the need for extensive hand labelling and augmentation of our experimental dataset that would be needed to train an end-to-end pipeline. Here we also fit background pixels using low-degree polynomials, and perform a statistical assessment of the background and noise properties over the entire image database. Our results provide a guide for choices made in statistical inference models using stagnation image data and can be applied in the generation of synthetic datasets with realistic choices of noise statistics and background models used for machine learning tasks in MagLIF data analysis. We anticipate that the method may be readily extended to automate other MagLIF stagnation imaging applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Oceanic Storm Characteristics Off the Kennedy Space Center Coast

Natural cloud-to-ground lightning may behave differently depending on the characteristics of the attachment mediums, including the peak current (inferred from radiation fields) and the number of ground strike locations per flash. Existing literature has raised issues over the yea"rs on the behavior of lightning over ocean terrain and these phenomena are not yet well understood. To investigate lightning characteristics over differing terrain we will obtain identical observations over adjacent land and ocean regions during both clear air and thunderstorm periods comparing the electric field behavior over these various terrains. For this, a 3-meter NOAA buoy moored 20NM off the coast of the Kennedy Space Center was instrumented with an electric field mill and New Mexico Tech's slow antenna to measure the electric fields aloft and compared to the existing on-shore electric field mill suite of 31 sensors and a coastal slow antenna. New Mexico Tech's Lightning Mapping Array and the Eastern Range Cloud-to-Ground Lightning Surveillance System, along with the network of high-speed cameras being used to capture cloud-to-ground lightning strikes over the terrain regions to identify a valid data set and verify the electric fields. This is an on-going project with the potential for significant impact on the determination of lightning risk to objects on the ground. This presentation will provide results and instrumentation progress to date.

Wilson, J.↗

MERIDIONAL FLOW IN THE SOLAR CONVECTION ZONE. II. HELIOSEISMIC INVERSIONS OF GONG DATA

Meridional flow is thought to play a very important role in the dynamics of the solar convection zone; however, because of its relatively small amplitude, precisely measuring it poses a significant challenge. Here we present a complete time–distance helioseismic analysis of about 2 years of ground-based Global Oscillation Network Group (GONG) Doppler data to retrieve the meridional circulation profile for modest latitudes in an attempt to corroborate results from other studies. We use an empirical correction to the travel times due to an unknown center-to-limb systematic effect. The helioseismic inversion procedure is first tested and reasonably validated on artificial data from a large-scale numerical simulation followed by a test to broadly recover the solar differential rotation found from global seismology. From GONG data, we measure poleward photospheric flows at all latitudes with properties that are comparable with earlier studies and a shallow equatorward flow about 65 Mm beneath the surface, in agreement with recent findings from Helioseismic and Magnetic Imager (HMI) data. No strong evidence of multiple circulation cells in depth or latitude is found, yet the whole phase space has not yet been explored. Tests of mass flux conservation are then carried out on the inferred GONG and HMI flows and compared to a fiducial numerical baseline from models, and we find that the continuity equation is poorly satisfied. While the two disparate data sets do give similar results for about the outer 15% of the interior radius, the total inverted circulation pattern appears to be unphysical in terms of mass conservation when interpreted over modest time scales. We can likely attribute this to both the influence of realization noise and subtle effects in the data and measurement procedure.

79 ASTRONOMY AND ASTROPHYSICS↗

Subgrid modelling for two-dimensional turbulence using neural networks

Here in this investigation, a data-driven turbulence closure framework is introduced and deployed for the subgrid modelling of Kraichnan turbulence. The novelty of the proposed method lies in the fact that snapshots from high-fidelity numerical data are used to inform artificial neural networks for predicting the turbulence source term through localized grid-resolved information. In particular, our proposed methodology successfully establishes a map between inputs given by stencils of the vorticity and the streamfunction along with information from two well-known eddy-viscosity kernels. Through this we predict the subgrid vorticity forcing in a temporally and spatially dynamic fashion. Our study is both a priori and a posteriori in nature. In the former, we present an extensive hyper-parameter optimization analysis in addition to learning quantification through probability-density-function-based validation of subgrid predictions. In the latter, we analyse the performance of our framework for flow evolution in a classical decaying two-dimensional turbulence test case in the presence of errors related to temporal and spatial discretization. Statistical assessments in the form of angle-averaged kinetic energy spectra demonstrate the promise of the proposed methodology for subgrid quantity inference. In addition, it is also observed that some measure of a posteriori error must be considered during optimal model selection for greater accuracy. The results in this article thus represent a promising development in the formalization of a framework for generation of heuristic-free turbulence closures from data.

42 ENGINEERING↗

Dixie Valley Engineered Geothermal System Exploration Methodology Project, Baseline Conceptual Model Report

The Engineered Geothermal System (EGS) Exploration Methodology Project is developing an exploration approach for EGS through the integration of geoscientific data. The Project chose the Dixie Valley Geothermal System in Nevada as a field laboratory site for methodology calibration purposes because, in the public domain, it is a highly characterized geothermal system in the Basin and Range with a considerable amount of geoscience and most importantly, well data. The overall project area is 2500km2 with the Calibration Area (Dixie Valley Geothermal Wellfield) being about 170km2. The project was subdivided into five tasks (1) collect and assess the existing public domain geoscience data; (2) design and populate a GIS database; (3) develop a baseline (existing data) geothermal conceptual model, evaluate geostatistical relationships, and generate baseline, coupled EGS favorability/trust maps from +1km above sea level (asl) to -4km asl for the Calibration Area at 0.5km intervals to identify EGS drilling targets at a scale of 5km x 5km; (4) collect new geophysical and geochemical data, and (5) repeat Task 3 for the enhanced (baseline + new ) data. Favorability maps were based on the integrated assessment of the three critical EGS exploration parameters of interest: rock type, temperature and stress. A complimentary trust map was generated to compliment the favorability maps to graphically illustrate the cumulative confidence in the data used in the favorability mapping. The Final Scientific Report (FSR) is submitted in two parts with Part I describing the results of project Tasks 1 through 3 and Part II covering the results of project Tasks 4 through 5 plus answering nine questions posed in the proposal for the overall project. FSR Part I presents (1) an assessment of the readily available public domain data and some proprietary data provided by Terra-Gen Power, LLC, (2) a re-interpretation of these data as required, (3) an exploratory geostatistical data analysis, (4) the baseline geothermal conceptual model, and (5) the EGS favorability/trust mapping. The conceptual model presented applies to both the hydrothermal system and EGS in the Dixie Valley region. FSR Part II presents (1) 278 new gravity stations; (2) enhanced gravity-magnetic modeling; (3) 42 new ambient seismic noise survey stations; (4) an integration of the new seismic noise data with a regional seismic network; (5) a new methodology and approach to interpret this data; (5) a novel method to predict rock type and temperature based on the newly interpreted data; (6) 70 new magnetotelluric (MT) stations; (7) an integrated interpretation of the enhanced MT data set; (8) the results of a 308 station soil CO2 gas survey; (9) new conductive thermal modeling in the project area; (10) new convective modeling in the Calibration Area; (11) pseudo-convective modeling in the Calibration Area; (12) enhanced data implications and qualitative geoscience correlations at three scales (a) Regional, (b) Project, and (c) Calibration Area; (13) quantitative geostatistical exploratory data analysis; and (14) responses to nine questions posed in the proposal for this investigation. Enhanced favorability/trust maps were not generated because there was not a sufficient amount of new, fully-vetted (see below) rock type, temperature, and stress data. The enhanced seismic data did generate a new method to infer rock type and temperature. However, in the opinion of the Principal Investigator for this project, this new methodology needs to be tested and evaluated at other sites in the Basin and Range before it is used to generate the referenced maps. As in the baseline conceptual model, the enhanced findings can be applied to both the hydrothermal system and EGS in the Dixie Valley region.

15 GEOTHERMAL ENERGY↗

Dixie Valley Engineered Geothermal System Exploration Methodology Project, Baseline Conceptual Model Report

The Engineered Geothermal System (EGS) Exploration Methodology Project is developing an exploration approach for EGS through the integration of geoscientific data. The overall project area is 2500km2 with the Calibration Area (Dixie Valley Geothermal Wellfield) being about 170km2. The Final Scientific Report (FSR) is submitted in two parts (I and II). FSR part I presents (1) an assessment of the readily available public domain data and some proprietary data provided by terra-gen power, llc, (2) a re-interpretation of these data as required, (3) an exploratory geostatistical data analysis, (4) the baseline geothermal conceptual model, and (5) the EGS favorability/trust mapping. The conceptual model presented applies to both the hydrothermal system and EGS in the Dixie Valley region. FSR Part II presents (1) 278 new gravity stations; (2) enhanced gravity-magnetic modeling; (3) 42 new ambient seismic noise survey stations; (4) an integration of the new seismic noise data with a regional seismic network; (5) a new methodology and approach to interpret this data; (5) a novel method to predict rock type and temperature based on the newly interpreted data; (6) 70 new magnetotelluric (MT) stations; (7) an integrated interpretation of the enhanced MT data set; (8) the results of a 308 station soil CO2 gas survey; (9) new conductive thermal modeling in the project area; (10) new convective modeling in the Calibration Area; (11) pseudo-convective modeling in the Calibration Area; (12) enhanced data implications and qualitative geoscience correlations at three scales (a) Regional, (b) Project, and (c) Calibration Area; (13) quantitative geostatistical exploratory data analysis; and (14) responses to nine questions posed in the proposal for this investigation. Enhanced favorability/trust maps were not generated because there was not a sufficient amount of new, fully-vetted (see below) rock type, temperature, and stress data. The enhanced seismic data did generate a new method to infer rock type and temperature (However, in the opinion of the Principal Investigator for this project, this new methodology needs to be tested and evaluated at other sites in the Basin and Range before it is used to generate the referenced maps. As in the baseline conceptual model, the enhanced findings can be applied to both the hydrothermal system and EGS in the Dixie Valley region).

15 GEOTHERMAL ENERGY↗

GW190814: Gravitational Waves from the Coalescence of a 23 Solar Mass Black Hole with a 2.6 Solar Mass Compact Object

We report the observation of a compact binary coalescence involving a 22.2–24.3 Me black hole and a compact object with a mass of 2.50–2.67Me (all measurements quoted at the 90% credible level). The gravitational-wavesignal, GW190814, was observed during LIGO’s and Virgo’s third observing run on 2019 August 14 at 21:10:39 UTC and has a signal-to-noise ratio of 25 in the three-detector network. The source was localized to18.5 deg2 at a distance of-+2414541Mpc; no electromagnetic counterpart has been confirmed to date. The source has the most unequal mass ratio yet measured with gravitational waves,-+0.1120.0090.008, and its secondary component is either the lightest black hole or the heaviest neutron star ever discovered in a double compact-object system. The dimensionless spin of the primary black hole is tightly constrained to 0.07. Tests of general relativity reveal no measurable deviations from the theory, and its prediction of higher-multipole emission is confirmed at high confidence. We estimate a merger rate density of 1–23 Gpc−3yr−1 for the new class of binary coalescence sources that GW190814 represents. Astrophysical models predict that binaries with mass ratios similar to this event can form through several channels, but are unlikely to have formed in globular clusters. However, the combination of mass ratio, component masses, and the inferred merger rate for this event challenges all current models of the formation and mass distribution of compact-object binaries.

Gravitational wave astronomy↗

Lithium-ion Battery Strain Gauge Monitoring and Depth of Discharge Estimation

Lithium-ion battery diagnostics and prognostics rely on measurements of electrical impedance, capacity, and voltage to infer the internal state of the battery. Mechanical changes to the cell structure represent an additional measure of the battery’s state because these changes are related to the overall battery health. As lithium-ion batteries are charged and discharged, lithium ions are inserted or removed from the anode and cathode, a process called intercalation and de-intercalation. As lithium ions intercalate and de-intercalate, they can cause changes to the lattice of the electrode particles, resulting in volumetric changes. These volumetric changes cause mechanical stresses and strains on the lithium-ion battery electrodes, and subsequently, the whole cell’s thickness varies as it is charged and discharged. This paper presents a study on the use of surface-mounted strain gauges for in-situ measurement of structural changes to lithium-ion batteries, along with a characterization of the unit-to-unit differences in strain response. A neural network modeling structure is then used to predict the battery’s depth of discharge under dynamic discharge conditions.

Lithium-ion batteries, depth of discharge, strain ↗

The Global Methane Budget 2000-2012

The global methane (CH4) budget is becoming an increasingly important component for managing realistic pathways to mitigate climate change. This relevance, due to a shorter atmospheric lifetime and a stronger warming potential than carbon dioxide, is challenged by the still unexplained changes of atmospheric CH4 over the past decade. Emissions and concentrations of CH4 are continuing to increase, making CH4 the second most important human-induced greenhouse gas after carbon dioxide. Two major difficulties in reducing uncertainties come from the large variety of diffusive CH4 sources that overlap geographically, and from the destruction of CH4 by the very short-lived hydroxyl radical (OH). To address these difficulties, we have established a consortium of multi-disciplinary scientists under the umbrella of the Global Carbon Project to synthesize and stimulate research on the methane cycle, and producing regular (approximately biennial) updates of the global methane budget. This consortium includes atmospheric physicists and chemists, biogeochemists of surface and marine emissions, and socio-economists who study anthropogenic emissions. Following Kirschke et al. (2013), we propose here the first version of a living review paper that integrates results of top-down studies (exploiting atmospheric observations within an atmospheric inverse-modeling framework) and bottom-up models, inventories and data-driven approaches (including process-based models for estimating land surface emissions and atmospheric chemistry, and inventories for anthropogenic emissions, data-driven extrapolations).For the 2003-2012 decade, global methane emissions are estimated by top-down inversions at 558 TgCH4 yr(exp -1), range 540-568. About 60 of global emissions are anthropogenic (range 50-65%). Since 2010, the bottom-up global emission inventories have been closer to methane emissions in the most carbon-intensive Representative Concentrations Pathway (RCP8.5) and higher than all other RCP scenarios. Bottom-up approaches suggest larger global emissions (736 TgCH4 yr(exp -1), range 596-884) mostly because of larger natural emissions from individual sources such as inland waters, natural wetlands and geological sources. Considering the atmospheric constraints on the top-down budget, it is likely that some of the individual emissions reported by the bottom-up approaches are overestimated, leading to too large global emissions. Latitudinal data from top-down emissions indicate a predominance of tropical emissions (approximately 64% of the global budget, less than 30deg N) as compared to mid (approximately 32%, 30-60deg N) and high northern latitudes (approximately 4%, 60-90deg N). Top-down inversions consistently infer lower emissions in China (approximately 58 TgCH4 yr(exp -1), range 51-72, minus14% ) and higher emissions in Africa (86 TgCH4 yr(exp -1), range 73-108, plus 19% ) than bottom-up values used as prior estimates. Overall, uncertainties for anthropogenic emissions appear smaller than those from natural sources, and the uncertainties on source categories appear larger for top-down inversions than for bottom-up inventories and models. The most important source of uncertainty on the methane budget is attributable to emissions from wetland and other inland waters. We show that the wetland extent could contribute 30-40% on the estimated range for wetland emissions. Other priorities for improving the methane budget include the following: (i) the development of process-based models for inland-water emissions, (ii) the intensification of methane observations at local scale (flux measurements) to constrain bottom-up land surface models, and at regional scale (surface networks and satellites) to constrain top-down inversions, (iii) improvements in the estimation of atmospheric loss by OH, and (iv) improvements of the transport models integrated in top-down inversions. The data presented here can be downloaded from the Carbon Dioxide Information Analysis Center (http://doi.org/10.3334/CDIAC/GLOBAL_ METHANE_BUDGET_2016_V1.1) and the Global Carbon Project.

Saunois, Marielle↗

Data-driven models of nonautonomous systems

Nonautonomous dynamical systems are characterized by time-dependent inputs, which complicates the discovery of predictive models describing the spatiotemporal evolution of the state variables of quantities of interest from their temporal snapshots. When dynamic mode decomposition (DMD) is used to infer a linear model, this difficulty manifests itself in the need to approximate the time-dependent Koopman operators. Our approach is to approximate the original nonautonomous system with a modified system derived via a local parameterization of the time-dependent inputs. The modified system comprises a sequence of local parametric systems, which are subsequently approximated by a parametric surrogate model using the DRIPS (dimension reduction and interpolation in parameter space) framework. The offline step of DRIPS relies on DMD to build a linear surrogate model, endowed with reduced-order bases for the observables mapped from training data. The online step interpolates on suitable manifolds to construct a sequence of iterative parametric surrogate models; the target/test parameter points on these manifolds are specified by a local parameterization of the test time-dependent inputs. Here, we use numerical experimentation to demonstrate the robustness of our method and compare its performance with that of deep neural networks.

97 MATHEMATICS AND COMPUTING↗

A data-driven global soil heterotrophic respiration dataset and the drivers of its inter-annual variability

Soil heterotrophic respiration (SHR), one of the primary carbon fluxes from terrestrial ecosystems to the atmosphere, is important for carbon-climate feedbacks because of its sensitivity to available litter and soil carbon, climatic conditions, and nutrient availability. However, until recently limited SHR data were available, and most published global SHR estimates have either a short time span, coarse spatial resolution, or reply on overly-simple model formulations. To better understand and quantify the global distribution of SHR and its sensitivity to climate variability, we produced a new global SHR dataset using Random Forest algorithms, up-scaling 455 point data from the Global Soil Respiration Database (SRDB 4.0) with gridded fields of climatic, edaphic and productivity as explanatory variables. We estimated a global total SHR of 46.8 Pg C yr-1 over 1985-2013 (95% confidence interval: 38.6-56.3 Pg C yr-1), with a significant increasing trend of 0.03 Pg C yr-2 during this period. We found that the choice of soil moisture datasets contributes more to the difference among these data-driven SHR members rather than that of productivity, temperature and precipitation data sources. We also analyzed the influence of climatic variables on the inter-annual variability (IAV) of our SHR product. Water availability was the dominant driver of IAV at global scales, although the inferred sensitivity depends on the choice of the soil moisture gridded dataset. At the ecosystem scale, temperature strongly controls the IAV of SHR in tropical forests, while water availability dominates in extra-tropical forest and semi-arid regions. Our machine-learning gridded SHR dataset and outputs from process-based land surface models (TRENDYv6) show agreement for a strong association between water variability and SHR IAV at the global scale, but the two approaches lead to different temporal trend globally and different controlling variables for IAV at the ecosystem scale. Our study provides evidence for the pervasive and important role of water availability in driving SHR, indicating both a direct effect limiting decomposition rates and an indirect effect through the amount of fresh organic matter made available to SHR from productivity. In consideration of potential limitations and uncertainties remaining in our data-driven SHR datasets, we call for a more scientifically designed observation network for SHR, more observation data compilation, and increased use of deep learning methods making maximum use of observation data in hand. This will benefit process-based models, and improve our understanding of SHR response to future anomalous environmental conditions.

Yao, Yitong↗

Recent Radar Observations of the Sub-Centimeter Orbital Debris Environment

The NASA Orbital Debris Program Office (ODPO) has conducted radar observations of the orbital debris environment since the early 1990’s to provide measurement data that supports orbital debris models and risk mitigation activities in support of NASA mission objectives. Orbital debris radar observations are a unique mode for radar operation, employing a fixed beam configuration to statistically sample the environment. An advantage of conducting operations in this fashion is that it enables observations of smaller classes of orbital debris than would otherwise be available from the same sensor operating in a traditional tracking mode. Orbital debris-mode radar observations are used to fill in the gaps, which exist in the currently available data from the Space Surveillance Network (SSN), on small size orbital debris populations that represent significant risk to NASA programs. These gaps have typically covered orbital debris with characteristic sizes less than approximately 10 cm down to approximately 3 mm in low Earth orbit (LEO) – depending upon the altitude and sensor configuration. The value of orbital debris radar measurements lies in the ability to extract partial orbital element information about orbital debris in the centimeter to several millimeter size regimes in low Earth orbit – which are not available from other measurement sources. This paper will discuss observations of this smaller class of orbital debris observed in recent years from the radars at the MIT Haystack Observatory in Westford, Massachusetts, and the Goldstone Solar System Radar near Barstow, California. The former radar is able to observe orbital debris down to approximately 5 mm, and the latter, orbital debris with characteristic sizes near 3 mm – at altitudes less than 1000 km. The characteristics and inferences about the current LEO orbital debris environment, and the different subpopulations that are identifiable in the observations are highlighted.

Kennedy, Timothy↗