Search NASA⌕ Search

SEARCH · Search NASA

Results for “global uncertainty 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 289 records · Page 16

Mitigating Satellite-Based Fire Sampling Limitations in Deriving Biomass Burning Emission Rates: Application to WRF-Chem Model over the Northern sub-Saharan African Region

Largely used in several independent estimates of fire emissions, fire products based on MODIS sensors aboard the Terra and Aqua polar-orbiting satellites have a number of inherent limitations, including (a) inability to detect fires below clouds, (b) significant decrease of detection sensitivity at the edge of scan where pixel sizes are much larger than at nadir, and (c) gaps between adjacent swaths in tropical regions. To remedy these limitations, an empirical method is developed here and applied to correct fire emission estimates based on MODIS pixel level fire radiative power measurements and emission coefficients from the Fire Energetics and Emissions Research (FEER) biomass burning emission inventory. The analysis was performed for January 2010 over the northern sub-Saharan African region. Simulations from WRF-Chem model using original and adjusted emissions are compared with the aerosol optical depth (AOD) products from MODIS and AERONET as well as aerosol vertical profile from CALIOP data. The comparison confirmed an 30-50% improvement in the model simulation performance (in terms of correlation, bias, and spatial pattern of AOD with respect to observations) by the adjusted emissions that not only increases the original emission amount by a factor of two but also results in the spatially continuous estimates of instantaneous fire emissions at daily time scales. Such improvement cannot be achieved by simply scaling the original emission across the study domain. Even with this improvement, a factor of two underestimations still exists in the modeled AOD, which is within the current global fire emissions uncertainty envelope. Plain Language Summary Polar-orbiting satellites sensors, such as MODIS, have limitations in detecting fires under clouds or when viewing angles are large or in the gaps among satellites' different ground swaths. Here we developed an empirical method to mitigate the effect of these limitations in fire emission estimate. The method is applied to a fire emission inventory (FEER) based on MODIS. We show that, with our method, the adjusted emission inventory improves WRF-Chem simulation of smoke transport and distribution.

emissions properties↗

Integrated Multi-Satellite Evaluation for the Global Precipitation Measurement: Impact of Precipitation Types on Spaceborne Precipitation Estimation

Integrated multi-sensor assessment is proposed as a novel approach to advance satellite precipitation validation in order to provide users and algorithm developers with an assessment adequately coping with the varying performances of merged satellite precipitation estimates. Gridded precipitation rates retrieved from space sensors with quasi-global coverage feed numerous applications ranging from water budget studies to forecasting natural hazards caused by extreme events. Characterizing the error structure of satellite precipitation products is recognized as a major issue for the usefulness of these estimates. The Global Precipitation Measurement (GPM) mission aims at unifying precipitation measurements from a constellation of low-earth orbiting (LEO) sensors with various capabilities to detect, classify and quantify precipitation. They are used in combination with geostationary observations to provide gridded precipitation accumulations. The GPM Core Observatory satellite serves as a calibration reference for consistent precipitation retrieval algorithms across the constellation. The propagation of QPE uncertainty from LEO active/passive microwave (PMW) precipitation estimates to gridded QPE is addressed in this study, by focusing on the impact of precipitation typology on QPE from the Level-2 GPM Core Observatory Dual-frequency Precipitation Radar (DPR) to the Microwave Imager (GMI) to Level-3 IMERG precipitation over the Conterminous U.S. A high-resolution surface precipitation used as a consistent reference across scales is derived from the ground radar-based Multi-Radar/Multi-Sensor. While the error structure of the DPR, GMI and subsequent IMERG is complex because of the interaction of various error factors, systematic biases related to precipitation typology are consistently quantified across products. These biases display similar features across Level-2 and Level-3, highlighting the need to better resolve precipitation typology from space and the room for improvement in global-scale precipitation estimates. The integrated analysis and framework proposed herein applies more generally to precipitation estimates from sensors and error sources affecting low-earth orbiting satellites and derived gridded products.

Rainfall↗

A photometric study of Enceladus

