Search NASA⌕ Search

SEARCH · Search NASA

Results for “Spatial patterns”

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 73 records · Page 4

Maps of growing season gross primary production and net ecosystem exchange for Council Road Mile Marker 71, Seward Peninsula, Alaska, [2017-2023]

This data archive is in support of the Next-Generation Ecosystem Experiments in the Arctic (NGEE Arctic) publication "Integrating Characteristic Arctic Vegetation in a Land Surface Model Improves Representation of Carbon Dynamics Across a Tundra Landscape", by Murphy et al. (2025a). Murphy et al. (2025a) evaluated whether incorporating observed Arctic vegetation heterogeneity into ELM, the land model of the Department of Energy’s Energy Exascale Earth System Model (E3SM), improved simulations of tundra carbon cycling. The associated model archive can be found at Murphy et al. (2025b). The study focused on the spatial patterns and net landscape-level growing season productivity and carbon uptake. As part of this evaluation, observationally derived maps of average growing season (June–August) net ecosystem exchange (NEE) and gross primary production (GPP) were developed for the same domain. These maps, which form the dataset described here, integrate eddy covariance flux tower, remote sensing, and vegetation community data to provide spatially explicit benchmarks for model evaluation. The maps provide spatially explicit estimates of average growing season NEE and GPP across 13 tundra vegetation communities within the study domain. By combining flux tower observations with Airborne Visible-Infrared Imaging Spectrometer-Next Generation (AVIRIS-NG) hyperspectral imagery and drone-based normalized difference vegetation index (NDVI), these maps capture the heterogeneity of carbon fluxes associated with different Arctic vegetation types. While they represent average seasonal conditions rather than interannual variability, the maps provide a unique dataset for evaluating model performance, comparing vegetation community contributions to landscape-scale carbon cycling, and supporting regional analyses of Arctic carbon dynamics. This data archive contains 5 m resolution maps of vegetation communities, vegetation community average growing season GPP, and vegetation community average growing season NEE (three *.tif files), a User’s Guide (*pdf file), and Table 1 of the User’s Guide displaying vegetation community coverage and average growing season NEE and GPP values (*.csv file).

Murphy, Bailey [ORNL] (ORCID:0000000203995221)↗

Poleward Migration of the Latitude of Maximum Tropical Cyclone Intensity—Forced or Natural?

Abstract Past studies have shown a significant observed poleward trend in the latitude at which tropical cyclones reach their lifetime maximum intensity (LMI), especially in the northwest Pacific basin. Given the brevity of the historical record, it remains difficult to separate the forced trend from internal variability of the climate system. A recently developed tropical cyclone downscaling model is used to downscale the Community Earth System Model, version 2 (CESM2), preindustrial control simulation. It is found that the observed trend in the latitude at which tropical cyclones reach their LMI in the northwest Pacific is very unlikely to be caused by internal variability. The same downscaling model is then used to downscale CESM2 simulations under historical forcing. The resulting trend distribution shows a significant poleward migration of tropical cyclone LMI even after regressing out both natural variability and the part of the forced warming pattern that projects onto natural variability. The results indicate that the observed poleward migration of the latitude at which tropical cyclones reach their LMI in the northwest Pacific basin is likely to be, at least in part, forced. However, the magnitude of the projected poleward trend in climate models can be significantly modulated by the simulated spatial pattern of ocean warming. This highlights how discrepancies between models and observations, with regard to projected changes to the equatorial zonal sea surface temperature gradient under anthropogenic forcing, can lead to large uncertainties in projected changes to the LMI latitude of tropical cyclones. Significance Statement Observations in the northwest Pacific basin show that the latitude at which tropical cyclones are at their most intense has been trending northward in the recent half century. These changes are important since tropical cyclones could bring hazardous weather to coastal areas that are poorly equipped to handle them. Here, we show that natural variations in Earth’s climate are very unlikely to explain the observed poleward trend in the latitude that tropical cyclone reach their maximum intensity. We find that it is much more likely that the observed trend is forced by human-related emissions, though the spatial pattern of warming in response to greenhouse emissions can have significant impacts on the magnitude of the trend.

