Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SLOPE”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,117 records · Page 62

Atmospheric Energy Spectra in Global Kilometre-Scale Models

Eleven 40-day long integrations of five different global models with horizontal resolutions of less than 9 km are compared in terms of their global energy spectra. The method of normal-mode function decomposition is used to distinguish between balanced (Rossby wave; RW) and unbalanced (inertia-gravity wave; IGW) circulation. The simulations produce the expected canonical shape of the spectra, but their spectral slopes at mesoscales, and the zonal scale at which RW and IGW spectra intersect differ significantly. The partitioning of total wave energies into RWs an IGWs is most sensitive to the turbulence closure scheme and this partitioning is what determines the spectral crossing scale in the simulations, which differs by a factor of up to two. It implies that care must be taken when using simple spatial filtering to compare gravity wave phenomena in storm-resolving simulations, even when the model horizontal resolutions are similar. In contrast to the energy partitioning between the RWs and IGWs, changes in turbulence closure schemes do not seem to strongly affect spectral slopes, which only exhibit major differences at mesoscales. Despite their minor contribution to the global (horizontal kinetic plus potential available) energy, small scales are important for driving the global mean circulation. Our results support the conclusions of previous studies that the strength of convection is a relevant factor for explaining discrepancies in the energies at small scales. The models studied here produce the major large-scale features of tropical precipitation patterns. However, particularly at large horizontal wavenumbers, the spectra of upper tropospheric vertical velocity, which is a good indicator for the strength of deep convection, differ by factors of three or more in energy. High vertical kinetic energies at small scales are mostly found in those models that do not use any convective parameterisation.

DYAMOND↗

The Influence of Land Use and Land Cover Change on Landslide Susceptibility in the Lower Mekong River Basin

The Lower Mekong River Basin in Southeast Asia experiences frequent rainfall-triggered landslides especially during the monsoon season. In this study, the influence of land use and land cover change and other causative factors on landslide susceptibility is evaluated in the Lower Mekong Basin. Frequency ratio analysis is performed to quantify the relationship between LULC change and susceptibility. Detailed landslide inventory maps are used for analysis with yearly LULC maps. The LULC change is used as a contributing variable in a logistic regression-based susceptibility model with other variables including distance to roads, slope, aspect, forest loss, and soil properties. The Receiver Operating Characteristic (ROC) curve and Area Under the Curve are estimated for the model trained by each landslide inventory. The models show good performance, with AUC values ranging from 0.697 to 0.958 and an average AUC equal to 0.820. Both the Frequency Ratio analysis and the Logistic Regression models indicate LULC change from agricultural land to forest has a positive correlation with landslide occurrence. The most significant factors in the models are found to be distance to roads, slope, and aspect. A better understanding of the effects of LULC on landslide susceptibility can be useful for local land and disaster management and for the implementation of LULC as a factor in future susceptibility models. Using datasets that are unique to the Lower Mekong region, this study provides additional insights into the relationship between causative factors and landslide activity to better inform regional and global landslide susceptibility modeling.

River Basin↗

Transition Region Contribution to AIA Observations in the Context of Coronal Heating

We investigate the ratio of coronal and transition region intensity in coronal loops observed by the AtmosphericImaging Assembly(AIA)on the Solar Dynamics Observatory(SDO). Using Enthalpy-based Thermal Evolution ofLoops(EBTEL)hydrodynamic simulations, we model loops with multiple lengths and energyfluxes heatedrandomly by events drawn from power-law distributions with different slopes and minimum delays between eventsto investigate how each of these parameters influences observable loop properties. We generate AIA intensitiesfrom the corona and transition region for each realization. The variations within and between models generatedwith these different parameters illustrate the sensitivity of narrowband imaging to the details of coronal heating.We then analyze the transition region and coronal emission from a number of observed active regions andfindbroad agreement with the trends in the models. In both models and observations, the transition region brightness issignificant, often greater than the coronal brightness in all six“coronal”AIA channels. We also identify an inverserelationship, consistent with heating theories, between the slope of the differential emission measure(DEM)coolward of the peak temperature and the observed ratio of coronal to transition region intensity. These resultshighlight the use of narrowband observations and the importance of properly considering the transition region ininvestigations of coronal heating

S J Schonfeld↗

