Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Uncertainties”

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,387 records · Page 77

ERRATUM: FERMI Large Area Telescope Study of Cosmic-Rays and the Interstellar Medium in Nearby Molecular Clouds

In the published version of the paper, errors were made in calculating the exposure time due to an analysis mistake. While they do not affect gas emissivities of the R CrA and Cepheus & Polaris flare regions significantly (the differences are within the systematic uncertainty), that of the Chamaeleon region is increased by approx.20%. Although we claimed a difference of ∼50% in gas emissivity among these molecular cloud regions in the original paper, it is decreased to ∼30% (comparable to the sum of the statistical and systematic uncertainties) in the revised analysis. Therefore, our conclusion of the original paper, that a small variation (approx. 20%) of the CR density in the solar neighborhood exists, is not supported by the data if we take these uncertainties into account. On the other hand, the obtained XCO and XAv values, and the masses of gas calculated from them are not changed significantly (the differences are within the statistical errors). Errors and corrections in the original paper are summarized below. 1. In the Abstract (lines 5-6) and Section 3 (lines 4-5 in the 3rd paragraph) in the original paper, the gamma -ray emissivity above 250 MeV for the Chamaeleon region should be (7.2 +/- 0.1stat +/- 1.0sys) × 10(exp −27) photons/s/sr/H-atom, not (5.9 +/-0.1stat +0.9−1.0sys) × 10(exp −27) photons/s/sr/H-atom. 2. In the Abstract (lines 8-10), "Whereas the energy dependences of the emissivities agree well with that predicted from direct CR observations at the Earth, the measured emissivities from 250 MeV to 10 GeV indicate a variation of the CR density by approx.20% in the neighborhood of the solar system, even if we consider the systematic uncertainties." should be changed to "The energy dependences of the emissivities agree well with that predicted from direct CR observations at the Earth. Although the measured emissivities from 250 MeV to 10 GeV differ by approx.30% among these molecular cloud regions, the difference is not significant if we take the systematic uncertainty into account." 3. Table 1 and Figure 13, which show gas emissivities and spectra for the Chamaeleon region in the original paper, should be changed to the Table 1 and Figure 1 as shown below. 4. Figure 16, which compares Hi gas emissivities among several regions in the original paper, should be changed to Figure 2 as shown below. 5. The text from the line 13 to the last one in the first paragraph of Section 4.1, "The spectral shapes for the three regions..., indicating a difference of the CR density between the Chamaeleon and the others as shown in Figure 16." should be changed to the paragraph that follows. "The shaded area of each spectrum indicates the systematic uncertainty as described in Section 3. We note that the systematic uncertainty of the LAT effective area (5% at 100 MeV and 20% at 10 GeV; Rando et al. 2009) does not affect the relative value of emissivities. The effect of unresolved point sources is small; we have verified that the obtained emissivities are almost unaffected by decreasing the threshold for point sources from TS = 100 to TS = 50. We also confirmed that the residual excess of photons around (l = 280deg to 288deg, b = −20deg to −12deg; see the bottom panel of Figure 8) in the Chamaeleon region does not affect the local Hi emissivity very much. Thus the total systematic uncertainty is reasonably expressed by the shaded area shown in Fig. 1.

emissivities↗

Attribution of Flux Partitioning Variations Between Land Surface Models over the Continental U.S.

Accurate quantification of the terrestrial evapotranspiration (ET) components of plant transpiration (T), soil evaporation (E) and evaporation of the intercepted water (I) is necessary for improving our understanding of the links between the carbon and water cycles. Recent studies have noted that, among the modeled estimates, large disagreements exist in the relative contributions of T, E and I to the total ET. As these models are often used in data assimilation environments for incorporating and extending ET relevant remote sensing measurements, understanding the sources of inter-model differences in ET components is also necessary for improving the utilization of such remote sensing measurements. This study quantifies the contributions of two key factors explaining inter-model disagreements to the uncertainty in total ET: (1) contribution of the local partitioning and (2) regional distribution of ET. The analysis is conducted by using outputs from a suite of land surface models in the North American Land Data Assimilation System (NLDAS) configuration. For most of these models, transpiration is the dominant component of the ET partition. The results indicate that the uncertainty in local partitioning dominates the inter-model spread in modeled soil evaporation E. The inter-model differences in T are dominated by the uncertainty in the distribution of ET over the Eastern U.S. and the local partitioning uncertainty in the Western U.S. The results also indicate that uncertainty in the T estimates is the primary driver of total ET errors. Over the majority of the U.S., the contribution of the two factors of uncertainty to the overall uncertainty is non-trivial.

