Search NASA⌕ Search

SEARCH · Search NASA

Results for “wave forecast”

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

A Baseline for Global Weather and Climate Simulations at 1 km Resolution

In an attempt to advance the understanding of the Earth's weather and climate by representing deep convection explicitly, we present a global, four-month simulation (November 2018 to February 2019) with ECMWF's hydrostatic Integrated Forecasting System (IFS) at an average grid spacing of 1.4 km. The impact of explicitly simulating deep convection on the atmospheric circulation and its variability is assessed by comparing the 1.4 km simulation to the equivalent well-tested and calibrated global simulations at 9 km grid spacing with and without parametrized deep convection. The explicit simulation of deep convection at 1.4 km results in a realistic large-scale circulation, better representation of convective storm activity, and stronger convective gravity wave activity when compared to the 9 km simulation with parametrized deep convection. Comparison of the 1.4 km simulation to the 9 km simulation without parametrized deep convection shows that switching off deep convection parametrization at a too coarse resolution (i.e., 9 km) generates too strong convective gravity waves. Based on the limited statistics available, improvements to the Madden-Julian Oscillation or tropical precipitation are not observed at 1.4 km, suggesting that other Earth system model components and/or their interaction are important for an accurate representation of these processes and may well need adjusting at deep convection resolving resolutions. Overall, the good agreement of the 1.4 km simulation with the 9 km simulation with parametrized deep convection is remarkable, despite one of the most fundamental parametrizations being turned off at 1.4 km resolution and despite no adjustments being made to the remaining parametrizations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Radiosonde Observations of Environments Supporting Deep Moist Convection Initiation during RELAMPAGO-CACTI

The Remote Sensing of Electrification, Lightning, and Mesoscale/Microscale Processes with Adaptive Ground Observations (RELAMPAGO) and Cloud, Aerosol, and Complex Terrain Interactions (CACTI) projects deployed a high-spatiotemporal-resolution radiosonde network to examine environments supporting deep convection in the complex terrain of central Argentina. This study aims to characterize atmospheric profiles most representative of the near-cloud environment (in time and space) to identify the mesoscale ingredients affecting storm initiation and growth. Spatiotemporal autocorrelation analysis of the soundings reveals that there is considerable environmental heterogeneity, with boundary layer thermodynamic and kinematic fields becoming statistically uncorrelated on scales of 1–2 h and 30 km. Using this as guidance, we examine a variety of environmental parameters derived from soundings collected within close proximity (30 km in space and 30 min in time) of 44 events over 9 days where the atmosphere either: 1) supported the initiation of sustained precipitating convection, 2) yielded weak and short-lived precipitating convection, or 3) produced no precipitating convection in disagreement with numerical forecasts from convection-allowing models (i.e., Null events). There are large statistical differences between the Null event environments and those supporting any convective precipitation. Null event profiles contained larger convective available potential energy, but had low free-tropospheric relative humidity, higher freezing levels, and evidence of limited horizontal convergence near the terrain at low levels that likely suppressed deep convective growth. Here, we also present evidence from the radiosonde and satellite measurements that flow–terrain interactions may yield gravity wave activity that affects CI outcome.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Diabatic Eddy Forcing Increases Persistence and Opposes Propagation of the Southern Annular Mode in MERRA-2

Abstract As a dominant mode of jet variability on subseasonal time scales, the Southern Annular Mode (SAM) provides a window into how the atmosphere can produce internal oscillations on longer-than-synoptic time scales. While SAM’s existence can be explained by dry, purely barotropic theories, the time scale for its persistence and propagation is set by a lagged interaction between barotropic and baroclinic mechanisms, making the exact physical mechanisms challenging to identify and to simulate, even in latest generation models. By partitioning the eddy momentum flux convergence in MERRA-2 using an eddy–mean flow interaction framework, we demonstrate that diabatic processes (condensation and radiative heating) are the main contributors to SAM’s persistence in its stationary regime, as well as the key for preventing propagation in this regime. In SAM’s propagating regime, baroclinic and diabatic feedbacks also dominate the eddy–jet feedback. However, propagation is initiated by barotropic shifts in upper-level wave breaking and then sustained by a baroclinic response, leading to a roughly 60-day oscillation period. This barotropic propagation mechanism has been identified in dry, idealized models, but here we show evidence of this mechanism for the first time in reanalysis. The diabatic feedbacks on SAM are consistent with modulation of the storm-track latitude by SAM, altering the emission temperature and cloud cover over individual waves. Therefore, future attempts to improve the SAM time scale in models should focus on the storm-track location, as well as the roles of the cloud and moisture parameterizations. Significance Statement As they circumnavigate the planet, the tropospheric jet streams slowly drift north and south over about 30 days, longer than the normal limit of weather prediction. Understanding the source of this “memory” could improve our knowledge of how the atmosphere organizes itself and our ability to make long-term forecasts. Current theories have identified several possible internal atmospheric interactions responsible for this memory. Yet most of the theories for understanding the jets’ behavior assume that this behavior is only weakly influenced by atmospheric water vapor. We show that this assumption is not enough to understand jet persistence. Instead, clouds and precipitation are more important contributors in reanalysis data than internal “dry” mechanisms to this memory of the Southern Hemisphere jet.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Climate Dynamics Preceding Summer Forest Fires in California and the Extreme Case of 2018