We have supplemented Voyager imaging data from Enceladus (limited to phase angles of 13 deg-43 deg) with recent Earth-based CCD observations to obtain an improved determination of the Bond albedo, to construct an albedo map of the satellite, and to constrain parameters in Hapke's (1986) photometric equation. A major result is evidence of regional variations in the physical properties of Enceladus' surface. The average global photometric properties are described by single scattering albedo omega(sub 0) average = 0.998 +/- 0.001, macroscopic roughness parameter theta average = 6 +/- 1 deg, and Henyey-Greenstein asymmetry parameter g = -0.399 +/- 0.005. The value of theta average is smaller than the 14 deg found by fitting whole-disk data, which include all terrains on Enceladus. The opposition surge amplitude B(sub 0) = 0.21 +/- 0.07 and regolith compaction parameter h = 0.014 +/- 0.02 are loosely constrained by the scarcity of and uncertainty in near-opposition observations. From the solar phase curve we determine the geometric albedo of Enceladus p(sub v) = 0.99 +/- 0.06 and phase integral q = 0.92 +/- 0.05, which corresponds to a spherical albedo A = p(sub v)q = 0.91 +/- 0.1. Since the spectrum of Enceladus is fairly flat, we can approximate the Bond albedo A(sub B) with the spherical albedo. Our photometric analysis is summarized in terms of an albedo map which generally reproduces the satellite's observed lightcurve and indicates that normal reflectances range from 0.9 on the leading hemisphere to 1.4 on the trailing one. The albedo map also revels an albedo variation of 15% from longitudes 170 deg to 200 deg, corresponding to the boundary between the leading and trailing hemispheres.

Verbiscer, Anne J.↗

Satellite-Based Evaluation of AeroCom Model Bias in Biomass Burning Regions

Global models are widely used to simulate biomass burning aerosol (BBA). Exhaustive evaluations on model representation of aerosol distributions and properties are fundamental to assess health and climate impacts of BBA. Here we conducted a comprehensive comparison of Aerosol Comparisons between Observations and Models (AeroCom) project model simulations with satellite observations. A total of 59 runs by 18 models from three AeroCom Phase-III experiments (i.e., biomass burning emissions, CTRL16, and CTRL19) and 14 satellite products of aerosols were used in the study. Aerosol optical depth (AOD) at 550 nm was investigated during the fire season over three key fire regions reflecting different fire dynamics (i.e., deforestation-dominated Amazon, Southern Hemisphere Africa where savannas are the key source of emissions, and boreal forest burning in boreal North America). The 14 satellite products were first evaluated against AErosol RObotic NETwork (AERONET) observations, with large uncertainties found. But these uncertainties had small impacts on the model evaluation that was dominated by modeling bias. Through a comparison with Polarization and Directionality of the Earth’s Reflectances measurements with the Generalized Retrieval of Aerosol and Surface Properties algorithm (POLDER-GRASP), we found that the modeled AOD values were biased by −93 % to 152 %, with most models showing significant underestimations even for the state-of-the-art aerosol modeling techniques (i.e., CTRL19). By scaling up BBA emissions, the negative biases in modeled AOD were significantly mitigated, although it yielded only negligible improvements in the correlation between models and observations, and the spatial and temporal variations in AOD biases did not change much. For models in CTRL16 and CTRL19, the large diversity in modeled AOD was in almost equal measures caused by diversity in emissions, lifetime, and the mass extinction coefficient (MEC). We found that in the AeroCom ensemble, BBA lifetime correlated significantly with particle deposition (as expected) and in turn correlated strongly with precipitation. Additional analysis based on Cloud-Aerosol LIdar with Orthogonal Polarization (CALIOP) aerosol profiles suggested that the altitude of the aerosol layer in the current models was generally too low, which also contributed to the bias in modeled lifetime. Modeled MECs exhibited significant correlations with the Ångström exponent (AE, an indicator of particle size). Comparisons with the POLDER-GRASP-observed AE suggested that the models tended to overestimate the AE (underestimated particle size), indicating a possible underestimation of MECs in models. The hygroscopic growth in most models generally agreed with observations and might not explain the overall underestimation of modeled AOD. Our results imply that current global models contain biases in important aerosol processes for BBA (e.g., emissions, removal, and optical properties) that remain to be addressed in future research.

biomass burning aerosols↗

MOBY, A Radiometric Buoy for Performance Monitoring and Vicarious Calibration of Satellite Ocean Color Sensors: Measurement and Data Analysis Protocols

