Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reanalysis”

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,477 records · Page 82

Improved Representation of the Global Dust Cycle Using Observational Constraints on Dust Properties and Abundance

Even though desert dust is the most abundant aerosol by mass in Earth’s atmosphere, atmospheric models 25 struggle to accurately represent its spatial and temporal distribution. These model errors are partially caused byfundamental difficulties in simulating dust emission in coarse-resolution models and in accurately representing dust microphysical properties. Here we mitigate these problems by developing a new methodology that yields an improved representation of the global dust cycle. We present an analytical framework that uses inverse modeling to integrate an ensemble of global model simulations with observational constraints on the dust size distribution, 30 extinction efficiency, and regional dust aerosol optical depth. We then compare the inverse model results against independent measurements of dust surface concentration and deposition flux and find that errors are reduced by approximately a factor of two relative to current model simulations of the Northern Hemisphere dust cycle. The inverse model results show smaller improvements in the less dusty Southern Hemisphere, most likely because both the model simulations and the observational constraints used in the inverse model are less accurate. On a global 35 basis, we find that the emission flux of dust with geometric diameter up to 20 μm (PM20) is approximately 5,000 Tg/year, which is greater than most models account for. This larger PM20 dust flux is needed to match observational constraints showing a large atmospheric loading of coarse dust. We obtain gridded data sets of dust emission, vertically integrated loading, dust aerosol optical depth, (surface) concentration, and wet and dry deposition fluxes that are resolved by season and particle size. As our results indicate that this data set is more accurate than current 40 model simulations and the MERRA-2 dust reanalysis product, it can be used to improve quantifications of dust impacts on the Earth system.

aerosols, dust↗

Large Contribution of Meteorological Factors to Inter-Decadal Changes in Regional Aerosol Optical Depth

Aerosol optical depth (AOD) has become a crucial metric for assessing global climate change. Although global and regional AOD trends have been studied extensively, it remains unclear what factors are driving the inter-decadal variations in regional AOD and how to quantify the relative contribution of each dominant factor. This study used a long-term (1980–2016) aerosol dataset from the Modern-Era Retrospective Analysis for Research and Applications, version 2 (MERRA-2) reanalysis, along with two satellite-based AOD datasets (MODIS/Terra and MISR) from 2001 to 2016, to investigate the long-term trends in global and regional aerosol loading. Statistical models based on emission factors and meteorological parameters were developed to identify the main factors driving the inter-decadal changes of regional AOD and to quantify their contribution. Evaluation of the MERRA-2 AOD with the ground-based measurements of AERONET indicated significant spatial agreement on the global scale (r= 0.85, root-mean-square error = 0.12, mean fractional error = 38.7 %, fractional gross error = 9.86 % and index of agreement = 0.94). However, when AOD observations from the China Aerosol Remote Sensing Network (CARSNET) were employed for independent verification, the results showed that MERRA-2 AODs generally underestimated CARSNET AODs in China (relative mean bias = 0.72 and fractional gross error =−34.3 %). In general, MERRA-2 was able to quantitatively reproduce the annual and seasonal AOD trends on both regional and global scales, as observed by MODIS/Terra, although some differences were found when compared to MISR. Over the 37-year period in this study, significant decreasing trends were observed over Europe and the eastern United States. In contrast, eastern China and southern Asia showed AOD increases, but the increasing trend of the former reversed sharply in the most recent decade. The statistical analyses suggested that the meteorological parameters explained a larger proportion of the AOD variability (20.4 %–72.8 %) over almost all regions of interest (ROIs) during 1980–2014 when compared with emission factors (0 %–56 %). Further analysis also showed that SO2 was the dominant emission factor, explaining 12.7 %–32.6 % of the variation in AOD over anthropogenic-aerosol-dominant regions, while black carbon or organic carbon was the leading factor over the biomass-burning-dominant (BBD) regions, contributing 24.0 %–27.7 % of the variation. Additionally, wind speed was found to be the leading meteorological parameter, explaining 11.8 %–30.3 % of the variance over the mineral-dust-dominant regions, while ambient humidity (including soil moisture and relative humidity) was the top meteorological parameter over the BBD regions, accounting for 11.7 %–35.5 % of the variation. The results of this study indicate that the variation in meteorological parameters is a key factor in determining the inter-decadal change in regional AOD.

