Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interval analysis”

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,153 records · Page 64

Eclipse-Free-Time Assessment Tool for IRIS

IRIS_EFT is a scientific simulation that can be used to perform an Eclipse-Free- Time (EFT) assessment of IRIS (Infrared Imaging Surveyor) mission orbits. EFT is defined to be those time intervals longer than one day during which the IRIS spacecraft is not in the Earth s shadow. Program IRIS_EFT implements a special perturbation of orbital motion to numerically integrate Cowell's form of the system of differential equations. Shadow conditions are predicted by embedding this integrator within Brent s method for finding the root of a nonlinear equation. The IRIS_EFT software models the effects of the following types of orbit perturbations on the long-term evolution and shadow characteristics of IRIS mission orbits. (1) Non-spherical Earth gravity, (2) Atmospheric drag, (3) Point-mass gravity of the Sun, and (4) Point-mass gravity of the Moon. The objective of this effort was to create an in-house computer program that would perform eclipse-free-time analysis. of candidate IRIS spacecraft mission orbits in an accurate and timely fashion. The software is a suite of Fortran subroutines and data files organized as a "computational" engine that is used to accurately predict the long-term orbit evolution of IRIS mission orbits while searching for Earth shadow conditions.

Eagle, David↗

Urban and regional land use analysis: CARETS and Census Cities experiment package

The author has identified the following significant results. The arrival of the so-called energy crisis makes the portion of this experiment dealing with land use climatology of more immediate significance than before, since in addition to helping to understand the processes of climatic change associated with urbanization, the knowledge obtained may be useful in assigning an energy balance impact factor to proposed changes in land use in and around cities. Thermal maps derived from S-192 data are to be used as a measure of the energy being radiated into space from the mosaic of different surfaces in and around the city. While presenting excellent spatial sampling potential for a metropolitan area tests site, the Skylab data permit a very poor temporal sampling opportunity, owing to the large number of factors beyond the investigator's control that determine when data will be taken over a given test site. The strategy is to augment the thermal maps derived from S-192 with a modeling technique which enables the simulation of a number of components of the surface energy balance, calculated at regular time intervals throughout the day or year. Preliminary tests on the performance of the model are still underway, using airborne MSS data from NASA aircraft flights. Results look extremely promising.

Alexander, R. H.↗

Global Premature Mortality By Dust and Pollution PM 2.5 Estimated From Aerosol Reanalysis of the Modern-Era Retrospective Analysis for Research and Applications, Version 2

This study quantifies global premature deaths attributable to long-term exposure of ambient PM 2.5 , or PM 2.5 -attributable mortality, by dust and pollution sources. We used NASA’s Modern-Era Retrospective Analysis for Research and Applications, Version 2 (MERRA-2) aerosol reanalysis product for PM 2.5 and the cause-specific relative risk (RR) from the integrated exposure-response (IER) model to estimate global PM2.5-attributable mortality for five causes of deaths, namely ischaemic heart disease (IHD), cerebrovascular disease (CEV) or stroke, lung cancer (LC), chronic obstructive pulmonary disease (COPD), and acute lower respiratory infection (ALRI). The estimated yearly global PM 2.5 -attributable mortality in 2019 amounts to 2.89 (1.38–4.48) millions, which is composed of 1.19 (0.73–1.84) million from IHD, 1.01 (0.35–1.55) million from CEV, 0.29 (0.11–0.48) million from COPD, 0.23 (0.14–0.33) million from ALRI, and 0.17 (0.04–0.28) million from LC (the numbers in parentheses represent the estimated mortality range due corresponding to RR spread at the 95% confidence interval). The mortality counts vary with geopolitical regions substantially, with the highest number of deaths occurring in Asia. China and India account for 40% and 23% of the global PM 2.5 -attributable deaths, respectively. In terms of sources of PM 2.5 , about 22% of the global all-cause PM 2.5 -attributable deaths are caused by desert dust. The largest dust attribution is 37% for ALRI. The relative contributions of dust and pollution sources vary with the causes of deaths and geographical regions. Enforcing air pollution regulations to transfer areas from PM 2.5 nonattainment to PM2.5 attainment can have great health benefits. Being attainable with the United States air quality standard (AQS) of 15 μg/m 3 globally would have avoided nearly 40% or 1.2 million premature deaths. The most recent update of PM 2.5 guideline from 10 to 5 μg/m 3 by the World Health Organization (WHO) would potentially save additional one million lives. Our study highlights the importance of distinguishing aerodynamic size from geometric size in accurately assessing the global health burden of PM 2.5 and particularly for dust. A use of geometric size in diagnosing dust PM 2.5 from the model simulation, a common approach in current health burden assessment, could overestimate the PM 2.5 level in the dust belt by 40–170%, leading to an overestimate of global all-cause mortality by 1 million or 32%.