Lin, Jonathan↗

Global Patterns of Lightning Properties Derived by OTD and LIS

The satellite instruments Optical Transient Detector (OTD) and Lightning Imaging Sensor (LIS) provide unique empirical data about the frequency of lightning flashes around the globe (OTD), and the tropics (LIS), which 5 has been used before to compile a well received global climatology of flash rate densities. Here we present a statistical analysis of various additional lightning properties derived from OTD/LIS, i.e. the number of so-called "events" and "groups" per flash, as well as 10 the mean flash duration, footprint and radiance. These normalized quantities, which can be associated with the flash "strength", show consistent spatial patterns; most strikingly, oceanic flashes show higher values than continental flashes for all properties. Over land, regions with high (Eastern US) 15 and low (India) flash strength can be clearly identified. We discuss possible causes and implications of the observed regional differences. Although a direct quantitative interpretation of the investigated flash properties is difficult, the observed spatial patterns provide valuable information for the 20 interpretation and application of climatological flash rates. Due to the systematic regional variations of physical flash characteristics, viewing conditions, and/or measurement sensitivities, parametrisations of lightning NOx based on total flash rate densities alone are probably affected by regional biases.

Beirle, Steffen↗

High-Resolution Mapping of Nitrogen Oxide Emissions in Large US Cities From TROPOMI Retrievals of Tropospheric Nitrogen Dioxide Columns

Satellite-derived spatiotemporal patterns of nitrogen oxide (NO x ) emissions can improve accuracy of emission inventories to better support air quality and climate research and policy studies. In this study, we develop a new method by coupling the chemical transport Model-Independent SATellite-derived Emission estimation Algorithm for Mixed-sources (MISATEAM) with a divergence method to map high-resolution NO x emissions across US cities using TROPOspheric Monitoring Instrument (TROPOMI) tropospheric nitrogen dioxide (NO 2 ) retrievals. The accuracy of the coupled method is validated through application to synthetic NO 2 observations from the NASA-Unified Weather Research and Forecasting (NU-WRF) model, with a horizontal spatial resolution of 4 km × 4 km for 33 large and mid-size US cities. Validation reveals excellent agreement between inferred and NU-WRF-provided emission magnitudes (R= 0.99, normalized mean bias, NMB = −0.01) and a consistent spatial pattern when comparing emissions for individual grid cells (R=0.88 ± 0.06). We then develop a TROPOMI-based database reporting annual emissions for 39 US cities at a horizontal spatial resolution of 0.05° × 0.05° from 2018 to 2021. This database demonstrates a strong correlation (R= 0.90) with the National Emission Inventory (NEI) but reveals some bias (NMB = −0.24). There are noticeable differences in the spatial patterns of emissions in some cities. Our analysis suggests that uncertainties in TROPOMI-based emissions and potential misallocation of emissions and/or missing sources in bottom-up emission inventories both contribute to these differences.

NOx emissions↗

Spaceborne Lidar Retrievals of PM2.5 for Air Quality Studies and Applications