Huizheng Che↗

Attribution of Chemistry-Climate Model Initiative (CCMI) Ozone Radiative Flux Bias from Satellites

The top-of-atmosphere (TOA) outgoing longwave flux over the 9.6-μm ozone band is a fundamental quantity for understanding chemistry-climate coupling. However, observed TOA fluxes are hard to estimate as they exhibit considerable variability in space and time that depend on the distributions of clouds, ozone (O3), water vapor (H2O), air temperature (Ta), and surface temperature (Ts). Benchmarking present day fluxes and quantifying the relative influence of their drivers is the first step for estimating climate feedbacks from ozone radiative forcing and predicting radiative forcing evolution. To that end, we constructed observational instantaneous radiative kernels (IRKs) under clear-sky conditions, representing the sensitivities of the TOA flux in the 9.6-μm ozone band to the vertical distribution of geophysical variables, including O3, H2O, Ta, and Ts based upon the Aura Tropospheric Emission Spectrometer (TES) measurements. Applying these kernels to present-day simulations from the Chemistry-Climate Model Initiative (CCMI) project as compared to a 2006 reanalysis assimilating satellite observations, we show that the models have large differences in TOA flux, attributable to different geophysical variables. In particular, model simulations continue to diverge from observations in the tropics, as reported in previous studies of the Atmospheric Chemistry Climate Model Inter-comparison Project (ACCMIP) simulations. The principal culprits are tropical mid and upper tropospheric ozone followed by tropical lower tropospheric H2O. Five models out of the eight studied here have TOA flux biases exceeding 100 mWm-2 attributable to tropospheric ozone bias. Another set of five models have flux biases over 50 mWm-2 due to H2O. On the other hand, Ta radiative bias is negligible in all models (no more than 30 mWm-2). We found that AM3 and CMAM have the lowest TOA flux biases globally but are a result of cancellation of opposite biases due to difference processes. Overall, the multi-model ensemble mean bias is –133±98 mWm-2, indicating that they are too atmospherically opaque due to trapping too much radiation in the atmosphere by overestimated tropical tropospheric O3 and H2O. Having too much O3 and H2O in the troposphere would have different impacts on the sensitivity of TOA flux to O3 and these competing effects add more uncertainties on the ozone radiative forcing. We find that the inter-model TOA outgoing longwave radiation (OLR) difference is well anti-correlated with their ozone band flux bias. This suggests that there is significant radiative compensation in the calculation of model outgoing longwave radiation.

Aura Tropospheric Emission Spectrometer (TES) meas↗

Impact of Regional Northern Hemisphere Mid-latitude Anthropogenic Sulfur Dioxide Emissions on Local and Remote Tropospheric Oxidants