Winds at the Mars 2020 Landing Site: 1. Near-Surface Wind Patterns at Jezero Crater

This is the first part of a two-part paper. NASA's Mars 2020 Perseverance rover measured winds on the Jezero crater floor close to the delta of an ancient river. A mostly repeatable diurnal cycle was observed and presented two regimes: (a) a convective regime, from dawn to sunset, with average easterly to southeasterly winds, during which maximum wind speeds were measured, and (b) a nighttime regime with westerly-northwesterly winds followed by a relatively calm period with highly variable wind directions as a function of sol and time of night. The timing and magnitude of the observed regimes are consistent with primary control by regional and local slope flows. Data suggest that the surface circulation at Jezero region in northern spring and summer is highly unaffected by large-scale circulation except during particular periods in the diurnal cycle or generally during dust storms, which is supported by MarsWRF model simulations. Consequently, the observed seasonal variability was weak. However, sol-to-sol and seasonal variability were measured, most of it during certain nighttime periods. Traveling waves consistent with baroclinic instability were clearly observed in surface winds at L_s ∼ 75°. The early MY36/2022A regional dust storm at L_s ∼ 153° disturbed the wind patterns with changes suggesting enhanced tidal flows. After sunset, the dust storm also produced detectable gravity wave activity, increasing the mixing in the nighttime planetary boundary layer during storm conditions. Inferred wind directions from dust devil movies strongly suggest that prevailing winds continued to be slope-driven during the late summer, fall and early winter seasons.

D. Viúdez-Moreiras↗

A Theoretical and Computational Revisit of Conversions Between Whitham’s F-Function and Equivalent Area

This paper proves mathematically that the integral transforms between Whitham’s F-function and equivalent area are the inverse transforms of each other if and only if the slope of the equivalent area at the origin is zero. This mathematical fact contradicts the accepted unconditional inverse relation between Whitham’s F-function and equivalent area in the sonic boom research literature. Piecewise linear approximations of an F-function and of the second derivative of an equivalent area are used to derive numerical formulas for conversions between Whitham’s F-function and equivalent area. Numerical results are included to show convergence of the numerical conversions as the maximum length of the segments for piecewise linear approximations goes to zero. These numerical conversions are approximately the inverse transforms of each other when the second derivative of an equivalent area is continuous and the slope of the equivalent area at the origin is zero.

Whitham's F-function↗

Analysis of the Smap Roughness Parameter and the Smap Vegetation Optical Depth

The SMAP product provides the soil moisture (SM) computed using three different retrieval algorithms: the single channel H and V algorithms (SCA-H and SCA-V), and the dual-channel algorithm (DCA) which in addition provides the vegetation optical depth (VOD). The roughness model and the roughness parameters play an important role in the determination of the soil moisture and the VOD. In this regard, the SMAP SCA and DCA utilize different approaches to incorporate the effect of roughness. In this work we will summarize those approaches and we will evaluate the effect of the DCA approach on the retrieval of VOD. We will compare the SMAP DCA roughness parameter h with topographic parameters such as DEM height, DEM slope, DEM height standard deviation and DEM slope standard deviation.

Mousavi, M.↗

Arctic Biogeochemical and Optical Properties of Dissolved Organic Matter Across River to Sea Gradients

Arctic landscapes are warming and becoming wetter due to changes in precipitation and the timing of snowmelt which consequently alters seasonal runoff and river discharge patterns. These changes in hydrology lead to increased mobilization and transport of terrestrial dissolved organic matter (DOM) to Arctic coastal seas where significant impacts on biogeochemical cycling can occur. Here, we present measurements of dissolved organic carbon (DOC) and chromophoric DOM (CDOM) in the Yukon River-to-Bering Sea system and two river plumes on the Alaska North Slope which flow into the Beaufort Sea. Our sampling characterized optical and biogeochemical properties of DOM during high and low river discharge periods for the Yukon River-Bering Sea system. The average DOC concentration at the multiple Yukon River mouths ranged from a high of 10.36 mg C L-1 during the ascending limb of the 2019 freshet (late May), 6.4 mg C L-1 during the descending limb of the 2019 freshet (late June), and a low of 3.86 mg C L-1 during low river discharge in August 2018. CDOM absorption coefficient at 412 nm (aCDOM(412)) averaged 8.23 m-1, 5.07 m-1, and 1.9 m-1, respectively. Several approaches to model DOC concentration based on its relationship with CDOM properties demonstrated cross-system seasonal and spatial robustness for these Arctic coastal systems despite spanning an order of magnitude decrease in DOC concentration from the lower Yukon River to the Northern Bering Sea as well as the North Slope systems. “Snapshot” fluxes of DOC and CDOM across the Yukon River Delta to Norton Sound were calculated from our measurements and modeled water fluxes forced with upstream USGS river gauge data. Our findings suggest that during high river flow, DOM reaches the delta largely unaltered by inputs or physical and biogeochemical processing and that the transformations of Yukon River DOM largely occur in the plume. However, during low summer discharge, multiple processes including local precipitation events, microbial decomposition, photochemistry, and likely others can alter the DOM properties within the lower Yukon River and Delta prior to flowing into Norton Sound.

