Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Statistical error”

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,459 records · Page 81
All sourcesNASA NTRSDOE OSTI
NASA NTRS2009

Predicting the Size of Sunspot Cycle 24 on the Basis of Single- and Bi-Variate Geomagnetic Precursor Methods

Examined are single- and bi-variate geomagnetic precursors for predicting the maximum amplitude (RM) of a sunspot cycle several years in advance. The best single-variate fit is one based on the average of the ap index 36 mo prior to cycle minimum occurrence (E(Rm)), having a coefficient of correlation (r) equal to 0.97 and a standard error of estimate (se) equal to 9.3. Presuming cycle 24 not to be a statistical outlier and its minimum in March 2008, the fit suggests cycle 24 s RM to be about 69 +/- 20 (the 90% prediction interval). The weighted mean prediction of 11 statistically important single-variate fits is 116 +/- 34. The best bi-variate fit is one based on the maximum and minimum values of the 12-mma of the ap index; i.e., APM# and APm*, where # means the value post-E(RM) for the preceding cycle and * means the value in the vicinity of cycle minimum, having r = 0.98 and se = 8.2. It predicts cycle 24 s RM to be about 92 +/- 27. The weighted mean prediction of 22 statistically important bi-variate fits is 112 32. Thus, cycle 24's RM is expected to lie somewhere within the range of about 82 to 144. Also examined are the late-cycle 23 behaviors of geomagnetic indices and solar wind velocity in comparison to the mean behaviors of cycles 2023 and the geomagnetic indices of cycle 14 (RM = 64.2), the weakest sunspot cycle of the modern era.

Wilson, Robert M.↗
NASA NTRS2014

Regional Characteristics of the Relationship Between Columnar AOD and Surface PM2.5: Application of Lidar Aerosol Extinction Profiles over Baltimore-Washington Corridor During DISCOVER-AQ

The first field campaign of DISCOVER-AQ (Deriving Information on Surface conditions from COlumn and VERtically resolved observations relevant to Air Quality) took place in July 2011 over Baltimore-Washington Corridor (BWC). A suite of airborne remote sensing and in-situ sensors was deployed along with ground networks for mapping vertical and horizontal distribution of aerosols. Previous researches were based on a single lidar station because of the lack of regional coverage. This study uses the unique airborne HSRL (High Spectral Resolution Lidar) data to baseline PM2.5 (particulate matter of aerodynamic diameter less than 2.5 μm) estimates and applies to regional air quality with satellite AOD (Aerosol Optical Depth) retrievals over BWC (∼6500 sq. km). The linear approximation takes into account aerosols aloft above AML (Aerosol Mixing Layer) by normalizing AOD with haze layer height (i.e., AOD/HLH). The estimated PM2.5 mass concentrations by HSRL AOD/HLH are shown within 2 RMSE (Root Mean Square Error ∼9.6 μg/cu. m) with correlation ∼0.88 with the observed over BWC. Similar statistics are shown when applying HLH data from a single location over the distance of 100 km. In other words, a single lidar is feasible to cover the range of 100 km with expected uncertainties. The employment of MPLNET-AERONET (MicroPulse Lidar NETwork - AErosol RObotic NETwork) measurements at NASA GSFC produces similar statistics of PM2.5 estimates as those derived by HSRL. The synergy of active and passive remote sensing aerosol measurements provides the foundation for satellite application of air quality on a daily basis. For the optimal range of 10 km, the MODIS-estimated PM2.5 values are found satisfactory at 27 (out of 36) sunphotometer locations with mean RMSE of 1.6-3.3 μg/cu. m relative to PM2.5 estimated by sunphotometers. The remaining 6 of 8 marginal sites are found in the coastal zone, for which associated large RMSE values ∼4.5-7.8 μg/cu. m are most likely due to overestimated AOD because of water-contaminated pixels.

DISCOVER-AQ↗
NASA NTRS2019

High Reliability Requires More than Providing Spares