uncertainty↗

Measurements for liquid rocket engine performance code verification

The goal of the rocket engine performance code verification tests is to obtain the I sub sp with an accuracy of 0.25% or less. This needs to be done during the sequence of four related tests (two reactive and two hot gas simulation) to best utilize the loss separation technique recommended in this study. In addition to I sub sp, the measurements of the input and output parameters for the codes are needed. This study has shown two things in regard to obtaining the I sub sp uncertainty within the 0.25% target. First, this target is generally not being realized at the present time, and second, the instrumentation and testing technology does exist to obtain this 0.25% uncertainty goal. However, to achieve this goal will require carefully planned, designed, and conducted testing. In addition, the test-stand (or system) dynamics must be evaluated in the pre-test and post-test phases of the design of the experiment and data analysis, respectively always keeping in mind that a .25% overall uncertainty in I sub sp is targeted. A table gives the maximum allowable uncertainty required for obtaining I sub sp with 0.25% uncertainty, the currently-quoted instrument specification, and present test uncertainty for the parameters. In general, it appears that measurement of the mass flow parameter within the required uncertainty may be the most difficult.

Praharaj, Sarat C.↗

NOVA 201 ultrasonic thickness gage (NOVA Gage)

The measurement integrity of the NOVA 201 digital ultrasonic thickness gage (NOVA gage) was demonstrated by comparing the NOVA gage measurements to the thickness gage measurements, and determining the bias and uncertainty of the NOVA gage when measuring redesigned solid rocket motor (RSRM) hardware per engineering test plans (ETP). The NOVA gage was tested by three different operators on steel and aluminum RSRM hardware for wall thickness. The results show that the measurement bias is not consistent. The uncertainty of the bias is caused by the heterogeneous material properties of the RSRM components that influence the time of flight of ultrasonic waves. The measurement uncertainty inherent to the design and operation of the NOVA gage is less in comparison to the uncertainty of the bias. The total measurement uncertainty cannot be substantially reduced by taking more than one measurement. There is no correlation between bias and the surface finish range of this test unless 3-in-One oil is used as a couplant, in which case there appears to be a slight trend. There is no correlation between uncertainty and the surface finish range. The measurement uncertainty of the NOVA gage can be reduced using 3-in-One oil as a couplant.

Garecht, Diane↗

Asteroid approach covariance analysis for the Clementine mission

The Clementine mission is designed to test Strategic Defense Initiative Organization (SDIO) technology, the Brilliant Pebbles and Brilliant Eyes sensors, by mapping the moon surface and flying by the asteroid Geographos. The capability of two of the instruments available on board the spacecraft, the lidar (laser radar) and the UV/Visible camera is used in the covariance analysis to obtain the spacecraft delivery uncertainties at the asteroid. These uncertainties are due primarily to asteroid ephemeris uncertainties. On board optical navigation reduces the uncertainty in the knowledge of the spacecraft position in the direction perpendicular to the incoming asymptote to a one-sigma value of under 1 km, at the closest approach distance of 100 km. The uncertainty in the knowledge of the encounter time is about 0.1 seconds for a flyby velocity of 10.85 km/s. The magnitude of these uncertainties is due largely to Center Finding Errors (CFE). These systematic errors represent the accuracy expected in locating the center of the asteroid in the optical navigation images, in the absence of a topographic model for the asteroid. The direction of the incoming asymptote cannot be estimated accurately until minutes before the asteroid flyby, and correcting for it would require autonomous navigation. Orbit determination errors dominate over maneuver execution errors, and the final delivery accuracy attained is basically the orbit determination uncertainty before the final maneuver.