The Marine Optical Buoy (MOBY) is the centerpiece of the primary ocean measurement site for calibration of satellite ocean color sensors based on independent in situ measurements. Since late 1996, the time series of normalized water-leaving radiances L(sub WN)(lambda) determined from the array of radiometric sensors attached to MOBY are the primary basis for the on-orbit calibrations of the USA Sea-viewing Wide Field-of-view Sensor (SeaWiFS), the Japanese Ocean Color and Temperature Sensor (OCTS), the French Polarization Detection Environmental Radiometer (POLDER), the German Modular Optoelectronic Scanner on the Indian Research Satellite (IRS1-MOS), and the USA Moderate Resolution Imaging Spectrometer (MODIS). The MOBY vicarious calibration L(sub WN)(lambda) reference is an essential element in the international effort to develop a global, multi-year time series of consistently calibrated ocean color products using data from a wide variety of independent satellite sensors. A longstanding goal of the SeaWiFS and MODIS (Ocean) Science Teams is to determine satellite-derived L(sub WN)(labda) with a relative combined standard uncertainty of 5 %. Other satellite ocean color projects and the Sensor Intercomparison for Marine Biology and Interdisciplinary Oceanic Studies (SIMBIOS) project have also adopted this goal, at least implicitly. Because water-leaving radiance contributes at most 10 % of the total radiance measured by a satellite sensor above the atmosphere, a 5 % uncertainty in L(sub WN)(lambda) implies a 0.5 % uncertainty in the above-atmosphere radiance measurements. This level of uncertainty can only be approached using vicarious-calibration approaches as described below. In practice, this means that the satellite radiance responsivity is adjusted to achieve the best agreement, in a least-squares sense, for the L(sub WN)(lambda) results determined using the satellite and the independent optical sensors (e.g. MOBY). The end result of this approach is to implicitly absorb unquantified, but systematic, errors in the atmospheric correction, incident solar flux, and satellite sensor calibration into a single correction factor to produce consistency with the in situ data.

Clark, Dennis K.↗

Total ozone-solar activity relationship

Approximately six years of global ozone monthly mean data from the Nimbus 4 BUV instrument are compared with monthly values of solar activity using 10.7-cm flux, F(10.7), as a parameter. Several techniques are explored in calculating the correlation between the two data sets, and all are found to yield relatively high correlations, ranging from R = 0.68 for 'raw' monthly means to 0.94 using a six-month running mean for each data set. It is shown, however, that the bulk of the correlation derives from the long-term decreasing trends in both data sets. When the long-term trends are removed, a cross-correlation analysis produces a maximum with no phase shift or with the ozone variations leading the solar variations by one month, thereby reducing the likelihood of a cause and effect relationship on time scales of this order. In view of the current uncertainty in the long-term stability of the BUV instrument and the resulting uncertainty in any long-term trend derived from its data, it is considered unrealistic to draw firm conclusions about a solar cycle influence on total ozone from this satellite data set.

Reber, C. A.↗

Section on Observed Impacts on El Nino