The unintended consequences of reductions in regional anthropogenic sulfur dioxide (SO2) emissions implemented to protect human health are poorly understood. SO2 decreases began in the 1970s in the US and Europe and are expected to continue into the future, while recent emissions decreases in China are also projected to continue. In addition to the well-documented climate effects (warming) from reducing aerosols, tropospheric oxidation is impacted via aerosol modification of photolysis rates and radical sinks. Impacts on the hydroxyl radical and other trace constituents directly affect climate and air quality metrics such as surface ozone levels. We use the Geophysical Fluid Dynamics Laboratory Atmospheric Model version 3 nudged towards National Centers for Environmental Prediction (NCEP) reanalysis wind velocities to estimate the impact of SO2 emissions from the US, Europe, and China by differencing a control simulation with an otherwise identical simulation in which 2015 anthropogenic SO2 emissions are set to zero over one of the regions. Springtime sulfate aerosol changes occur both locally to the emission region and also throughout the Northern Hemispheric troposphere, including remote oceanic regions and the Arctic. Hydroperoxy (HO2) radicals are directly removed via heterogeneous chemistry on aerosol surfaces, including sulfate, in the model, and we find that sulfate aerosol produced by SO2 emissions from the three individual northern mid-latitude regions strongly reduces both HO2 and hydroxyl (OH) by up to 10 % year-round throughout most of the troposphere north of 30∘ N latitude. Regional SO2 emissions significantly increase nitrogen oxides (NOx) by about 5 %–8 % throughout most of the free troposphere in the Northern Hemisphere by increasing the NOx lifetime as the heterogeneous sink of HO2 on sulfate aerosol declines. Despite the NOx increases, tropospheric ozone decreases at northern mid-latitudes by 1 %–4 % zonally averaged and by up to 5 ppbv in summertime surface air over China, where the decreases in HO2 and OH suppress O3 production. Since 2015 anthropogenic SO2 emissions in China exceed those in the US or Europe, the oxidative response is greatest for the China perturbation simulation. Chemical effects of aerosols on oxidation (reactive uptake) dominate over radiative effects (photolysis rates), the latter of which are only statistically significant locally for the large perturbation over China. We find that the SO2 emissions decrease in China, which has yet to be fully realized, will have the largest impact on oxidants and related species in the Northern Hemisphere free troposphere compared to future decreases in Europe or the US. Our results bolster previous calls for a multipollutant strategy for air pollution mitigation to avoid the unintended consequence of aerosol removal leading to surface ozone increases that offset or mask surface ozone gains achieved by regulation of other pollutants, especially in countries where current usage of high-sulfur emitting fuels may be phased out in the future.

anthropogenic sulfur dioxide emissions↗

On Assessing ERA5 and MERRA2 Representations of Cold-Air Outbreaks Across the Gulf Stream

The warm Gulf Stream sea surface temperatures (SSTs) strongly impact the evolution of winter clouds behind atmospheric cold fronts. Such cloud evolution remains challenging to model. The Gulf Stream is too wide within the ERA5 and MERRA2 reanalyses, affecting the turbulent surface fluxes. Known problems within the ERA5 boundary layer (too-dry and too-cool with too strong westerlies), ascertained primarily from ACTIVATE 2020 campaign aircraft dropsondes and secondarily from older buoy measurements, reinforce surface flux biases. In contrast, MERRA2 winter surface winds and air-sea temperature/humidity differences are slightly too weak, producing surface fluxes that are too low. Reanalyses boundary layer heights in the strongly-forced winter cold-air-outbreak regime are realistic, whereas late-summer quiescent stable boundary layers are too shallow. Nevertheless, the reanalysis biases are small, and reanalyses adequately support their use for initializing higher-resolution cloud process modeling studies of cold-air outbreaks.

Gulf stream↗

A Brief Assessment of the Impact of Nearly 40 Years of Assimilated Observations Over the Amazon Basin

Adding a Forecast Sensitivity-based Observation Impact component to Version 2 of the Modern Era Retrospective-analysis for Research and Applications, the present study provides an assessment of the impact of nearly 40 years of observations on short-range (24-hr) forecasts over the Amazon basin. Under self-verification, forecast errors are found to slightly increase from the early data-sparse days to the more recent years, when data dramatically increase. Throughout the reanalysis, satellite radiances dominate in volume, but only before 1999 they dominate the impacts. Beyond 1999, over 50% of forecast error reduction is associated with conventional observations (radiosondes). Atmospheric Motion Vectors are also found to be large contributors to error reduction, but their contribution reduces in dry periods. In opposition to Atmospheric Motion Vectors, satellite radiances tend to contribute more in the dry season. Results provide motivation for additional conventional observations and the use of all-sky treatment of radiances.