Arctic↗

The Warm Neptune GJ 3470b Has a Polar Orbit

The warm Neptune GJ 3470b transits a nearby (d = 29 pc) bright slowly rotating M1.5-dwarf star. Using spectroscopic observations during two transits with the newly commissioned NEID spectrometer on the WIYN 3.5 m Telescope at Kitt Peak Observatory, we model the classical Rossiter–McLaughlin effect, yielding a sky-projected obliquity of λ=98 +15 −12 ˚and a v sin i = 0.85 +0.27 −0.33 kms -1 . Leveraging information about the rotation period and size of the host star, our analysis yields a true obliquity of ψ=95 +9 −8 ◦ , revealing that GJ 3470b is on a polar orbit. Using radial velocities from HIRES, HARPS, and the Habitable-zone Planet Finder, we show that the data are compatible with a long-term radial velocity (RV) slope of 𝛾̀=-0.0022±0.0011 ms -1 day -1 over a baseline of 12.9 yr. If the RV slope is due to acceleration from another companion in the system, we show that such a companion is capable of explaining the polar and mildly eccentric orbit of GJ 3470b using two different secular excitation models. The existence of an outer companion can be further constrained with additional RV observations, Gaia astrometry, and future high-contrast imaging observations. Lastly, we show that tidal heating from GJ 3470b’s mild eccentricity has most likely inflated the radius of GJ 3470b by a factor of ∼1.5–1.7, which could help account for its evaporating atmosphere.

Exoplanet astronomy↗

Evaluating the Feasibility of Data Simulation for Radiation Research: Estimating RBE Values for HZE Radiation

As the drive for deep space exploration intensifies, a comprehensive understanding of the health effects of radiation exposure becomes paramount to the future of human space flight. However, epidemiological data for radiation exposure, particularly to high-energy (HZE) ions, is limited, partially due to the financial and logistical costs of radiation studies. As an alternative, this study aims to assess the viability of data simulation strategies to accurately model potential study parameters prior to utilizing laboratory conditions. This study estimates a relative biological effectiveness (RBE) factor based on the solid tumor data for outbred mice provided by Edmundson et al. 2020. Excess relative risk (ERR) models for HZE-Fe ions and gamma radiation were estimated using Poisson regression with Weibull models to represent the background solid tumor hazard without radiation. RBE values were calculated from the ratio of the heavy ion linear slope to the gamma linear slope. The parameters from these models were then applied to simulate iterations of 300 datasets across HZE-Fe doses of 0.05, 0.2, 0.4, and 0.75 Gy; gamma radiation doses of 0.75, 2, and 3 Gy, as well as an unirradiated control group. Sample size per dose varied from 100-500 mice across simulations (800-4000 mice total per dataset). 1500 datasets were generated total (300 each for sample sizes 800, 1600, 3200, and 4000). For each dataset, ERR per radiation type and RBE for HZE-Fe were calculated. The RBE from the Edmondson data was calculated to be 5.55. The RBEs from the simulated data converged around this value across the increasing sample sizes. For n = 4000, the mean RBE across the 300 datasets was 5.71 (95% CI: 5.47-5.95). These results suggest that RBEs calculated from simulated data are credible. Based on this exploration, data simulation is a viable method of testing radiation studies. It provides a method of testing study assumptions and refining research questions prior to embarking on costly laboratory experiments.

