Search NASA⌕ Search

SEARCH · Search NASA

Results for “atmospheric reconstruction”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Improvements to the Simulated CubeSpark Satellite Constellation and Their Effects on Lightning Geolocation Accuracy from Orbit

The CubeSpark mission concept is being developed as a constellation of up to six satellites in low-Earth orbit (LEO) for high resolution 3D lightning mapping using optical and radio frequency (RF) sensors. Individual lightning VHF signals are simulated from Earth’s atmosphere through the ionosphere to each satellite, using their arrival times to reconstruct source locations. Here we present recent updates to these simulations based on improved ionospheric modeling, with a focus on the expected three-dimensional resolution. These studies include testing from the equator up to high latitudes, with varying vertical total electron content (vTEC), and using between one and six orbiting stations. In addition to the more robust ionosphere model, the constellation formation has also been updated to reduce its resulting errors and increase the effective range of VHF geolocation from space. The goals of CubeSpark include mapping thundercloud charge structure as well as lightning channel lengths relevant to climatology, meteorology, and more. These applications require location uncertainty less than 1-2 km in each dimension. This improved algorithm shows sufficient resolution up to high latitudes, including significantly larger areas having 3D resolution less than 1 km. Analysis of the distributions of biases in simulated arrival times has also revealed the unexpected relationship between the shape of those distributions and the resulting uncertainties. This work helps to refine our understanding of the sources of error in lightning geolocation and reinforces the potential for post-processing improvements in this and other similar systems.

Lightning↗
NASA NTRS2012

Calibration and Image Reconstruction for the Hurricane Imaging Radiometer (HIRAD)

The Hurricane Imaging Radiometer (HIRAD) is a new airborne passive microwave synthetic aperture radiometer designed to provide wide swath images of ocean surface wind speed under heavy precipitation and, in particular, in tropical cyclones. It operates at 4, 5, 6 and 6.6 GHz and uses interferometric signal processing to synthesize a pushbroom imager in software from a low profile planar antenna with no mechanical scanning. HIRAD participated in NASA s Genesis and Rapid Intensification Processes (GRIP) mission during Fall 2010 as its first science field campaign. HIRAD produced images of upwelling brightness temperature over a aprox 70 km swath width with approx 3 km spatial resolution. From this, ocean surface wind speed and column averaged atmospheric liquid water content can be retrieved across the swath. The calibration and image reconstruction algorithms that were used to verify HIRAD functional performance during and immediately after GRIP were only preliminary and used a number of simplifying assumptions and approximations about the instrument design and performance. The development and performance of a more detailed and complete set of algorithms are reported here.

Ruf, Christopher↗
NASA NTRS2005

Orbit determination results and trajectory reconstruction for the Cassini/Huygens Mission

During Cassini's third orbit around Saturn, the Huygens Probe was successfully released on a trajectory that resulted in the probe entering Titan's atmosphere on January 14, 2005, making it both the most distant spacecraft landing and the first spacecraft to successfully land on the moon of another planet. This paper documents the reconstruction of both the orbiter and probe trajectoriespanning the Titan-B and Titan-C encounters.

reconstruction↗
DOE OSTI2021

Deep-Learning based Reconstruction of the Shower Maximum $X_{\mathrm{max}}$ using the Water-Cherenkov Detectors of the Pierre Auger Observatory

The atmospheric depth of the air shower maximum X max is an observable commonly used for the determination of the nuclear mass composition of ultra-high energy cosmic rays. Direct measurements of X max are performed using observations of the longitudinal shower development with fluorescence telescopes. At the same time, several methods have been proposed for an indirect estimation of X max from the characteristics of the shower particles registered with surface detector arrays. In this paper, we present a deep neural network (DNN) for the estimation of X max . The reconstruction relies on the signals induced by shower particles in the ground based water-Cherenkov detectors of the Pierre Auger Observatory. The network architecture features recurrent long short-term memory layers to process the temporal structure of signals and hexagonal convolutions to exploit the symmetry of the surface detector array. We evaluate the performance of the network using air showers simulated with three different hadronic interaction models. Thereafter, we account for long-term detector effects and calibrate the reconstructed X max using fluorescence measurements. Finally, we show that the event-by-event resolution in the reconstruction of the shower maximum improves with increasing shower energy and reaches less than 25 g/cm 2 at energies above 2 × 10 19 eV.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Reconstruction of events recorded with the surface detector of the Pierre Auger Observatory

