Search NASA⌕ Search

SEARCH · Search NASA

Results for “mass balance errors”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Mass balance assessment using GPS

Mass balance is an integral part of any comprehensive glaciological investigation. Unfortunately, it is hard to determine at remote locations where there is no fixed reference. The Global Positioning System (GPS) offers a solution. Simultaneous GPS observations at a known location and the remote field site, processed differentially, will accurately position the camp site. From there, a monument planted in the firn atop the ice can also be accurately positioned. Change in the monument's vertical position is a direct indicator of ice thickness change. Because the monument is not connected to the ice, its motion is due to both mass balance change and to the settling of firn as it densifies into ice. Observations of relative position change between the monument and anchors at various depths within the firn are used to remove the settling effect. An experiment to test this method has begun at Byrd Station on the West Antarctic Ice Sheet and the first epoch of observations was made. Analysis indicates that positioning errors will be very small. It appears likely that the largest errors involved with this technique will arise from ancillary data needed to determine firn settling.

Hulbe, Christina L.↗
NASA NTRS2000

Troughs on Martian Ice Sheets: Analysis of Their Closure and Mass Balance

At the Copenhagen workshop on Martian polar processes, Ralf Greve commented that the flow regime surrounding scarps and troughs of the Martian polar ice sheets cannot be modeled using traditional "plan view" ice-sheet models. Such models are inadequate because they typically use reduced equations that embody certain simplifications applicable only to terrestrial ice sheets where the upper ice sheet surface is smooth. In response to this suggestion, we have constructed a 2-dimensional, time dependent "side view" (two spatial dimensions: one horizontal, one vertical) model of scarp closure that is designed to overcome the difficulties described by Greve. The purpose of the model is to evaluate the scales of stress variation and styles of flow closure so as to estimate errors that may be encountered by "plan view" models. We show that there may be avenues whereby the complications associated with scarp closure can be overcome in "plan view" models through appropriate parameterizations of 3-dimensional effects. Following this, we apply the flow model to simulate the evolution of a typical scarp on the North Polar Cap of Mars. Our simulations investigate: (a) the role of "radiation trapping" (see our companion abstract) in creating and maintaining "spiral-like" scarps on the ice sheet, (b) the consequences of different flowlaws and ice compositions on scarp evolution and, in particular, scarp age, and (c) the role of dust and debris in scarp evolution.

Fountain, A.↗
NASA NTRS2001

Diagnosing Tropical Divergent Circulations From Satellite-Derived Diabatic Forcing and Variational Integral Constraints

Large-scale divergent circulations are part of the atmospheric dynamic response to diabatic heating from condensation, radiative processes, and surface energy fluxes. Vertical motion and the associated divergent wind is thus intimately tied to the hydrologic cycle and the global beat balance. Despite its importance, the divergent circulation is too small in comparison to the rotational flow to measure directly with any accuracy. Vertical motions are recovered diagnostically from reanalyses and, as such, are subject to shortcomings in model physics, numerics, and data availability. While reanalysis estimates of tropical divergent circulations are much improved over those from the Global Weather Experiment, there are still substantial differences between products from operational centers. This is because these circulations are still forced largely by model physics and only secondarily by observations. In order to produce a refined estimate of tropical divergence and its interannual variability we have used a number of remotely-sensed data sets along with variational constraints to improve upon reanalysis estimates. Among these are: precipitation from SSM/I and GPCP, TOVS Path-A vertical cloud distributions; ISCCP radiative cooling rates; TOA radiative fluxes from ERBS, surface radiative fluxes from the SRB project, and surface latent and sensible flux estimates from SSM/I. The TOVS Path-A data constrain the divergent outflow in precipitating regions to have the same vertical structure as observed cloudiness. Using integral constraints for moisture, heat, and mass balance, we retrieve consistent divergent wind flows. We examine the ability of this type analysis to capture regional details of ENSO related perturbations to the divergent wind and associated tropical energy balance. Precipitation from satellite is found to be the major constraint in supplying this horizontal structure. We also consider the ability of this analysis to quantify integrated (land vs. ocean) fluctuations in the global monsoonal flow and energy balance. This depends, of course, on the error characteristics of the data constraints. Some estimate of these properties and their relative importance is provided.