It is sometimes optimistically hoped that a space life support system can be kept working throughout a long duration mission by repairing failed components, as long as sufficient spares are flown. It is usually assumed that the components have constant known failure rates. Then the needed numbers of spares can be computed to have any particular probability that all failed components can be replaced by available spares. This approach can provide high reliability if its favorable assumptions, including constant known failure rates, are satisfied. Other favorable assumptions are that the failures are statistically independent, repair will be successful without causing further failures, and all failures are due to internal component failures. These assumptions are not usually justified. The failure rates may be estimates that are inadequately verified because of insufficient testing. Failure rates may change due to materials substitutions, manufacturing changes, redesigns to fix failures, and new failures caused by redesigns. Failures that are not statistically independent may result from one common cause, such as a design or manufacturing error or a cascade of cause and effect, possibly caused by an external event such as a power outage. Repair may be unsuccessful or cause damage. Many failures occur at component interfaces or at the overall systems level, not within isolated components. Other failures causes are completely external to the system, due to assembly, maintenance, and operational errors or to unexpected environmental challenges. Replacement with sufficient spares can compensate for expected internal component failures but may not be able to cope with unpredictable design and manufacturing flaws, human errors, and environmental impacts. Reliability estimates based on providing sufficient spares to compensate for expected failures may be far too high. They are essentially upper bounds on reliability that might be approached if many frequent but often unconsidered failure causes can be eliminated.

spares↗
NASA NTRS2021

Evaluation of VIIRS and MODIS snow cover fraction in High-Mountain Asia using Landsat 8 OLI

We present thefirst application of the Snow Covered Area and Grain size model (SCAG) tothe Visible Infrared imaging Radiometer Suite (VIIRS) and assess these retrievals withfiner-resolution fractional snow cover maps from Landsat 8 Operational Land Imager (OLI).Because Landsat 8 OLI avoids saturation issues common to Landsat 1–7 in the visiblewavelengths, we re-assess the accuracy of the SCAG fractional snow cover maps fromModerate Resolution Imaging Spectroradiometer (MODIS) that were previously evaluatedusing data from earlier Landsat sensors. Use of the fractional snow cover maps fromLandsat 8 OLI shows a negative bias of−0.5% for MODSCAG and−1.3% for VIIRSCAG,whereas previous MODSCAG evaluations found a bias of−7.6% in the Himalaya. Wefindsimilar root mean squared error (RMSE) values of 0.133 and 0.125 for MODIS and VIIRS,respectively. The Recall statistic (probability of detection) for cells with more than 15%snow cover in this challenging steep topography was found to be 0.90 for both MODSCAGand VIIRSCAG, significantly higher than previous evaluations based on Landsat 5Thematic Mapper (TM) and 7 Enhanced Thematic Mapper Plus (ETM+). In addition,daily retrievals from MODIS and VIIRS are consistent across gradients of elevation, slope,and aspect. Different native resolutions of the gridded products at 1 km and 500 m forVIIRS and MODIS, respectively, result in snow cover maps showing a slightly differentdistribution of values with VIIRS having more mixed pixels and MODIS having 7% morepure snow pixels. Despite the resolution differences, the snow maps from both sensorsproduce similar total snow-covered areas and snow-line elevations in this region, withR2values of 0.98 and 0.88, respectively. Wefind that the SCAG algorithm performsconsistently across various spatial resolutions and that fractional snow cover mapsfrom the VIIRS instruments aboard Suomi NPP, JPPS–1, and JPPS–2 can be asuitable replacement as MODIS sensors reach their ends of life.

Karl Rittger↗
DOE OSTI2019

Application of Markov Chain Monte Carlo for Uncertainty Quantification in Quantitative Imaging Problems

The Differential Evolution Adaptive Metropolis (DREAM) method, a Markov chain Monte Carlo approach, is used for optimization and uncertainty analysis of a radioactive source/shield system using pixelated measured data generated by radiation imagers, where each pixel in the image contains a gamma-ray spectrum with statistical uncertainty. DREAM uses this spectral information to determine source thickness and source strength while also propagating uncertainty from the measured data to the solution. Using measurements simulated with GEANT4, successful parameter reconstruction is demonstrated on a numerical test case in a slab/shield geometry. In the presence of statistical noise of 5-9%, parameters are calculated to better than 95% with 2σ error bars that generally encompass the actual values. In a test problem with real-world measurements, source strength per energy line is calculated to within 78-97% of the actual value.

Bledsoe, Keith↗
DOE OSTI2020

Contaminant Transport Modeling for Technology Evaluation and Long-Term Monitoring in the Tims Branch Testbed, SC - 20343