PM2.5↗

Initial Validation of NDVI time seriesfrom AVHRR, VEGETATION, and MODIS

The paper will address Theme 7: Multi-sensor opportunities for VEGETATION. We present analysis of a long-term vegetation record derived from three moderate resolution sensors: AVHRR, VEGETATION, and MODIS. While empirically based manipulation can ensure agreement between the three data sets, there is a need to validate the series. This paper uses atmospherically corrected ETM+ data available over the EOS Land Validation Core Sites as an independent data set with which to compare the time series. We use ETM+ data from 15 globally distributed sites, 7 of which contain repeat coverage in time. These high-resolution data are compared to the values of each sensor by spatially aggregating the ETM+ to each specific sensors' spatial coverage. The aggregated ETM+ value provides a point estimate for a specific site on a specific date. The standard deviation of that point estimate is used to construct a confidence interval for that point estimate. The values from each moderate resolution sensor are then evaluated with respect to that confident interval. Result show that AVHRR, VEGETATION, and MODIS data can be combined to assess temporal uncertainties and address data continuity issues and that the atmospherically corrected ETM+ data provide an independent source with which to compare that record. The final product is a consistent time series climate record that links historical observations to current and future measurements.

Morisette, Jeffrey T.↗

Method and apparatus for I-V data acquisition from solar cells

A method and apparatus for logging current-voltage (I-V) characteristic d of a solar cell module (10) in two modes using a portable instrument. One mode controls the load current through a circuit (36) in 256 equal intervals while voltage is measured from open circuit to at least halfway into the knee of the curve and the other mode controls the load voltage through a circuit (34) in 256 equal intervals from the lowest voltage measurement possible (short circuit) to at least halfway into the knee of the curve, under control of a microcomputer (12). All measurements are packed by discarding each measurement that is within 0.5% of the value predicted from two previous measurements, except every ninth (9th) measurement which is retained. The remaining data is further packed into a memory block of a detachable storage medium (14) by recording the data points in sequence following a header containing data common to all points, with each point having the value of the controlled parameter recorded as the number of increments from the previous point recorded followed by the measured value. The detachable storage medium is preferably a solid state device for reliability, and is transferable to a playback terminal which unpacks the data for analysis and display.

Cole, Steven W.↗

Analysis of Pull Force Test Results for Crimped Connections

Crimped electrical contacts are reliable when strict process controls are followed during manufacturing and accompanied by continuous process verification through pull force testing. Cable and wire harness assemblies’ standards are developed and refined over time to provide the minimum pull force that a crimp contact must meet before it breaks from the wire. However, in practice the failures occur at a much higher tensile strength than the minimum required. The first section of this paper reviewed 780 pull force test results provided by NASA Centers that were collected and analyzed to determine how the data compare to NASA’s pre-existing requirements from the cable/harness standards NASA-STD-8739.4 and IPC/WHMA-A-620B-S. The measured tensile strength of most of the contact/conductor pairs (i.e., Contact/Wires or C/W) exceeded the minimum pull force values of NASA-STD-8739.4 and IPC/WHMA-A-620 standards by at least 100%. The C/W pair samples’ tensile strength values followed a normal distribution with an average tensile strength value that was at least 182% of the minimum requirement, and all the samples analyzed passed the pull force testing. In addition, the 95% confidence interval of the average tensile strength distributions for several C/W pairs was determined and plotted as error bars to show that the C/W pairs will meet and surpass the requirements. The frequency of pull force testing can be problematic for projects because of the cost and availability of spare contacts for the destructive test. It is possible to reduce the frequency of pull force testing if at the beginning of the production run, the conditions of the crimp tool and materials are verified, and the settings of the tool remain unchanged throughout the process. However, the project needs to understand and evaluate the impact to risk from reducing the frequency of testing prior to implementing process changes.

Alejandra Constante↗

Analysis of Pull Force Test Results for Crimped Connections