Sneha Mehta↗

Preferential Acceleration of Heavy Ions in a Spontaneously Fragmenting Flare Current Sheet

We study the ion acceleration in a mesoscale, spontaneously fragmenting flare current sheet (SFCS) characterized by the presence of a plasmoid cascade. The main subject of our investigation is to determine whether and how plasmoid cascades at intermediate scales in a fragmented current sheet of a solar flare can impact the (preferential) acceleration of specific ions. The time evolution of the SFCS is obtained from high-resolution 2.5D MHD simulations. The ion trajectories (in the background fields resulting from the MHD model), energies, and pitch angles are calculated using a relativistic test-particle code based on the half-acceleration–rotation–half-acceleration method. For light ions, the main acceleration effects of electromagnetic fields within the SFCS are analyzed using the guiding center approximation. We identify regions with the most-efficient ion acceleration within the SFCS, the accelerator efficiency, and spectra of the accelerated ions. The influence of the charge-to-mass ratio on ion behavior is also studied and resulting ion abundances are compared with observational data. The main ion acceleration takes place in the regions with a strong polarization term, which is part of the first-order Fermi acceleration. Because the term is mass dependent, heavier ions undergo preferential acceleration. The ion energy spectra, abundance-enhancement factors, and differential fluxes, obtained from the model, exhibit power-law profiles, in agreement with observed solar energetic particle events. Nonetheless, the obtained slopes for the abundance-enhancement factor do not exactly match the observed data. The computed slopes and profiles are not sensitive to changes in the initial plasma temperature.

David Kramoliš↗

Ancient Martian aeolian sand dune deposits recorded in the stratigraphy of Valles Marineris and implications for past climates

Aeolian sediment transport, deposition, and erosion have been ongoing throughout Mars’ history. This record of widespread aeolian processes is preserved in landforms and geologic units that retain important clues about past environmental conditions including wind patterns. In this study we describe landforms within Melas Chasma, Valles Marineris, that occur in distinct groups with linear to crescentic shapes, arranged with a characteristic wavelength; some possess slope profiles analogous to modern sand dunes, yet show evidence for lithification. Based on the features’ dimensions, asymmetry, and spatial patterns relative to modern equivalents, we interpret these landforms to be two classes of aeolian bedforms: decameter-scale mega-ripples and sand dunes. The presence of superposed erosional features and depositional units indicates these landforms were cemented and likely ancient. Melas paleo-dunes are found atop Hesperian-aged layered deposits, but we estimate them to be younger, likely lithified in the Amazonian period. Although a range of degradation was observed, some paleo-dunes are >10 m tall and maintain steep lee sides (>25°), an uncommon scenario for terrestrial examples as other geologic processes lead to dune obliteration. The preserved paleo-bedform geometries are largely consistent with those of modern aeolian indicators, suggesting no major shifts in wind regime or contributing boundary conditions. Finally, we propose their appearance and context require sequential periods of dune migration, stabilization following catastrophic burial, cementation, differential erosion, exposure, and burial. The presence of wholly preserved duneforms appears to be more common on Mars compared to the Earth and may signal something important about Martian landscape evolution. Plain Language Summary Wind-driven sand dunes are common on modern Mars and the presence of certain sedimentary rock layers indicates these landforms occurred there in the past. Here, we explore features in the canyons of Valles Marineris which show characteristics frequently attributed to dunes, yet evidence for their lithification and burial is clear. Their horizontal extent, height, shape, slopes, and collective spatial patterns led to our interpretation of them as lithified dune fields. The heavily eroded appearance along with superposed elements such as craters and boulders led to the conclusion that these are relatively ancient landforms. Despite this erosion This level of preservation is rare for terrestrial sand dunes due to erosion and tectonics, thus providing an opportunity to reconstruct the various factors contributing to their history. Understanding the circumstances that led to the preservation of these ancient duneforms will yield crucial information regarding planetary sedimentary processes and the geologic history of the region.

Matthew Chojnacki↗

Relationship Between the Sub-Micron Fraction (SMF) and Fine-Mode Fraction (FMF) in the Context of AERONET Retrievals