Ionasescu, Rodica↗

SeaWiFs Technical Report Series. Volume 34: The Third SeaWiFS Intercalibration Round-Robin Experiment (SIRREX-3), 19-30 September 1994

This report presents results of the third Sea-viewing Wide Field-of-view Sensor (SeaWiFS) Intercalibration Round- Robin Experiment (SIRREX-3), which was held at the San Diego State University (SDSU) Center for Hydro-Optics and Remote Sensing (CHORS) on 19-30 September 1994. Spectral irradiances of FEL lamps belonging to each participant were intercompared by reference to the National Institute of Standards and Technology (NIST) scale of spectral irradiance using secondary standard lamps F268, F269, and F182, with a Type A uncertainty between 1.1-1.5%. This level of uncertainty was achieved despite difficulties with lamp F269. The average spectral irradiances of FEL lamps, compared in both SIRREX-2 and SIRREX-3, differed between the two experiments by 1.5%, which probably indicates that the values assigned to the secondary standard lamp at the time of SIRREX-2 were in error. With two exceptions, spectral radiance values of integrating sphere sources were measured during SIRREX-3 with uncertainties in temporal stability of less than 0.3% and absolute uncertainties of 1.5-2.0%. This is a significant improvement over similar intercomparisons in SIRREX- I and SIRREX-2. Plaque reflectances were intercompared with an uncertainty of about 1-2%, but the absolute uncertainty is undefined. Although this is an improvement over results of previous SIRREXS, the sources and magnitude of uncertainty associated with transfers of spectral radiance using plaques requires further evaluation in future experiments.

Hooker, Stanford B.↗

Space Radiation Cancer Risks

Space radiation presents major challenges to astronauts on the International Space Station and for future missions to the Earth s moon or Mars. Methods used to project risks on Earth need to be modified because of the large uncertainties in projecting cancer risks from space radiation, and thus impact safety factors. We describe NASA s unique approach to radiation safety that applies uncertainty based criteria within the occupational health program for astronauts: The two terrestrial criteria of a point estimate of maximum acceptable level of risk and application of the principle of As Low As Reasonably Achievable (ALARA) are supplemented by a third requirement that protects against risk projection uncertainties using the upper 95% confidence level (CL) in the radiation cancer projection model. NASA s acceptable level of risk for ISS and their new lunar program have been set at the point-estimate of a 3-percent risk of exposure induced death (REID). Tissue-averaged organ dose-equivalents are combined with age at exposure and gender-dependent risk coefficients to project the cumulative occupational radiation risks incurred by astronauts. The 95% CL criteria in practice is a stronger criterion than ALARA, but not an absolute cut-off as is applied to a point projection of a 3% REID. We describe the most recent astronaut dose limits, and present a historical review of astronaut organ doses estimates from the Mercury through the current ISS program, and future projections for lunar and Mars missions. NASA s 95% CL criteria is linked to a vibrant ground based radiobiology program investigating the radiobiology of high-energy protons and heavy ions. The near-term goal of research is new knowledge leading to the reduction of uncertainties in projection models. Risk projections involve a product of many biological and physical factors, each of which has a differential range of uncertainty due to lack of data and knowledge. The current model for projecting space radiation cancer risk relies on the three assumptions of linearity, additivity, and scaling along with the use of population averages. We describe uncertainty estimates for this model, and new experimental data that sheds light on the accuracy of the underlying assumptions. These methods make it possible to express risk management objectives in terms of quantitative metrics, i.e., the number of days in space without exceeding a given risk level within well defined confidence limits. The resulting methodology is applied to several human space exploration mission scenarios including lunar station, deep space outpost, and a Mars mission. Factors that dominate risk projection uncertainties and application of this approach to assess candidate mitigation approaches are described.

Cucinotta, Francis A.↗

Adapting Covariance Propagation to Account for the Presence of Modeled and Unmodeled Maneuvers

