Search NASA⌕ Search

SEARCH · Search NASA

Results for “lower bounds”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Solving the duck curve in a smart grid environment using a non-cooperative game theory and dynamic pricing profiles

With the intermittency that comes with electricity generation from renewables, utilizing dynamic pricing will encourage the demand-side to respond in a smart way that would minimize the electricity costs and flatten the net electricity demand curve. Determining the optimal dynamic pricing profile that would leverage distributed storage to flatten the curve is a novel idea that needs to be studied. Moreover, the economic feasibility of utilizing distributed electrical energy storage is still not given in the literature. Therefore, in this paper, a novel way of solving a citywide dynamic model using a bilevel programming algorithm is introduced. The problem is developed as a novel non-cooperative Stackelberg game that utilizes air-conditioning systems and electrical storage through the end-users to determine the optimal dynamic pricing profile. The results show that the combined effect of utilizing demand-side air-conditioning systems and distributed storage together can flatten the curve while employing the optimal dynamic pricing profile. An economic study is performed to determine the economic feasibility of 20 different cases with different battery designs and the level of solar penetration. Three metrics were used to evaluate the economic performance of each case: the levelized cost of storage, the levelized cost of energy, and the simple payback period. Most cases had levelized cost of storage values lower than 0.457 $/kWh, which is the lower bound available in the literature. Seven out of 16 cases have a simple payback period shorter than the lifetime of the system (25 years). The case with a 100 MW PV power plant and a battery storage of size 597 MWh, was found to be the most promising case with a simple payback period of 12.71 years for the photovoltaic plant and 19.86 years for the demand-side investments.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2023

Granger Causality for prediction in Dynamic Mode Decomposition: Application to power systems

Here, the dynamic mode decomposition (DMD) technique extracts the dominant modes characterizing the innate dynamical behavior of the system within the measurement data. For appropriate identification of dominant modes from the measurement data, the DMD algorithm necessitates ensuring the quality of the input measurement data sequences. On that account, for validating the usability of the dataset for the DMD algorithm, the paper proposed two conditions: Persistence of excitation (PE) and the Granger Causality Test (GCT). The virtual data sequences are designed with the hankel matrix representation such that the dimensions of the subspace spanning the essential system modes are increased with the addition of new state variables. The PE condition provides the lower bound for the trajectory length, and the GCT provides the order of the model. Satisfying the PE condition enables estimating an approximate linear model, but the predictability with the identified model is only assured with the temporal causation among data searched with GCT. The proposed methodology is validated with the application for coherency identification (CI) in a multi-machine power system (MMPS), an essential phenomenon in transient stability analysis. The significance of PE condition and GCT is demonstrated through various case studies implemented on 22 bus six generator system.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Multi-fidelity Bayesian neural networks: Algorithms and applications

Here we propose a new class of Bayesian neural networks (BNNs) that can be trained using noisy data of variable fidelity, and we apply them to learn function approximations as well as to solve inverse problems based on partial differential equations (PDEs). These multi-fidelity BNNs consist of three neural networks: The first is a fully connected neural network, which is trained following the maximum a posteriori probability (MAP) method to fit the low-fidelity data; the second is a Bayesian neural network employed to capture the cross-correlation with uncertainty quantification between the low- and high-fidelity data; and the last one is the physics-informed neural network, which encodes the physical laws described by PDEs. For the training of the last two neural networks, we first employ the mean-field variational inference (VI) to maximize the evidence lower bound (ELBO) to obtain informative prior distributions for the hyperparameters in the BNNs, and subsequently we use the Hamiltonian Monte Carlo (HMC) method to estimate accurately the posterior distributions for the corresponding hyperparameters. We demonstrate the accuracy of the present method using synthetic data as well as real measurements. Specifically, we first approximate a one- and four-dimensional function, and then infer the reaction rates in one- and two-dimensional diffusion-reaction systems. Moreover, we infer the sea surface temperature (SST) in the Massachusetts and Cape Cod Bays using satellite images and in-situ measurements. Taken together, our results demonstrate that the present method can capture both linear and nonlinear correlation between the low- and high-fidelity data adaptively, identify unknown parameters in PDEs, and quantify uncertainties in predictions, given a few scattered noisy high-fidelity data. Finally, we demonstrate that we can effectively and efficiently reduce the uncertainties and hence enhance the prediction accuracy with an active learning approach, using as examples a specific one-dimensional function approximation and an inverse PDE problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Scalable semidefinite programming approach to variational embedding for quantum many-body problems

