Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Surface estimation”

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 595 records · Page 33

Surface radiation in the tropical Pacific

Monthly surface radiative fluxes have been calculated for the tropical Pacific between January 1970 and February 1978, using a radiative transfer parameterization. The radiative transfer parameterization included detailed treatments of the molecular and droplet absorptions, and surface and cloud reflections. The input data used in the calculations were obtained from the National Climatic Center (NCC), the National Center for Atmospheric Research (NCAR), and from the University of Hawaii. The results show that the distribution of surface radiation closely follows the distribution of cloudiness, and, to a lesser degree, humidity. The rms net error in the surface radiation estimates was about 15 W per sq m, with the largest contribution from uncertainties in the cloud cover and humidity data. The sensitivity of surface radiation parameterizations to input data errors is discussed, and some accuracy requirements for satellite retrievals of atmospheric and cloud parameters are proposed. The calculations are presented in the annually-averaged maps of surface radiation variations.

Chou, M.-D.↗

A climatology of stratospheric aerosol

A global climatology of stratospheric aerosol is created by combining nearly a decade (1979-1981 and 1984-1990) of contemporaneous observations from the Stratospheric Aerosol and Gas Experiment (SAGE I and II) and Stratospheric Aerosol Measurement (SAM II) instruments. One goal of this work is to provide a representative distribution of the aerosol layer for use in radiative and chemical modeling. A table of decadal average 1 micron extinction values is included, extending from the tropopause to 35 km and 80 deg S to 85 deg N, which allows estimation of surface area density. We find that the aerosol layer is distinctly volcanic in nature and suggest that the decadal average is a more useful estimate of future aerosol loading than a 'background' loading, which is never clearly achieved during the data record. This climatology lends insight into the general circulation of the stratosphere. Latitude - altitude sections of extinction radio at 1 micron are shown, average by decade, season, and phase of the quasi-biennial oscillation (QBO). A tropical reservoir region is diagnosed, with an 'upper' and a 'lower' transport regime. In the tropics above 22 km (upper regime), enhanced lofting occurs in the summer, with suppressed lofting or eddy dilution in the winter. In the extratropics within two scale heights of the tropopause (lower regime), poleward and downward transport is most robust during winter, especially in the northern hemisphere. The transport patterns persist into the subsequent equinoctial season. Ascent associated with QBO easterly shear favors detrainment in the upper regime, while relative descent and poleward spreading during QBO westerly shear favors detrainment in the lower regime. Extinction radio differences between the winter-spring and summer-fall hemispheres, and differences between the two phases of the QBO, are typically 20-50%. Dynamical implications of the aerosol distributions are explored, with focus on interhemispheric differences, strong subtropical gradients, and the pronounced annual cycle.

Hitchman, Matthew H.↗

Joint Retrieval of Surface BRDF from Geostationary and Polar-Orbiting Satellite Sensors

The latest geostationary sensors like GOES 16/17 ABI and Himawari 8/9 AHI provide high frequent observations of the Earth surface with continuously changing solar illumination geometries, which allow us to retrieve the surface Bidirectional Reflectance Distribution Function (BRDF) with dedicated atmospheric correction algorithms like MAIAC (Multi-Angle Implementation of Atmospheric Correction). However, because the viewing geometry of a specific location from the geostationary satellites are fixed, the angular sampling of surface BRDF by GEO (Geostationary Earth Orbit) sensors is far from comprehensive. This study tries to address this issue by exploring a GEO-LEO (Low-Earth-Orbit) synergy, in particular, jointly retrieving surface BRDF parameters with concurrent ABI/AHI and VIIRS top-of-atmosphere (TOA) reflectance for the near-infrared (NIR) band. The NIR band is chosen because the ABI, AHI, and VIIRS instruments have very similar spectral response functions in this band and therefore simplifies the requirements for cross-sensor radiometric calibration. We compile ABI/AHI and VIIRS TOA data with the best quality (cloud-free and low aerosol loading) at chosen AERONET sites. We then run the GeoNEX MAIAC algorithm to retrieve the Ross-Thick-Li-Sparse (RTLS) surface BRDF parameters with or without the AEORNET measured atmospheric aerosol optical depth (AOD) as inputs. The joint retrieval results are considered the best estimate of surface BRDF. We compare the joint BRDF retrievals with the corresponding MAIAC BRDF products, retrieved with ABI/AHI or VIIRS separately, to evaluate their differences. We expect that the jointly retrieved BRDF data are more robust than the standard products and may help us reduce uncertainties in higher-level earth observation satellite products.