This paper explores techniques that can be used to adapt the standard linearized propagation of an orbital covariance matrix to the case where there is a maneuver and an associated execution uncertainty. A Monte Carlo technique is used to construct a final orbital covariance matrix for a 'prop-burn-prop' process that takes into account initial state uncertainty and execution uncertainties in the maneuver magnitude. This final orbital covariance matrix is regarded as 'truth' and comparisons are made with three methods using modified linearized covariance propagation. The first method accounts for the maneuver by modeling its nominal effect within the state transition matrix but excludes the execution uncertainty by omitting a process noise matrix from the computation. The second method does not model the maneuver but includes a process noise matrix to account for the uncertainty in its magnitude. The third method, which is essentially a hybrid of the first two, includes the nominal portion of the maneuver via the state transition matrix and uses a process noise matrix to account for the magnitude uncertainty. The first method is unable to produce the final orbit covariance except in the case of zero maneuver uncertainty. The second method yields good accuracy for the final covariance matrix but fails to model the final orbital state accurately. Agreement between the simulated covariance data produced by this method and the Monte Carlo truth data fell within 0.5-2.5 percent over a range of maneuver sizes that span two orders of magnitude (0.1-20 m/s). The third method, which yields a combination of good accuracy in the computation of the final covariance matrix and correct accounting for the presence of the maneuver in the nominal orbit, is the best method for applications involving the computation of times of closest approach and the corresponding probability of collision, PC. However, applications for the two other methods exist and are briefly discussed. Although the process model ("prop-burn-prop") that was studied is very simple - point-mass gravitational effects due to the Earth combined with an impulsive delta-V in the velocity direction for the maneuver - generalizations to more complex scenarios, including high fidelity force models, finite duration maneuvers, and maneuver pointing errors, are straightforward and are discussed in the conclusion.

Schiff, Conrad↗

Orbit Determination of Spacecraft in Earth-Moon L1 and L2 Libration Point Orbits

The ARTEMIS mission, part of the THEMIS extended mission, is the first to fly spacecraft in the Earth-Moon Lissajous regions. In 2009, two of the five THEMIS spacecraft were redeployed from Earth-centered orbits to arrive in Earth-Moon Lissajous orbits in late 2010. Starting in August 2010, the ARTEMIS P1 spacecraft executed numerous stationkeeping maneuvers, initially maintaining a lunar L2 Lissajous orbit before transitioning into a lunar L1 orbit. The ARTEMIS P2 spacecraft entered a L1 Lissajous orbit in October 2010. In April 2011, both ARTEMIS spacecraft will suspend Lissajous stationkeeping and will be maneuvered into lunar orbits. The success of the ARTEMIS mission has allowed the science team to gather unprecedented magnetospheric measurements in the lunar Lissajous regions. In order to effectively perform lunar Lissajous stationkeeping maneuvers, the ARTEMIS operations team has provided orbit determination solutions with typical accuracies on the order of 0.1 km in position and 0.1 cm/s in velocity. The ARTEMIS team utilizes the Goddard Trajectory Determination System (GTDS), using a batch least squares method, to process range and Doppler tracking measurements from the NASA Deep Space Network (DSN), Berkeley Ground Station (BGS), Merritt Island (MILA) station, and United Space Network (USN). The team has also investigated processing of the same tracking data measurements using the Orbit Determination Tool Kit (ODTK) software, which uses an extended Kalman filter and recursive smoother to estimate the orbit. The orbit determination results from each of these methods will be presented and we will discuss the advantages and disadvantages associated with using each method in the lunar Lissajous regions. Orbit determination accuracy is dependent on both the quality and quantity of tracking measurements, fidelity of the orbit force models, and the estimation techniques used. Prior to Lissajous operations, the team determined the appropriate quantity of tracking measurements that would be needed to meet the required orbit determination accuracies. Analysts used the Orbit Determination Error Analysis System (ODEAS) to perform covariance analyses using various tracking data schedules. From this analysis, it was determined that 3.5 hours of DSN TRK-2-34 range and Doppler tracking data every other day would suffice to meet the predictive orbit knowledge accuracies in the Lissajous region. The results of this analysis are presented. Both GTDS and ODTK have high-fidelity environmental orbit force models that allow for very accurate orbit estimation in the lunar Lissajous regime. These models include solar radiation pressure, Earth and Moon gravity models, third body gravitational effects from the Sun, and to a lesser extent third body gravitational effects from Jupiter, Venus, Saturn, and Mars. Increased position and velocity uncertainties following each maneuver, due to small execution performance errors, requires that several days of post-maneuver tracking data be processed to converge on an accurate post-maneuver orbit solution. The effects of maneuvers on orbit determination accuracy will be presented, including a comparison of the batch least squares technique to the extended Kalman filter/smoother technique. We will present the maneuver calibration results derived from processing post-maneuver tracking data. A dominant error in the orbit estimation process is the uncertainty in solar radiation pressure and the resultant force on the spacecraft. An estimation of this value can include many related factors, such as the uncertainty in spacecraft reflectivity and surface area which is a function of spacecraft orientation (spin-axis attitude), uncertainty in spacecraft wet mass, and potential seasonal variability due to the changing direction of the Sun line relative to the Earth-Moon Lissajous reference frame. In addition, each spacecraft occasionally enters into Earth or Moon penumbra or umbra and these shadow crossings reduche solar radiation force for several hours. The effects of these events on orbit determination accuracy will be presented. In order to plan for upcoming stationkeeping maneuvers, the maneuver planning team must take the current orbit estimate, propagate it forward to the planned maneuver time, and determine the optimal maneuver to maintain the Lissajous orbit for one or more revolutions. The propagation is performed using a Runge-Kutta 7/8 integrator and typically the position and velocity uncertainty increases with propagation time, increasing the overall uncertainty of the orbit state at the maneuver execution time. The effect of orbit knowledge uncertainty on stationkeeping operations will be presented.