The sub-micron (SM) aerosol optical depth (AOD) is an optical separation based on the fraction of particles below a specified cutoff radius of the particle size distribution (PSD) at a given particle radius. It is fundamentally different from spectrally separated FM (fine-mode) AOD. We present a simple (AOD-normalized) SM fraction versus FM fraction (SMF vs. FMF) linear equation that explains the well-recognized empirical result of SMF generally being greater than the FMF. The AERONET inversion (AERinv) products (combined inputs of spectral AOD and sky radiance) and the spectral deconvolution algorithm (SDA) products (input of AOD spectra) enable, respectively, an empirical SMF vs. FMF comparison at similar (columnar) remote sensing scales across a variety of aerosol types. SMF (AERinv-derived) vs. FMF (SDA-derived) behavior is primarily dependent on the relative truncated portion (ε c ) of the coarse-mode (CM) AOD associated with the cutoff portion of the CM PSD and, to a second order, the cutoff FM PSD and FM AOD (ε f ). The SMF vs. FMF equation largely explains the SMF vs. FMF behavior of the AERinv vs. SDA products as a function of PSD cutoff radius (“inflection point”) across an ensemble of AERONET sites and aerosol types (urban-industrial, biomass burning, dust, maritime and a mixed class of Arctic aerosols). The overarching dynamic was that the linear SMF vs. FMF relation pivots clockwise about the approximate (SMF, FMF) singularity of (1, 1) in a “linearly inverse” fashion (slope and intercept of approximately 1−ε c and ε c ) with increasing cutoff radius. SMF vs. FMF slopes and intercepts derived from AERinv and SDA retrievals confirmed the general domination of ε c over ε f in controlling that dynamic. A more general conclusion is the apparent confirmation that the optical impact of truncating modal (whole) PSD features can be detected by an SMF vs. FMF analysis.

Norman T. O'Neill↗

A computationally efficient statistically downscaled 100 m resolution Greenland product from the regional climate model MAR

The Greenland Ice Sheet (GrIS) has been contributing directly to sea level rise, and this contribution is projected to accelerate over the next decades. A crucial tool for studying the evolution of surface mass loss (e.g., surface mass balance, SMB) consists of regional climate models (RCMs), which can provide current estimates and future projections of sea level rise associated with such losses. However, one of the main limitations of RCMs is the relatively coarse horizontal spatial resolution at which outputs are currently generated. Here, we report results concerning the statistical downscaling of the SMB modeled by the Modèle Atmosphérique Régional (MAR) RCM from the original spatial resolution of 6 km to 100 m building on the relationship between elevation and mass losses in Greenland. To this goal, we developed a geospatial framework that allows the parallelization of the downscaling process, a crucial aspect to increase the computational efficiency of the algorithm. Using the results obtained in the case of the SMB, surface and air temperature are assessed through the comparison of the modeled outputs with in situ and satellite measurement. The downscaled products show a considerable improvement in the case of the downscaled product with respect to the original coarse output, with the coefficient of determination (R 2 ) increasing from 0.868 for the original MAR output to 0.935 for the SMB downscaled product. Moreover, the value of the slope and intercept of the linear regression fitting modeled and measured SMB values shifts from 0.865 for the original MAR to 1.015 for the downscaled product in the case of the slope and from the value −235 mm w.e. yr -1 (original) to −57 mm w.e. yr -1 (downscaled) in the case of the intercept, considerably improving upon results previously published in the literature.

Greenland Ice Sheet↗

Morphometric Characterization of Lunar Landing Sites

As ambitious surface exploration of the Moon commences in the 2020s, it is important to develop reliable and objective metrics for understanding the quality of future landing sites. A key prerequisite for any exploration and utilization of the lunar surface is a safe landing. One of the most important methods to provide understanding of potential metrics for landing site safety in the lunar context is the systematic comparison of candidate and historic landing sites. The goal of this project is to determine the morphometric parameters of the lunar surface at 16 successful lunar landing sites where adequate data exists to execute a quantitative comparison using various parameters. Quantitative analysis of Lunar Reconnaissance Orbiter (LRO) data should inform mission planning activities by providing morphologic metrics for landing sites including slope, Terrain Ruggedness Index (TRI), and rock abundance. These comparisons will assist mission planners and exploration scientists by providing “calibration points” for using Lunar Reconnaissance Orbiter (LRO) data to successfully plan and execute lunar powered descents in the future. The metrics included in this analysis describe the form of the terrain and can be calculated for any potential future landing site. Another objective of this project is to determine how the derived morphologic parameters for the same landing site change between common pixel scales. In ideal lighting conditions, 2m/px Narrow Angle Camera Digital Terrain Models (NAC DTMs) can be assembled using stereoscopic imagery from two or more concentric orbits of the LRO. However, due to low light conditions near the lunar poles, the highest quality data for many potential lunar landing sites comes from the Lunar Orbiter Laser Altimeter (LOLA). Therefore, we seek to assess how the derived slope and TRI values change as a function of changes in the DTM postings from 2m/px to 5m/px.