Fine particulate matter (PM2.5) substantially contributes to air pollution and negatively affects human health. While many studies have investigated the use of passive column-integrated aerosol optical depth to infer surface PM2.5, the use of lidar observations for air quality characterization is not nearly as extensive. Lidar measurements are critical, however, due to the vertical aerosol information they provide, including near the surface. In this presentation, we first provide an overview of various lidar-based approaches for estimating PM2.5 concentrations and then discuss how lidar measurements can assist other air quality applications. For example, estimates of PM2.5 have been obtained in a physics-based approach through CALIOP near-surface aerosol extinction retrievals, assumptions on the mass extinction efficiency, and incorporating other parameters (an aerosol hygroscopic growth factor and PM2.5/PM10 ratio). Application of this algorithm over the contiguous United States (CONUS) from 2006 to 2018 yielded larger PM2.5 values over the eastern and western CONUS (~10-15 μg/m³) and lower PM2.5 levels in the central CONUS (~5 μg/m³). These spatial patterns were similar to those from gridded PM2.5 concentrations obtained through in situ measurements at ground stations operated by the US Environmental Protection Agency. In another approach, the Cloud Aerosol Transport System (CATS) lidar was used with the Goddard Earth Observing System (GEOS) model in a 1D ensemble-based variational technique to obtain PM2.5 over the US and Europe, and the spatial patterns of the CATS/GEOS based PM2.5 concentrations generally captured those from surface stations (with corresponding hourly EPA PM2.5 vs CATS PM2.5 statistics of R=0.4 and bias=1.5 μg/m³). In our recent work, as part of the Models, In situ, and Remote sensing of Aerosols (MIRA) Working Group, we have applied both the CALIOP and CATS/GEOS based approaches over the highly polluted country of India during the post-monsoon season (September-October 2016). We derived elevated levels of two-month mean PM2.5 (~100 μg/m³) in northern India, especially near New Delhi. These high PM2.5 concentrations in the Indo-Gangetic plain are driven in large part from the seasonal burning of crop residue and meteorological conditions typical at this time of the year, such as low wind speeds and a shallow boundary layer. While the satellite-derived PM2.5 moderately replicates (R = ~0.7-0.9) the spatial variability in the two-month mean of surface in situ PM2.5 from monitoring sites operated by the Central and State Pollution Control Boards, we show results from specific scenes for which there are large deviations between the satellite-derived PM2.5 and in situ measurements. Other current work on this topic focuses on developing PM2.5 estimates using airborne high spectral resolution lidar measurements through machine learning regression algorithms and involves several parameters (e.g., aerosol extinction, color ratio, lidar ratio). Application of this method over major metropolitan areas in the US and Asia have resulted in high correlations (R = 0.93) with surface measurements. This airborne lidar approach can be adapted to spaceborne lidar measurements, and all three of these approaches can be applied to ESA’s EarthCARE Atmospheric Lidar instrument, setting the stage for the future Cloud Aerosol Lidar for Global Scale Observations of the Ocean-Land Atmosphere System (CALIGOLA) mission. Ultimately, beyond estimates of PM2.5, the aerosol vertical distribution from lidars can benefit studies involving passive sensor approaches for PM2.5 proxies (including from geostationary satellites), wildfire smoke plume injection heights, volcanic emissions (e.g., ash height retrievals), and aerosol/air quality model assimilation, evaluation, and forecasts.

Travis D Toth↗

A Rotation Procedure to Improve Seasonally Varying Empirical Orthogonal Function Bases for MJO Indices

Various indices have been defined to characterize the phase and amplitude of the Madden-Julian oscillation (MJO). One widely used index is the Outgoing Longwave Radiation (OLR) based MJO index (OMI), which is calculated using the spatial pattern of 30-96-day eastward-filtered OLR. The EOFs used to calculate the OMI in observations are prone to degeneracy and exhibit oscillations on the order of 10-20 days, despite initial filtering of the OLR. We propose a simple modification to the OMI that involves aligning the EOFs between neighboring days while retaining the spatial pattern described by the EOFs. This rotation method is implemented as a postprocessing procedure of the current OMI calculation and cleanly removes the spurious oscillations and degeneracy issues seen in the standard method. A similar rotation procedure can be implemented in calculations of other MJO indices.

Sarah Weidman↗

Impacts of Climate Change on Surface Ozone and Intercontinental Ozone Pollution: A Multi-Model Study