In quantum embedding theories, a quantum many-body system is divided into localized clusters of sites which are treated with an accurate ‘high-level’ theory and glued together self-consistently by a less accurate ‘low-level’ theory at the global scale. The recently introduced variational embedding approach for quantum many-body problems combines the insights of semidefinite relaxation and quantum embedding theory to provide a lower bound on the ground-state energy that improves as the cluster size is increased. The variational embedding method is formulated as a semidefinite program (SDP), which can suffer from poor computational scaling when treated with black-box solvers. Here, we exploit the interpretation of this SDP as an embedding method to develop an algorithm which alternates parallelizable local updates of the high-level quantities with updates that enforce the low-level global constraints. Moreover, we show how translation invariance in lattice systems can be exploited to reduce the complexity of projecting a key matrix to the positive semidefinite cone.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Determination of single crystal thermal expansion in Uranium-6wt%Niobium shape memory alloy using in-situ diffraction and modeling of textured polycrystalline samples

In-situ neutron diffraction and ex-situ dilatometry experiments were performed on a chemically banded, quenched uranium-6wt% niobium shape memory alloy to study the impact of deformation-(detwinning-) induced texture on the evolution of thermally induced strains (and associated stresses). Thermal heating and cooling cycles between 5K and 473K were performed in-situ with neutron diffraction, and the lattice strain evolution is reported for the observed monoclinic α" phase. Comparisons between the measured diffraction strains and macroscopic dilatometry experiments reveal relationships between micro- and macro-level thermal expansion. Softening of the texture during heating suggests that twin boundary motion can accommodate the large internal thermal strains which are approximately 2% greater than that typically observed during the 175K heating interval used to age more randomly oriented polycrystalline material. Assuming weak constraint of neighboring grains in randomly textured polycrystals, the single crystal thermal expansion tensor is extracted from measurements of lattice strains over the range 5K to 473K. Predictions of polycrystalline thermal expansion, based upon this single crystal thermal expansion tensor, are shown to compare favorably with bulk thermal expansion observations of the as quenched microstructure. However, such a lower-bound estimate is insufficient to explain all aspects of the behavior of the textured material, where the matrix is not isotropic. In conclusion, it is hypothesized that relaxation processes which occur within the quenched microstructure during heating are responsible for the distinct thermal expansion behavior observed during the initial heating cycle compared to cooling and subsequent cycling.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Mapping strain and structural heterogeneities around bubbles in amorphous ionically conductive Bi 2 O 3

While amorphous materials are often approximated to have a statistically homogeneous atomic structure, they frequently exhibit localized structural heterogeneity that challenges simplified models. This study uses 4D scanning transmission electron microscopy to investigate the strain and structural modifications around gas bubbles in amorphous Bi 2 O 3 induced by argon irradiation. We present a method for determining strain fields surrounding bubbles that can be used to measure the internal pressure of the gas. Compressive strain is observed around the cavities, with higher-order crystalline symmetries emerging near the cavity interfaces, suggesting paracrystalline ordering as a result of bubble coarsening. This ordering, along with a compressive strain gradient, indicates that gas bubbles induce significant localized changes in atomic packing. By analyzing strain fields with maximum compressive strains of 3%, we estimate a lower bound on the internal pressure of the bubbles at 2.5 GPa. These findings provide insight into the complex structural behavior of amorphous materials under stress, particularly in systems with gas inclusions, and offer new methods for probing the local atomic structure in disordered materials. Although considering structural heterogeneity in amorphous systems is non-trivial, these features have crucial impacts on material functionalities, such as mechanical strength, ionic conductivity, and electronic mobility.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Assessing high fidelity multi-component models to facilitate safeguards at Gas Centrifuge Enrichment Plants