Agricultural applications of El Nino forecasts are already underway in some countries and need to be evaluated or re-evaluated. For example, in Peru, El Nino forecasts have been incorporated into national planning for the agricultural sector, and areas planted with rice and cotton (cotton being the more drought-tolerant crop) are adjusted accordingly. How well are this and other such programs working? Such evaluations will contribute to the governmental and intergovernmental institutions, including the Inter-American Institute for Global Change Research and the US National Ocean and Atmospheric Agency that are fostering programs to aid the effective use of forecasts. As El Nino climate forecasting grows out of the research mode into operational mode, the research focus shifts to include the design of appropriate modes of utilization. Awareness of and sensitivity to the costs of prediction errors also grow. For example, one major forecasting model failed to predict the very large El Nino event of 1997, when Pacific sea-surface temperatures were the highest on record. Although simple correlations between El Nino events and crop yields may be suggestive, more sophisticated work is needed to understand the subtleties of the interplay among the global climate system, regional climate patterns, and local agricultural systems. Honesty about the limitations of an forecast is essential, especially when human livelihoods are at stake. An end-to-end analysis links tools and expertise from the full sequence of ENSO cause-and-effect processes. Representatives from many disciplines are needed to achieve insights, e.g, oceanographers and atmospheric scientists who predict El Nino events, climatologists who drive global climate models with sea-surface temperature predictions, agronomists who translate regional climate connections in to crop yield forecasts, and economists who analyze market adjustments to the vagaries of climate and determine the value of climate forecasts. Methods include historical studies to understand past patterns and to test hindcasts of the prediction tools, crop modeling, spatial analysis and remote sensing. This research involves expanding, deepening, and applying the understanding of physical climate to the fields of agronomy and social science; and the reciprocal understanding of crop growth and farm economics to climatology. Delivery of a regional climate forecast with no information about how the climate forecast was derived limits its effectiveness. Explanation of a region's major climate driving forces helps to place a seasonal forecast in context. Then, a useful approach is to show historical responses to previous El Nino events, and projections, with uncertainty intervals, of crop response from dynamic process crop growth models. Regional ID forecasts should be updated with real-time weather conditions. Since every El Nino event is different, it is important to track, report and advise on each new event as it unfolds. The stability of human enterprises depends on understanding both the potentialities and the limits of predictability. Farmers rely on past experience to anticipate and respond to fluctuations in the biophysical systems on which their livelihoods depend. Now scientists are improving their ability to predict some major elements of climate variability. The improvements in the reliability of El Nino forecasts are encouraging, but seasonal forecasts for agriculture are not, and will probably never be completely infallible, due to the chaotic nature of the climate system. Uncertainties proliferate as we extend beyond Pacific sea-surface temperatures to climate teleconnections and agricultural outcomes. The goal of this research is to shed as a clear light as possible on these inherent uncertainties and thus to contribute to the development of appropriate responses to El Nino and other seasonal forecasts for a range of stakeholders, which, ultimately, includes food consumers everywhere.

Rosenzweig, Cynthia↗

IMERG and GPCP Seasonality and Response to Climate and Weather Variability

The Integrated Multi-satellitE Retrievals for GPM (IMERG) and the The Global Precipitation Climatology Project (GPCP) are two of the most popular precipitation products. IMERG is a relatively new dataset that targets the needs primarily of the hydrological community by resolving the hydrological cycle of precipitation at fine temporal (30-minutes) and spatial (10-km) scales. IMERG only recently exceeded the user base of the highly successful, but discontinued in 2019, TRMM Multi-satellite Precipitation Analysis (TMPA). GPCP, on the other hand, has traditionally been strong in the climate research community, and recently has been revised under the framework of NASA's Making Earth System Data records for Use in Research Environments (MEaSUREs) program. Both IMERG and GPCP are similar in the underlying approaches to achieve global coverage, in particular using satellite microwave and infrared observations, and adjusting the precipitation retrieval with rain gauge information. While there is a tendency to use both datasets interchangeably, differences remain and some of them limit the usage of IMERG as a climate data record at this point. By applying Principal Component Analysis, we identify the differences and similarities mode-by-mode, and further the guidance on suggested usage of IMERG as a climate data record. As an example in the attached figure, IMERG vs GPSP differences in the explained variance by the two leading seasonal modes can be identified, mainly in the extreme southern latitudes, and over western boundary currents (Gulfstream and Kuroshio). At the preparation time for this presentation, the new version "07"of IMERG was in the works that may resolve the issues presented here. Nevertheless, our analysis can help to gauge the uncertainties of the studies already done using the currently existing IMERG version "06", and evaluate the improvements in the upcoming version "07".

Andrey Savtchenko↗

Using Global Market Demand Analysis to Guide Conceptual Design of Low-Boom Supersonic Transports

This paper uses a mathematical framework to identify the interdependencies of key variables in supersonic transportation demand analysis. The existing quantitative models for supersonic transportation demand analysis are compared for consistency in modeling the interdependencies. Unlike other quantitative models, the Low-Boom Systems Analysis Model (LBSAM2) can propagate important design characteristics of a supersonic transport aircraft concept to the level of economic metrics (such as the number of future supersonic passengers), with a consistent coupling of the market demand analysis, detailed mission analysis, and low-boom constraint. This allows the use of the detailed demand analysis results from LBSAM2 to maximize the economic viability of a supersonic transport aircraft by finding favorable system-level trades between weight, range, fuel burn, and assumed sonic boom ground noise limit for supersonic overland flight. In this paper, LBSAM2 is integrated with conceptual low-boom design to improve the economic viability of low-boom supersonic transport aircraft. A brief discussion of uncertainties in the LBSAM2 analysis is also included, focusing on their impacts on the relative economic advantages between low-boom concepts.