Woodard, Mark↗

Chance-Constrained Guidance With Non-Convex Constraints

Missions to small bodies, such as comets or asteroids, require autonomous guidance for descent to these small bodies. Such guidance is made challenging by uncertainty in the position and velocity of the spacecraft, as well as the uncertainty in the gravitational field around the small body. In addition, the requirement to avoid collision with the asteroid represents a non-convex constraint that means finding the optimal guidance trajectory, in general, is intractable. In this innovation, a new approach is proposed for chance-constrained optimal guidance with non-convex constraints. Chance-constrained guidance takes into account uncertainty so that the probability of collision is below a specified threshold. In this approach, a new bounding method has been developed to obtain a set of decomposed chance constraints that is a sufficient condition of the original chance constraint. The decomposition of the chance constraint enables its efficient evaluation, as well as the application of the branch and bound method. Branch and bound enables non-convex problems to be solved efficiently to global optimality. Considering the problem of finite-horizon robust optimal control of dynamic systems under Gaussian-distributed stochastic uncertainty, with state and control constraints, a discrete-time, continuous-state linear dynamics model is assumed. Gaussian-distributed stochastic uncertainty is a more natural model for exogenous disturbances such as wind gusts and turbulence than the previously studied set-bounded models. However, with stochastic uncertainty, it is often impossible to guarantee that state constraints are satisfied, because there is typically a non-zero probability of having a disturbance that is large enough to push the state out of the feasible region. An effective framework to address robustness with stochastic uncertainty is optimization with chance constraints. These require that the probability of violating the state constraints (i.e., the probability of failure) is below a user-specified bound known as the risk bound. An example problem is to drive a car to a destination as fast as possible while limiting the probability of an accident to 10(exp -7). This framework allows users to trade conservatism against performance by choosing the risk bound. The more risk the user accepts, the better performance they can expect.

FROM↗

Robust Control Design for Uncertain Nonlinear Dynamic Systems