Robertson, Franklin R.↗
DOE OSTI2024

Development of Steady-State and Dynamic Mass and Energy Constrained Neural Networks for Distributed Chemical Systems Using Noisy Transient Data

The paper presents the development of algorithms for mass and energy constrained neural network models that can exactly conserve the overall mass and energy of distributed chemical process systems, even though the noisy transient data used for optimal model training violate the same. In contrast to approximately satisfying mass and energy balance constraints of a system by soft penalization of objective function, algorithms have been developed for solving equality-constrained nonlinear optimization problems, thus providing the guarantee of exactly satisfying the system mass and energy conservation laws. For developing dynamic mass-energy constrained network models for distributed systems, hybrid series and parallel dynamic-static neural networks have been leveraged. The developed algorithms for solving both the training and forward problems are validated using both steady-state and dynamic data in the presence of various noise characteristics. The developed data-driven algorithms are flexible to exactly satisfy mass and energy balance constraints for dynamic chemical processes if the system holdup information is available. The proposed network structures and algorithms are applied to the development of data-driven lumped and distributed models of an adiabatic superheater/reheater system, a nonisothermal continuous stirred tank reactor, as well as an electrically heated plug-flow reactor system where one form of energy gets transformed to another. It has been observed that the mass-energy constrained neural networks yield a root mean squared error of <1% with respect to the system truth for the case studies evaluated in this work.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2010

Application of a Physics-Based Stabilization Criterion to Flight System Thermal Testing

The theory shown here can provide thermal stability criteria based on physics and a goal steady state error rather than on an arbitrary "X% Q/mC(sub P)" method. The ability to accurately predict steady-state temperatures well before thermal balance is reached could be very useful during testing. This holds true for systems where components are changing temperature at different rates, although it works better for the components closest to the sink. However, the application to these test cases shows some significant limitations: This theory quickly falls apart if the thermal control system in question is tightly coupled to a large mass not accounted for in the calculations, so it is more useful in subsystem-level testing than full orbiter tests. Tight couplings to a fluctuating sink causes noise in the steady state temperature predictions.

Baker, Charles↗
NASA NTRS2018

High-resolution Wave Propagation Method for Stratified Flows

The implementation of the multidimensional f-waves Riemann solver for the time-dependent, three-dimensional, nonhydrostatic, meso- and microscale atmospheric flows is described in detail. The Riemann solver employs flux-based wave decomposition (f-waves) for the calculation of Godunov fluxes in which the flux differences are written directly as the linear combination of the right eigenvectors of the hyperbolic system. The scheme incorporates the source term due to gravity without introducing discretization errors which is an important property in the context of atmospheric flows. The resulting flow solver is conservative, accurate, stable, and well-balanced. The implementation of the solver is evaluated using benchmark test cases for atmospheric dynamics.

Riemann problem↗
NASA NTRS2021

Contribution of the World's Main Dust Source Regions to the Global Cycle of Desert Dust

Even though desert dust is the most abundant aerosol by mass in Earth's atmosphere, the relative contributions of the world's major source regions to the global dust cycle remain poorly constrained. This problem hinders accounting for the potentially large impact of regional differences in dust properties on clouds, the Earth's energy balance, and terrestrial and marine biogeochemical cycles. Here, we constrain the contribution of each of the world's main dust source regions to the global dust cycle. We use an analytical framework that integrates an ensemble of global aerosol model simulations with observationally informed constraints on the dust size distribution, extinction efficiency, and regional dust aerosol optical depth (DAOD). We obtain a dataset that constrains the relative contribution of nine major source regions to size-resolved dust emission, atmospheric loading, DAOD, concentration, and deposition flux. We find that the 22–29 Tg (1 standard error range) global loading of dust with a geometric diameter up to 20 µm is partitioned as follows: North African source regions contribute ∼ 50 % (11–15 Tg), Asian source regions contribute ∼ 40 % (8–13 Tg), and North American and Southern Hemisphere regions contribute ∼ 10 % (1.8–3.2 Tg). These results suggest that current models on average overestimate the contribution of North African sources to atmospheric dust loading at ∼ 65 %, while underestimating the contribution of Asian dust at ∼ 30 %. Our results further show that each source region's dust loading peaks in local spring and summer, which is partially driven by increased dust lifetime in those seasons. We also quantify the dust deposition flux to the Amazon rainforest to be ∼ 10 Tg yr−1, which is a factor of 2–3 less than inferred from satellite data by previous work that likely overestimated dust deposition by underestimating the dust mass extinction efficiency. The data obtained in this paper can be used to obtain improved constraints on dust impacts on clouds, climate, biogeochemical cycles, and other parts of the Earth system.