Cosmic rays arriving at Earth collide with the upper parts of the atmosphere, thereby inducing extensive air showers. When secondary particles from the cascade arrive at the ground, they are measured by surface detector arrays. We describe the methods applied to the measurements of the surface detector of the Pierre Auger Observatory to reconstruct events with zenith angles less than $60^\circ$ using the timing and signal information recorded using the water-Cherenkov detector stations. In addition, we assess the accuracy of these methods in reconstructing the arrival directions of the primary cosmic ray particles and the sizes of the induced showers.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS1992

Fractionation of carbon isotopes by phytoplankton and estimates of ancient CO2 levels

Reports of the 13C content of marine particulate organic carbon are compiled and on the basis of GEOSECS data and temperatures, concentrations, and isotopic compositions of dissolved CO2 in the waters in which the related phytoplankton grew are estimated. In this way, the fractionation of carbon isotopes during photosynthetic fixation of CO2 is found to be significantly correlated with concentrations of dissolved CO2. Because ancient carbon isotopic fractionations have been determined from analyses of sedimentary porphyrins [Popp et al., 1989], the relationship between isotopic fractionation and concentrations of dissolved CO2 developed here can be employed to estimate concentrations of CO2 dissolved in ancient oceans and, in turn, partial pressures of CO2 in ancient atmospheres. The calculations take into account the temperature dependence of chemical and isotopic equilibria in the dissolved-inorganic-carbon system and of air-sea equilibria. Paleoenvironmental temperatures for each sample are estimated from reconstructions of paleogeography, latitudinal temperature gradients, and secular changes in low-latitude sea surface temperature. It is estimated that atmospheric partial pressures of CO2 were over 1000 micro atm 160 - 100 Ma ago, then declined to values near 300 micro atm during the next 100 Ma. Analysis of a high-resolution record of carbon isotopic fractionation at the Cenomanian-Turonian boundary suggests that the partial pressure of CO2 in the atmosphere was drawn down from values near 840 micro atm to values near 700 micro atm during the anoxic event.

NASA Discipline Exobiology↗
DOE OSTI2026

Neutron reconstruction via blips in liquid argon time projection chambers

Neutrons are important final-state particles in neutrino interactions, yet they are not considered or reconstructed in most current neutrino liquid argon time projection chamber (LArTPC) physics analyses. In this paper, we present a simulation-based proof-of-concept study of neutron reconstruction in a generic LArTPC detector. Leveraging isolated, MeV-scale energy deposits, or blips, from neutron inelastic scattering, and using realistic blip response from published experimental results, we demonstrate the capability to identify neutrons and to reconstruct the direction and energy of the final-state neutron system in sub-GeV neutrino interactions. We then explore how neutron-related blip attributes can be used to improve physics studies of neutrino interactions, such as enhancing neutrino-antineutrino separation in atmospheric neutrinos and reverse-horn-current beam neutrinos. This study provides an initial quantification of LArTPC neutron reconstruction capabilities, which we expect to improve with future advancements in blip reconstruction, identification, and classification algorithms, as well as the modeling of neutrons.