F. L. R. Diniz↗

Effect of volcanic emissions on clouds during the 2008 and 2018 Kilauea degassing events

Volcanic eruptions in otherwise clean environments are “natural experiments” wherein the effects of aerosol emissions on clouds and climate can be partitioned from meteorological variability and anthropogenic activities. In this work, we combined satellite retrievals, reanalysis products, and atmospheric modeling to analyze the mechanisms of aerosol–cloud interactions during two degassing events at the Kilauea volcano in 2008 and 2018. The eruptive nature of the 2008 and 2018 degassing events was distinct from long-term volcanic activity for Kilauea. Although previous studies assessed the modulation of cloud properties from the 2008 event, this is the first time such an analysis has been reported for the 2018 event and that multiple degassing events have been analyzed and compared at this location. Both events resulted in significant changes in cloud effective radius and cloud droplet number concentration that were decoupled from local meteorology and in line with an enhanced cloud albedo. However, it is likely that the effects of volcanic emissions on liquid water path and cloud fraction were largely offset by meteorological variability. Comparison of cloud anomalies between the two events suggested a threshold response of aerosol–cloud interactions to overcome meteorological effects, largely controlled by aerosol loading. In both events, the ingestion of aerosols within convective parcels enhanced the detrainment of condensate in the upper troposphere, resulting in deeper clouds than observed under pristine conditions. Accounting for ice nucleation on ash particles led to enhanced ice crystal concentrations at cirrus levels and a slight decrease in ice water content, improving the correlation of the model results with the satellite retrievals. Overall, aerosol loading, plume characteristics, and meteorology contributed to changes in cloud properties during the Kilauea degassing events.

Katherine H. Breen↗

Improved representation of the global dust cycle using observational constraints on dust properties and abundance

Even though desert dust is the most abundant aerosol by mass in Earth's atmosphere, atmospheric models struggle to accurately represent its spatial and temporal distribution. These model errors are partially caused by fundamental difficulties in simulating dust emission in coarse-resolution models and in accurately representing dust microphysical properties. Here we mitigate these problems by developing a new methodology that yields an improved representation of the global dust cycle. We present an analytical framework that uses inverse modeling to integrate an ensemble of global model simulations with observational constraints on the dust size distribution, extinction efficiency, and regional dust aerosol optical depth. We then compare the inverse model results against independent measurements of dust surface concentration and deposition flux and find that errors are reduced by approximately a factor of 2 relative to current model simulations of the Northern Hemisphere dust cycle. The inverse model results show smaller improvements in the less dusty Southern Hemisphere, most likely because both the model simulations and the observational constraints used in the inverse model are less accurate. On a global basis, we find that the emission flux of dust with a geometric diameter up to 20 µm (PM20) is approximately 5000 Tg yr−1, which is greater than most models account for. This larger PM20 dust flux is needed to match observational constraints showing a large atmospheric loading of coarse dust. We obtain gridded datasets of dust emission, vertically integrated loading, dust aerosol optical depth, (surface) concentration, and wet and dry deposition fluxes that are resolved by season and particle size. As our results indicate that this dataset is more accurate than current model simulations and the MERRA-2 dust reanalysis product, it can be used to improve quantifications of dust impacts on the Earth system.

Dust cycle↗

FRB131104 Swift/BAT Data Revisited: No Evidence of a Gamma-Ray Counterpart

We present a reanalysis of the Swift Burst Alert Telescope (BAT) data around the radio detection of FRB131104. Possible evidence of a gamma-ray counterpart was presented by DeLaunay et al. However, based on our analysis using all the available BAT data, no significant emission is found in either the temporal or the image domain. We place a 5σ fluence upper limit of 3.3 × 10−6 erg cm−2 and 2.7 × 10−6 erg cm−2 (15–150 keV) with an integration time of 300 s assuming a simple power-law spectrum with photon index of −1.2 and −2.0, respectively. Our result does not support the association of this FRB with a high-energy counterpart, in agreement with growing observational evidence that most FRBs are not associated with catastrophic events such as gamma-ray bursts.