Dust↗
DOE OSTI2024

Dynamic population balance in molecular-level simulations of hypersonic flows

This report summarizes the work towards developing stochastic weighted particle methods (SWPM) for future application in hypersonic flows. Extensive changes to Sandia’s direct simulation Monte Carlo (DSMC) solver, SPARTA (Stochastic Particle Real Time Analyzer), were made to enable the necessary particle splitting and reduction capabilities for SWPM. The results from one-dimensional Couette and Fourier flows suggest that SWPM can reproduce the correct transport for a large range of Knudsen numbers with adequate accuracy. The associated velocity and temperature profiles are in good agreement with DSMC. An issue with particle placement during particle number reduction, is identified, to which, a simple but effective solution based on minimizing the center of mass error is proposed. High Mach wheel flows are simulated using the SWPM and DSMC methods. SWPM is capable of providing nearly an order of magnitude increase in efficiency over DSMC while retaining high accuracy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1987

Estimation of blade airloads from rotor blade bending moments

A method is developed to estimate the blade normal airloads by using measured flap bending moments; that is, the rotor blade is used as a force balance. The blade's rotation is calculated in vacuum modes and the airloads are then expressed as an algebraic sum of the mode shapes, modal amplitudes, mass distribution, and frequency properties. The modal amplitudes are identified from the blade bending moments using the Strain Pattern Analysis Method. The application of the method is examined using simulated flap bending moment data that have been calculated for measured airloads for a full-scale rotor in a wind tunnel. The estimated airloads are compared with the wind tunnel measurements. The effects of the number of measurements, the number of modes, and errors in the measurements and the blade properties are examined, and the method is shown to be robust.

Bousman, William G.↗
NASA NTRS2020

The Effect of Europa and Enceladus Analog Seawater Composition on Isotopic Measurements of Volatile CO2

Science goals for icy ocean worlds missions include characterizing the chemical compositionof the surface and interior using remote sensing and in situ techniques. The next class of flightmass spectrometers for these missions will obtain compositional identifications and isotope ratiomeasurements of volatiles evolved from the ice surface (exosphere) and plumes. These massspectra will be combined with data from other flight instruments to infer the composition of theinterior from volatiles. To ensure accurate interpretation of these measurements, it is critical toverify whether the fundamental assumption that volatiles observed in icy ocean world exospheresor plumes will be a direct reflection of the sub-ice ocean. The present study evaluates whether isotopologues from an initial CO2 gas fractionate by interacting with seawater (brine) of varyingsalt composition and concentration. δ13CCO2 are affected by the pH of the brine and subsequentspeciation of {CO2}, where {CO2} represents the combination of CO2, H2CO3, HCO3-, and CO32-.{CO2} for low pH brines hypothesized for Europa will preferentially speciate as CO2. Analyzedδ13CCO2 for low pH brines are within error of the original δ13CCO2, demonstrating that volatile CO2in a low pH system will be a direct reflection of the original CO2. However, Europa’s radiativeenvironment and rapid depressurization due to plume ejection may impose fractionation effects.In contrast, high pH systems relevant to Enceladus or a more alkaline Europa are expected toform all carbonate species, while favoring speciation as HCO3-. High pH brines in theseexperiments include both an original CO2 gas and an isotopically distinct HCO3- (from NaHCO3salt). These alkaline experiments demonstrate that δ13CCO2 values are highly variable, and dependon the concentration of the dominant carbonate species, (Na)HCO3-. Mass balance estimatesindicate that measured δ13CCO2 is a thermodynamically predictable mixture of both carbonsources, suggesting that measurements of δ13C at Enceladus would directly reflect of the sourcesof CO2 and carbonate buffering in the ocean. δ18O measurements for CO2 interacting with KCland MgCl2 follow established models for δ18O-ionic strength. CO2 interacting with MgSO4,Na2SO4, and NaCl demonstrate an offset from established δ18O-ionic strength models dependingon the concentration of initial CO2. These results suggest that current predictive models for δ18Oin brines need to be resolved for changing concentrations of CO2