Recent record-breaking wildfire seasons in California prompt an investigation into the climate patterns that typically precede anomalous summer burned forest area. Using burned-area data from the U.S. Forest Service’s Monitoring Trends in Burn Severity (MTBS) product and climate data from the fifth major global reanalysis produced by the European Centre for Medium-Range Weather Forecasts (ERA5) over 1984–2018, relationships between the interannual variability of antecedent climate anomalies and July California burned area are spatially and temporally characterized. Lag correlations show that antecedent high vapor pressure deficit (VPD), high temperatures, frequent extreme high temperature days, low precipitation, high subsidence, high geopotential height, low soil moisture, and low snowpack and snowmelt anomalies all correlate significantly with July California burned area as far back as the January before the fire season. Seasonal regression maps indicate that a global midlatitude atmospheric wave train in late winter is associated with anomalous July California burned area. July 2018, a year with especially high burned area, was to some extent consistent with the general patterns revealed by the regressions: low winter precipitation and high spring VPD preceded the extreme burned area. However, geopotential height anomaly patterns were distinct from those in the regressions. Extreme July heat likely contributed to the extent of the fires ignited that month, even though extreme July temperatures do not historically significantly correlate with July burned area. While the 2018 antecedent climate conditions were typical of a high-burned-area year, they were not extreme, demonstrating the likely limits of statistical prediction of extreme fire seasons and the need for individual case studies of extreme years.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

LISA Galactic binaries in the Roman Galactic Bulge Time-Domain Survey

ABSTRACT Short-period Galactic white dwarf binaries detectable by Laser Interferometer Space Antenna are the only guaranteed persistent sources for multimessenger gravitational-wave astronomy. Large-scale surveys in the 2020s present an opportunity to conduct preparatory science campaigns to maximize the science yield from future multimessenger targets. The Nancy Grace Roman Space Telescope Galactic Bulge Time-Domain Survey will (in its Reference Survey design) image seven fields in the Galactic Bulge approximately 40 000 times each. Although the Reference Survey cadence is optimized for detecting exoplanets via microlensing, it is also capable of detecting short-period white dwarf binaries. In this paper, we present forecasts for the number of detached short-period binaries the Roman Galactic Bulge Time-Domain Survey will discover and the implications for the design of electromagnetic surveys. Although population models are highly uncertain, we find a high probability that the baseline survey will detect of the order of ∼5 detached white dwarf binaries. The Reference Survey would also have a ${\gtrsim} 20\,{\rm per\,cent}$ chance of detecting several known benchmark white dwarf binaries at the distance of the Galactic Bulge.

Astronomy & Astrophysics↗
DOE OSTI2023

Assessing Latent and Kinetic Energy Trend Changes in Extratropical Cyclones From 1940 to 2020: Results From ERA‐5 Reanalysis

Abstract Baroclinic or extratropical cyclones (ETCs) transport heat and moisture to higher latitudes, making it fundamentally important to understand how their influence changes as Earth's climate evolves. A 2–8‐day Lanzcos bandpass filter is applied to European Center for Medium Range Weather Forecasting 5th Generation Reanalysis latent energy (LE) and kinetic energy (KE) data to assess how ETCs have changed from 1940 to 2020 relative to full‐scale changes in LE and KE. Full‐scale KE trends are more positive at high latitudes relative to mid‐latitudes, confirming several previous studies that ETCs have shifted poleward. LE increases have occurred globally, and trends in both full‐scale LE and KE are statistically significant in the southern high latitudes. The high relative fractional contribution of 2–8‐day LE wave power and trend clearly suggest that ETCs have an increasingly important role in poleward moisture transport but are not solely responsible for the observed statistically significant increases.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Integrating Reanalysis and Satellite Cloud Information to Estimate Surface Downward Long-Wave Radiation