Remote Sensing↗

The Global Methane Budget 2000-2012

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

Saunois, Marielle↗

Cooling of pulsars

Cooling rates are calculated for superfluid neutron stars of about one solar mass and 10 km radius, with magnetic fields from zero to about 10 to the 14th power Gauss, when possible internal friction effects are neglected. The results show that most old pulsars are so cold that thermal ionization of surface atoms would be negligible. At an age of a million years and with canonical magnetic fields of 10 to the 12th power Gauss, the estimated stellar surface temperature is several thousand to a hundred thousand degrees. However, if we neglect magnetic fields and superfluid states of nucleons, the same surfaces would be about a million degrees.

Tsuruta, S.↗

Toward proof of satellite anemometer concept

The S-193 oceanographic program is designed to demonstrate the operational feasibility of measuring ocean surface winds on a global basis with a composite radar-radiometer sensor. Near-simultaneous radiometric/scatterometric measurements offer a better opportunity to make improved estimates of surface winds through a variety of atmospheric conditions. Support is presented in developing interpretational techniques for the joint measurements, designing experiments to verify the joint interpretation techniques, predicting the observations from theory, evaluating system effects on the measurements, and developing system design criteria.

Classen, J.↗

Paleoatmospheric temperature structure

Radiative equilibrium and radiative convective temperature profiles for the earth's evolving atmosphere have been calculated. If the atmosphere evolved from one rich in carbon dioxide, and deficient in oxygen, to its present composition, the temperature structure showed considerable change. The models of 3 to 4 billion years ago display steadily decreasing temperatures with altitude, being 185 K at pressures associated with the present-day upper stratosphere. A lapse rate feature similar to the present-day tropopause is not indicated until about 1 billion years ago; but the stratospheric region is approximately 15 K colder than presently found at comparable pressures. Surface temperatures approximately 10 K warmer than at present existed until nearly 1 billion years ago. When the oxygen content exceeded roughly 0.1 times the present level, surface temperatures began to decrease. If biological processes are important to carbon dioxide - ozone variations, such as has been suggested during the Ice Ages, then estimates of surface temperature should include the effects of both gases.

Morss, D. A.↗

Polar motion-induced gravity

Variations in the geocentric position of the earth's rotation axis (polar motion) cause deformation within the earth. The effects of this deformation on surface gravity and on radial and horizontal positions of points on the earth's surface are estimated. The effects of the oceans and of the earth's anelasticity on this deformation are found to be negligible. Peak-to-peak variations in surface gravity of 10 microgals or more, and in radial motion of 1-2 cm are possible over six months or so. These numbers are small enough that they can probably not be used to learn about the earth; however, they are large enough to noticeably affect present high-quality geodetic observations.

Wahr, John M.↗

Systematic Errors in Simulated L-Band Brightness Temperature in the SMAP Level-4 Soil Moisture Analysis