We report that the International Atomic Energy Agency (IAEA) inspectors routinely carry out environmental sampling (ES) as a verification method. Collection of environmental swipe samples at various locations in Gas Centrifuge Enrichment Facilities (GCEPs) is an important process in detecting misuse of a declared facility and possibly the existence of undeclared nuclear material. These samples are measured for isotopic composition in uranium containing particles by Thermal Ionization Mass Spectrometry (TIMS) or Inductively Coupled Plasma Mass Spectrometry (ICP-MS). Even though ES is highly effective in detecting the absolute value of enrichments and their deviations from the declared values, it cannot explain the cause of those changes. Several potential explanations can serve as possibilities for particles detected above or other than the declared enrichment. These include normal and non-malicious events such as the design of the enrichment cascades, unintentional failures of the machines, or overshoot during startup of a cascade. It can also be the result of deliberate misuse by the facility operators. The primary objective of this work is to understand how these factors affect the enrichments produced by a cascade and quantify anticipated multi-isotopic concentrations for each case. The following methodology is employed to determine signatures at a particular facility. 1) Utilize a new two-dimensional multi-component diffusion code to obtain centrifuge performance data and use that information to design and perform cascade analysis. Compare and contrast the results with previous 1-D radially averaged solutions from the Pancake code. 2) Design a GCEP cascade with the production goal of 19.75% 235 U for each set of machine data above. Investigate two cascade scenarios that include enrichment of natural uranium (NU) feed to 19.75% 235U in a single cascade compared to a two-step process of NU to 5% and then 5% to 19.75%. 3) Simulate the intentional vs. unintentional off-normal scenarios in the cascades to assess the differences in isotopic concentrations. A non-ideal squared-off cascade model developed at the University of Virginia is used to calculate flow rates and isotopic concentrations of the process gas. The analysis is performed using the Rome machine model operated at 600 m/s rotor speed. The upper and lower bounds of normal and abnormal enrichments in a typical facility are used in conjunction with ES results to understand the root causes of such observations.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2025

Evaluating design safety margins in the American Society of Mechanical Engineers graphite core components design-by-analysis assessments

Graphite is an important material being used for core components in next-generation high-temperature gas-cooled nuclear reactors. The selection of graphite grade for a specific Designer is a complex task, dependent on reactor conditions, component functionality, and required reliability. The American Society of Mechanical Engineers (ASME) provides two semi-probabilistic design-by-analysis assessments to evaluate graphite core components against design reliability targets. The simplified assessment uses a 2-parameter Weibull distribution to describe the graphite grade’s tensile-strength distribution to establish component stress limits. The full assessment uses the 3-parameter Weibull distribution and a modified Weakest-Link Theory approach to calculate a component design probability of failure. The paper defines recommended assessment rules, which are the as-written simplified assessment and the full assessment with parameter lower bounds, the modulus update with threshold reduction, and the 2027 grouping rules. Code rules are applied to three grades: 2114, IG-110, and NBG-18. The baseline margin calculation is developed using the experimental tensile dogbone specimen. Percent margin is defined as the percent reduction in the median experimental load to obtain the allowable load per ASME assessments. Under the recommended rules, the SRC–1 margin in the simplified assessment ranged from 40.2 % to 52.7 % among the grades in this study and from 36.1 % to 49.8 % in the full assessment. The full assessment only decreases the margins by 2.5–4.5 % for the SRC-1 components and 0–1.5 % for the SRC-2 components for this baseline case. Margin is inversely related to material median strength (i.e., the strongest grade, 2114, has the lowest margin).

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Exact hopfion vortices in a 3D Heisenberg ferromagnet

Here, we find exact static soliton solutions for the unit spin vector field of an inhomogeneous, anisotropic three-dimensional Heisenberg ferromagnet. Each soliton is labeled by two integers n and m. It is a (modified) skyrmion in the z = 0 plane with winding number n, which twists out of the plane m times in the z-direction to become a 3D soliton. Here m arises due to the periodic boundary condition at the z-boundaries. We use Whitehead’s integral expression to find that the Hopf invariant of the soliton is an integer H = nm. It represents a hopfion vortex. Plots of the preimages of this topological soliton show that they are either unknots or nontrivial knots, depending on n and m. Any pair of preimage curves links H times, corroborating the interpretation of H as a linking number. We also calculate the exact energy of the hopfion vortex, and show that its topological lower bound has a sublinear dependence on H. Using Derrick’s scaling analysis, we demonstrate that the presence of a spatial inhomogeneity in the anisotropic interaction, which in turn introduces a characteristic length scale in the system, leads to the stability of the hopfion vortex.