Ocean Worlds↗
NASA NTRS2004

Turbofan Engine Simulated in a Graphical Simulation Environment

Recently, there has been an increase in the development of intelligent engine technology with advanced active component control. The computer engine models used in these control studies are component-level models (CLM), models that link individual component models of state space and nonlinear algebraic equations, written in a computer language such as Fortran. The difficulty faced in performing control studies on Fortran-based models is that Fortran is not supported with control design and analysis tools, so there is no means for implementing real-time control. It is desirable to have a simulation environment that is straightforward, has modular graphical components, and allows easy access to health, control, and engine parameters through a graphical user interface. Such a tool should also provide the ability to convert a control design into real-time code, helping to make it an extremely powerful tool in control and diagnostic system development. Simulation time management is shown: Mach number versus time, power level angle versus time, altitude versus time, ambient temperature change versus time, afterburner fuel flow versus time, controller and actuator dynamics, collect initial conditions, CAD output, and component-level model: CLM sensor, CAD input, and model output. The Controls and Dynamics Technologies Branch at the NASA Glenn Research Center has developed and demonstrated a flexible, generic turbofan engine simulation platform that can meet these objectives, known as the Modular Aero-Propulsion System Simulation (MAPSS). MAPSS is a Simulink-based implementation of a Fortran-based, modern high pressure ratio, dual-spool, low-bypass, military-type variable-cycle engine with a digital controller. Simulink (The Mathworks, Natick, MA) is a computer-aided control design and simulation package allows the graphical representation of dynamic systems in a block diagram form. MAPSS is a nonlinear, non-real-time system composed of controller and actuator dynamics (CAD) and component-level model (CLM) modules. The controller in the CAD module emulates the functionality of a digital controller, which has a typical update rate of 50 Hz. The CLM module simulates the dynamics of the engine components and uses an update rate of 2500 Hz, which is needed to iterate to balance mass and energy among system components. The actuators in the CAD module use the same sampling rate as those in the CLM. Two graphs of normalized spool speed versus time in seconds and one graph of normalized average metal temperature versus time in seconds is shown. MAPSS was validated via open-loop and closed-loop comparisons with the Fortran simulation. The preceding plots show the normalized results of a closed-loop comparison looking at three states of the model: low-pressure spool speed, high-pressure spool speed, and the average metal temperature measured from the combustor to the high-pressure turbine. In steady state, the error between the simulations is less than 1 percent. During a transient, the difference between the simulations is due to a correction in MAPSS that prevents the gas flow in the bypass duct inlet from flowing forward instead of toward the aft end, which occurs in the Fortran simulation. A comparison between MAPSS and the Fortran model of the bypass duct inlet flow for power lever angles greater than 35 degrees is shown.

Parker, Khary I.↗
NASA NTRS2021

Data Assimilation of Terrestrial Water Storage Observations to Estimate Precipitation Fluxes: A Synthetic Experiment

The Gravity Recovery and Climate Experiment (GRACE) mission and its Follow-On (GRACE-FO) mission provide unprecedented observations of terrestrial water storage (TWS) dynamics at basin to continental scales. Established GRACE data assimilation techniques directly adjust the simulated water storage components to improve the estimation of groundwater, streamflow, and snow water equivalent. Such techniques artificially add/subtract water to/from prognostic variables, thus upsetting the simulated water balance. To overcome this limitation, we propose and test an alternative assimilation scheme in which precipitation fluxes are adjusted to achieve the desired changes in simulated TWS. Using a synthetic data assimilation experiment, we show that the scheme improves performance skill in precipitation estimates in general, but that it is more robust for snowfall than for rainfall, and it fails in certain regions with strong horizontal gradients in precipitation. The results demonstrate that assimilation of TWS observations can help correct (adjust) the model’s precipitation forcing and, in turn, enhance model estimates of TWS, snow mass, soil moisture, runoff, and evaporation. A key limitation of the approach is the assumption that all errors in TWS originate from errors in precipitation. Nevertheless, the proposed approach produces more consistent improvements in simulated runoff than the established GRACE data assimilation techniques.