Morquecho, Miguel Hernandez [Univ. of Minnesota, M↗
NASA NTRS2021

Reduced El NiñO Variability in the Mid-Pliocene According to the Pliomip2 Ensemble

The mid-Pliocene warm period (3.264–3.025 Ma) is the most recent geological period during which atmospheric CO 2 levels were similar to recent historical values (∼400 ppm). Several proxy reconstructions for the mid-Pliocene show highly reduced zonal sea surface temperature (SST) gradients in the tropical Pacific Ocean, indicating an El Niño-like mean state. However, past modelling studies do not show these highly reduced gradients. Efforts to understand mid-Pliocene climate dynamics have led to the Pliocene Model Intercomparison Project (PlioMIP). Results from the first phase (PlioMIP1) showed clear El Niño variability (albeit significantly reduced) and did not show the greatly reduced time-mean zonal SST gradient suggested by some of the proxies. In this work, we study El Niño–Southern Oscillation (ENSO) variability in the PlioMIP2 ensemble, which consists of additional global coupled climate models and updated boundary conditions compared to PlioMIP1. We quantify ENSO amplitude, period, spatial structure and “flavour”, as well as the tropical Pacific annual mean state in mid-Pliocene and pre-industrial simulations. Results show a reduced ENSO amplitude in the model-ensemble mean (−24 %) with respect to the pre-industrial, with 15 out of 17 individual models showing such a reduction. Furthermore, the spectral power of this variability considerably decreases in the 3–4-year band. The spatial structure of the dominant empirical orthogonal function shows no particular change in the patterns of tropical Pacific variability in the model-ensemble mean, compared to the pre-industrial. Although the time-mean zonal SST gradient in the equatorial Pacific decreases for 14 out of 17 models (0.2 ∘C reduction in the ensemble mean), there does not seem to be a correlation with the decrease in ENSO amplitude. The models showing the most “El Niño-like” mean state changes show a similar ENSO amplitude to that in the pre-industrial reference, while models showing more “La Niña-like” mean state changes generally show a large reduction in ENSO variability. The PlioMIP2 results show a reasonable agreement with both time-mean proxies indicating a reduced zonal SST gradient and reconstructions indicating a reduced, or similar, ENSO variability.

mid-Pliocene warm period↗
NASA NTRS1972

Profile inversion in presence of ray bending

Inversion of radio occultation data for planetary atmospheres and ionospheres has been performed using the seismological Herglotz-Wiechert method, as adapted by Phinney and Anderson to the radio-occultation case. Profile reconstruction performed in computer simulated experiments with this approach have been compared with the ones obtained with the straight-ray Abel transform. For a thin atmosphere and ionosphere, like the ones encountered on Mars, microwave occultation data can be inverted accurately with both methods. For a dense ionosphere like the sun's corona, ray bending of microwaves is severe, and recovered refractivity by the Herglotz-Wiechert method provides significant improvement over the straight-ray Abel transform: the error reduces from more than 60% to less than 20% at a height of 60,000 km above the base of the corona.

Wallio, H. A.↗
NASA NTRSDate not supplied

NASA Space Launch System Artemis I & II Post Flight Ascent Aerothermal Environments Overview

Since 2011 the Aerosciences Branch/EV33 at NASA Marshall Space Flight Center has been involved with the development of ascent external aerothermal environments for the NASA Space Launch System (SLS) Block 1 launch vehicle for the purposes of supporting thermal analysis and the design of thermal protection systems. The SLS Block 1 Artemis I and II launch vehicles successfully launched from Pad39B at NASA Kennedy Space Center on November 16th, 2022 and April 1st, 2026, respectively. Over 70 aerothermal islands, consisting of over 265 operational instruments captured aerodynamic heating and plume induced environments throughout the launch vehicles. Gauges consisted of calorimeters, radiometers, gas temperature probes, pressure transducers, bi-directional pressure probes and thermocouples. Prior to launch, aerothermal design environment models were generated to predict ascent aerodynamic heating and plume induced environments over a design space that covered a range of vehicle trajectories that varied atmospheric, vehicle performance, and off-nominal, engine-out conditions. Post flight reconstruction models were developed for each flight island using the Day-of-Launch (DOL) Best Equivalent Trajectory (BET) that provided freestream conditions and propulsion system boundary conditions. This paper discusses a summary of the ascent aerothermal environments observed during the flights and the respective modelling approaches and the performance of them through comparisons of flight data and predictions.

aerothermodynamics↗
DOE OSTI2020

Studies on the response of a water-Cherenkov detector of the Pierre Auger Observatory to atmospheric muons using an RPC hodoscope

Extensive air showers, originating from ultra-high energy cosmic rays, have been successfully measured through the use of arrays of water-Cherenkov detectors (WCDs). Sophisticated analyses exploiting WCD data have made it possible to demonstrate that shower simulations, based on different hadronic-interaction models, cannot reproduce the observed number of muons at the ground. The accurate knowledge of the WCD response to muons is paramount in establishing the exact level of this discrepancy. In this work, we report on a study of the response of a WCD of the Pierre Auger Observatory to atmospheric muons performed with a hodoscope made of resistive plate chambers (RPCs), enabling us to select and reconstruct nearly 600 thousand single muon trajectories with zenith angles ranging from 0$^\circ$ to 55$^\circ$. Comparison of distributions of key observables between the hodoscope data and the predictions of dedicated simulations allows us to demonstrate the accuracy of the latter at a level of 2%. As the WCD calibration is based on its response to atmospheric muons, the hodoscope data are also exploited to show the long-term stability of the procedure.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Paleoclimatic implications of glacial fluctuations in the Sierra Nevada del Cocuy, northern Andes, Colombia, during the Lateglacial and Holocene

The reconstruction of former mountain glaciers from geomorphic mapping and cosmogenic-nuclide surface-exposure dating provides a unique opportunity to infer patterns of past terrestrial climate variability. Tropical mountain glaciers are particularly valuable as there are comparatively few terrestrial climate proxies at equatorial latitudes relative to higher latitudes. As the single largest climate zone on Earth, the tropics play an outsized role in mediating global climate via the ocean-atmosphere transfer of latent heat and water vapor. Nonetheless, there remains a persistent gap in our understanding of how the tropics influenced – or were influenced by – the high-magnitude climate shifts of the Late Pleistocene, and whether this high-energy region simply responded to extratropical forcing or was itself a driver of global climatic change. To help address this knowledge gap, we analyzed geologic evidence for past glacial fluctuations in three adjacent valleys in the Sierra Nevada del Cocuy, the highest subrange of the Eastern Cordillera in the Colombian Andes, to provide a terrestrial record of atmospheric temperature during the latter part of Termination 1. Coupled with geomorphic mapping and paleo-snowline reconstructions, our beryllium-10 glacial chronology indicates that glaciers in the humid inner tropics underwent pronounced growth and gradual decay during the Antarctic Cold Reversal (14.5–12.8 ka) and Younger Dryas (12.8–11.7 ka) periods, respectively, following a trend that, according to directly dated moraine records from throughout both polar hemispheres, appears to have been global. While the specific mechanism(s) behind this large-scale behavior remains to be corroborated, we revisit the hypothesis that ocean-atmosphere heat transfer and water vapor flux are key drivers of abrupt Lateglacial temperature fluctuations. Subsequent to the Lateglacial, deglaciation of the Sierra Nevada del Cocuy accelerated during the Early Holocene, a pattern also observed in other tropical glacier records. More recently, the magnitude of snowline rise and glacier retreat over the last two centuries supports the view that modern tropospheric warming is anomalously strong at least relative to the last ∼16,000 years.

Andes↗
NASA NTRS1992

Modelling the Phanerozoic carbon cycle and climate: constraints from the 87Sr/86Sr isotopic ratio of seawater

A numerical model describing the coupled evolution of the biogeochemical cycles of carbon, sulfur, calcium, magnesium, phosphorus, and strontium has been developed to describe the long-term changes of atmospheric carbon dioxide and climate during the Phanerozoic. The emphasis is on the effects of coupling the cycles of carbon and strontium. Various interpretations of the observed Phanerozoic history of the seawater 87Sr/86Sr ratio are investigated with the model. More specifically, the abilities of continental weathering, volcanism, and surface lithology in generating that signal are tested and compared. It is suggested that the observed fluctuations are mostly due to a changing weatherability over time. It is shown that such a conclusion is very important for the modelling of the carbon cycle. Indeed, it implies that the conventional belief that the evolution of atmospheric carbon dioxide and climate on a long time scale is governed by the balance between the volcanic input of CO2 and the rate of silicate weathering is not true. Rather carbon exchanges between the mantle and the exogenic system are likely to have played a key role too. Further, the increase of the global weathering rates with increasing surface temperature and/or atmospheric CO2 pressure usually postulated in long-term carbon cycle and climate modelling is also inconsistent with the new model. Other factors appear to have modulated the weatherability of the continents through time, such as mountain building and the existence of glaciers and ice sheets. Based on these observations, a history of atmospheric carbon dioxide and climate during Phanerozoic time, consistent with the strontium isotopic data, is reconstructed with the model and is shown to be compatible with paleoclimatic indicators, such as the timing of glaciation and the estimates of Cretaceous paleotemperatures.

NASA Discipline Number 52-30↗
NASA NTRSDate not supplied

Mars 2020 Reconstructed Aerothermal Environments and Design Margins

The Mars 2020 aeroshell's thermal protection system was nearly identical to the Mars Science Laboratory system that successfully completed its mission for the Curiosity rover's landing in 2012. It was predicted that, like Mars Science Laboratory, the Mars 2020 heatshield would experience boundary layer transition and that the thermal protection system would provide sufficient material thickness margins against the aerothermal environments, even after adding radiative heating that was not included for Mars Science Laboratory. The Mars 2020 flight instrumentation suite included heatshield sensors similar to Mars Science Laboratory and new backshell measurements. A full set of surface pressure and in-depth temperature data were collected during atmospheric entry. This paper provides an initial analysis of the flight data and focuses on the reconstructed total surface heat flux inferred from the measured temperatures. Turbulent boundary layer conditions again were observed on the Mars 2020 heatshield, and the temperatures at all eleven heatshield and six backshell thermocouple in-depth sensor locations were well within system capabilities due to favorable entry conditions and conservative uncertainties. New computational fluid dynamics results on the reconstructed entry trajectory are compared to the measured surface pressures and reconstructed total heat fluxes. The predicted heatshield surface pressures at six locations match the data qualitatively and quantitatively well, as expected. Smooth-wall laminar heating predictions prior to boundary layer transition fall above the reconstructed heating on the heatshield. After the observed boundary layer transition time, total heat flux based on algebraic turbulence model calculations generally match the heat flux trends relative to one another. The convective heat flux was predicted to constitute the majority of the total heating. On the backshell, smooth-wall laminar total heat flux predictions generally exceed the reconstructed total heating at six locations. At each of these locations, it is estimated that radiative heating provided almost all of the total heating. The Mars 2020 as-flown aerothermal environments were well below the design levels for all measurement locations due to conservative design assumptions and a stressing design trajectory.

Mars 2020↗
NASA NTRSDate not supplied

Mars 2020 Reconstructed Aerothermal Environments and Design Margins

The Mars 2020 aeroshell's thermal protection system was nearly identical to the Mars Science Laboratory system that successfully completed its mission for the Curiosity rover's landing in 2012. It was predicted that, like Mars Science Laboratory, the Mars 2020 heatshield would experience boundary layer transition and that the thermal protection system would provide sufficient material thickness margins against the aerothermal environments, even after adding radiative heating that was not included for Mars Science Laboratory. The Mars 2020 flight instrumentation suite included heatshield sensors similar to Mars Science Laboratory and new backshell measurements. A full set of surface pressure and in-depth temperature data were collected during atmospheric entry. This paper provides an initial analysis of the flight data and focuses on the reconstructed total surface heat flux inferred from the measured temperatures. Turbulent boundary layer conditions again were observed on the Mars 2020 heatshield, and the temperatures at all eleven heatshield and six backshell thermocouple in-depth sensor locations were well within system capabilities due to favorable entry conditions and conservative uncertainties. New computational fluid dynamics results on the reconstructed entry trajectory are compared to the measured surface pressures and reconstructed total heat fluxes. The predicted heatshield surface pressures at six locations match the data qualitatively and quantitatively well, as expected. Smooth-wall laminar heating predictions prior to boundary layer transition fall above the reconstructed heating on the heatshield. After the observed boundary layer transition time, total heat flux based on algebraic turbulence model calculations generally match the heat flux trends relative to one another. The convective heat flux was predicted to constitute the majority of the total heating. On the backshell, smooth-wall laminar total heat flux predictions generally exceed the reconstructed total heating at six locations. At each of these locations, it is estimated that radiative heating provided almost all of the total heating. The Mars 2020 as-flown aerothermal environments were well below the design levels for all measurement locations due to conservative design assumptions and a stressing design trajectory.

Mars 2020↗
NASA NTRS2021

Evaluation of CloudSat Radiative Kernels Using ARM and CERES Observations and ERA5 Reanalysis

Despite the widespread use of the radiative kernel technique for studying radiative feedbacks and radiative forcings, there has not been any systematic, observation-based validation of the radiative kernel method. Here, we utilize observed and reanalyzed radiative fluxes and atmospheric profiles from the Atmospheric Radiation Measurement (ARM) program and ERA5 reanalysis to assess a set of observation-based radiative kernels from CloudSat for six ARM sites. The CloudSat radiative kernels, convoluted with the ERA5 state variables, can almost perfectly reconstruct the monthly anomalies of shortwave (SW) and longwave (LW) radiative fluxes in ERA5 at the surface (SFC) and top-of atmosphere (TOA) with correlations significantly being greater than 0.95. The biases of kernel-estimated flux anomalies calculated using the ARM-observed state variables can be more than twice as large when compared with the ARM-observed surface flux anomalies and Clouds and Earth’s Radiant Energy System (CERES) observed anomalies at the TOA. Generally, clouds contribute to most (>60%) of the variance of flux anomalies at Southern Great Plain (SGP), Tropical Western Pacific (TWP), and Eastern North Atlantic (ENA), and surface albedo dominates (>69%) the variance of SW flux anomalies at North Slope of Alaska (NSA). The radiative kernels exhibit the lowest correlation (r~[0.55,0.85]) when reconstructing SFC LW flux anomalies at SGP, TWP, and ENA, whose biases are related to the possibility that the kernels may not fully capture the characteristics 47 associated with MJO and ENSO at TWP and the presence of clouds at SGP and ENA.

Radiative kernels↗
DOE OSTI2021

Evaluation of CloudSat Radiative Kernels Using ARM and CERES Observations and ERA5 Reanalysis

Despite the widespread use of the radiative kernel technique for studying radiative feedbacks and radiative forcings, there has not been any systematic, observation-based validation of the radiative kernel method. Here, we utilize observed and reanalyzed radiative fluxes and atmospheric profiles from the Atmospheric Radiation Measurement (ARM) program and ERA5 reanalysis to assess a set of observation-based radiative kernels from CloudSat for six ARM sites. The CloudSat radiative kernels, convoluted with the ERA5 state variables, can almost perfectly reconstruct the monthly anomalies of shortwave (SW) and longwave (LW) radiative fluxes in ERA5 at the surface (SFC) and top-of-atmosphere (TOA) with correlations significantly being greater than 0.95. The biases of kernel-estimated flux anomalies calculated using the ARM-observed state variables can be more than twice as large when compared with the ARM-observed surface flux anomalies and Clouds and Earth’s Radiant Energy System (CERES) observed anomalies at the TOA. Generally, clouds contribute to most (>60%) of the variance of flux anomalies at Southern Great Plain (SGP), Tropical Western Pacific (TWP), and Eastern North Atlantic (ENA), and surface albedo dominates (>69%) the variance of SW flux anomalies at North Slope of Alaska (NSA). Furthermore, the radiative kernels exhibit the lowest correlation (r ~ [0.55,0.85]) when reconstructing SFC LW flux anomalies at SGP, TWP, and ENA, whose biases are related to the possibility that the kernels may not fully capture the characteristics associated with MJO and ENSO at TWP and the presence of clouds at SGP and ENA.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1990

Atmospheric drag perturbations on the Mars Observer Orbiter

The nominal launch date of the Mars Observer spacecraft is 09 - 16 - 02, with Mars encounter on 08 - 19 - 93 and the mapping orbit established on 12 - 06 - 93. Because this orbit is nearly circular at low altitude (361 km), both gravity and atmospheric drag will dominate the orbital evolution and are also the major error sources. This paper provides: (1) a summary of atmospheric density and errors for the Mars Observer orbit and time frame, (2) an assessment of spacecraft-position errors using current density estimates, and (3) an estimate of atmospheric density modeling improvement directly from the Doppler data acquired during flight operations. With an improved gravity field model, reconstructed spacecraft-position errors are at the level of 0.5 km (1sigma). However, when predicting the spacecraft's position over 14 days, the density error (at the 90-percent confidence level) dominates, yielding a down-track error of almost 70 km. By modeling the mean atmospheric density, this error can be reduced by at least a factor of two.

Esposito, P. B.↗