The NASA Soil Moisture Active Passive (SMAP) mission has been providing L-band (1.4 GHz) brightness temperature (Tb) observations since April 2015. By assimilating the Tb observations into the NASA Catchment land surface model using a spatially distributed ensemble Kalman filter (EnKF), the SMAP Level-4 Soil Moisture (L4_SM) product provides global, 3-hourly, 9-km resolution estimates of surface (0-5 cm) and root-zone (0-100 cm) soil moisture with ~2.5-day latency for use in research and applications. The EnKF-based L4_SM analysis assumes unbiased forecast errors. Consequently, the seasonally varying bias between the model forecast Tb and the observed values is removed prior to the assimilation of the SMAP Tb observations. The L4_SM system is thus designed to only correct errors in synoptic-scale and interannual variations from the long-term mean seasonal cycle while maintaining the model’s (possibly wrong) climatology. In this paper, we examine the Tb observation-minus-forecast (O-F) residuals from the L4_SM Version 7 product (computed after rescaling the Tb observations to the mean seasonal cycle of the simulated Tb). The long-term average of the Tb O-F residuals has a global mean of only 0.13 K and locally small values, ranging from -1 to 3 K. The model forecast Tb, however, still exhibits undesirable systematic errors relative to the (rescaled) Tb observations. At some locations, the time-average Tb O-F values strongly depend on surface soil moisture (SM). At the Yanco SMAP core validation site, for example, the Tb O-F residuals typically range from 5 to 15 K under dry soil moisture conditions (SM < 0.15 m3 m-3) yet are predominantly negative under wet soil moisture conditions (SM > 0.25 m3 m-3), with values ranging from 0 to -40 K. This results in soil moisture analysis increments that persistently make the soil drier under dry SM conditions and persistently make the soil wetter under wet SM conditions, suggesting an error in the dynamic range of the simulated Tb, soil moisture or soil temperature. In this paper, we describe the higher-order systematic Tb forecast errors in more detail, examine their impact on the L4_SM product quality, and explore potential avenues to improve the L4_SM algorithm.

Rolf Reichle↗

Consequences of Different Air-Sea Feedbacks on Ocean Using MITgcm and MERRA-2 Forcing: Implications for Coupled Data Assimilation Systems

Ocean surface flux estimates from atmospheric and oceanic reanalyses contain errors that compensate for inaccuracies in the respective atmosphere and ocean models used to generate these reanalyses. A conundrum for climate studies is the discrepancy between surface fluxes that minimize model-data differences for an atmosphere-only model vs surface fluxes that minimize model-data differences for an ocean model. As a first step towards a consistent coupled ocean-atmosphere data-assimilation (DA) system, we compare surface net heat flux from a state-of-the-art atmospheric reanalysis, the Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2), to net heat flux from a state-of-the-art ocean state estimate, the Estimating the Circulation and Climate of the Ocean Version 4 (ECCO-v4). The possible impacts of the MERRA-2 and ECCO-v4 air-sea net heat flux difference in a coupled DA system were assessed using a set of experiments designed to imitate different “flavors” of a coupled DA system in an ocean-only setup. This was done by forcing the ECCO-v4 underlying ocean model - the Massachusetts Institute of Technology general circulation model (MITgcm) - with different sets of MERRA-2 fields and utilizing different forcing methods. By doing so we were able to turn off different air-sea feedbacks which, in a coupled DA setup, are partially muted by the constraining observations. The set of experiments, therefore, represents a range of active feedbacks in different “flavors” of coupled data-assimilation systems. For the period 1992–2011, MERRA-2 net heat flux has a global mean difference of -4.9 Wm(exp -2) relative to ECCO-v4. When MERRA-2 surface fields are used to force MITgcm, imbalances in the energy and the hydrological cycles of MERRA-2, which are directly related to the fact that MERRA-2 was created without an interactive ocean, propagate to the ocean. The experiment in which MITgcm is forced with MERRA-2 fluxes (MERRA-2-flux experiment) results in a 2.5°C global mean Sea Surface Temperature (SST) cooling, a 1m reduction in global mean sea level, and other drastic changes in the large scale ocean circulation relative to those resulting when the MITgcm is forced with the optimized ECCO-v4 net heat flux (the ECCO-v4 experiment itself). When MITgcm is forced with MERRA-2 state variables (MERRA-2-state experiment), the SST is somewhat restored to the observed SST, but the errors are shifted to the water cycle, resulting in a global mean sea level increase of 2.7 m. To further explore the pros and cons of these two approaches, we introduce a new intermediate forcing method in which the ocean is forced with turbulent fluxes but has a long wave feedback. This method, unlike MERRA-2 state, preserves the MERRA-2 water and salinity cycles, and it reduces the SST error compared to the MERRA-2-flux experiment, but the SST is not as good as that in the MERRA-2-state experiment. Our results have implications for ocean-model forcing recipes and clearly reveal the undesirable consequences of limiting the feedbacks in either these types of experiments or in coupled DA.