Crimped electrical contacts are reliable when strict process controls are followed during manufacturing and accompanied by continuous process verification through pull force testing. Cable and wire harness assemblies’ standards are developed and refined over time to provide the minimum pull force that a crimp contact must meet before it breaks from the wire. However, in practice the failures occur at a much higher tensile strength than the minimum required. The first section of this paper reviewed 780 pull force test results provided by NASA Centers were collected and analyzed to determine how the data compares to NASA’s pre-existing requirements from the cable/harness standards NASA-STD-8739.4 and IPC/WHMA-A-620B-S. The measured tensile strength of most of the contact/conductor pairs (i.e. Contact/Wires or C/W) exceeded the minimum pull force values of NASA-STD-8739.4 and IPC/WHMA-A-620 standards by at least 100 %. The tensile strength samples of the C/W pairs values followed a normal distribution with an average tensile strength value that was at least 182 % of the minimum requirement, and all the samples analyzed passed the pull force testing. In addition, the 95 % confidence interval of the average tensile strength distributions for several C/W pairs was determined and plotted as error bars to show that there is good likelihood to observe C/W pair normal distribution with a 95 % confidence interval that will meet and surpass the standard requirements. The frequency of pull force testing can be problematic for projects because of the cost and availability of spare contacts for the destructive test. It could be possible to reduce the frequency of pull force testing if at the beginning of the production run, the conditions of the crimp tool and materials are verified, and the settings of the tool remain unchanged throughout the process. However, the project must evaluate the impact to risk from reducing the frequency of testing.

Crimped connections↗

Markov reward processes

Numerous applications in the area of computer system analysis can be effectively studied with Markov reward models. These models describe the behavior of the system with a continuous-time Markov chain, where a reward rate is associated with each state. In a reliability/availability model, upstates may have reward rate 1 and down states may have reward rate zero associated with them. In a queueing model, the number of jobs of certain type in a given state may be the reward rate attached to that state. In a combined model of performance and reliability, the reward rate of a state may be the computational capacity, or a related performance measure. Expected steady-state reward rate and expected instantaneous reward rate are clearly useful measures of the Markov reward model. More generally, the distribution of accumulated reward or time-averaged reward over a finite time interval may be determined from the solution of the Markov reward model. This information is of great practical significance in situations where the workload can be well characterized (deterministically, or by continuous functions e.g., distributions). The design process in the development of a computer system is an expensive and long term endeavor. For aerospace applications the reliability of the computer system is essential, as is the ability to complete critical workloads in a well defined real time interval. Consequently, effective modeling of such systems must take into account both performance and reliability. This fact motivates our use of Markov reward models to aid in the development and evaluation of fault tolerant computer systems.

Smith, R. M.↗

Tolerance and UQ4SIM: Nimble Uncertainty Documentation and Analysis Software

Ultimately, scientific numerical models need quantified output uncertainties so that modeling can evolve to better match reality. Documenting model input uncertainties and variabilities is a necessary first step toward that goal. Without known input parameter uncertainties, model sensitivities are all one can determine, and without code verification, output uncertainties are simply not reliable. The basic premise of uncertainty markup is to craft a tolerance and tagging mini-language that offers a natural, unobtrusive presentation and does not depend on parsing each type of input file format. Each file is marked up with tolerances and optionally, associated tags that serve to label the parameters and their uncertainties. The evolution of such a language, often called a Domain Specific Language or DSL, is given in [1], but in final form it parallels tolerances specified on an engineering drawing, e.g., 1 +/- 0.5, 5 +/- 10%, 2 +/- 10 where % signifies percent and o signifies order of magnitude. Tags, necessary for error propagation, can be added by placing a quotation-mark-delimited tag after the tolerance, e.g., 0.7 +/- 20% 'T_effective'. In addition, tolerances might have different underlying distributions, e.g., Uniform, Normal, or Triangular, or the tolerances may merely be intervals due to lack of knowledge (uncertainty). Finally, to address pragmatic considerations such as older models that require specific number-field formats, C-style format specifiers can be appended to the tolerance like so, 1.35 +/- 10U_3.2f. As an example of use, consider figure 1, where a chemical reaction input file is has been marked up to include tolerances and tags per table 1. Not only does the technique provide a natural method of specifying tolerances, but it also servers as in situ documentation of model uncertainties. This tolerance language comes with a utility to strip the tolerances (and tags), to provide a path to the nominal model parameter file. And, as shown in [1], having the ability to quickly mark and identify model parameter uncertainties facilitates error propagation, which in turn yield output uncertainties.

Kleb, Bil↗

Development and testing of a portable wind sensitive directional air sampler