36 MATERIALS SCIENCE↗
DOE OSTI2023

System-size dependence of the charged-particle pseudorapidity density at $\sqrt{s_{NN}}$ = 5.02 TeV for pp, p–Pb, and Pb–Pb collisions

We present the first systematic comparison of the charged-particle pseudorapidity densities for three widely different collision systems, pp, p–Pb, and Pb–Pb, at the top energy of the Large Hadron Collider ($\sqrt{s_{NN}}$ = 5.02 TeV) measured over a wide pseudorapidity range (-3.5 < η < 5), the widest possible among the four experiments at that facility. The systematic uncertainties are minimised since the measurements are recorded by the same experimental apparatus (ALICE). The distributions for p–Pb and Pb–Pb collisions are determined as a function of the centrality of the collisions, while results from pp collisions are reported for inelastic events with at least one charged particle at midrapidity. The charged-particle pseudorapidity densities are, under simple and robust assumptions, transformed to charged-particle rapidity densities. This allows for the calculation and the presentation of the evolution of the width of the rapidity distributions and of a lower bound on the Bjorken energy density, as a function of the number of participants in all three collision systems. We find a decreasing width of the particle production, and roughly a smooth ten fold increase in the energy density, as the system size grows, which is consistent with a gradually higher dense phase of matter

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Combination of searches for pair-produced leptoquarks at s = 13 TeV with the ATLAS detector

A statistical combination of various searches for pair-produced leptoquarks is presented, using the full LHC Run 2 (2015–2018) data set of 139 fb -1 collected with the ATLAS detector from proton–proton collisions at a centre-of-mass energy of $\sqrt{s}$ = 13 TeV. All possible decays of the leptoquarks into quarks of the third generation and charged or neutral leptons of any generation are investigated. Since no significant deviations from the Standard Model expectation are observed in any of the individual analyses, combined exclusion limits are set on the production cross-sections for scalar and vector leptoquarks. The resulting lower bounds on leptoquark masses exceed those from the individual analyses by up to 100 GeV, depending on the signal hypothesis.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Graphene Oxide Nanoribbons for High Early-Strength Cement Concrete

Longitudinal oxidative unzipping of the outer walls of multiwalled carbon nanotubes (MWCNTs) yields graphene oxide nanoribbons (GONRs), which exhibit greater open surface area and functional edge content than MWCNTs. This paper presents a study of the nano-amendment of Portland cement concrete with GONRs in concentrations between 0.05% (in weight of cement, wt%) and 0.0005 wt%, thus up to two orders of magnitude lower than that reported as lower-bound in the archival literature. The dispersibility in aqueous solution as a function of GONR concentration and oxygen weight content (O%) was assessed through dynamic light scattering (DLS) and zeta potential analysis. The results indicated that less effective suspensions were obtained for 0.05 wt% of GONRs and 22.7 O%. Therefore, GONR water suspensions with 30-40% O% were used to manufacture 50 mm × 100 mm cylindrical concrete specimens. After 7 days of curing, results from uniaxial compression tests using four specimens per configuration (MWCNT concentration and O%) showed that the incorporation of GONRs resulted in an average increase in compressive strength up to 45%. Consistent with the DLS and compression test results, SEM micrographs showed well-dispersed GONRs together with accelerated and preferential formation of calcium silicate hydrates (C-S-H) for all GONR concentrations. The results indicate, for the first time, that the incorporation of very small concentrations (as low as 0.0005 wt%) of well-dispersed GONR amendments can significantly enhance the early-age concrete strength. However, such enhancement became insignificant after 28 days of curing.

cement↗
DOE OSTI2023

Cloud condensation nuclei activity of internally mixed particle populations at a remote marine free troposphere site in the North Atlantic Ocean