J M McCallion↗

Shape from spectra

We introduce a new unified atmospheric–topographic correction approach that estimates surface geometry directly from the radiance measurement. Surface topography influences the at-sensor radiance measurement, making precise topography modeling critical in applications like vegetation or snow studies in mountainous terrain. Currently, elevation maps are used to derive topographic variables such as the slope and sky-view factor. This process is error-prone since static global digital elevation models do not generally achieve the accuracy required, and even minor mismatches in spatial resolution can introduce significant artifacts in downstream processing. Here we demonstrate that it is possible to estimate topographic parameters directly from spectral data, ensuring perfect physical consistency, temporal coincidence, and spatial alignment. We present experiments estimating topographic slope in two scenes in Southern California, with data from NASA’s Next Generation Airborne Visible/Near Infrared Imaging Spectrometer (AVIRIS-NG). We compared our radiance-based estimates against high-resolution lidar datasets. Our initial validation result showed a correlation of R 2 = 0.864 (n = 160) over the homogeneous surface of Beckman Auditorium’s cone-shaped roof on the Caltech campus in Pasadena, California. We then validate the model over a larger study site near Santa Clarita, California, finding R 2 = 0.923 (n = 40,000) in a 350 x 350m area. The accuracy of our model estimates, combined with its systematic advantages over the alternative, show the potential of the approach for use in both airborne campaigns and orbital missions.

Nimrod Carmon↗

The Spectral Characteristics of Lunar Agglutinates: Visible-Near-Infrared Spectroscopy of Apollo Soil Separates

The lunar surface evolves over time due to space weathering, and the visible–near-infrared spectra of more mature (i.e., heavily weathered) soils are lower in reflectance and steeper in spectral slope (i.e., darker and redder) than their immature counterparts. These spectral changes have traditionally been attributed to the space-weathered rims of soil grains (and particularly nanophase iron therein). However, understudied thus far is the spectral role of agglutinates—the agglomerates of mineral and lithic fragments, nanophase iron, and glass that are formed by micrometeoroid impacts and are ubiquitous in mature lunar soils. We separated agglutinates and non-agglutinates from six lunar soils of varying maturity and composition, primarily from the 125–250 μm size fraction, and measured their visible–near-infrared reflectance spectra. For each soil, the agglutinate spectra are darker, redder, and have weaker absorption bands than the corresponding non-agglutinate and unsorted soil spectra. Moreover, greater soil maturity corresponds to darker agglutinate spectra with weaker absorption bands. These findings suggest that agglutinates (rather than solely the space-weathered rims) play an important role in both the darkening and reddening of mature soils—at least for the size fractions examined here. Comparisons with analog soils suggest that high nanophase iron abundance in agglutinates is likely responsible for their low reflectance and spectrally red slope. Additional studies of agglutinates are needed both to more comprehensively characterize their spectral properties (across size fractions and in mixing with non-agglutinates) and to assess the relative roles of agglutinates and rims in weathering-associated spectral changes. Plain Language Summary - In scientific study of the Moon, one key focus is surface processes: how do physical and chemical properties of the Moon’s surface change over time due to weathering (e.g., bombardment by micrometeoroids and by particles from the Sun)? Such investigations provide valuable insights into the Moon’s history (such as the ages of impact craters) that are often deduced from measurements of reflected light; as a soil is weathered it reflects light differently, which manifests visually as a progressive darkening of the soil. This phenomenon had primarily been attributed to weathering-associated development of rims on individual soil grains, but in this work we explored an alternative cause: soil particles known as agglutinates (misshapen, vesicular agglomerates of mineral fragments, iron, and glass that form due to weathering processes). We isolated agglutinates of six soil samples from the Moon and measured how they reflect light. We find that they reflect light in patterns reminiscent of how the Moon’s surface does when weathered. These findings suggest that agglutinates play a more important role than previously thought in determining the light-reflecting properties of the Moon’s surface, thus warranting greater and more nuanced consideration in future studies of how the Moon’s surface changes over time.