A portable wind sensitive directional air sampler was developed as part of an air pollution source identification system. The system is designed to identify sources of air pollution based on the directional collection of field air samples and their analysis for TSP and trace element characteristics. Sources can be identified by analyzing the data on the basis of pattern recognition concepts. The unit, designated Air Scout, receives wind direction signals from an associated wind vane. Air samples are collected on filter slides using a standard high volume air sampler drawing air through a porting arrangement which tracks the wind direction and permits collection of discrete samples. A preset timer controls the length of time each filter is in the sampling position. At the conclusion of the sampling period a new filter is automatically moved into sampling position displacing the previous filter to a storage compartment. Thus the Air Scout may be set up at a field location, loaded with up to 12 filter slides, and left to acquire air samples automatically, according to the wind, at any timer interval desired from 1 to 30 hours.

Deyo, J.↗

Uncertainty Analysis of Slug Calorimeters in the HyMETS Arc-Jet Facility

The objective of this work is to perform an uncertainty analysis of the deduced stagnation heat flux environment on a slug calorimeter for conditions that span the performance envelope of the Hypersonic Materials Environmental Test System arc-jet facility located at NASA Langley Research Center. Analytical solutions are developed for boundary-value problems on the slug element accounting for non-ideal effects, including spatial variation in the slug heat flux, multi-dimensional thermal conduction, and back-face losses, which departs from the state-of-the-art method derived from the American Society of Testing and Materials. Boundary-value problem definitions are informed by preliminary finite element thermal analysis of the slug calorimeter assembly (including both slug and housing) and just the slug element. The analytical solutions are presented in a general sense and in a truncated form from error analysis. Results are shown in optimizing and validating the analytical models against available slug back-face thermal data. The optimization results indicate that the appropriate epistemic uncertainty of the deduced stagnation heat flux on the slug calorimeter is at most±2.5% for both a high-and low-enthalpy test condition. In addition, a numerical approach is used to determine the aleatory (probabilistic) uncertainty component in the slug stagnation heat flux by applying a marching least-squares slope routine through the steady-state portion of the slug back-face thermal response. Results indicate a compromise between the number of samples and the filter frequency of slug back-face thermal data points when evaluating the standard deviation of the deduced stagnation heat flux statistics. When combining the mixed uncertainty, both aleatory and epistemic, the interval of uncertainty in the deduced stagnation heat flux is determined to be up to ±4%, which is at least a 60% reduction from the standard uncertainty used in the state-of-the-art method.

uncertainty↗

Comparison of mode estimation methods and application in molecular clock analysis

BACKGROUND: Distributions of time estimates in molecular clock studies are sometimes skewed or contain outliers. In those cases, the mode is a better estimator of the overall time of divergence than the mean or median. However, different methods are available for estimating the mode. We compared these methods in simulations to determine their strengths and weaknesses and further assessed their performance when applied to real data sets from a molecular clock study. RESULTS: We found that the half-range mode and robust parametric mode methods have a lower bias than other mode methods under a diversity of conditions. However, the half-range mode suffers from a relatively high variance and the robust parametric mode is more susceptible to bias by outliers. We determined that bootstrapping reduces the variance of both mode estimators. Application of the different methods to real data sets yielded results that were concordant with the simulations. CONCLUSION: Because the half-range mode is a simple and fast method, and produced less bias overall in our simulations, we recommend the bootstrapped version of it as a general-purpose mode estimator and suggest a bootstrap method for obtaining the standard error and 95% confidence interval of the mode.

Evolution, Molecular↗

Analysis of large solid propellant rocket engine exhaust plumes using the direct simulation Monte Carlo method

A new solution procedure has been developed to analyze the flowfield properties in the vicinity of the Inertial Upper Stage/Spacecraft during the 1st stage (SRMI) burn. Continuum methods are used to compute the nozzle flow and the exhaust plume flowfield as far as the boundary where the breakdown of translational equilibrium leaves these methods invalid. The Direct Simulation Monte Carlo (DSMC) method is applied everywhere beyond this breakdown boundary. The flowfield distributions of density, velocity, temperature, relative abundance, surface flux density, and pressure are discussed for each species for 2 sets of boundary conditions: vacuum and freestream. The interaction of the exhaust plume and the freestream with the spacecraft and the 2-stream direct interaction are discussed. The results show that the low density, high velocity, counter flowing free-stream substantially modifies the flowfield properties and the flux density incident on the spacecraft. A freestream bow shock is observed in the data, located forward of the high density region of the exhaust plume into which the freestream gas does not penetrate. The total flux density incident on the spacecraft, integrated over the SRM1 burn interval is estimated to be of the order of 10 to the 22nd per sq m (about 1000 atomic layers).