The estimation of downward long-wave radiation (DLR) at the surface is very important for the understanding of the Earth’s radiative budget with implications in surface–atmosphere exchanges, climate variability, and global warming. Theoretical radiative transfer and observationally based studies identify the crucial role of clouds in modulating the temporal and spatial variability of DLR. In this study, a new machine learning algorithm that uses multivariate adaptive regression splines (MARS) and the combination of near-surface meteorological data with satellite cloud information is proposed. The new algorithm is compared with the current operational formulation used by the European Organization for the Exploitation of Meteorological Satellites (EUMETSAT) Satellite Application Facility on Land Surface Analysis (LSA-SAF). Both algorithms use near-surface temperature and dewpoint temperature along with total column water vapor from the latest European Centre for Medium-range Weather Forecasts (ECMWF) reanalysis ERA5 and satellite cloud information from the Meteosat Second Generation. The algorithms are trained and validated using both ECMWF-ERA5 and DLR acquired from 23 ground stations as part of the Baseline Surface Radiation Network (BSRN) and the Atmospheric Radiation Measurement (ARM) user facility. Results show that the MARS algorithm generally improves DLR estimation in comparison with other model estimates, particularly when trained with observations. When considering all the validation data, root mean square errors (RMSEs) of 18.76, 23.55, and 22.08 W·m –2 are obtained for MARS, operational LSA-SAF, and ERA5, respectively. The added value of using the satellite cloud information is accessed by comparing with estimates driven by ERA5 total cloud cover, showing an increase of 17% of the RMSE. The consistency of MARS estimate is also tested against an independent dataset of 52 ground stations (from FLUXNET2015), further supporting the good performance of the proposed model.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Effect of Rocky Mountains and Tibetan Plateau 1998 Spring Land Temperature on N. American and East Asian Summer Precipitation Anomalies

This work follows up on the GEWEX/LS4P Phase I (LS4P-I) experiments, a community effort highlighting the spring land surface temperature anomalies in the Tibetan Plateau (TP) as a useful source for subseasonal to seasonal (S2S) prediction of summer precipitation in global hot spot regions, particularly in East Asia and North America. This paper extends the investigation to both the US Rocky Mountain (RM) region and the TP, considering the 1998 summer drought/flood event in North America/East Asia, respectively, as a case study. A previously developed initialization method for land surface temperature/subsurface temperature (LST/SUBT) is used in the NCEP Global Forecast System, coupled with a land model, SSiB2 (GFS/SSiB2), to produce observed RM cold May temperature anomaly. Forward simulation yields June precipitation anomalies at five remote locations. Likewise, the TP warm May temperature anomaly also produces June precipitation anomalies at these five locations. The effects of RM (cold) and TP (warm) temperature anomalies are consistent in the US South Coastal regions and the south Yangtze River Basin, yielding 49% (42%) of observed drought and 34% (44%) of observed flood, respectively. These LST/SUBT effects in RM and TP induce a global large-scale wave train linking North America with the TP, affecting the subtropical westerly jet and thereby modulating summer precipitation. Global SST effect is examined for comparison but does not yield statistically significant June precipitation anomalies in GFS/SSiB2. Furthermore, this study adds to evidence that high-mountain LST effects in the RM and TP are first-order sources of S2S precipitation predictability in summer months.