Studies conducted at U.S. DOE sites have shown the presence of heavy metals, radionuclides, and volatile organic compounds in surface water, groundwater, and soil as a result of nuclear activity in the Cold War era. Since the 1990's, innovative cleanup methods have been implemented in the Tims Branch watershed at Savannah River Site (SRS) to limit the contaminant flux to the stream that have reduced the contaminant concentrations to acceptable regulatory levels in the dissolved phase. A tin-based treatment which effectively eliminated all local anthropogenic mercury inputs to this ecosystem resulted in a known step function addition of inert tin oxide particles which now serve as a potential tracer for sedimentation and particle transport processes in the stream. The long-term effectiveness of this and other remediation techniques and the potential for remobilization of adsorbed contaminant in sediment during extreme hydrologic conditions however remains unclear. It is therefore important to understand not only the fate and transport of dissolved contaminants, but also the movement of sediment and the relevant interactions with dissolved contaminant. To narrow this knowledge gap, a study is being conducted using the Tims Branch watershed as a stream-scale ecosystem test-bed to identify the primary transport processes of major contaminants of concern (such as mercury, nickel and uranium) with an emphasis on interactions with sediment transport. This involves the development of a fully distributed hydrologic watershed model of the Tims Branch watershed to predict streamflow under extreme weather conditions, as well as the development of a comprehensive contaminant transport model that can properly account for coupled contaminant and sediment transport. Review of relevant research reports and peer-reviewed journals revealed that aside from advection-dispersion transport of dissolved contaminants, adsorption and desorption with suspended solids and bed sediment also play an important role in the transport of those contaminants of concern. To develop the fully distributed hydrologic watershed model, the MIKE SHE 2-dimensional (2D) land surface/3D groundwater model that simulates surface/subsurface hydrologic processes (such as overland flow, evapotranspiration, and infiltration) was coupled with a 1D streamflow model that accounts for stream water hydraulics (such as hydraulic structures, cross-sections, and network). To model contaminant transport, the MIKE 11 streamflow component was coupled with the MIKE 11 AD module that simulates solute transport through advection and dispersion, and MIKE ECO Lab module that accounts for both sediment transport and interactions with dissolved contaminant. At this stage, the development and optimization of the fully distributed hydrologic model has been completed, achieving satisfactory statistical results between observed and predicted discharge as indicated by a root mean square error (RMSE) of 0.039 cms and a Nash-Sutcliffe efficiency coefficient (NSE) of 0.764. The ongoing development of the contaminant transport model has also yielded realistic results from preliminary tests. Results from this study are a key to evaluating the effectiveness of tin (II)-based mercury treatment of wetlands at the SRS site, and are also relevant to evaluating the potential of using this type of novel remediation technology in other mercury-contaminated stream systems. Knowledge acquired from this research will also support interpretation of historical data on the trends of contaminant concentration distribution in Tims Branch, particularly considering the effect of extreme hydrological events on the stream flow and pollutant transport. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2014

The R-factor gap in macromolecular crystallography: an untapped potential for insights on accurate structures

In macromolecular crystallography, the agreement between observed and predicted structure factors (Rcryst and Rfree) is seldom better than 20%. This is much larger than the estimate of experimental error (Rmerge). The difference between Rcryst and Rmerge is the R-factor gap. There is no such gap in small-molecule crystallography, for which calculated structure factors are generally considered more accurate than the experimental measurements. Perhaps the true noise level of macromolecular data is higher than expected? Or is the gap caused by inaccurate phases that trap refined models in local minima? By generating simulated diffraction patterns using the program MLFSOM, and including every conceivable source of experimental error, we show that neither is the case. Processing our simulated data yielded values that were indistinguishable from those of real data for all crystallographic statistics except the final Rcryst and Rfree. These values decreased to 3.8% and 5.5% for simulated data, suggesting that the reason for high R-factors in macromolecular crystallography is neither experimental error nor phase bias, but rather an underlying inadequacy in the models used to explain our observations. The present inability to accurately represent the entire macromolecule with both its flexibility and its proteinsolvent interface may be improved by synergies between small-angle X-ray scattering, computational chemistry and crystallography. The exciting implication of our finding is that macromolecular data contain substantial hidden and untapped potential to resolve ambiguities in the true nature of the nanoscale, a task that the second century of crystallography promises to fulfill.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1985

An error-dependent model of instrument-scanning behavior in commercial airline pilots