Here, this study reports results from research conducted at the Observatory of Mount Pico (OMP), 2225 m above mean sea level on Pico Island in the Azores archipelago in June and July 2017. We investigated the chemical composition, mixing state, and cloud condensation nuclei (CCN) activities of long-range transported free tropospheric (FT) particles. FLEXible PARTicle Lagrangian particle dispersion model (FLEXPART) simulations reveal that most air masses that arrived at the OMP during the sampling period originated in North America and were highly aged (average plume age > 10 days). We probed size-resolved chemical composition, mixing state, and hygroscopicity parameter ($\mathcal{κ}$) of individual particles using computer-controlled scanning electron microscopy with an energy-dispersive X-ray spectrometer (CCSEM-EDX). Based on the estimated individual particle mass from elemental composition, we calculated the mixing state index, $χ$. During our study, FT particle populations were internally mixed ($χ$ of samples are between 53% and 87%), owing to the long atmospheric aging time. We used data from a miniature Cloud Condensation Nucleus Counter (miniCCNC) to derive the hygroscopicity parameter, $\mathcal{κ}$ CCNC . Combining $\mathcal{κ}$ CCNC and FLEXPART, we found that air masses recirculated above the North Atlantic Ocean with lower mean altitude had higher $\mathcal{κ}$ CCNC due to the higher contribution of sea salt particles. We used CCSEM-EDX and phase state measurements to predict single-particle $\mathcal{κ}$ ($\mathcal{κ}$ CCSEM-EDX ) values, which overlap with the lower range of $\mathcal{κ}$ CCNC measured below 0.15 % SS. Therefore, CCSEM-EDX measurements can be useful in predicting the lower bound of $\mathcal{κ}$, which can be used in climate models to predict CCN activities, especially in remote locations where online CCN measurements are unavailable.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Sensitivity-based voltage constraints for optimal power flow in low-voltage distribution feeders

The optimal power flow (OPF) problem for distribution systems can include network details down to the low-voltage (LV) points of interconnection of individual customers. This paper addresses the implementation of voltage magnitude constraints, and sets forth a practicable approach for capturing the effects on voltage from the switching behavior of loads (e.g., heat pumps, air conditioners, water heaters, or pool pumps) and from the variability of renewable generation (e.g., rooftop solar). The proposed method adjusts the OPF voltage constraints based on forecasts of load and generation upper and lower bounds, in conjunction with sensitivity factors derived from the power flow equations. An illustrative OPF formulation is also provided, which incorporates transformer models that include core loss. We demonstrate that accurate modeling of these LV network components is critical to avoid voltage violations at customer points of interconnection. Furthermore, the ideas are validated through numerical case studies on a realistic distribution feeder.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

Thermomechanical analysis and modeling of a high-temperature light trapping planar cavity solar receiver

The development of durable particle-based high-temperature solar receivers is critical for advancing concentrating solar-thermal (CST) technologies to enable high-efficiency power generation and industrial process heat. Here, this study presents a computational framework to evaluate the thermomechanical performance of a proposed enclosed light-trapping planar cavity receiver designed for particle-based thermal energy systems. The receiver incorporates absorptive cavities and fluidized particle-bed channels to enhance heat capture and reduce thermal losses. Finite element analysis (FEA) is employed to assess stress, strain, and creep-fatigue behavior under concentrated solar flux using realistic thermal boundary conditions derived from coupled system models and experimental assembly parameters. The analysis investigates the influence of particle-to-wall heat transfer coefficients (HTC) ranging from 800 to 1800 W/m 2 .K on the thermomechanical response of six candidate high-temperature alloys: Alloy 740H, Alloy 282, Alloy 617, 316H, Alloy 230, and 800H. Results show that increasing HTC reduces thermal gradients, leading to lower stresses and strains and extended minimum predicted creep life. While all materials satisfy fatigue life requirements under the investigated conditions, significant differences in creep resistance are observed. Alloy 740H consistently exhibits the longest minimum predicted creep life and the most favorable durability margins, followed by Alloy 282, with the remaining materials showing reduced creep resistance under identical loading. The reported creep lives are conservative lower-bound estimates intended for comparative material evaluation. This framework highlights the critical roles of material selection and geometry optimization in improving mechanical durability and reliability of solar-thermal receivers, forming a foundation for future experimental validation and design optimization.

14 SOLAR ENERGY↗
DOE OSTI2022

Joint routing of conventional and range-extended electric vehicles in a large metropolitan network