Nayak, Hara Prasad [University of California, Los ↗
DOE OSTI2021

Experimental Investigation of Elastodynamic Nonlinear Response of Dry Intact, Fractured and Saturated Rock

Nonlinear elastodynamic response of fractured rocks carries crucial information on fracture features that can be exploited to forecast flow properties, friction constitutive behavior and poromechanical response. Well-controlled laboratory experiments are designed to measure the nonlinear elastodynamic response of Westerly granite in three states: dry intact, dry fractured and saturated fractured. We study the effect of fracturing and saturation in modifying the elastodynamic response of the rock. Each sample is tested at a normal stress level of 15 MPa. We measure the elastodynamic response of an intact L-shaped sample of Westerly Granite subjected to normal stress oscillations of prescribed amplitudes (0.2–1 MPa) and frequencies (0.1, 1, 10 Hz). Ultrasonic waves transmitted across the sample are used to monitor the evolution of wave velocity before, during and after dynamic stressing. Furthermore, the nonlinearity of the elastodynamic response is measured in terms of: (1) the offset in normalized wave velocity; (2) the amplitude of wave velocity fluctuation during the oscillations; and (3) recovery rate of the wave velocity post-oscillation. We observe that the three nonlinearity parameters show a similar trend. Irrespective of the parameter, the nonlinearity measures higher for sample in dry-intact condition than that for dry-fractured and the saturated-fractured sample exhibits smaller nonlinearity than the dry-fractured sample. As expected, the saturated sample exhibits less nonlinearity than the dry intact and fractured samples due to the presence of interstitial fluid and the resulting increased interface stiffness. Conversely, the dry intact rock shows a higher nonlinearity than the dry fractured. We use numerical simulations to show that the presence of fracture significantly alters the strain distribution across the bulk of the sample and only the contacting asperities are highly strained, thus resulting in a decrease in the measured elastodynamic nonlinearity.

58 GEOSCIENCES↗
DOE OSTI2022

The ICON-A model for direct QBO simulations on GPUs (version icon-cscs:baf28a514)

Abstract. Classical numerical models for the global atmosphere, as used for numerical weather forecasting or climate research, have been developed for conventional central processing unit (CPU) architectures. This hinders the employment of such models on current top-performing supercomputers, which achieve their computing power with hybrid architectures, mostly using graphics processing units (GPUs). Thus also scientific applications of such models are restricted to the lesser computer power of CPUs. Here we present the development of a GPU-enabled version of the ICON atmosphere model (ICON-A), motivated by a research project on the quasi-biennial oscillation (QBO), a global-scale wind oscillation in the equatorial stratosphere that depends on a broad spectrum of atmospheric waves, which originates from tropical deep convection. Resolving the relevant scales, from a few kilometers to the size of the globe, is a formidable computational problem, which can only be realized now on top-performing supercomputers. This motivated porting ICON-A, in the specific configuration needed for the research project, in a first step to the GPU architecture of the Piz Daint computer at the Swiss National Supercomputing Centre and in a second step to the JUWELS Booster computer at the Forschungszentrum Jülich. On Piz Daint, the ported code achieves a single-node GPU vs. CPU speedup factor of 6.4 and allows for global experiments at a horizontal resolution of 5 km on 1024 computing nodes with 1 GPU per node with a turnover of 48 simulated days per day. On JUWELS Booster, the more modern hardware in combination with an upgraded code base allows for simulations at the same resolution on 128 computing nodes with 4 GPUs per node and a turnover of 133 simulated days per day. Additionally, the code still remains functional on CPUs, as is demonstrated by additional experiments on the Levante compute system at the German Climate Computing Center. While the application shows good weak scaling over the tested 16-fold increase in grid size and node count, making also higher resolved global simulations possible, the strong scaling on GPUs is relatively poor, which limits the options to increase turnover with more nodes. Initial experiments demonstrate that the ICON-A model can simulate downward-propagating QBO jets, which are driven by wave–mean flow interaction.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Sensitivity of mesoscale modeling to urban morphological feature inputs and implications for characterizing urban sustainability

We examine the differences in meteorological output from the Weather Research and Forecasting (WRF) model run at 270 m horizontal resolution using 10 m, 100 m and 1 km resolution 3D neighborhood morphological inputs and with no morphological inputs. We find that the spatial variability in temperature, humidity, and other meteorological variables across the city can vary with the resolution and the coverage of the 3D urban morphological input, and that larger differences occur between simulations run without 3D morphological input and those run with some type of 3D morphology. We also find that the inclusion of input-building-defined roughness length calculations would improve simulation results further. We show that these inputs produce different patterns of heat wave spatial heterogeneity across the city of Washington, DC. These findings suggest that understanding neighborhood level urban sustainability under extreme heat waves, especially for vulnerable neighborhoods, requires attention to the representation of surface terrain in numerical weather models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Ab initio calculations for void swelling bias in α - and δ -plutonium

Void swelling can develop in materials under persistent irradiation when nonequilibrium vacancy and self-interstitial populations migrate under sufficiently asymmetric interaction biases. In conventional metals, the propensity is determined to the first approximation by comparing point-defect relaxation strains. We thus present DFT-based calculations of structures and formation energies and volumes of point defects in the α and the δ phases of plutonium. Here we discuss the pros and cons of various levels of electronic structure theory: spin polarization, spin-orbit coupling, and orbital polarization. Our results show that lattice defects in δ-Pu, in contrast to most fcc metals, have surprisingly small formation volumes. Equally unexpected are the large defect formation volumes found in the low-symmetry α-Pu phase. Both these unusual properties can be satisfactorily explained from defect-induced spin/orbital moment formation and destruction in the Pu phases. Surprisingly, the point defects in α-Pu are found to induce far larger transformation of the local electronic structure than in δ-Pu. When we use the calculated defect properties to estimate the classic void swelling bias in each of the phases, we find it to be unusually small in δ-Pu but likely much larger in α-Pu. Hence, swelling rates and mechanisms can diverge dramatically between the different phases of Pu. Especially in the transient regime before the formation of large defect clusters, the swelling rate of α-Pu can reliably be expected to be much larger than δ-Pu. However, accurate forecasts over longer times will require the conventional void-swelling theory to be modified to handle the complexities presented by the different Pu phases. As a case in point, we show the possible anomalous temperature dependence of vacancy properties in δ-Pu, caused by entropic contributions from defect-induced spin-lattice fluctuations. Such complications may affect defect-defect interactions and thus alter the void swelling bias.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Probing physics beyond the standard model: limits from BBN and the CMB independently and combined

We present new Big Bang Nucleosynthesis (BBN) limits on the cosmic expansion rate or relativistic energy density, quantified via the number N ν of equivalent neutrino species. We use the latest light element observations, neutron mean lifetime, and update our evaluation for the nuclear rates d + d → 3 He + n and d + d → 3 H+p. Combining this result with the independent constraints from the cosmic microwave background (CMB) yields tight limits on new physics that perturbs N ν and η prior to cosmic nucleosynthesis: a joint BBN+CMB analysis gives N ν = 2.898 ± 0.141, resulting in N ν < 3.180 at 2σ. We apply these limits to a wide variety of new physics scenarios including right-handed neutrinos, dark radiation, and a stochastic gravitational wave background. The strength of the independent BBN and CMB constraints now opens a new window: we can search for limits on potential changes in N ν and/or the baryon-to-photon ratio η between the two epochs. The present data place strong constraints on the allowed changes in N ν between BBN and CMB decoupling; for example, we find -0.708 < N ν CMB - N ν BBN < 0.328 in the case where η and the primordial helium mass fraction Yp are unchanged between the two epochs; we also give limits on the allowed variations in η or in (η, N ν ) jointly. We discuss scenarios in which such changes could occur, and show that BBN+CMB results combine to place important constraints on some early dark energy models to explain the H 0 tension. Looking to the future, we forecast the tightened precision for N ν arising from both CMB Stage 4 measurements as well as improvements in astronomical 4 He measurements. Here, we find that CMB-S4 combined with present BBN and light element observation precision can give σ(N ν ) ≃ 0.03. Such future precision would reveal the expected effect of neutrino heating (N eff -3 = 0.044) of the CMB during BBN, and would be near the level to reveal any particle species ever in thermal equilibrium with the standard model. Improved Y p measurements can push this precision even further.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Dissipation Scaled Internal Wave Drag in a Global Heterogeneously Coupled Internal/External Mode Total Water Level Model

This study showcases a global, heterogeneously coupled total water level system wherein salinity and temperature outputs from a coarser-resolution (~12 km) ocean general circulation model are used to calculate density-driven terms within a global, higher-resolution (~2.5 km) depth-averaged total water level model. We demonstrate that the inclusion of baroclinic forcing in the barotropic model requires modification of the internal wave drag term to prevent excess degradation of tidal results compared to the barotropic model. By scaling the internal tide dissipation by an easy to calculate dissipation ratio, the resulting heterogeneously coupled model has complex root mean square errors (RMSE) of 2.27 cm in the deep ocean and 12.16 cm in shallow waters for the M 2 tidal constituent. While this represents a 10%–20% deterioration as compared to the barotropic model, the improvements in total water level prediction more than offset this degradation. Global median RMSE compared to observations of total water levels, 30-day sea levels, and non-tidal residuals improve by 1.86 (18.5%), 2.55 (42.5%), and 0.36 (5.3%) cm respectively. The drastic improvement in model performance highlights the importance of including density-driven effects within global hydrodynamic models and will help to improve the results of both hindcasts and forecasts in modeling extreme and nuisance flooding. With only an 11% increase in model run time compared to the fully barotropic total water level model, this approach paves the way for high resolution coastal water level and flood models to be used alongside climate models, improving operational forecasting of total water levels.

Blakely, Coleman Peter [University of Notre Dame, ↗
DOE OSTI2025

Rapid wavefield forecasting for earthquake early warning via deep sequence to sequence learning

We propose a deep learning model, WaveCastNet, to forecast high-dimensional wavefields. WaveCastNet integrates a convolutional long expressive memory architecture into a sequence-to-sequence forecasting framework, enabling it to model long-term dependencies and multiscale patterns in both space and time. By sharing weights across spatial and temporal dimensions, WaveCastNet requires significantly fewer parameters than more resource-intensive models such as transformers, resulting in faster inference times. Crucially, WaveCastNet also generalizes better than transformers to rare and critical seismic scenarios, such as high-magnitude earthquakes. Here, we show the ability of the model to predict the intensity and timing of destructive ground motions in real time, using simulated data from the San Francisco Bay Area. Furthermore, we demonstrate its zero-shot capabilities by evaluating WaveCastNet on real earthquake data. Our approach does not require estimating earthquake magnitudes and epicenters, steps that are prone to error in conventional methods, nor does it rely on empirical ground-motion models, which often fail to capture strongly heterogeneous wave propagation effects.

Geophysics↗
DOE OSTI2026

Wind and Temperature Consensus at Horn Point, HU-Beltsville, Piney Run (Maryland) in support of CoURAGE

The Maryland Department of the Environment (MDE) operates a ground-based atmospheric profiling network consisting of collocated radar wind profilers (RWP) and radio acoustic sounding systems (RASS) as part of its Ambient Air Monitoring Program. This network provides continuous observations of wind and temperature structure in the lower troposphere to support air quality forecasting, regulatory analysis, and atmospheric research. The network currently includes three fixed sites across Maryland: Horn Point (HP, lower eastern shore) [38.587525°,-76.141006°], Howard University-Beltsville (HUB, central Maryland) [39.055277°, -76.878632°], and Piney Run (PR, western Maryland) [39.705950°, -79.012000°] The network is designed to capture regional variability in atmospheric transport and boundary-layer processes. These systems measure vertical profiles of horizontal wind speed and direction using Doppler radar techniques, with observations typically spanning from ~100 m above ground level up to approximately 2.5–4 km. Measurements are derived from the Doppler shift of backscattered electromagnetic signals, enabling retrieval of wind vectors at multiple altitudes with high temporal resolution (e.g., 30-minute averages reported every 6 minutes). Each radar wind profiler is paired with a Radio Acoustic Sounding System (RASS) to provide profiles of virtual temperature in the lower atmosphere (~100–200 m AGL) by measuring the propagation speed of acoustic waves. Together, the RWP/RASS system yields a coupled data set of thermodynamic and kinematic atmospheric structure, including additional parameters such as vertical velocity, radial velocity, signal-to-noise ratio, and spectral width for advanced analysis. There are two types of files for each station: wind data (files with a "w" prefix) and virtual temperature RASS data (files with a "t" prefix). The wind data files are in the format wYYDDD.cns, where YY is the 2-digit year and DDD is the day of the year. The RASS virtual temperature data files are in the format tYYDDD.cns. Each record has the following header structure: Line 1 : Station Name RASS files Line 2 : RASS rev DeTect_2.0, WINDS files Line 2 : WINDS rev ATI 5.1 Line 3 : N latitude, W longitude, and site elevation (m) Line 4 : Date and begin time of consensus: yy mm dd hh mn ss plus # minutes to add to get UTC Line 5 : Consensus averaging time (minutes); number of beams; number of range gates Line 6 : Number of records required to make consensus (num) total number of records (tot) and the consensus window size (m/s) in the format: num:tot (window) RASS files Line 7 : no. of coded cells, no. of spec, pulse width (ns), and inter-pulse period (µs), WINDS files Line 7 : No. of coded cells, no. of spectra, pulse width (ns), and inter-pulse period (µs), each with a pair of values: first value is for oblique beams, second for vertical RASS files Line 8 : Full scale Doppler value (m/s) Delay to first gate (ns) Number of gates Spacing of gates (ns), WINDS files Line 8 : Full scale Doppler velocity (m/s), oblique and vertical Vertical correction applied to oblique beams? (0 = no, 1 = yes) Delay to first gate (ns), oblique and vertical Number of gates, oblique and vertical Spacing of gates (ns), oblique and vertical Line 9 : Azimuth and elevation (9s indicate vertical beam not used) RASS files Line 10, values : HT = Height above ground (km), T = Uncorrected virtual temperature consensus (deg C), Tc = Corrected virtual temperature consensus (deg C), W = Vertical wind consensus (9s indicate vertical beam not used, w-component, positive upward, m/s), CNT = Number of records that made consensus (for the 3 values in same order), SNR = Average signal to noise ratio (dB) of records in consensus (same order) WINDS files Line 10, values : HT = Height above ground (km), SPD = Wind speed (m/s), DIR = Wind direction (deg E of N from N), RAD = Radial velocities for each beam (m/s) in order given in azimuth and elevation line (positive toward radar; 9s indicate vertical beam not used, CNT = Number of records that made consensus, SNR = Average signal to noise ratio (dB) of records in consensus

{"wind speed and direction",temperature}↗
DOE OSTI2022

Simons Observatory: Constraining inflationary gravitational waves with multitracer B-mode delensing

We introduce and validate a delensing framework for the Simons Observatory (SO), which will be used to improve constraints on inflationary gravitational waves by reducing the lensing noise in measurements of the B modes in CMB polarization. SO will initially observe CMB by using three small aperture telescopes and one large-aperture telescope. While polarization maps from small-aperture telescopes will be used to constrain inflationary gravitational waves, the internal CMB lensing maps used to delens will be reconstructed from data from the large-aperture telescope. Since lensing maps obtained from the SO data will be noise dominated on subdegree scales, the SO lensing framework constructs a template for lensing-induced B modes by combining internal CMB lensing maps with maps of the cosmic infrared background from Planck as well as galaxy density maps from the LSST survey. We construct a likelihood for constraining the tensor-to-scalar ratio r that contains auto and cross spectra between observed B modes and the lensing B-mode template. We test our delensing analysis pipeline on map-based simulations containing survey nonidealities, but that, for this initial exploration, does not include contamination from Galactic and extragalactic foregrounds. We find that the SO survey masking and inhomogeneous and atmospheric noise have very little impact on the delensing performance, and the r constraint becomes σ(r)≈0.0015 which is close to that obtained from the idealized forecasts in the absence of the Galactic foreground and is nearly a factor of 2 tighter than without delensing. We also find that uncertainties in the external large-scale structure tracers used in our multitracer delensing pipeline lead to bias much smaller than the 1σ statistical uncertainties.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Simulations of axion minihalos

The axion, motivated as a solution to the strong CP problem, is also a viable dark matter candidate. The axion field takes random values in causally disconnected regions if the symmetry breaking that establishes the particle occurs after inflation, leading to white-noise density fluctuations at low wave numbers and forming dense minihalos with sub-planetary masses subsequently. There have been two recent proposals that appear capable of testing this scenario, namely using pulsar timing arrays and studying cosmological microlensing caustics. Motivated by these proposals, we use N-body simulations to study the formation of substructures from white-noise density fluctuations. The density profiles of our relaxed axion minihalos can be described by the Navarro–Frenk–White profile, and the minihalos’ concentration number agrees well with a simple, physically-motivated model. We develop a semi-analytic formula to fit the mass function from our simulation, which agrees broadly at different redshifts and only differs at factor of two level from classic halo mass functions. This analytic mass function allows us to consider uncertainties in the post-inflation axion scenario, as well as extrapolate our high-redshift simulations results to the present. Our work estimates the present-day abundance of axion substructures, as is necessary for predicting their effect on cosmological microlensing caustics and pulsar timing. Our calculations suggest that if pulsar timing and microlensing probes can reach recent sensitivity forecasts, they may be sensitive to the post-inflation axion dark matter scenario, even when accounting for uncertainties pertaining to axion strings. For pulsar timing, the most significant caveat is whether axion minihalos are disrupted by stars, which our estimates show is mildly important at the most relevant masses. In conclusion, as our gravitational simulations are scale invariant, the results can be extended to models where the dark matter is comprised of other axion-like particles and even clusters of primordial black holes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