Robustness to parametric uncertainty is fundamental to successful control system design and as such it has been at the core of many design methods developed over the decades. Despite its prominence, most of the work on robust control design has focused on linear models and uncertainties that are non-probabilistic in nature. Recently, researchers have acknowledged this disparity and have been developing theory to address a broader class of uncertainties. This paper presents an experimental application of robust control design for a hybrid class of probabilistic and non-probabilistic parametric uncertainties. The experimental apparatus is based upon the classic inverted pendulum on a cart. The physical uncertainty is realized by a known additional lumped mass at an unknown location on the pendulum. This unknown location has the effect of substantially altering the nominal frequency and controllability of the nonlinear system, and in the limit has the capability to make the system neutrally stable and uncontrollable. Another uncertainty to be considered is a direct current motor parameter. The control design objective is to design a controller that satisfies stability, tracking error, control power, and transient behavior requirements for the largest range of parametric uncertainties. This paper presents an overview of the theory behind the robust control design methodology and the experimental results.

Kenny, Sean P.↗

Robustness Analysis and Optimally Robust Control Design via Sum-of-Squares

A control analysis and design framework is proposed for systems subject to parametric uncertainty. The underlying strategies are based on sum-of-squares (SOS) polynomial analysis and nonlinear optimization to design an optimally robust controller. The approach determines a maximum uncertainty range for which the closed-loop system satisfies a set of stability and performance requirements. These requirements, de ned as inequality constraints on several metrics, are restricted to polynomial functions of the uncertainty. To quantify robustness, SOS analysis is used to prove that the closed-loop system complies with the requirements for a given uncertainty range. The maximum uncertainty range, calculated by assessing a sequence of increasingly larger ranges, serves as a robustness metric for the closed-loop system. To optimize the control design, nonlinear optimization is used to enlarge the maximum uncertainty range by tuning the controller gains. Hence, the resulting controller is optimally robust to parametric uncertainty. This approach balances the robustness margins corresponding to each requirement in order to maximize the aggregate system robustness. The proposed framework is applied to a simple linear short-period aircraft model with uncertain aerodynamic coefficients.

Dorobantu, Andrei↗

A Statistics-Based Material Property Analysis to Support TPS Characterization

Accurate characterization of entry capsule heat shield material properties is a critical component in modeling and simulating Thermal Protection System (TPS) response in a prescribed aerothermal environment. The thermal decomposition of the TPS material during the pyrolysis and charring processes is poorly characterized and typically results in large uncertainties in material properties as inputs for ablation models. These material property uncertainties contribute to large design margins on flight systems and cloud re- construction efforts for data collected during flight and ground testing, making revision to existing models for entry systems more challenging. The analysis presented in this work quantifies how material property uncertainties propagate through an ablation model and guides an experimental test regimen aimed at reducing these uncertainties and characterizing the dependencies between properties in the virgin and charred states for a Phenolic Impregnated Carbon Ablator (PICA) based TPS. A sensitivity analysis identifies how the high-fidelity model behaves in the expected flight environment, while a Monte Carlo based uncertainty propagation strategy is used to quantify the expected spread in the in-depth temperature response of the TPS. An examination of how perturbations to the input probability density functions affect output temperature statistics is accomplished using a Kriging response surface of the high-fidelity model. Simulations are based on capsule configuration and aerothermal environments expected during the Mars Science Laboratory (MSL) entry sequence. We identify and rank primary sources of uncertainty from material properties in a flight-relevant environment, show the dependence on spatial orientation and in-depth location on those uncertainty contributors, and quantify how sensitive the expected results are.

Copeland, Sean R.↗

A Regional CO2 Observing System Simulation Experiment for the ASCENDS Satellite Mission