Ehud Strobach↗

Towards an Indian Land Data Assimilation System (ILDAS): A Coupled Hydrologic-Hydraulic System for Water Balance Assessments

Effective management of water resources requires reliable estimates of land surface states and fluxes, including water balance components. But most land surface models run in uncoupled mode and do not produce river discharge at catchment scales to be useful for water resources management applications. Such integrated systems are also rare over India where hydrometeorological extremes have wreaked havoc on the economy and people. So, an Indian Land Data Assimilation System (ILDAS) with a coupled land surface and a hydrodynamic model has been developed and driven by multiple meteorological forcings (0.1°, daily) to estimate land surface states, channel discharge, and floodplain inundation. ILDAS benefits from an integrated framework as well as the largest suite of observation records collected over India and has been used to produce a reanalysis product for 1981–2021 using four forcing datasets, namely, Modern-Era Retrospective Analysis for Research and Applications, Version 2 (MERRA-2), Climate Hazards Group InfraRed Precipitation with Station data (CHIRPS), ECMWF’s ERA-5, and Indian Meteorological Department (IMD) gridded precipitation. We assessed the uncertainty and bias in these precipitation datasets and validated all major components of the terrestrial water balance, i.e., surface runoff, soil moisture, terrestrial water storage anomalies, evapotranspiration, and streamflow, against a combination of satellite and in situ observation datasets. Our assessment shows that ILDAS can represent the hydrological processes reasonably well over the Indian landmass with IMD precipitation showing the best relative performance. Evaluation against ESA-CCI soil moisture shows that MERRA-2 based estimates outperform the others, whereas ERA-5 performs best in simulating evapotranspiration when evaluated against MODIS ET. Evaluations against observed records show that CHIRPS-based estimates have the highest performance in reconstructing surface runoff and streamflow. Once operational, this system will be useful for supporting transboundary water management decision making in the region.

Indian Land Data Assimilation System (ILDAS)↗

Probabilistic inversion of circular phase spectra: application to two-station phase-velocity dispersion estimation in western Canada

SUMMARY Periodic directional and temporal measurements are common in seismology, and necessitate specific statistical analyses that are appropriate for circular quantities. In this work, we explore the use of a von Mises distribution as a representation of errors on circular seismological observations. Specifically, we automate the estimation of surface-wave phase-velocity dispersion for the teleseismic two-station method, which generally suffers from a 2π phase ambiguity. The use of Bayesian inverse techniques, which aim to rigorously quantify model parameter uncertainty, have become widespread throughout seismology over the last decade. Here, we apply Bayesian inversion to measurements of surface-wave phase spectra in order to estimate 1-D, path-averaged Earth structure between station pairs. The dispersion curve and associated uncertainties are additional results of the inversion, which can then be used as input for subsequent analyses (e.g. tomography). We demonstrate this technique through application to surface-wave recordings from long-running seismic stations throughout western Canada. Our results for over 10 000 station pairs reveal first-order tectonic features consistent with previous studies, which provides confidence in our approach as well as an initial step towards resolving a full 3-D seismic velocity model for the region. This work also presents a foundation for the inversion of surface-wave phase spectra to estimate 3-D Earth structure directly. Finally, the ideas presented in this work are not limited to the inversion of surface-wave phase spectra, but can also be considered for Bayesian geophysical inversion of any circular quantities.

Gosselin, Jeremy M. (ORCID:0000000203754102)↗

Venus gravity - A high-resolution map