Data Assimilation↗
NASA NTRS1999

Elevation Change of the Southern Greenland Ice Sheet from Satellite Radar Altimeter Data

Long-term changes in the thickness of the polar ice sheets are important indicators of climate change. Understanding the contributions to the global water mass balance from the accumulation or ablation of grounded ice in Greenland and Antarctica is considered crucial for determining the source of the about 2 mm/yr sea-level rise in the last century. Though the Antarctic ice sheet is much larger than its northern counterpart, the Greenland ice sheet is more likely to undergo dramatic changes in response to a warming trend. This can be attributed to the warmer Greenland climate, as well as a potential for amplification of a global warming trend in the polar regions of the Northern Hemisphere. In collaboration with Drs. Curt Davis and Craig Kluever of the University of Missouri, we are using data from satellite radar altimeters to measure changes in the elevation of the Southern Greenland ice sheet from 1978 to the present. Difficulties with systematic altimeter measurement errors, particularly in intersatellite comparisons, beset earlier studies of the Greenland ice sheet thickness. We use altimeter data collected contemporaneously over the global ocean to establish a reference for correcting ice-sheet data. In addition, the waveform data from the ice-sheet radar returns are reprocessed to better determine the range from the satellite to the ice surface. At JPL, we are focusing our efforts principally on the reduction of orbit errors and range biases in the measurement systems on the various altimeter missions. Our approach emphasizes global characterization and reduction of the long-period orbit errors and range biases using altimeter data from NASA's Ocean Pathfinder program. Along-track sea-height residuals are sequentially filtered and backwards smoothed, and the radial orbit errors are modeled as sinusoids with a wavelength equal to one revolution of the satellite. The amplitudes of the sinusoids are treated as exponentially-correlated noise processes with a time-constant of six days. Measurement errors (e.g., altimeter range bias) are simultaneously recovered as constant parameters. The corrections derived from the global ocean analysis are then applied over the Greenland ice sheet. The orbit error and measurement bias corrections for different missions are developed in a single framework to enable robust linkage of ice-sheet measurements from 1978 to the present. In 1998, we completed our re-evaluation of the 1978 Seasat and 1985-1989 Geosat Exact Repeat Mission data. The estimates of ice thickness over Southern Greenland (south of 72N and above 2000 m) from 1978 to 1988 show large regional variations (+/-18 cm/yr), but yield an overall rate of +1.5 +/- 0.5 cm/yr (one standard error). Accounting for systematic errors, the estimate may not be significantly different from the null growth rate. The average elevation change from 1978 to 1988 is too small to assess whether the Greenland ice sheet is undergoing a long-term change.

Haines, Bruce J.↗
NASA NTRS2000

A Cost Effective System Design Approach for Critical Space Systems

NASA-JSC required an avionics platform capable of serving a wide range of applications in a cost-effective manner. In part, making the avionics platform cost effective means adhering to open standards and supporting the integration of COTS products with custom products. Inherently, operation in space requires low power, mass, and volume while retaining high performance, reconfigurability, scalability, and upgradability. The Universal Mini-Controller project is based on a modified PC/104-Plus architecture while maintaining full compatibility with standard COTS PC/104 products. The architecture consists of a library of building block modules, which can be mixed and matched to meet a specific application. A set of NASA developed core building blocks, processor card, analog input/output card, and a Mil-Std-1553 card, have been constructed to meet critical functions and unique interfaces. The design for the processor card is based on the PowerPC architecture. This architecture provides an excellent balance between power consumption and performance, and has an upgrade path to the forthcoming radiation hardened PowerPC processor. The processor card, which makes extensive use of surface mount technology, has a 166 MHz PowerPC 603e processor, 32 Mbytes of error detected and corrected RAM, 8 Mbytes of Flash, and I Mbytes of EPROM, on a single PC/104-Plus card. Similar densities have been achieved with the quad channel Mil-Std-1553 card and the analog input/output cards. The power management built into the processor and its peripheral chip allows the power and performance of the system to be adjusted to meet the requirements of the application, allowing another dimension to the flexibility of the Universal Mini-Controller. Unique mechanical packaging allows the Universal Mini-Controller to accommodate standard COTS and custom oversized PC/104-Plus cards. This mechanical packaging also provides thermal management via conductive cooling of COTS boards, which are typically designed for convection cooling methods.