Consistent coupling of disciplinary analyses↗

Using Global Market Demand Analysis to Guide Conceptual Design of Low-Boom Supersonic Transports

This paper uses a mathematical framework to identify the interdependencies of key variables in supersonic transportation demand analysis. The existing quantitative models for supersonic transportation demand analysis are compared for consistency in modeling the interdependencies. Unlike other quantitative models, the Low-Boom Systems Analysis Model (LBSAM2) can propagate important design characteristics of a supersonic transport aircraft concept to the level of economic metrics (such as the number of future supersonic passengers), with a consistent coupling of the market demand analysis, detailed mission analysis, and low-boom constraint. This allows the use of the detailed demand analysis results from LBSAM2 to maximize the economic viability of a supersonic transport aircraft by finding favorable system-level trades between weight, range, fuel burn, and assumed sonic boom ground noise limit for supersonic overland flight. In this paper, LBSAM2 is integrated with conceptual low-boom design to improve the economic viability of low-boom supersonic transport aircraft. A brief discussion of uncertainties in the LBSAM2 analysis is also included, focusing on their impacts on the relative economic advantages between low-boom concepts.

Equilibrium state of supersonic market↗

Comparison of a 2D Photochemical Model to Data Using Statistical Trend Analysis

We have analyzed our 23-year merged ozone data set for variability and trends with a statistical time-series model. To assist in that analysis, we have analyzed the Goddard ED photochemical model for the same time period with the same time-series model. Multiple runs of the photochemical model allow us to separate the effects of various terms on ozone, such as solar cycle and volcanic eruptions. We use this to help us separate those signals from each other in the data. We also used a smoothed version of the photochemical model's prediction of global ozone change in place of a simple linear trend. We find a number of interesting results. This particular photochemical model is more sensitive to chlorine perturbations than the atmosphere appears to be. It is less sensitive to solar cycle. It predicts an effect from the Pinatubo eruption that is nearly symmetric in the two hemispheres, but the data appears to have not responded to Pinatubo in the southern mid-latitudes. These results and their uncertainties will be discussed.

Stolarski, Richard S.↗

Comparison of a 2D Photochemical Model to Data Using Statistical Trend Analysis

We have analyzed our 23-year merged ozone data set for variability and trends with a statistical time-series model. To assist in that analysis, we have analyzed the Goddard 2D photochemical model for the same time period with the same time-series model. Multiple runs of the photochemical model allow us to separate the effects of various terms on ozone, such as solar cycle and volcanic eruptions. We use this to help us separate those signals from each other in the data. We also used a smoothed version of the photochemical model's prediction of global ozone change in place of a simple linear trend. We find a number of interesting results. This particular photochemical model is more sensitive to chlorine perturbations than the atmosphere appears to be. It is less sensitive to solar cycle. It predicts an effect from the Pinatubo eruption that is nearly symmetric in the two hemispheres, but the data appears to have not responded to Pinatubo in the southern mid-latitudes. These results and their uncertainties will be discussed.

Bhartia, P. K.↗

Variation in Coronal Activity from Solar Cycle 24 Minimum to Maximum Using Three-Dimensional Reconstructions of the Coronal Electron Density from STEREO/COR1

Three-dimensional electron density distributions in the solar corona are reconstructed for 100 Carrington rotations (CR 2054 - 2153) during 2007/03 - 2014/08 using the spherically symmetric method from polarized white-light observations with the inner coronagraph (COR1) onboard the twin Solar Terrestrial Relations Observatory (STEREO). These three-dimensional electron density distributions are validated by comparison with similar density models derived using other methods such as tomography and a magnetohydrodynamics (MHD) model as well as using data from the Solar and Heliospheric Observatory (SOHO)/Large Angle and Spectrometric Coronagraph (LASCO)-C2. Uncertainties in the estimated total mass of the global corona are analyzed based on differences between the density distributions for COR1-A and -B. Long-term variations of coronal activity in terms of the global and hemispheric average electron densities (equivalent to the total coronal mass) reveal a hemispheric asymmetry during the rising phase of Solar Cycle 24, with the northern hemisphere leading the southern hemisphere by a phase shift of 7 - 9 months. Using 14 CR (approx. equal to 13-month) running averages, the amplitudes of the variation in average electron density between Cycle 24 maximum and Cycle 23/24 minimum (called the modulation factors) are found to be in the range of 1.6 - 4.3. These modulation factors are latitudinally dependent, being largest in polar regions and smallest in the equatorial region. These modulation factors also show a hemispheric asymmetry: they are somewhat larger in the southern hemisphere. The wavelet analysis shows that the short-term quasi-periodic oscillations during the rising and maximum phases of Cycle 24 have a dominant period of 7 - 8 months. In addition, it is found that the radial distribution of the mean electron density for streamers at Cycle 24 maximum is only slightly larger (by approx. equal to 30%) than at cycle minimum.