Top-down estimates of the spatiotemporal variations in emissions and uptake of CO2 will benefit from the increasing measurement density brought by recent and future additions to the suite of in situ and remote CO2 measurement platforms. In particular, the planned NASA Active Sensing of CO2 Emissions over Nights, Days, and Seasons (ASCENDS) satellite mission will provide greater coverage in cloudy regions, at high latitudes, and at night than passive satellite systems, as well as high precision and accuracy. In a novel approach to quantifying the ability of satellite column measurements to constrain CO2 fluxes, we use a portable library of footprints (surface influence functions) generated by the WRF-STILT Lagrangian transport model in a regional Bayesian synthesis inversion. The regional Lagrangian framework is well suited to make use of ASCENDS observations to constrain fluxes at high resolution, in this case at 1 degree latitude x 1 degree longitude and weekly for North America. We consider random measurement errors only, modeled as a function of mission and instrument design specifications along with realistic atmospheric and surface conditions. We find that the ASCENDS observations could potentially reduce flux uncertainties substantially at biome and finer scales. At the 1 degree x 1 degree, weekly scale, the largest uncertainty reductions, on the order of 50 percent, occur where and when there is good coverage by observations with low measurement errors and the a priori uncertainties are large. Uncertainty reductions are smaller for a 1.57 micron candidate wavelength than for a 2.05 micron wavelength, and are smaller for the higher of the two measurement error levels that we consider (1.0 ppm vs. 0.5 ppm clear-sky error at Railroad Valley, Nevada). Uncertainty reductions at the annual, biome scale range from 40 percent to 75 percent across our four instrument design cases, and from 65 percent to 85 percent for the continent as a whole. Our uncertainty reductions at various scales are substantially smaller than those from a global ASCENDS inversion on a coarser grid, demonstrating how quantitative results can depend on inversion methodology. The a posteriori flux uncertainties we obtain, ranging from 0.01 to 0.06 Pg C yr-1 across the biomes, would meet requirements for improved understanding of long-term carbon sinks suggested by a previous study.

remote sensing↗

Model Sensitivity Studies of the Decrease in Atmospheric Carbon Tetrachloride

Carbon tetrachloride (CCl4) is an ozone-depleting substance, which is controlled by the Montreal Protocol and for which the atmospheric abundance is decreasing. However, the current observed rate of this decrease is known to be slower than expected based on reported CCl4 emissions and its estimated overall atmospheric lifetime. Here we use a three-dimensional (3-D) chemical transport model to investigate the impact on its predicted decay of uncertainties in the rates at which CCl4 is removed from the atmosphere by photolysis, by ocean uptake and by degradation in soils. The largest sink is atmospheric photolysis (74% of total), but a reported 10% uncertainty in its combined photolysis cross section and quantum yield has only a modest impact on the modelled rate of CCl4 decay. This is partly due to the limiting effect of the rate of transport of CCl4 from the main tropospheric reservoir to the stratosphere, where photolytic loss occurs. The model suggests large interannual variability in the magnitude of this stratospheric photolysis sink caused by variations in transport. The impact of uncertainty in the minor soil sink (9%of total) is also relatively small. In contrast, the model shows that uncertainty in ocean loss (17%of total) has the largest impact on modelled CCl4 decay due to its sizeable contribution to CCl4 loss and large lifetime uncertainty range (147 to 241 years). With an assumed CCl4 emission rate of 39 Gg year(exp -1), the reference simulation with the best estimate of loss processes still underestimates the observed CCl4 (overestimates the decay) over the past 2 decades but to a smaller extent than previous studies. Changes to the rate of CCl4 loss processes, in line with known uncertainties, could bring the model into agreement with in situ surface and remote-sensing measurements, as could an increase in emissions to around 47 Gg year(exp -1). Further progress in constraining the CCl4 budget is partly limited by systematic biases between observational datasets. For example, surface observations from the National Oceanic and Atmospheric Administration (NOAA) network are larger than from the Advanced Global Atmospheric Gases Experiment (AGAGE) network but have shown a steeper decreasing trend over the past 2 decades. These differences imply a difference in emissions which is significant relative to uncertainties in the magnitudes of the CCl4 sinks.

Carbon Trechloride↗

Attribution of Flux Partitioning Variations Between Land Surface Models Over the Continental U.S.