Hueser, J. E.↗

Components Refurbishment and Chemical Analysis Facility, SWMU #041 - Per- and Polyfluoroalkyl Substances Site Assessment Report Kennedy Space Center, Florida

This PFAS Site Assessment (SA) Report presents the activities and results associated with PFAS investigations at the Components Refurbishment and Chemical Analysis (CRCA) facility located at Kennedy Space Center (KSC), Florida. In 2022, CRCA was identified as an Area of Potential Concern because the facility stores several potential PFAS-containing chemicals. A groundwater sample collected from an onsite monitoring well detected PFOA and PFOS at concentrations greater than State of Florida provisional Groundwater Cleanup Target Levels (pGCTLs). PFAS SA activities were conducted between March 2022 and February 2024. During this timeframe, a total of 16 direct-push technology (DPT) locations and 82 discrete samples were collected from these locations, along with 90 monitoring well samples. The analytical data screening process focused on State of Florida pGCTLs and United States Environmental Protection Agency (USEPA) Regional Screening Level (RSLs) from November 2023 to evaluate the data. Analytical results identified PFOS, PFOA, and PFBA at concentrations exceeding their respective RSLs at each depth interval (shallow, intermediate, and deep). PFOS and PFOA had the largest footprint of RSL exceedances in each depth interval, but pGCTL exceedances were only observed at limited locations in the shallow and intermediate intervals. PFBA had the highest detections of any PFAS compound, with concentrations exceeding 180,000 nanograms per liter (ng/L), which is 100-times the RSL of 1,800 ng/L in the shallow and intermediate intervals and 10-times the RSL in the deep interval. The maximum PFBA detection was 841,000 ng/L at shallow monitoring well, MW0006. The PFAS SA also included samples collected from the onsite hydraulic containment system (HCS) which was installed to control and treat the onsite chlorinated volatile organic compound (CVOC) plume. Influent and effluent aqueous samples were collected monthly from the system. PFAS concentrations were relatively the same for both influent and effluent samples, indicating that while the HCS has been effective for CVOC treatment, it does not provide any additional treatment for PFAS compounds. Since the HCS has achieved its objectives, the system was shut down in December 2024. PFAS data gaps still exist, to include surface water and soil, which were not sampled during this SA. Soil sampling near the Chemical Process Area is recommended. Surface water and soil sampling at select stormwater outfalls is also recommended. Additionally, further groundwater sampling is recommended (DPT and monitoring well) in all depth intervals to delineate the extent of PFAS impacts at CRCA.

K Alex Murphy↗

Roles of divergent and rotational winds in the kinetic energy balance during intense convective activity

Contributions of divergent and rotational wind components to the synoptic-scale kinetic energy balance are described using rawinsonde data at 3 and 6 h intervals from NASA's fourth Atmospheric Variability experiment. Two intense thunderstorm complexes occurred during the period. Energy budgets are described for the entire computational region and for limited volumes that enclosed storm-induced, upper level wind maxima located poleward of convection. Although small in magnitude, the divergent wind component played an important role in the cross-contour generation and horizontal flux divergence of kinetic energy. The importance of V(D) appears directly related to the presence and intensity of convection. Although K(D) usually comprised less than 10 percent of the total kinetic energy content, generation of kinetic energy by V(D) was a major factor in the creation of upper-level wind maxima to the north of the storm complexes. Omission of the divergent wind apparently would lead to serious misrepresentations of the energy balance. A random error analysis is presented to assess confidence limits in the various energy parameters.

Fuelberg, H. E.↗

Space Station Furnace Facility Management Information System (SSFF-MIS) Development