Wang, Tongjiang↗

Global Comparison of CFD and Wind-Tunnel-Derived Force and Moment Databases for the Space Launch System

Recently a very large (739 runs) collection of high-fidelity RANS CFD solutions was obtained for Space Launch System ascent aerodynamics for the vehicle to be used for the first exploratory (unmanned) mission (EM-1). The extensive computations, at full-scale conditions, were originally developed to obtain detailed line and protuberance loads and surface pressures for venting analyses. The line loads were eventually integrated for comparison of the resulting forces and moments to the database that was derived from wind tunnel tests conducted at sub-scale conditions. The comparisons presented herein cover the ranges 0.5 < or = M(infinity) < or = 5, −6deg < or = alpha < or = 6deg, and −6deg < or = beta < or = 6deg. For detailed comparisons, slender-body-theory-based component build-up aero models from missile aerodynamics are used. The differences in the model fit coefficients are shown to be relatively small except for the low supersonic Mach number range, 1.1 < or = M(infinity) < or = 2.0. The analysis is intended to support process improvement and development of uncertainty models.

Hemsch, Michael J.↗

Assessment of Clear Sky Radiative Forcing in the Caribbean Region Using an Aerosol Dispersion Model and Ground Radiometry During Puerto Rico Dust Experiment

This study investigates the surface and top of the atmosphere solar radiative forcing by long-range transport of Saharan dust. The calculations of radiative forcing are based on measurements collected in the Puerto Rico Dust Experiment (PRIDE) carried out during July, 2000. The purpose of the experiment was the characterization of the Saharan dust plume, which frequently reaches the Caribbean region during the summer. The experiment involved the use of three approaches to study the plume: space and ground based remote sensing, airborne and ground based in-situ measurements and aerosol dispersion modeling. The diversity of measuring platforms provides an excellent opportunity for determination of the direct effect of dust on the clear sky radiative forcing. Specifically, comparisons of heating rates, surface and TOA fluxes derived from the Navy global aerosol dispersion model NAAPS (NRL Aerosol Analysis and Prediction System) and actual measurements of fluxes from ground and space based platforms are shown. In addition, the direct effect of dust on the clear sky radiative forcing is modeled. The extent and time of evolution of the radiative properties of the plume are computed with the aerosol concentrations modeled by NAAPS. Standard aerosol parameterizations, as well as in-situ composition and size distributions measured during PRIDE, are utilized to compute the aerosol optical depth, single scattering albedo and asymmetry factor. Radiative transfer computations are done with an in-house modified spectral radiative transfer code (Fu-Liou). The code includes gas absorption and cloud particles (ice and liquid phase) and it allows the input of meteorological data. The code was modified to include modules for the aerosols contribution to the calculated fluxes. This comparison study helps to narrow the current uncertainty in the dust direct radiative forcing, as recently reported in the 2001 IPCC assessment.

Gasso, Santiago↗

Regional and Model-Specific Response Types in A Global Gridded Crop Model Ensemble