T Sakamoto↗

Relationship between circum-Arctic atmospheric wave patterns and large-scale wildfires in boreal summer

Long-term assessment of severe wildfires and associated air pollution and related climate patterns in and around the Arctic is essential for assessing healthy human life status. To examine the relationships, we analyzed the National Aeronautics and Space Administration (NASA) modern-era retrospective analysis for research and applications, version 2 (MERRA-2). Our investigation based on this state-of-the-art atmospheric reanalysis data reveals that 13 out of the 20 months with the highest PM2.5 (corresponding to the highly elevated organic carbon in the particulate organic matter [POM] form) monthly mean mass concentration over the Arctic for 2003–2017 were all in summer (July and August), during which POM of ⩾0.5 μg/cu. m and PM2.5 were positively correlated. This correlation suggests that high PM2.5 in the Arctic is linked to large wildfire contributions and characterized by significant anticyclonic anomalies (i.e. clockwise atmospheric circulation) with anomalous surface warmth and drier conditions over Siberia and subpolar North America, in addition to Europe. A similar climate pattern was also identified through an independent regression analysis for the July and August mean data between the same atmospheric variables and the sign-reversed Scandinavian pattern index. We named this pattern of recent atmospheric circulation anomalies the circum-Arctic wave (CAW) pattern as a manifestation of eastward group-velocity propagation of stationary Rossby waves (i.e. large-scale atmospheric waves). The CAW induces concomitant development of warm anticyclonic anomalies over Europe, Siberia, Alaska, and Canada, as observed in late June 2019. Surprisingly, the extended regression analysis of the 1980–2017 period revealed that the CAW pattern was not prominent before 2003. Understanding the CAW pattern under future climate change and global warming would lead to better prediction of co-occurrences of European heatwaves and large-scale wildfires with air pollution over Siberia, Alaska, and Canada in and around the Arctic in summer.

wildfire↗

Present and future aerosol impacts on Arctic climate change in the GISS-E2.1 Earth system model

The Arctic is warming 2 to 3 times faster than the global average, partly due to changes in short-lived climate forcers (SLCFs) including aerosols. In order to study the effects of atmospheric aerosols in this warming, recent past (1990–2014) and future (2015–2050) simulations have been carried out using the GISS-E2.1 Earth system model to study the aerosol burdens and their radiative and climate impacts over the Arctic (>60°N), using anthropogenic emissions from the Eclipse V6b and the Coupled Model Intercomparison Project Phase 6 (CMIP6) databases, while global annual mean greenhouse gas concentrations were prescribed and kept fixed in all simulations. Results showed that the simulations have underestimated observed surface aerosol levels, in particular black carbon (BC) and sulfate (SO2−4), by more than 50 %, with the smallest biases calculated for the atmosphere-only simulations, where winds are nudged to reanalysis data. CMIP6 simulations performed slightly better in reproducing the observed surface aerosol concentrations and climate parameters, compared to the Eclipse simulations. In addition, simulations where atmosphere and ocean are fully coupled had slightly smaller biases in aerosol levels compared to atmosphere-only simulations without nudging. Arctic BC, organic aerosol (OA), and SO2−4 burdens decrease significantly in all simulations by 10 %–60 % following the reductions of 7 %–78 % in emission projections, with the Eclipse ensemble showing larger reductions in Arctic aerosol burdens compared to the CMIP6 ensemble. For the 2030–2050 period, the Eclipse ensemble simulated a radiative forcing due to aerosol–radiation interactions (RFARI) of −0.39 ± 0.01 W/sq. m, which is −0.08 W/sq. m larger than the 1990–2010 mean forcing (−0.32 W/sq. m), of which −0.24 ± 0.01 W/sq. m was attributed to the anthropogenic aerosols. The CMIP6 ensemble simulated a RFARI of −0.35 to −0.40 W/sq. m for the same period, which is −0.01 to −0.06 W/sq. m larger than the 1990–2010 mean forcing of −0.35 W/sq. m. The scenarios with little to no mitigation (worst-case scenarios) led to very small changes in the RFARI, while scenarios with medium to large emission mitigations led to increases in the negative RFARI, mainly due to the decrease in the positive BC forcing and the decrease in the negative SO2−4 forcing. The anthropogenic aerosols accounted for −0.24 to −0.26 W/sq. m of the net RFARI in 2030–2050 period, in Eclipse and CMIP6 ensembles, respectively. Finally, all simulations showed an increase in the Arctic surface air temperatures throughout the simulation period. By 2050, surface air temperatures are projected to increase by 2.4 to 2.6 °C in the Eclipse ensemble and 1.9 to 2.6 °C in the CMIP6 ensemble, compared to the 1990–2010 mean. Overall, results show that even the scenarios with largest emission reductions leads to similar impact on the future Arctic surface air temperatures and sea-ice extent compared to scenarios with smaller emission reductions, implying reductions of greenhouse emissions are still necessary to mitigate climate change.