Abbott, Larry Wayne↗
NASA NTRS1998

Laser pointing determination for the geoscience laser altimeter system

The Geoscience Laser Altimeter System (GLAS) is a space-based lidar being developed to monitor changes in the mass balance of the Earth's polar ice sheets (Thomas et al. 1985). GLAS is part of NASA's Earth Observing System (Schutz 1995), and is being designed to launch into a 600 km circular polar orbit in the year 2001, for continuous operation over 3 to 5 years. The orbit's 94 degree inclination has been selected to allow good coverage and profile patterns over the ice sheets of Greenland and Antarctica. The GLAS mission uses a small dedicated spacecraft provided by Ball Aerospace, which is required to have a very stable nadir and zenith pointing platform which points to within approximately 100 urad (20 arcseconds) of Nadir. Accurate knowledge of the laser beam's pointing angle (in the far field) is critical since pointing the laser beam away from nadir biases the altimetry measurements (Gardner 1992, Bufton et al. 1991). This error is a function of the distance of the laser centroid off nadir multiplied by the orbit altitude and the tangent of the slope angle of the terrain. Most of the ice sheet surface slopes are less than 1? resulting in pointing knowledge bias of only 7.6 cm with 7.3 urad accuracy, and overall single shot height accuracy of approximately 15 cm. However, over a 3 deg surface slope pointing knowledge to approximately 7.3 urad is the largest error source (23 cm) in achieving 26 cm height accuracy. The GLAS design incorporates a stellar reference system (SRS) to relate the laser beam pointing angle to the star field to an accuracy of 7.3 urad. The stellar reference system combines an attitude determination system (ADS) operating from 4 to 10 Hz coupled to a 40 Hz laser reference system (LRS) to perform this task.

Miller, Pamela S.↗
NASA NTRS1998

Ice Surface Elevation Changes in East Antarctica from Satellite Altimetry

Estimates of the overall mass balance and seasonal and inter-annual variations in the surface mass balance are obtainable from time-series of ice surface elevations measured by satellite altimetry. Beginning in 2001, NASA's ICESat laser altimeter and lidar mission will significantly improve the range accuracy, the orbit accuracy, and the spatial coverage for measurement of ice sheet elevations (to 86 S) , as compared to previous radar altimeters designed for ocean measurements The radar altimeters on Seasat and Geosat provided ice sheet measurements to 72 S, and on ERS-1 and ERS-2 to 81 S. Although radar altimetry has significant limitations in coverage (due to loss of tracking) and accuracy over sloping surfaces, information on ice-sheet surface-elevation changes has been derived for parts of Antarctica. Recently, the accuracy of the ice measurements by Seasat (3 months of 1978) and Geosat (1985 to 1989) have been improved by new calculations of the satellite orbit heights and other altimeter corrections. Residual orbit errors and inter-satellite biases are evaluated by crossover analysis and by global adjustments to an ocean surface derived from altimeter data. The standard deviation of the orbit error is less than 9 cm, and the long-term trend in the error appears to be less than 1 cm/yr. Orbit errors can be further reduced by adjustment to the ocean surface, but false signals of several cm/yr may be also introduced by the adjustments. These false signals are caused mainly by residual errors in the altimeter corrections over the ocean, and secondary by real changes in the ocean surface elevation. Maps of ice sheet elevation changes north of 72 S are derived from Seasat-Geosat crossovers and from 4.5 years of Geosat crossovers. A notable ice thinning rate of about 50 cm/yr is found at elevations below 2200 meters between 70 and 72 S to the East of the Amery ice shelf, in both the Seasat-Geosat and Geosat-Geosat time intervals Above 2200 meters, to the ridge around 95 E, the thinning is about 10 to 20 cm/yr. For the Seasat/Geosat interval, a thickening is noted on the area East of the ridge between about 95 E and 115 E. Preliminary results from analysis of ERS-1 and ERS-2 radar altimeter data are also described.