spectroscopy↗

Shock Capturing via Limiting for High-Order Methods Including Discontinuous Galerkin

As is well-known, popular high-order methods such as discontinuous Galerkin (DG) and flux reconstruction (FR) tend to generate oscillations near shocks and discontinuities, which can lead to negative pressure or density, ultimately causing code breakdown. For standard second-order methods, one approach to mitigating oscillations is to impose constraints on the calculated slopes, ensuring they do not become excessively steep. This limiting process is rooted in the idea of preserving monotonicity introduced by Van Leer (1974): when the data are monotone, limiting the slopes result in a monotone piecewise linear solution. However, a drawback of this approach is the loss of accuracy near extrema, where the non-monotone solution is somewhat flattened and achieves only first-order accuracy.

numerical methods↗

Satellite Remote-Sensing Capability to Assess Tropospheric-Column Ratios of Formaldehyde and Nitrogen Dioxide: Case Study During the Long Island Sound Tropospheric Ozone Study 2018 (LISTOS 2018) Field Campaign

Satellite retrievals of tropospheric-column formaldehyde (HCHO) and nitrogen dioxide (NO 2 ) are frequently used to investigate the sensitivity of ozone (O 3 ) production to emissions of nitrogen oxides and volatile organic carbon compounds. This study inter-compared the systematic biases and uncertainties in retrievals of NO 2 and HCHO, as well as resulting HCHO–NO 2 ratios (FNRs), from two commonly applied satellite sensors to investigate O 3 production sensitivities (Ozone Monitoring Instrument, OMI, and TROPOspheric Monitoring Instrument, TROPOMI) using airborne remote-sensing data taken during the Long Island Sound Tropospheric Ozone Study 2018 between 25 June and 6 September 2018. Compared to aircraft-based HCHO and NO 2 observations, the accuracy of OMI and TROPOMI were magnitude-dependent with high biases in clean environments and a tendency towards more accurate comparisons to even low biases in moderately polluted to polluted regions. OMI and TROPOMI NO 2 systematic biases were similar in magnitude (normalized median bias, NMB = 5 %–6 %; linear regression slope ≈ 0.5–0.6), with OMI having a high median bias and TROPOMI resulting in small low biases. Campaign-averaged uncertainties in the three satellite retrievals (NASA OMI; Quality Assurance for Essential Climate Variables, QA4ECV OMI; and TROPOMI) of NO 2 were generally similar, with TROPOMI retrievals having slightly less spread in the data compared to OMI. The three satellite products differed more when evaluating HCHO retrievals. Campaign-averaged tropospheric HCHO retrievals all had linear regression slopes ∼0.5 and NMBs of 39 %, 17 %, 13 %, and 23 % for NASA OMI, QA4ECV OMI, and TROPOMI at finer (0.05° x 0.05°) and coarser (0.15° x 0.15°) spatial resolution, respectively. Campaign-averaged uncertainty values (root mean square error, RMSE) in NASA and QA4ECV OMI HCHO retrievals were ~9.0 x 10 15 molecules cm –2 (∼ 50 %–55 % of mean column abundance), and the higher-spatial-resolution retrievals from TROPOMI resulted in RMSE values ∼30 % lower. Spatially averaging TROPOMI tropospheric-column HCHO, along with NO 2 and FNRs, to resolutions similar to the OMI reduced the uncertainty in these retrievals. Systematic biases in OMI and TROPOMI NO 2 and HCHO retrievals tended to cancel out, resulting in all three satellite products comparing well to observed FNRs. However, while satellite-derived FNRs had minimal campaign-averaged median biases, unresolved errors in the indicator species did not cancel out in FNR calculations, resulting in large RMSE values compared to observations. Uncertainties in HCHO retrievals were determined to drive the unresolved biases in FNR retrievals.

Matthew S. Johnson↗