Aerosols↗

Response of the Quasi-Biennial Oscillation to historical volcanic eruptions

The impact of volcanic eruptions on surface climate is well-appreciated, but their in situ impact on the Quasi-Biennial Oscillation (QBO) has received comparatively little attention. This study examines the QBO responses to Krakatoa and Pinatubo using five configurations of the NASA Goddard Institute for Space Studies (GISS) Model E2.2 and MERRA-2 reanalysis. A dynamically consistent response is found in terms of static stability, zonal wind, and upwelling. Eruptions are found to bias the QBO towards a westerly state, such that the QBO period response depends upon the phase at the time of eruption. The QBO does not have a clear amplitude response to an eruption, based on these simulations. The underlying mechanisms appear not to be influenced to fi rst-order by interactive composition, sea surface temperatures, or long-term trends in CO2 and ozone-depleting substances.

volcanic eruptions↗

Preconditioning of overcast-to-broken cloud transitions by riming in marine cold air outbreaks

Marine cold air outbreaks (CAOs) commonly form overcast cloud decks that transition into broken cloud fields downwind, dramatically altering the local radiation budget. In this study, we investigate the impact of frozen hydrometeors on these transitions. We focus on a CAO case in the NW Atlantic, the location of the multi-year flight campaign ACTIVATE (Aerosol Cloud meTeorology Interactions oVer the western ATlantic Experiment). We use MERRA-2 (Modern-Era Retrospective analysis for Research and Applications, version 2) reanalysis fields to drive large eddy simulations with mixed-phase two-moment microphysics in a Lagrangian framework. We find that transitions are triggered by substantial rain (rainwater paths >25 g m−2), and only simulations that allow for aerosol depletion result in sustained breakups, as observed. Using a range of diagnostic ice nucleating particle concentrations, Ninp, we find that increasing ice progressively accelerates transitions, thus abbreviating the overcast state. Ice particles affect the cloud-topped boundary layer evolution, primarily through riming-related processes prior to substantial rain, leading to (1) a reduction in cloud liquid water, (2) early consumption of cloud condensation nuclei, and (3) early and light precipitation cooling and moistening below cloud. We refer to these three effects collectively as “preconditioning by riming”. Greater boundary layer aerosol concentrations available as cloud condensation nuclei (CCN) delay the onset of substantial rain. However, cloud breakup and low CCN concentration final stages are found to be inevitable in this case, due, primarily, to liquid water path buildup. An ice-modulated cloud transition speed suggests the possibility of a negative cloud–climate feedback. To address prevailing uncertainties in the model representation of mixed-phase processes, the magnitude of ice formation and riming impacts and, thereby, the strength of an associated negative cloud–climate feedback process, requires further observational evaluation by targeting riming hot spots with in situ imaging probes that allow for both the characterization of ice particles and abundance of supercooled droplets.