The impact of climate change between 2000 and 2095 SRES A2 climates on surface ozone (O)3 and on O3 source-receptor (S-R) relationships is quantified using three coupled climate-chemistry models (CCMs). The CCMs exhibit considerable variability in the spatial extent and location of surface O3 increases that occur within parts of high NOx emission source regions (up to 6 ppbv in the annual average and up to 14 ppbv in the season of maximum O3). In these source regions, all three CCMs show a positive relationship between surface O3 change and temperature change. Sensitivity simulations show that a combination of three individual chemical processes-(i) enhanced PAN decomposition, (ii) higher water vapor concentrations, and (iii) enhanced isoprene emission-largely reproduces the global spatial pattern of annual-mean surface O3 response due to climate change (R2 = 0.52). Changes in climate are found to exert a stronger control on the annual-mean surface O3 response through changes in climate-sensitive O3 chemistry than through changes in transport as evaluated from idealized CO-like tracer concentrations. All three CCMs exhibit a similar spatial pattern of annual-mean surface O3 change to 20% regional O3 precursor emission reductions under future climate compared to the same emission reductions applied under present-day climate. The surface O3 response to emission reductions is larger over the source region and smaller downwind in the future than under present-day conditions. All three CCMs show areas within Europe where regional emission reductions larger than 20% are required to compensate climate change impacts on annual-mean surface O3.

Troposphere: constituent transport and chemistry↗

A Numerical Study of the Liquid Motion in Titan’s Subsurface Ocean