Zwally, H. Jay↗
DOE OSTI2025

Baryon fraction from the BAO amplitude: a consistent approach to parameterizing perturbation growth

Galaxy clustering constrains the baryon fraction Omega_b/Omega_m through the amplitude of baryon acoustic oscillations and the suppression of perturbations entering the horizon before recombination. This produces a different pre-recombination distribution of baryons and dark matter. After recombination, the gravitational potential responds to both components in proportion to their mass, allowing robust measurement of the baryon fraction. This is independent of new-physics scenarios altering the recombination background (e.g. Early Dark Energy). The accuracy of such measurements does, however, depend on how baryons and CDM are modeled in the power spectrum. Previous template-based splitting relied on approximate transfer functions that neglected part of information. We present a new method that embeds an extra parameter controlling the balance between baryons and dark matter in the growth terms of the perturbation equations in the CAMB Boltzmann solver. This approach captures the baryonic suppression of CDM prior to recombination, avoids inconsistencies, and yields a clean parametrization of the baryon fraction in the linear power spectrum, separating out the simple physics of growth due to the combined matter potential. We implement this framework in an analysis pipeline using Effective Field Theory of Large-Scale Structure with HOD-informed priors and validate it against noiseless LCDM and EDE cosmologies with DESI-like errors. The new scheme achieves comparable precision to previous splitting while reducing systematic biases, providing a more robust way to baryon-fraction measurements. In combination with BBN constraints on the baryon density and Alcock-Paczynski estimates of the matter density, these results strengthen the use of baryon fraction measurements to derive a Hubble constant from energy densities, with future DESI and Euclid data expected to deliver competitive constraints.

Crespi, Andrea [U. Waterloo (main); Waterloo U., I↗
NASA NTRS2023

Sensitivity of the MAR Regional Climate Model Snowpack to the Parameterization of the Assimilation of Satellite-Derived Wet-Snow Masks on the Antarctic Peninsula

Both regional climate models (RCMs) and remote sensing (RS) data are essential tools in understanding the response of polar regions to climate change. RCMs can simulate how certain climate variables, such as surface melt, runoff and snowfall, are likely to change in response to different climate scenarios but are subject to biases and errors. RS data can assist in reducing and quantifying model uncertainties by providing indirect observations of the modeled variables on the present climate. In this work, we improve on an existing scheme to assimilate RS wet snow occurrence data with the “Modèle Atmosphérique Régional” (MAR) RCM and investigate the sensitivity of the RCM to the parameters of the scheme. The assimilation is performed by nudging the MAR snowpack temperature to match the presence of liquid water observed by satellites. The sensitivity of the assimilation method is tested by modifying parameters such as the depth to which the MAR snowpack is warmed or cooled, the quantity of water required to qualify a MAR pixel as “wet” (0.1 % or 0.2 % of the snowpack mass being water), and assimilating different RS datasets. Data assimilation is carried out on the Antarctic Peninsula for the 2019–2021 period. The results show an increase in meltwater production (+66.7 % on average, or +95 Gt), along with a small decrease in surface mass balance (SMB) (−4.5 % on average, or −20 Gt) for the 2019–2020 melt season after assimilation. The model is sensitive to the tested parameters, albeit with varying orders of magnitude. The prescribed warming depth has a larger impact on the resulting surface melt production than the liquid water content (LWC) threshold due to strong refreezing occurring within the top layers of the snowpack. The values tested for the LWC threshold are lower than the LWC for typical melt days (approximately 1.2 %) and impact results mainly at the beginning and end of the melting period. The assimilation method will allow for the estimation of uncertainty in MAR meltwater production and will enable the identification of potential issues in modeling near-surface snowpack processes, paving the way for more accurate simulations of snow processes in model projections.

regional climate models↗