Range-extended electric vehicles combine the higher efficiency and environmental benefits of battery-powered electric motors with the longer mileage and autonomy of conventional internal combustion engines. This combination is particularly advantageous for time-constrained delivery routing in dense urban areas, where battery recharging along routes can be too time-consuming to economically justify the use of all-electric vehicles. However, switching from electric to conventional fossil fuel modes also results in higher costs and emissions and lower efficiency. This paper analyzes this heterogeneous vehicle routing problem and describes two solution methods: an exact branch-price-and-cut algorithm and an iterated tabu search metaheuristic. From a methodological perspective, we find that the exact algorithm consistently obtains tight lower bounds that also serve to certify the metaheuristic solutions as near-optimal. From a policy standpoint, we examine a large-scale real-world case study concerning parcel deliveries in the Chicago metropolitan area and quantify various operational metrics including energy costs and vehicle miles traveled. We find that by deploying roughly 20% of range -extended vehicles with a modest all-electric range of 33 miles, parcel distributors can save energy costs by up to 17% while incurring less than 0.5% increase in vehicle miles traveled. Increasing the range to 60 miles further reduces costs by only 4%, which can alternatively be achieved by decreasing the average service time by 1 minute or increasing driver working time by 1 hour. In conclusion, our study reveals several key areas of improvement on which vehicle manufacturers, distributors, and policy makers can focus their attention.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

Effects of anisotropic shale creep on the stress and permeability evolution of a geological nuclear waste repository

To ensure long-term safety and performance, geological nuclear waste repositories require low-permeability barriers such as bentonite buffers and/or shale host rock. Shale is not only known for its low permeability but also for its trend to undergo time-dependent deformation (i.e., creep), which could heal damage, but the effects of shale creep on the long-term performance of nuclear waste repositories have not been clearly understood. In particular, the anisotropic nature of shale (i.e., bedding) could have a significant effect on its creep behavior, and consequently, on the long-term performance of nuclear waste repositories. In this research, numerical simulations were carried out with the objective of showing the effects of anisotropic shale creep on the stress and permeability evolution of a generic geological nuclear waste repository in shale. The TOUGH-FLAC simulator, a thermo-hydromechanically (THM) coupled numerical code, was used for the simulations. To achieve the objective, comparisons were performed between the results of anisotropic shale creep simulations and those of different simulation cases, namely, no creep (i.e., elastic), isotropic creep, and long-term creep shale cases. Results of the comparisons show that the elastic and isotropic creep shale cases respectively led to the overestimation and underestimation of stress and permeability in the repository, whereas the long-term creep shale case, which accumulated greater creep in later periods than in earlier periods, helped to keep large shear and tensile stresses from developing while maintaining compressive spherical stress, resulting in consistently low permeability levels. These results indicate that performance assessments with elastic and isotropic creep formation models will provide the upper and lower bound estimates of stress and permeability, while more reasonable estimates will be given by an anisotropic creep formation model, and that shale with long-term creep characteristics will be beneficial in many aspects of the safety and performance of nuclear waste repositories.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Deriving iceberg ablation rates using an on-iceberg autonomous phase-sensitive radar (ApRES)

Abstract The increase in iceberg discharge into the polar oceans highlights the importance of understanding how quickly icebergs are deteriorating and where the resulting freshwater injection is occurring. Recent advances in quantifying iceberg deterioration through combinations of modeling, remote sensing and direct in situ measurements have successfully calculated overall ablation rates, and surface and sidewall ablation; however, in situ measurements of basal melt rates have been difficult to obtain. Radar has successfully measured iceberg thickness, but repeat measurements, which would capture a change in iceberg thickness with time, have not yet been collected. Here we test the applicability of using an on-iceberg autonomous phase-sensitive radar (ApRES) to quantify basal ablation rates of a large (~800 m long) non-tabular Arctic iceberg during an intensive 2019 summer field campaign in Sermilik Fjord, southeast Greenland. We find that ApRES can be used to measure basal ablation even over a short deployment period (10 d), and also provide a lower bound on sidewall melt. This study fills a critical gap in iceberg research and pushes the limits of field instrumentation.

54 ENVIRONMENTAL SCIENCES↗