Accurate quantification of the terrestrial evapotranspiration (ET) components of plant transpiration (T), soil evaporation (E) and evaporation of the intercepted water (I) is necessary for improving our understanding of the links between the carbon and water cycles. Recent studies have noted that, among the modeled estimates, large disagreements exist in the relative contributions of T, E and I to the total ET. As these models are often used in data assimilation environments for incorporating and extending ET relevant remote sensing measurements, understanding the sources of inter-model differences in ET components is also necessary for improving the utilization of such remote sensing measurements. This study quantifies the contributions of two key factors explaining inter-model disagreements to the uncertainty in total ET: (1) contribution of the local partitioning and (2) regional distribution of ET. The analysis is conducted by using outputs from a suite of land surface models in the North American Land Data Assimilation System (NLDAS) configuration. For most of these models, transpiration is the dominant component of the ET partition. The results indicate that the uncertainty in local partitioning dominates the inter-model spread in modeled soil evaporation E. The inter-model differences in T are dominated by the uncertainty in the distribution of ET over the Eastern U.S. and the local partitioning uncertainty in the Western U.S. The results also indicate that uncertainty in the T estimates is the primary driver of total ET errors. Over the majority of the U.S., the contribution of the two factors of uncertainty to the overall uncertainty is non-trivial.

evapotranspiration partitioning; transpiration; so↗

HZETRN Radiation Transport Validation Using Balloon-Based Experimental Data

The deterministic radiation transport code HZETRN (High charge (Z) and Energy TRaNsport) was developed by NASA to study the effects of cosmic radiation on astronauts and instrumentation shielded by various materials. This work presents an analysis of computed differential flux from HZETRN compared with measurement data from three balloon-based experiments over a range of atmospheric depths, particle types, and energies. Model uncertainties were quantified using an interval-based validation metric that takes into account measurement uncertainty both in the flux and the energy at which it was measured. Average uncertainty metrics were computed for the entire dataset as well as subsets of the measurements (by experiment, particle type, energy, etc.) to reveal any specific trends of systematic over- or under-prediction by HZETRN. The distribution of individual model uncertainties was also investigated to study the range and dispersion of errors beyond just single scalar and interval metrics. The differential fluxes from HZETRN were generally well-correlated with balloon-based measurements; the median relative model difference across the entire dataset was determined to be 30%. The distribution of model uncertainties, however, revealed that the range of errors was relatively broad, with approximately 30% of the uncertainties exceeding ± 40%. The distribution also indicated that HZETRN systematically under-predicts the measurement dataset as a whole, with approximately 80% of the relative uncertainties having negative values. Instances of systematic bias for subsets of the data were also observed, including a significant underestimation of alpha particles and protons for energies below 2.5 GeV/u. Muons were found to be systematically over-predicted at atmospheric depths deeper than 50 g/cm(sup 2) but under-predicted for shallower depths. Furthermore, a systematic under-prediction of alpha particles and protons was observed below the geomagnetic cutoff, suggesting that improvements to the light ion production cross sections in HZETRN should be investigated.

James E Warner↗

Terrestrial Single-Station Analog for Constraining the Martian Core and Deep Interior: Implications for InSight

We used a terrestrial single-station seismometer to quantify the uncertainty of InSight (INterior explorations using Seismic Investigations, Geodesy and Heat Transport) data for determining Martian core size. To mimic Martian seismicity, we formed a catalog using 917 terrestrial earthquakes, from which we randomly selected events. We stacked ScS amplitudes on modeled arrival times and searched for where ScS produced coherent seismic amplitudes. A core detection was defined by a coherent peak with small offset between predicted and user-selected arrival times. Iterating the detection algorithm with varying signal-to-noise (SNR) ranges and quantity of events determined the selection frequency of each model and quantified core depth uncertainty. Increasing the quantity of events reduced core depth uncertainty while increasing the recovery rate, while increasing event SNR had little effect. Including ScS2 multiples increased the recovery rate and reduced core depth uncertainty when we used low quantities of events. The most-frequent core depths varied by back azimuth, suggesting our method is sensitive to the presence of mantle heterogeneities. When we added 1° in source distance errors, core depth uncertainty increased by up to 11 km and recovery rates decreased by <5%. Altering epicentral distances by 25% added ~35 km of uncertainty and reduced recovery rates to <50% in some cases. From these experiments, we estimate that if InSight can detect five events with high location precision (<10% epicentral distance errors), that there is at least an 88% chance of core depth recovery using ScS alone with uncertainty in core depth approaching 18 km and decreasing as more events are located.

Marusiak, Angela G.↗