A new flexible model of pilot instrument scanning behavior is presented which assumes that the pilot uses a set of deterministic scanning patterns on the pilot's perception of error in the state of the aircraft, and the pilot's knowledge of the interactive nature of the aircraft's systems. Statistical analyses revealed that a three stage Markov process composed of the pilot's three predicted lookpoints (LP), occurring 1/30, 2/30, and 3/30 of a second prior to each LP, accurately modelled the scanning behavior of 14 commercial airline pilots while flying steep turn maneuvers in a Boeing 737 flight simulator. The modelled scanning data for each pilot were not statistically different from the observed scanning data in comparisons of mean dwell time, entropy, and entropy rate. These findings represent the first direct evidence that pilots are using deterministic scanning patterns during instrument flight. The results are interpreted as direct support for the error dependent model and suggestions are made for further research that could allow for identification of the specific scanning patterns suggested by the model.

Jones, D. H.↗
NASA NTRS2014

On the Relationship Between Global Land-Ocean Temperature and Various Descriptors of Solar-Geomagnetic Activity and Climate

Examined are sunspot cycle- (SC-) length averages of the annual January-December values of the Global Land-Ocean Temperature Index ( ) in relation to SC-length averages of annual values of various descriptors of solar-geomagnetic activity and climate, incorporating lags of 0-5 yr. For the overall interval SC12-SC23, the is inferred to correlate best against the parameter incorporating lag = 5 yr, where the parameter refers to the resultant aa value having removed that portion of the annual aa average value due to the yearly variation of sunspot number (SSN). The inferred correlation between the and is statistically important at confidence level cl > 99.9%, having a coefficient of linear correlation r = 0.865 and standard error of estimate se = 0.149 degC. Excluding the most recent cycles SC22 and SC23, the inferred correlation is stronger, having r = 0.969 and se = 0.048 degC. With respect to the overall trend in the , which has been upwards towards warmer temperatures since SC12 (1878-1888), solar-geomagnetic activity parameters are now trending downwards (since SC19). For SC20-SC23, in contrast, comparison of the against SC-length averages of the annual value of the Mauna Loa carbon dioxide ( ) index is found to be highly statistically important (cl >> 99.9%), having r = 0.9994 and se = 0.012 degC for lag = 2 yr. On the basis of the inferred preferential linear correlation between the and , the current ongoing SC24 is inferred to have warmer than was seen in SC23 (i.e., >0.526 degC), probably in excess of 0.68 degC (relative to the 1951-1980 base period).

Wilson, Robert M.↗
NASA NTRS2010

Maximum Likelihood Time-of-Arrival Estimation of Optical Pulses via Photon-Counting Photodetectors

Many optical imaging, ranging, and communications systems rely on the estimation of the arrival time of an optical pulse. Recently, such systems have been increasingly employing photon-counting photodetector technology, which changes the statistics of the observed photocurrent. This requires time-of-arrival estimators to be developed and their performances characterized. The statistics of the output of an ideal photodetector, which are well modeled as a Poisson point process, were considered. An analytical model was developed for the mean-square error of the maximum likelihood (ML) estimator, demonstrating two phenomena that cause deviations from the minimum achievable error at low signal power. An approximation was derived to the threshold at which the ML estimator essentially fails to provide better than a random guess of the pulse arrival time. Comparing the analytic model performance predictions to those obtained via simulations, it was verified that the model accurately predicts the ML performance over all regimes considered. There is little prior art that attempts to understand the fundamental limitations to time-of-arrival estimation from Poisson statistics. This work establishes both a simple mathematical description of the error behavior, and the associated physical processes that yield this behavior. Previous work on mean-square error characterization for ML estimators has predominantly focused on additive Gaussian noise. This work demonstrates that the discrete nature of the Poisson noise process leads to a distinctly different error behavior.

Erkmen, Baris I.↗
NASA NTRS2021

Probabilistic Processing of Asynchronous Geiger-Mode Avalanche Photodiode Arrays for Background Rate Measurements