Marine cold air outbreaks↗

Tropical peatland hydrology simulated with a global land surface model

Tropical peatlands are among the most carbon-dense ecosystems on Earth, and their water storage dynamics strongly control these carbon stocks. The hydrological functioning of tropical peatlands differs from that of northern peatlands, which has not yet been accounted for in global land surface models (LSMs). Here, we integrated tropical peat-specific hydrology modules into a global LSM for the first time, by utilizing the peatland-specific model structure adaptation (PEATCLSM) of the NASA Catchment Land Surface Model (CLSM). We developed literature-based parameter sets for natural (PEATCLSM(Trop,Nat)) and drained (PEATCLSM(Trop,Drain)) tropical peatlands. Simulations with PEATCLSM(Trop,Nat) were compared against those with the default CLSM version and the northern version of PEATCLSM (PEATCLSM(North,Nat)) with tropical vegetation input. All simulations were forced with global meteorological reanalysis input data for the major tropical peatland regions in Central and South America, the Congo Basin, and Southeast Asia. The evaluation against a unique and extensive data set of in situ water level and eddy covariance-derived evapotranspiration showed an overall improvement in bias and correlation compared to the default CLSM version. Over Southeast Asia, an additional simulation with PEATCLSM(Trop,Drain) was run to address the large fraction of drained tropical peatlands in this region. PEATCLSM(Trop,Drain) outperformed CLSM, PEATCLSM(North,Nat) and PEATCLSM(Trop,Nat) over drained sites. Despite the overall improvements of PEATCLSM(Trop,Nat) over CLSM, there are strong differences in performance between the three study regions. We attribute these performance differences to regional differences in accuracy of meteorological forcing data, and differences in peatland hydrologic response that are not yet captured by our model.

Peatland↗

PBL Height from AIRS, GPS RO, and MERRA-2 Products in NASA GES DISC and Their 10 Year Seasonal Mean Intercomparison

Within the planetary boundary layer (PBL), surface forcing response, drag, turbulence, and vertical mixing are important processes and play a more critical role here than in the overlying “free atmosphere”. The PBL Height (PBLH) is an important parameter in climate models, weather forecasts, and air quality prediction. The NASA Goddard Earth Sciences Data and Information Services Center (GES DISC) provides data processing, archiving, and distribution services for numerous Earth science products. PBLH is a parameter in three products served by GES DISC, which are from the Atmospheric Infrared Sounder (AIRS), the Global Positioning System (GPS) radio occultation (RO) experiment, and the NASA reanalysis product Modern-Era Retrospective analysis for Research and Applications – 2 (MERRA-2). These products have different spatial and temporal resolutions and coverages, and their PBLH definitions are also different. To better serve the PBL research community, we have summarized the specifications of these products. A ten-year seasonal mean intercomparison is also conducted to provide further guidance to users. The intercomparison results show that MERRA-2 has a much shallower PBL than AIRS and GPS RO. An experimental study indicates the different PBLH definition in MERRA-2 caused smaller values of PBLH. The improvement of the water vapor retrieval in AIRS version 7 over version 6 results in the version 7 PBLH agreeing better with GPS RO and MERRA-2 than version 6, especially near the equator and low latitudes.

Feng Ding↗

The Arya Crop Yield Forecasting Algorithm: Application to the Main Wheat Exporting Countries