Crop models are often employed to project crop yields under changing conditions such as global warming and associated management change for adaptation. Multi-model ensembles are promoted to enhance the robustness of projections, but questions remain on what causes often large differences between projections of individual models. Global Gridded Crop Models (GGCMs) are especially exposed to this question when applied for assessing climate change impacts, adaptation, environmental impacts of agricultural production, because their results are used in downstream analyses, such as in integrated assessment or economic modeling for projecting future land-use change. Even though global gridded crop models are often based on detailed field-scale models or have implemented similar modeling principles in other ecosystem models, global-scale models are subject to substantial uncertainties from both model structure and parametrization as well as from calibration and input data quality. AgMIP’s Global Gridded Crop Model Intercomparison (GGCMI) has thus set out to intercompare GGCMs in order to evaluate model performance, describe model uncertainties, identify inconsistencies within the ensemble and underlying reasons, and to ultimately improve models and modeling capacities. In phase 2 of the GGCMI activities, 12 modeling groups followed a modeling protocol that asked for up to 1404 31-year global simulations at 0.5 arc-degree spatial resolution to assess models’ sensitivities to changes in carbon dioxide (C; 4 different levels) temperature (T; 7 different offset levels), water supply (W; 9 levels), and nitrogen (N; 3 levels), the so-called CTWN experiment (Franke et al. 2020; http://dx.doi.org/10.5194/gmd-13-2315-2020). We here present analyses of model response types using impact response surfaces along the C, T, W, and N dimensions, respectively and collectively. Doing so, we can understand differences in simulated responses per driver rather than aggregated changes in yields. We find that models’ sensitivities to the individual driver dimensions are substantially different and often more different across models than across regions. A cluster analysis finds regional and model-specific patterns. There is some agreement across models with respect to the spatial patterns of response types but strong differences in the distribution of response type clusters across models suggests that models need to undergo further scrutiny. We suggest establishing standards in model process evaluation not only against historical dynamics but also against dedicated experiments across the CTWN dimensions.

crop models↗

Clouds in the Arctic Climate System: A Force for Change?

Clouds play a central role in the global climate system through the modulation of Earth’s energy flows and as a mediator of precipitation. In the Arctic, clouds are also a major player; however, the processes that govern cloud evolution in the Arctic differ from most of the other regions of the globe. Thus, clouds are a key “wildcard” within the Arctic climate system that could have a substantial influence on the Arctic climate system response to anthropogenic forcing. What is the role of clouds within the phenomenon known as Arctic Amplification? This is question does not have a clear answer. Clouds seem to be a center to the important processes driving Arctic Amplification (e.g., the atmospheric response to sea ice loss and airmass transformation), however feedback analysis studies indicate that net cloud feedback in the Arctic is small. This seminar discusses the role of clouds within the Arctic Amplification processes reviewing aspects of what we know about Arctic clouds and the uncertainties that limit our ability to model them. Results from recently published and ongoing work are presented that provide an observationally-based estimate of the cloud-sea ice feedback and evaluate cloud properties within models. The goal of this presentation is to ignite a discussion and new collaborations around the best approaches to resolving uncertainties related to the role of clouds with the Arctic system and how to better represent Arctic clouds in models.

Patrick C Taylor↗

Sensitivity Studies for Assimilated Ozone Profiles

An ozone data assimilation system at the NASA/Goddard Data Assimilation Office (DAO) produces three-dimensional global ozone fields. They are obtained by assimilating ozone retrieved from the Solar Backscatter UltraViolet/2 (SBUV/2) instrument and the Earth Probe Total Ozone Mapping Spectrometer (EP TOMS) measurements into an off-line parameterized chemistry and transport model. In this talk we focus on the quality of lower stratospheric assimilated ozone profiles. Ozone in the lower stratosphere plays a key role in the forcing of climate. A biased ozone field in this region will adversely impact calculations of the stratosphere-troposphere exchange and, when used as a first guess in retrievals, the values determined from satellite observations. The SBUV/2 ozone data have a coarse vertical resolution with increased uncertainty below the ozone maximum, and TOMS provides only total ozone columns. Thus, the assimilated ozone profiles in the lower stratosphere are only weakly constrained by the incoming SBUV and TOMS data. Consequently, the assimilated ozone distribution should be sensitive to changes in inputs to the statistical analysis scheme. We investigate the sensitivity of assimilated ozone profiles to changes in a variety of system inputs: TOMS and SBUV/2 data selection, forecast and observations error covariance models, inclusion or omission of a parameterized chemistry model, and different versions of DAO assimilated wind fields used to drive the transport model. Comparisons of assimilated ozone fields with independent observations, primarily ozone sondes, are used to determine the impact of each of these changes.

Stajner, Ivanka↗