This report summarizes the chronology, results, and lessons learned from the development of the SSFF-MIS. This system has been nearly two years in development and has yielded some valuable insights into specialized MIS development. General: In December of 1994, the Camber Corporation and Science Applications International Corporation (SAIC) were contracted to design, develop, and implement a MIS for Marshall Space Flight Center's Space Station Furnace Facility Project. The system was to be accessible from both EBM-Compatible PC and Macintosh platforms. The system was required to contain data manually entered into the MIS as well as data imported from other MSFC sources. Electronic interfaces were established for each data source and retrieval was to be performed at prescribed time intervals. The SOW requirement that predominantly drove the development software selection was the dual-platform (IBM-PC and Macintosh) requirement. The requirement that the system would be maintained by Government personnel influenced the selection of Commercial Off-the-shelf software because of its inherent stability and readily available documentation and support. Microsoft FoxPro Professional 2.6 for Windows and Macintosh was selected as the development tool. This is a software development tool that has been in use for many years. It is stable and powerful. Microsoft has since released the replacement for this product, Microsoft Visual FoxPro, but at the time of this development, it was only available on the Windows platform. The initial contract included included the requirement for capabilities relating to the Work- and Organizational Breakdown Structures, cost (plan and actuals), workforce (plan and actuals), critical path scheduling, trend analysis, procurements and contracts, interface to manufacturing, Safety and Mission Assurance, risk analysis, and technical performance indicators. It also required full documentation of the system and training of users. During the course of the contract, the requirements for Safety and Mission Assurance interface, risk analysis, and technical performance indicators were deleted. Additional capabilities were added as reflected in the Contract Chronology below. Modification 4 added the requirement for Support Contractor manpower data, the ability to manually input data not imported from non-nal sources, a general 'health' indicator screen, and remote usage. Mod 6 included the ability to change the level of planning of Civil Service Manpower at any time and the ability to manually enter Op Codes in the manufacturing data where such codes were not provided by the EMPACS database. Modification 9 included a number of changes to report contents and formats. Modification 11 required the preparation of a detailed System Design Document.

Meade, Robert M.↗

The soft x ray telescope for Solar-A

The Solar-A satellite being prepared by the Institute for Sapce and Astronautical Sciences (ISAS) in Japan is dedicated to high energy observations of solar flares. The Soft X Ray Telescope (SXT) is being prepared to provide filtered images in the 2 to 60 A interval. The flight model is now undergoing tests in the 1000 foot tunnel at MSFC. Launch will be in September 1991. Earlier resolution and efficiency tests on the grazing incidence mirror have established its performance in soft x rays. The one-piece, two mirror grazing incidence telescope is supported in a strain free mount separated from the focal plane assembly by a carbon-epoxy metering tube whose windings and filler are chosen to minimize thermal and hygroscopic effects. The CCD detector images both the x ray and the concentric visible light aspect telescope. Optical filters provide images at 4308 and 4700 A. The SXT will be capable of producing over 8000 of the smallest partial frame images per day, or fewer but larger images, up to 1024 x 1024 pixel images. Image sequence with two or more of the five x ray analysis filters, with automatic exposure compensation to optimize the charge collection by the CCD detector, will be used to provide plasma diagnostics. Calculations using a differential emission measure code were used to optimize filter selection over the range of emission measure variations and to avoid redundancy, but the filters were chosen primarily to give ratios that are monotonic in plasma temperature.

Brown, W. A.↗

An Inversion Analysis of Recent Variability in Natural CO2 Fluxes Using GOSAT and In Situ Observations

About one-half of the global CO2 emissions from fossil fuel combustion and deforestation accumulates in the atmosphere, where it contributes to global warming. The rest is taken up by vegetation and the ocean. The precise contribution of the two sinks, and their location and year-to-year variability are, however, not well understood. We use two different approaches, batch Bayesian synthesis inversion and variational data assimilation, to deduce the global spatiotemporal distributions of CO2 fluxes during 2009-2010. One of our objectives is to assess different sources of uncertainties in inferred fluxes, including uncertainties in prior flux estimates and observations, and differences in inversion techniques. For prior constraints, we utilize fluxes and uncertainties from the CASA-GFED model of the terrestrial biosphere and biomass burning driven by satellite observations and interannually varying meteorology. We also use measurement-based ocean flux estimates and two sets of fixed fossil CO2 emissions. Here, our inversions incorporate column CO2 measurements from the GOSAT satellite (ACOS retrieval, filtered and bias-corrected) and in situ observations (individual flask and afternoon-average continuous observations) to estimate fluxes in 108 regions over 8-day intervals for the batch inversion and at 3 x 3.75 weekly for the variational system. Relationships between fluxes and atmospheric concentrations are derived consistently for the two inversion systems using the PCTM atmospheric transport model driven by meteorology from the MERRA reanalysis. We compare the posterior fluxes and uncertainties derived using different data sets and the two inversion approaches, and evaluate the posterior atmospheric concentrations against independent data including aircraft measurements. The optimized fluxes generally resemble those from other studies. For example, the results indicate that the terrestrial biosphere is a net CO2 sink, and a GOSAT-only inversion suggests a shift in the global sink from the tropics south to the north relative to the prior and to an in-situ-only inversion. We also find a smaller terrestrial sink in higher-latitude northern regions in boreal summer of 2010 relative to 2009.

Carbon dioxide↗