An ocean filled with liquid water lies beneath the icy surface of several Jovian and Saturnian moons. In such an ocean, the currents are driven by various phenomena such as the tidal forcing, the deformation of the ice shell lying at its top, the temperature gradient resulting from the surface and bottom heat fluxes…The flow induced by the first two forcings can be modelled by means of a 2D depth-averaged model, while the third one generates horizontal and vertical density variations whose effects can only be captured by a 3D baroclinic model. We study the tides of Titan’s subsurface ocean and the impact of the ice shell on the liquid motion by means of the Second-generation Louvain-la-Neuve Ice–ocean Model, SLIM (https://www.climate.be/slim). The impact of the ice shell lying at the top of the ocean is modelled by a surface friction term and surface pressure terms. The latter are a function of the difference between the ocean elevation and the vertical displacement of the shell and the time derivative of this difference. Because of Titan’s appreciable obliquity (0.306°), the tidal motion expected (and found) is similar to the Europa tidal scenario described by Tyler (2008): the surface elevation consists of two bulges rotating around Titan and the associated depth-averaged velocity field consists of two gyres, separated by an area of high speed flow, whose centre follows a sinusoidal path centred on the equator. The ice shell damps the surface motion, thus slowing down the flow, without significantly modifying the spatial patterns of these fields. The depth of the ocean and the mechanical characteristics of the ice shell being poorly constrained, a sensitivity analysis is conducted. The depth-averaged flow slows down when the depth is increased and a lag appears in the tidal phase but the tidal range remains similar. The ice shell mechanical characteristics influences both the elevation and depth-averaged velocity fields in terms of magnitude but does not modify the spatial patterns of these fields. The influence of the surface heat flux is studied by means of the 3D baroclinic version of SLIM. The heat flux derived from Titan’s topography by Kvorka et al. (2018) is used as surface boundary condition for the temperature equation while a uniform bottom heat flux is implemented. Its value is computed assuming that the heat budget of the ocean is at equilibrium. These boundary conditions cause density variations, which impact the hydrodynamics of the ocean. While the flow velocity induced by these variations is two orders of magnitude smaller than the tidal flow, its orientation is time-independent, hence impacting the orientation of the velocity field. Although the variations of ocean surface elevation and speed with respect to the shell mechanical properties can be larger than those induced by the surface heat flux, taking into account the latter results in large variations of the velocity field global patterns, which was not observed when modifying the shell mechanical properties. Future studies should therefore focus on modelling the surface and bottom heat fluxes while uncertainties about the mechanical characteristics of the shell can be tolerated.

Titan↗

High-Resolution Subtropical Summer Precipitation Derived from Dynamical Downscaling of the NCEP-DOE Reanalysis: How Much Small-Scale Information Is Added by a Regional Model?

This study assesses the regional-scale summer precipitation produced by the dynamical downscaling of analyzed large-scale fields. The main goal of this study is to investigate how much the regional model adds smaller scale precipitation information that the large-scale fields do not resolve. The modeling region for this study covers the southeastern United States (Florida, Georgia, Alabama, South Carolina, and North Carolina) where the summer climate is subtropical in nature, with a heavy influence of regional-scale convection. The coarse resolution (2.5deg latitude/longitude) large-scale atmospheric variables from the National Center for Environmental Prediction (NCEP)/DOE reanalysis (R2) are downscaled using the NCEP Environmental Climate Prediction Center regional spectral model (RSM) to produce precipitation at 20 km resolution for 16 summer seasons (19902005). The RSM produces realistic details in the regional summer precipitation at 20 km resolution. Compared to R2, the RSM-produced monthly precipitation shows better agreement with observations. There is a reduced wet bias and a more realistic spatial pattern of the precipitation climatology compared with the interpolated R2 values. The root mean square errors of the monthly R2 precipitation are reduced over 93 (1,697) of all the grid points in the five states (1,821). The temporal correlation also improves over 92 (1,675) of all grid points such that the domain-averaged correlation increases from 0.38 (R2) to 0.55 (RSM). The RSM accurately reproduces the first two observed eigenmodes, compared with the R2 product for which the second mode is not properly reproduced. The spatial patterns for wet versus dry summer years are also successfully simulated in RSM. For shorter time scales, the RSM resolves heavy rainfall events and their frequency better than R2. Correlation and categorical classification (above/near/below average) for the monthly frequency of heavy precipitation days is also significantly improved by the RSM.

NARR↗

Future Spatially Explicit Patterns of Land Transitions in the United States With Multiple Stressors

Climate change, income and population growth, and changing diets are major drivers of the global food system with implications for land use change. Land use in the U.S. will be affected directly by local and regional forces and indirectly through international trade. In order to investigate the effects of several potential forces on land use changes in the U.S., we advanced capabilities in representing the interactions between natural and human systems by linking a multisectoral and multiregional socio-economic model of the world economy to a model that downscales land use to a 0.5°grid scale. This enables us to translate regional projections of future land use into higher-resolution representations of time-evolving land cover (effectively spatially explicit land use transitions). We applied the framework over the U.S., with a particular interest in the Mississippi River Basin and its four sub-basins, to consider how a range of global drivers affect land use and cover in the target regions. Our results show that under scenarios of high pressure on the world food system a comparative advantage in livestock production amplifies the recent trend toward less cropland and more pastures in the U.S. Under low pressures on the world food system agricultural land is used less intensively. However, there can be key differences among the various land-use transitions at the sub- basin scale. Overall, these results highlighted the need for high resolution details to explicitly understand the implications of land use change on environmental impacts such as carbon storage, soil erosion, chemical use, hydrology, and water quality.

54 ENVIRONMENTAL SCIENCES↗

Average and Marginal Capacity Credit Values of Renewable Energy and Battery Storage in the United States Power System

As deployment of renewable resources and storage continue to significantly grow in the coming decades, these technologies will play increasingly important roles in maintaining power systems' resource adequacy. Few analyses so far offer comprehensive comparisons of forward-looking average and marginal capacity credits of variable renewable energy and storage in the U.S. interconnections across a wide range of possible futures. To fill this research gap, we quantify the average and marginal capacity credits of solar PV, onshore and offshore wind, and batteries between 2026 and 2050 across the U.S power systems to examine the temporal trends, spatial patterns, and trade-offs between these two capacity accreditation approaches. Across technologies, capacity credits of solar PV most clearly follow downward trends over time, reflecting the significant rise in solar PV generation share as the grid decarbonizes. While battery storages' generation shares also rise significantly over time, their capacity credits always remain stably high due to their capabilities to be dispatched strategically during critical periods to maintain reliability. On the other hand, capacity credits of wind technologies in general follow slight upward trends as their generation shares level off. There are strong spatial variabilities of both average and marginal capacity credits across technologies, but capacity credits of solar PV displaying the most obvious spatial patterns with high capacity credits concentrating in wind-rich, solar-poor regions in SPP, PJM, and MISO, suggesting potential reliability benefits of interconnection-wide planning for renewable energy deployments. Additionally, except for offshore wind, average capacity credits of all other renewable technologies tend to be higher than their marginal capacity credits, indicating that existing renewable resources tend to be accredited higher than new resources at almost any time.

25 ENERGY STORAGE↗

Generation of Continental Scale Percent Tree Cover Product Using Deep-learning and Multi-scale Remote Sensing Data

Spatially explicit percent tree cover (TC) estimation is critical for mapping forest aboveground biomass and its dynamics. While various TC products have been developed, there has not been a generalized framework that can be applied to diverse terrestrial ecosystems due to underlain extreme complexities. Deep learning algorithms can learn a spatial pattern and radiometric characteristics of tree canopy as a robust approximation of physical or empirical models, and thus have emerged as promising and efficient tools for large-scale TC mapping. In this study, we synergistically use very high-resolution aerial imageries (National Agriculture Imagery Program, NAIP) and medium resolution Landsat data to map continental-scale TC (CONUS and Mexico) through a hierarchical deep learning approach (Convolutional Neural Network), i.e., NAIP TC generated from a NAIP model is utilized to train a Landsat model. The produced TC product (hereafter, NEX-TC) is able to capture the spatial pattern of TC distribution and its changes driven by natural disturbance and human land management. We further explore and analyze the reliability and potential uncertainty of the NEX-TC by comparing it to lidar- (lidar-TC), National Land Cover Database (NLCD-TC), and MODIS Vegetation Continuous Field (MODIS-TC). This evaluation practice reveals that TC products based on passive optical sensors tend to underestimate TC across all land cover types while Landsat-based TCs (i.e., NEX-TC & NLCD-TC) perform better than the coarser MODIS TC estimate. Our results show that the NEX-TC is generally comparable to NLCD-TC but it particularly outperforms NLCD-TC and MODIS-TC over the dense forests where lidar-TC indicates >80% TC. These results indicate that our hierarchical deep learning approach and TC product will be effective and useful for characterizing large-scale tree cover and possibly associated carbon dynamics.

Landsat↗

Contrasting conditions of surface water balance in wet years and dry years as a possible land surface-atmosphere feedback mechanism in the West African Sahel

The climate of West Africa, in particular the Sahel, is characterized by multiyear persistence of anomalously wet or dry conditions. Its Southern Hemisphere counterpart, the Kalahari, lacks the persistence that is evident in the Sahel even though both regions are subject to similar large-scale forcing. It has been suggested that land surface-atmosphere feedback contributes to this persistence and to the severity of drought. In this study, surface energy and water balance are quantified for nine stations along a latitudinal transect that extends from the Sahara to the Guinea coast. In the wetter regions of West Africa, the difference between wet and dry years is primarily reflected in the magnitude of runoff. For the Sahel and drier locations, evapotranspiration and soil moisture are more sensitive to rainfall anomalies. The increase in evapotranspiration, and hence latent heating, over the Sahel in wet years alters the thermal structure and gradients of the overlying atmosphere and thus the strength of the African easterly jet (AEJ) at 700 mb. The difference between dry and wet Augusts corresponds to a decrease in magnitude of the AEJ at 15 deg N on the order of 2.6 m/s, which is consistent with previous studies of observed winds. Spatial patterns were also developed for surface water balance parameters for both West Africa and southern Africa. Over southern Africa, the patterns are not as spatially homogeneous as those over West Africa and are lower in magnitude, thus supporting the suggestion that the persistence of rainfall anomalies in the Sahel might be due, at least in part, to land-atmosphere feedback, and that the absence of such persistence in the Kalahari is a consequence of less significant changes in surface water and energy balance.

Lare, A. R.↗

Dynamics and pattern of a managed coniferous forest landscape in Oregon

We examined the process of fragmentation in a managed forest landscape by comparing rates and patterns of disturbance (primarily clear-cutting) and regrowth between 1972 and 1988 using Landsat imagery. A 2589-km(exp 2) managed forest landscape in western Oregon was classified into two forest types, closed-canopy conifer forest (CF) (typically, greater than 60% conifer cover) and other forest and nonforest types (OT) (typically, less than 40 yr old or deciduous forest). The percentage of CF declined from 71 to 58% between 1972 and 1988. Declines were greatest on private land, least in wilderness, and intermediate in public nonwilderness. High elevations (greater than 914 m) maintained a greater percentage of CF than lower elevations (less than 914 m). The percentage of the area at the edge of the two cover types increased on all ownerships and in both elevational zones, whereas the amount of interior habitat (defined as CF at least 100 m from OT) decreased on all ownerships and elevational zones. By 1988 public lands contained approximately 45% interior habitat while private lands had 12% interior habitat. Mean interior patch area declined from 160 to 62 ha. The annual rate of disturbance (primarily clear-cutting) for the entire area including the wilderness was 1.19%, which corresponds to a cutting rotation of 84 yr. The forest landscape was not in a steady state or regulated condition which is not projected to occur for at least 40 yr under current forest plans. Variability in cutting rates within ownerships was higher on private land than on nonreserve public land. However, despite the use of dispersed cutting patterns on public land, spatial patterns of cutting and remnant forest patches were nonuniform across the entire public ownership. Large remaining patches (less than 5000 ha) of contiguous interior forest were restricted to public lands designated for uses other than timber production such as wilderness areas and research natural areas.

Spies, Thomas A.↗

Global Gridded Crop Model Evaluation: Benchmarking, Skills, Deficiencies and Implications.

Crop models are increasingly used to simulate crop yields at the global scale, but so far there is no general framework on how to assess model performance. Here we evaluate the simulation results of 14 global gridded crop modeling groups that have contributed historic crop yield simulations for maize, wheat, rice and soybean to the Global Gridded Crop Model Intercomparison (GGCMI) of the Agricultural Model Intercomparison and Improvement Project (AgMIP). Simulation results are compared to reference data at global, national and grid cell scales and we evaluate model performance with respect to time series correlation, spatial correlation and mean bias. We find that global gridded crop models (GGCMs) show mixed skill in reproducing time series correlations or spatial patterns at the different spatial scales. Generally, maize, wheat and soybean simulations of many GGCMs are capable of reproducing larger parts of observed temporal variability (time series correlation coefficients (r) of up to 0.888 for maize, 0.673 for wheat and 0.643 for soybean at the global scale) but rice yield variability cannot be well reproduced by most models. Yield variability can be well reproduced for most major producing countries by many GGCMs and for all countries by at least some. A comparison with gridded yield data and a statistical analysis of the effects of weather variability on yield variability shows that the ensemble of GGCMs can explain more of the yield variability than an ensemble of regression models for maize and soybean, but not for wheat and rice. We identify future research needs in global gridded crop modeling and for all individual crop modeling groups. In the absence of a purely observation-based benchmark for model evaluation, we propose that the best performing crop model per crop and region establishes the benchmark for all others, and modelers are encouraged to investigate how crop model performance can be increased. We make our evaluation system accessible to all crop modelers so that other modeling groups can also test their model performance against the reference data and the GGCMI benchmark.

wheat↗

Evaluation of a GCM cirrus parameterization using satellite observations

This study applies a simple yet effective methodology to validate a general circulation model parameterization of cirrus ice water path. The methodology combines large-scale dynamic and thermodynamic fields from operational analyses with prescribed occurrence of cirrus clouds from satellite observations to simulate a global distribution of ice water path. The predicted cloud properties are then compared with the corresponding satellite measurements of visible optical depth and infrared cloud emissivity to evaluate the reliability of the parameterization. This methodology enables the validation to focus strictly on the water loading side of the parameterization by eliminating uncertainties involved in predicting the occurrence of cirrus internally within the parameterization. Overall the parameterization performs remarkably well in capturing the observed spatial patterns of cirrus optical properties. Spatial correlations between the observed and the predicted optical depths are typically greater than 0.7 for the tropics and northern hemisphere midlatitudes. The good spatial agreement largely stems from the strong dependence of the ice water path upon the temperature of the environment in which the clouds form. Poorer correlations (r approximately 0.3) are noted over the southern hemisphere midlatitudes, suggesting that additional processes not accounted for by the parameterization may be important there. Quantitative evaluation of the parameterization is hindered by the present uncertainty in the size distribution of cirrus ice particles. Consequently, it is difficult to determine if discrepancies between the observed and the predicted optical properties are attributable to errors in the parameterized ice water path or to geographic variations in effective radii.

Soden, B. J.↗

Variation and evolution of atmospheric structure over the SEA BB season as seen from aircraft and reanalysis

The presence of absorbing aerosols and atmospheric water vapor will impact clouds both radiatively and dynamically. The atmosphere over the southeast Atlantic Ocean (SEA) experiences consistent springtime biomass burning (BB) smoke from widespread agricultural fires on the African continent. This smoke layer is initially lofted high in a continental mixed layer (~5-6km) and is then transported westward into the free troposphere where it overlies and ultimately mixes into the underlying stratocumulus-topped oceanic boundary layer. Coincident with this smoke is an elevated humidity signal which is present from the time a given airmass is over the continental source region. ORACLES (ObseRvations of Aerosols above CLouds and their intEractionS) was a NASA Earth Venture Suborbital mission with the goal of measuring aerosol, cloud, and atmospheric properties over this region over three deployments in September 2016, August 2017, and October 2018. Measurements included in situ profile measurements of aerosol physical and optical properties, trace gas concentrations, and cloud properties, as well as column measurements of aerosol optical depth and trace gas concentrations from an airborne sun photometer, allowing for calculation of radiative heating at different altitudes. Here we present results from the three ORACLES deployments over the southeast Atlantic Ocean. The airborne observations capture both seasonal patterns (August through October) and spatial patterns (north/south and east/west) in the vertical smoke distribution and its correlation to atmospheric specific humidity over the SEA basin. We describe the observed variations and show that the variability in both meteorological and emissions is largely captured in the European Centre for Medium-range Weather Forecasting Copernicus Atmosphere Monitoring Service (ECMWF CAMS) reanalysis, which allows for a more systematic examination of the relative frequency of observed conditions and their radiative impact over the region. Finally, we discuss the broader radiative and dynamical implications of these results for conditions of biomass burning aerosols overlying stratocumulus clouds.

Kristina Marie Myers Pistone↗

Large-Scale Circulation and Climate Variability

The causes of regional climate trends cannot be understood without considering the impact of variations in large-scale atmospheric circulation and an assessment of the role of internally generated climate variability. There are contributions to regional climate trends from changes in large-scale latitudinal circulation, which is generally organized into three cells in each hemisphere-Hadley cell, Ferrell cell and Polar cell-and which determines the location of subtropical dry zones and midlatitude jet streams. These circulation cells are expected to shift poleward during warmer periods, which could result in poleward shifts in precipitation patterns, affecting natural ecosystems, agriculture, and water resources. In addition, regional climate can be strongly affected by non-local responses to recurring patterns (or modes) of variability of the atmospheric circulation or the coupled atmosphere-ocean system. These modes of variability represent preferred spatial patterns and their temporal variation. They account for gross features in variance and for teleconnections which describe climate links between geographically separated regions. Modes of variability are often described as a product of a spatial climate pattern and an associated climate index time series that are identified based on statistical methods like Principal Component Analysis (PC analysis), which is also called Empirical Orthogonal Function Analysis (EOF analysis), and cluster analysis.

Perlwitz, J.↗