The Doppler data from the radio tracking of the Pioneer Venus Orbiter (PVO) have been used in a two-stage analysis to develop a high-resolution map of the gravitational potential of Venus, represented by a central mass and a surface mass density. The two-stage procedure invokes a Kalman filter-smoother to determine the orbit of the spacecraft, and a stabilized linear inverter to estimate the surface mass density. The resultant gravity map is highly correlated with the topographic map derived from the PVO radar altimeter data. However, the magnitudes of the gravity variations are smaller than would be expected if the topography were uncompensated, indicating that at least partial compensation has taken place.

Reasenberg, R. D.↗

Tropical Rainfall Distributions Determined Using TRMM Combined with other Satellite and Raingauge Information

Abstract A technique is described to use Tropical Rain Measuring Mission (TRMM) combined radar/radiometer information to adjust geosynchronous infrared satellite data (the TRMM Adjusted GOES Precipitation Index, or TRMM AGPI). The AGPI is then merged with rain gauge information (mostly over land; the TRMM merged product) to provide fine- scale (1 deg latitude/longitude) pentad and monthly analyses, respectively. The TRMM merged estimates are 10% higher than those from the Global Precipitation Climatology Project (GPCP) when integrated over the tropical oceans (37 deg N-S) for 1998, with 20% differences noted in the most heavily raining areas. In the dry subtropics the TRMM values are smaller than the GPCP estimates. The TRMM merged-product tropical-mean estimates for 1998 are 3.3 mm/ day over ocean and 3.1 mm/ day over land and ocean combined. Regional differences are noted between the western and eastern Pacific Ocean maxima when TRMM and GPCP are compared. In the eastern Pacific rain maximum the TRMM and GPCP mean values are nearly equal, very different from the other tropical rainy areas where TRMM merged-product estimates are higher. This regional difference may indicate that TRMM is better at taking in to account the vertical structure of the rain systems and the difference in structure between the western and eastern (shallower) Pacific convection. Comparisons of these TRMM merged analysis estimates with surface data sets shows varied results; the bias is near zero when compared to western Pacific Ocean atoll raingauge data, but significantly positive compared to Kwajalein radar estimates (adjusted by rain gauges). Over land the TRMM estimates also show a significant positive bias. The inclusion of gauge information in the final merged product significantly reduces the bias over land, as expected. The monthly precipitation patterns produced by the TRMM merged data process clearly show the evolution of the ENSO tropical precipitation pattern from early 1998 (El Nino) through early 1999 (La Nina) and beyond. The El Nino minus La Nina difference map shows the eastern Pacific maximum, the maritime continent minima and other tropical and mid-latitude features. The differences in the Pacific are very similar to those detected by the GPCP analyses. However, summing the El Nino minus La Nina differences over the global tropical oceans yields divergent answers from TRMM, GPCP and other estimates. This emphasizes the need for additional validation and analysis before it is feasible to understand the relations between global precipitation anomalies and Pacific Ocean ENSO temperature changes.

Adler, Robert F.↗