A system-level performance evaluation of Geiger-mode avalanche photodiode (GmAPD) arrays requires accurate measurement and prediction of the background rate of the device due to dark counts and other spurious detection events. Since a GmAPD detector reports only a binary value and timestamp associated with an avalanche event, dark count rates are typically measured by averaging thousands of frames to support a statistically significant measurement. For both synchronous and asynchronous detector, the Poisson distributed background rates are referenced to the time each pixel is armed. Unlike for synchronous GmAPD imagers where all the pixels are armed to an array-wide arm signal, an asynchronous pixel operates independently from its neighboring pixels; requiring the background rates to be calculated using an interarrival histogram. For both types of imagers, the background rate is typically evaluated by fitting an exponential distribution to a fixed window within a measured histogram of time intervals between detection events However, if the statistics of the background rate are insufficient – whether that is due to low population sizes, saturation, or a large dynamic range of population size across the array, the pixel, or array-wide, performance metrics may report results with varying accuracy. This paper reports on an implementation of an algorithm that evaluates GmAPD background rates based on statistical metrics rather than fixed windows. The algorithm functions by determining the appropriate integration window within the interarrival time histogram based on a per-pixel count rate set by a predetermined tolerable measurement error. The implementation of the algorithm allows us to characterize GmAPD arrays with orders of magnitude spread in background rates across the detector using common statistical parameters.

Mumolo, J.↗
DOE OSTI2021

Modeling protein structures from predicted contacts with modern molecular dynamics potentials: accuracy, sensitivity, and refinement

Protein structure prediction has become increasingly popular and successful in recent years. An essential step for fragment-free, template-free methods is the generation of a final three-dimensional protein model from a set of predicted amino acid contacts that are often described by interresidue pairwise atomic distances. Here we explore the use of modern, open-source molecular dynamics (MD) engines, which have been continually developed over the last three decades with all-atom Hamiltonians to model biomolecular structure and dynamics, to generate accurate protein structures starting from a set of inferred pairwise distances. Additionally, the ability of MD empirical physical potentials to correct inaccuracies in the predicted geometries is tested. We rigorously characterize the effect of modeling parameters on results, the effect of different amounts of error in the predicted distances on the final structures, and test the ability of post-processing analysis to sort the best models out of a set of statistical replicas. We find that with exact distances and with noisy distances, the method can produce excellent structural models, and that the molecular dynamics force field seems to help correct errors in distance predictions, resisting the effects of applied noise.

Davidson, Russ↗
NASA NTRS2002

Evaluation of the DAO Retrospective Data Assimilation System

We have developed and implemented a retrospective data assimilation system (RDAS) as an upgrade to the operational DAO/Terra data assimilation system. This formulation aims at improving analysis over filter analysis by the dynamically consistent incorporation of observation information past a given analysis time. The current implementation of the RDAS uses the adjoint of the tangent linear model of a simplified version of the Terra general circulation model and extensions to the physical-space statistical analysis system to propagate observation information back in time. The RDAS adopts the same assumptions of the regular data assimilation system, particularly, no explicit propagation of error covariances are involved therefore rendering a procedure that is computationally affordable. In this study, we show results of experiments conducted to investigate the performance of the 6-hour (lag-1) RDAS. Statistical results obtained over one month during a winter season indicate that the RDAS represents considerable improvement over the regular assimilation. Plans for implementation of the RDAS capability in our new finite-volume data assimilation system will also be presented at the time of the conference.

Zhu, Yanqiu↗
NASA NTRS1980

Backward deletion to minimize prediction errors in models from factorial experiments with zero to six center points

Population model coefficients were chosen to simulate a saturated 2 to the 4th fixed-effects experiment having an unfavorable distribution of relative values. Using random number studies, deletion strategies were compared that were based on the F-distribution, on an order statistics distribution of Cochran's, and on a combination of the two. The strategies were compared under the criterion of minimizing the maximum prediction error, wherever it occurred, among the two-level factorial points. The strategies were evaluated for each of the conditions of 0, 1, 2, 3, 4, 5, or 6 center points. Three classes of strategies were identified as being appropriate, depending on the extent of the experimenter's prior knowledge. In almost every case the best strategy was found to be unique according to the number of center points. Among the three classes of strategies, a security regret class of strategy was demonstrated as being widely useful in that over a range of coefficients of variation from 4 to 65%, the maximum predictive error was never increased by more than 12% over what it would have been if the best strategy had been used for the particular coefficient of variation. The relative efficiency of the experiment, when using the security regret strategy, was examined as a function of the number of center points, and was found to be best when the design used one center point.

Holms, A. G.↗
NASA NTRS2001

DAO's Next Generation Physical-Space/Finite-Volume Data Assimilation System: Formulation and Initial Evaluation