Wheat is the most important commodity traded in the international food market. Thus, accurate and timely information on wheat production can help mitigate food price fluctuations. Within the existing operational regional and global scale agricultural monitoring systems that provide information on global crop yield and area forecasts, there are still fundamental gaps: #1. Lack of quantitative Earth Observation (EO) derived crop information, #2. Lack of global but detailed (national or subnational level) and timely crop production forecasts and #3. Lack of information on forecast uncertainties. In this study we present the Agriculture Remotely-sensed Yield Algorithm (ARYA) an EO-based method, advancing the state of EO data application and usage (addressing gap #1) to forecast wheat yield. The algorithm is based on the evolution of the Difference Vegetation Index (DVI) using MODIS data at 1km resolution and the Growing Degree Days (GDD) from reanalysis data. Additionally, we explore how Land Surface Temperature (LST) can be included into the model and whether this parameter adds any value to the model performance when combined with the optical information. ARYA is implemented at the national and subnational level to forecast winter wheat yield in the main wheat exporting countries of US, Russia, Ukraine, France, Germany, Australia and Argentina from 2001 to 2019 (covering over 70% of wheat exports globally) in a timely manner by providing daily forecasts (addressing gap #2). The results show that ARYA provides yield estimations with RMSE’s within 0.3 ± 0.1 t/ha at national level and 0.6 ± 0,1 t/ha at subnational level after Day Of the Year (DOY) 140 (mid May) in the Northern Hemisphere and DOY 280 (beginning of October) in the Southern Hemisphere. This means that ARYA can provide crop yield estimates of wheat yield with 5-15 % error at national and 7-20 % error at subnational level starting from 2 to 2.5 months prior to harvest.

Agriculture↗

Evaluating Twenty-Year Trends in Earth’s Energy Flows from Observations

Recent satellite and ocean in-situ observations indicate a doubling of Earth’s Energy Imbalance (EEI) since 2005 that is a consequence of anthropogenic forcing, internal variability, and climate feedbacks. In this study, we follow-up with an assessment of how other components of Earth’s energy budget have changed since 2000 and how well these can be tracked. Specifically, we use top-of-atmosphere and surface observations from the Clouds and the Earth’s Radiant Energy System (CERES) and calculations of atmospheric divergence of energy transport from ERA-5 to determine the change in total surface flux along with its radiative and turbulent heat flux components. We also assess trends in ocean heat uptake for different ocean layers using various ocean reanalysis datasets. We compare trends derived using various observational and modeling approaches to assess how robust the derived trends are. Our analysis is restricted to global, hemispheric and regional spatial scales and annual mean and seasonal timescales.

CERES↗

Polar Vortex Outbreak Air Transport: Observation using Satellite IR Sounder Derived Ozone Product and Comparison with Model

The Single Field of View (SFOV) Sounder Atmospheric Products (SiFSAP) derived from Cross-track Infrared Sounder (CrIS) on SNPP and JPSS have a spatial resolution ( ~14.5 km) better than most global weather and climate models. Most recent significant improvement in the quality of these products enables us to use these high-resolution observation-based sounding data for weather studies and model evaluation. The aim of this study is to explore the value to use these SFOV products for studying the dynamic transport associated with polar vortex outbreak. In a few cases studies, we used SiFSAP products, especially ozone, as well as the data of wind, geopotential height (GPH) and potential velocity (PV) from the fifth-generation ECMWF reanalysis (ERA5) data, to analyze the variation of total ozone, relative humidity and ozone in different layers, as well as their correlation with GPH and PV. Some comparison of the retrieved temperature and water vapor with ERA5 data, and ozone with the Ozone Mapping and Profiler Suite (OMPS) measurements have also been made. It is found (1) the transport of O3-riched polar air accompanying the polar vortex breakup to lower latitude and lower atmosphere, and (2) the transport path of O3 agreeing well with the path of polar cold air (illustrated using the retrieved RH, wind fields and GPH). These results demonstrate the 3-D structure of O3 distribution as derived from CrIS measurements provides some insights of the cold air transport, and can be used to track the dynamic transport of polar cold air following the outbreak of polar vortex. The observed enhancement of O3 following the polar vortex outbreak also suggests some possible linkage between the tropospheric cold air advected with stratospheric air source.

Xiaozhen Xiong↗