Photogrammetry: Develop methods for detecting deviations from expected geometry of components in a glovebox environment (Lawrence Livermore National Laboratory Aging and Lifetimes Program FY23, Milestone 8651, Grading Criterion #4)

Measurements in glovebox environments are challenging to conduct. Currently deployed techniques are time consuming and provide limited information. These limitations make current disposition and future assessments of products challenging. In FY23 we developed methods for detecting deviations of components in a glovebox environment that address these shortcomings by rapidly collecting information-dense measurements. Conducted tests demonstrate that camera pose (position and orientation) can be derived for images utilized in the photogrammetry process. This information, in conjunction with common image processing techniques, enables automated detection and size estimation of surface features. First, image processing is used to identify marks, or localized regions, that standout from the surrounding area. Then, the physical size of the region can be estimated because scale can be determined for photogrammetry images. Development of this methods shows that photogrammetry can identify measurable deviations from expected geometries. Although measurement error when using this technique through a glovebox window still needs to be reduced, this capability shows promise for eliminating or reducing the use of tedious manual measurements.

42 ENGINEERING↗

Assessment of the Impact of Discontinuity in Satellite Instruments and Retrievals on Global PM2.5

Exposure to fine particulate matter (PM2.5) is the leading environmental risk factor for mortality globally. Satellite-derived estimates of surface PM2.5 developed from a combination of satellites, simulations, and ground monitor data are relied upon for health impact studies. The ability to develop satellite-derived PM2.5 estimates requires the continued availability of aerosol optical depth (AOD) sources. This work examines the impact of the addition or loss of satellite AOD data sources on global PM2.5 estimation and the impact of continuing the longterm record with AOD from the Visible Infrared Imaging Radiometer Suite (VIIRS) on the Suomi-National Polar orbiting Partnership (S-NPP) satellite after the loss of the MODIS (MODerate resolution Imaging Spectroradiometer) and MISR (Multi-angle Imaging Spectroradiometer) instruments on board the Terra and Aqua satellites. We find that the addition of VIIRS S-NPP AOD products to geophysical PM2.5 estimates from satellites and simulations causes regional differences that correspond to differences in the VIIRS and MODIS Deep Blue AOD algorithms and sampling. Changes in long-term trends and timeseries due to the addition or loss of AOD data sources are generally within their uncertainties. Statistical fusion with ground monitor data partially corrects for changes due to sampling differences when introducing the VIIRS AOD products, but uncertainty remains over desert regions where ground monitor coverage is sparse. This work provides promise for the sustained development of global satellite-derived PM2.5 estimates, despite discontinuities in instruments and retrieval methods.

Remote sensing↗

Crater Population on Asteroid (101955) Bennu Indicates Impact Armouring and a Young Surface

The impactor-to-crater size scaling relationships that enable estimates of planetary surface ages rely on an accurate formulation of impactor–target physics. An armouring regime, specific to rubble-pile surfaces, has been proposed to occur when an impactor is comparable in diameter to a target surface particle (for example, a boulder). Armouring is proposed to reduce crater diameter, or prevent crater formation in the asteroid surface, at small crater diameters. Here, using measurements of 1,560 craters on the rubble-pile asteroid (101955) Bennu, we show that the boulder population controls a transition from crater formation to armouring at crater diameters ~2–3 m, below which crater formation in the bulk surface is increasingly rare. By combining estimates of impactor flux with the armouring scaling relationship, we find that Bennu’s crater retention age (surface age derived from crater abundance) spans from 1.6–2.2 Myr for craters less than a few meters to ~10–65 Myr for craters >100 m in diameter, reducing the maximum surface age by a factor of >15 relative to previous estimates. The range of crater retention ages, together with latitudinal variations in large-crater spatial density, indicate that ongoing resurfacing processes render the surface many times younger than the bulk asteroid.

Asteroids↗

Crater population on asteroid (101955) Bennu indicates impact armouring and a young surface

The impactor-to-crater size scaling relationships that enable estimates of planetary surface ages rely on an accurate formulation of impactor–target physics. An armouring regime, specific to rubble-pile surfaces, has been proposed to occur when an impactor is comparable in diameter to a target surface particle (for example, a boulder). Armouring is proposed to reduce crater diameter, or prevent crater formation in the asteroid surface, at small crater diameters. Here, using measurements of 1,560 craters on the rubble-pile asteroid (101955) Bennu, we show that the boulder population controls a transition from crater formation to armouring at crater diameters ~2–3m, below which crater formation in the bulk surface is increasingly rare. By combining estimates of impactor flux with the armouring scaling relationship, we find that Bennu’s crater retention age (surface age derived from crater abundance) spans from 1.6–2.2Myr for craters less than a few meters to ~10–65Myr for craters >100m in diameter, reducing the maximum surface age by a factor of >15 relative to previous estimates. The range of crater retention ages, together with latitudinal variations in large-crater spatial density, indicate that ongoing resurfacing processes render the surface many times younger than the bulk asteroid.

E. B. Bierhaus↗