The Physical-space/Finite-volume Data Assimilation System (fvDAS) is the next generation global atmospheric data assimilation system in development at the Data Assimilation Office at NASA's Goddard Space Flight Center. It is based on a new finite-volume general circulation model jointly developed by NASA and NCAR, and on the Physical-Space Statistical Analysis System (PSAS) developed at the DAO. In this talk we will describe the general system formulation, the adaptive quality control and general aspects of the error covariance modeling. The NASA-NCAR GCM is a completely new model which replaces the CEOs GCM used in the previous GEOS-1/2/3 Data Assimilation systems. A particular configuration of adaptive Statistical Quality Control and the Physical-space Statistical Analysis System (PSAS) are currently implemented in DAO's operational Data Assimilation System. However, the unique finite-volume formulation of the NASA-NCAR GCM, combined with the generality of the observation-space formulation of PSAS, provides for a very simple and accurate model-analysis interface. The system assimilates a variety of conventional and satellite observations. In particular, TOVS Level 1B radiances are assimilated using a 1-D variational scheme, both in clear sky and cloudy conditions. Computationally, the fvDAS runs approximately 10 times faster than the operational GEOS-Terra system. We will show that the next-generation fvDAS has much improved observation-minus-6hr forecast (O-F) statistics, as well as 5-day forecast skills. Top of the atmosphere radiation fields are in closer agreement with CERES measurements, with realistic precipitation and moisture fields. We will also show that the finite-volume formulation of the fvDAS produce assimilated fields which are more suitable for driving constituent transport models.

daSilva, A.↗
NASA NTRS2023

The Statistical Emulators of GGCMI Phase 2: Responses of Year-to-Year Variation of Crop Yield to CO 2 , Temperature, Water and Nitrogen Perturbations

Understanding the impact of climate change on year-to-year variation of crop yield is critical to global food stability and security. While crop model emulators are believed to be lightweight tools to replace the models, few emulators have been developed to capture such interannual variation of crop yield in response to climate variability. In this study, we developed a statistical emulator with a machine learning algorithm to reproduce the response of year-to-year variation of four crop yields to CO 2 (C), temperature (T), water (W), and nitrogen (N) perturbations defined in the Global Gridded Crop Model Intercomparison Project (GGCMI) phase 2. The emulators were able to explain more than 52 % of the variance of simulated yield and performed well in capturing the year-to-year variation of global average and gridded crop yield over current croplands in the baseline. With the changes in CO 2 –temperature–water–nitrogen (CTWN) perturbations, the emulators could reproduce the year-to-year variation of crop yield well over most current cropland. The variation of R and the mean absolute error was small under the single CTWN perturbations and dual-factor perturbations. These emulators thus provide statistical response surfaces of yield, including both its mean and interannual variability, to climate factors. They could facilitate spatiotemporal downscaling of crop model simulation, projecting the changes in crop yield variability in the future and serving as a lightweight tool for multi-model ensemble simulation. The emulators enhanced the flexibility of crop yield estimates and expanded the application of large-ensemble simulations of crop yield under climate change.

Climate change↗
DOE OSTI2016

A statistical approach for isolating fossil fuel emissions in atmospheric inverse problems

We study independent verification and quantification of fossil fuel (FF) emissions that constitutes a considerable scientific challenge. By coupling atmospheric observations of CO 2 with models of atmospheric transport, inverse models offer the possibility of overcoming this challenge. However, disaggregating the biospheric and FF flux components of terrestrial fluxes from CO 2 concentration measurements has proven to be difficult, due to observational and modeling limitations. In this study, we propose a statistical inverse modeling scheme for disaggregating winter time fluxes on the basis of their unique error covariances and covariates, where these covariances and covariates are representative of the underlying processes affecting FF and biospheric fluxes. The application of the method is demonstrated with one synthetic and two real data prototypical inversions by using in situ CO 2 measurements over North America. Also, inversions are performed only for the month of January, as predominance of biospheric CO 2 signal relative to FF CO 2 signal and observational limitations preclude disaggregation of the fluxes in other months. The quality of disaggregation is assessed primarily through examination of a posteriori covariance between disaggregated FF and biospheric fluxes at regional scales. Findings indicate that the proposed method is able to robustly disaggregate fluxes regionally at monthly temporal resolution with a posteriori cross covariance lower than 0.15 µmol m -2 s -1 between FF and biospheric fluxes. Error covariance models and covariates based on temporally varying FF inventory data provide a more robust disaggregation over static proxies (e.g., nightlight intensity and population density). However, the synthetic data case study shows that disaggregation is possible even in absence of detailed temporally varying FF inventory data.

54 ENVIRONMENTAL SCIENCES↗