Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reanalysis data”

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 631 records · Page 35

Spatial Extent of a Deep Moonquake Nest: A Preliminary Report of Reexamination

Deep moonquakes, occurring at depths about halfway to the center of the Moon, were discovered during the Apollo lunar landing missions, 1969-1972. Their near-monthly occurrence with nearly identical waveforms at any given seismic station suggests that they are strongly influenced by tides, caused by the Earth and the Sun, acting on certain limited regions of the deep lunar interior. However, much about them is unknown, why they are restricted to certain depths and to limited source regions (nests), and what they tell us about the nature of the material and dynamics of the interior of the Moon. A piece of information helpful to decipher their true nature is the spatial extent and distribution of their hypocenters. The occurrence of nearly identical waveforms suggests groups of hypocenters appear in close proximity to one another, but can we tell how closely they are located and how they are distributed? The Apollo PSE (Passive Seismic Experiment) data from all stations of the seismic network were received in real time at a common receiving station on Earth in digital form. This provided extremely high inter-station timing accuracy not achievable for most earthquake data on Earth at that time. We took advantage of this to compute relative locations of waveform-matched events, and concluded that deep moonquake foci in the A1 moonquake nest were concentrated on a nearly horizontal plane of less than 1 km in diameter. This earlier study was based on deep moonquake events visually identified on seismograms. A recent reanalysis of earlier unidentified seismic events fully utilizing the high capability of present-day computers expanded the list of positively identified deep moonquakes by more than a factor of five. The new list contains many events that are not visually matched in waveforms, yet correlated at a significant level when cross-correlated with a computer. Thus it became imperative to reexamine the spatial distribution of deep moonquake hypocenters including the newly identified events. Additional information is included in the original extended abstract.

Nakamura, Yosio↗

Evaporation from the Southern Ocean Estimated on the Basis of AIRS Satellite Data

Evaporation plays an important role in the global water and energy cycles and, hence, in climate change. Evaporation over the Southern Ocean, where the Antarctic sea ice coverage has a large annual cycle, is poorly quantified. In this study, daily evaporation is estimated for the Southern Ocean with a sea‐ice‐specific algorithm, using surface temperature and air humidity from National Aeronautics and Space Administration's Atmospheric Infrared Sounder (AIRS), and wind speeds from Modern‐Era Retrospective Analysis for Research and Applications, Version 2 (MERRA‐2), reanalysis during 2003–2016. An uncertainty of 34% was found in the evaporation product. The results indicate that annual evaporation has considerable interannual and regional variability, but with a decreasing trend during the study period over most of the Southern Ocean. There are, however, areas where evaporation has increased, specifically in the Ross Sea in winter and summer, with smaller positive trends in spring and fall. Overall, the changes in the difference between the surface specific humidity and the air specific humidity, and to a much lesser extent in the wind speed, are the main drivers for the changes in evaporation throughout the year. During spring and fall months, changes to the sea ice cover, which alter the surface specific humidity, are the main drivers for the change, but in summer and winter the main driver is the air‐specific humidity. Air masses originating from the Antarctic continent (south) are associated with cold and dry conditions, which increase evaporation, whereas air masses from lower latitudes in the Southern Ocean (north) are associated with warm and moist conditions, decreasing evaporation. Comparisons with other reanalysis evaporation products produce similar trends, although annual averages differ.

Boisvert, Linette↗

A Path to Improving Simulated Properties of Low Clouds over the Beaufort Sea using Airborne In Situ Observations of Subgrid-Scale Variability

Arctic low clouds influence the evolution of the Arctic system through their effects on radiative fluxes, boundary layer mixing, stability, turbulence, humidity, and precipitation. Unfortunately, atmospheric models and retrospective analysis (reanalysis) products struggle to accurately simulate the occurrence and properties of low clouds in the Arctic. One of the main reasons for this problem are the possible unrealistic assumptions that models/reanalyses make about the subgrid-scale (SGS) variability of meteorological properties, as well as the relationship between SGS variability and grid-scale (GS) cloud properties. We utilize cloud and thermodynamic data of low level (primarily) liquid clouds collected from two aircraft campaigns conducted over the Beaufort Sea to better understand and characterize this problem. Examining data from the September 2014 Arctic Radiation-IceBridge Sea and Ice Experiment (ARISE) airborne campaign reveals that GS cloud water variability is closely related with SGS distribution of total water (i.e. water vapor + cloud water). Multiple models/reanalyses use SGS supersaturated total water (relative to GS saturation) as a threshold for partitioning available water for condensation, and indeed we find significant correlation between GS cloud water and SGS supersaturation. However, we also find that the assumption of a static threshold of 100% saturation to be unrealistic. Empirical calculations from the ARISE data show a large sensitivity of this threshold to GS relative humidity, and so a microphysical parameterization allowing the threshold to vary according to GS thermodynamic properties may result in more realistic GS cloud water values. Finally, to determine how sensitive the ARISE-derived results are to that particular campaign, we include additional data from the First International Satellite Cloud Climatology Project (ISCCP) Regional Experiment (FIRE)–Arctic Cloud Experiment (ACE) conducted in 1998. The inclusion of the second dataset will help with demonstrating the robustness of the results and their utility in improving the representation of Arctic clouds in models and reanalyses.

J Brant Dodson↗

A Path to Improving Simulated Properties of Low Clouds over the Beaufort Sea using Airborne In Situ Observations of Subgrid-Scale Variability

Arctic low clouds influence the evolution of the Arctic system through their effects on radiative fluxes, boundary layer mixing, stability, turbulence, humidity, and precipitation. Unfortunately, atmospheric models and retrospective analysis (reanalysis) products struggle to accurately simulate the occurrence and properties of low clouds in the Arctic. One of the main reasons for this problem are the possible unrealistic assumptions that models/reanalyses make about the subgrid-scale (SGS) variability of meteorological properties, as well as the relationship between SGS variability and grid-scale (GS) cloud properties. We utilize cloud and thermodynamic data of low level (primarily) liquid clouds collected from two aircraft campaigns conducted over the Beaufort Sea to better understand and characterize this problem. Examining data from the September 2014 Arctic Radiation-IceBridge Sea and Ice Experiment (ARISE) airborne campaign and the 1998 First International Satellite Cloud Climatology Project (ISCCP) Regional Experiment (FIRE)–Arctic Cloud Experiment (ACE) reveals that GS cloud water variability is closely related with SGS distribution of total water (i.e. water vapor + cloud water). We investigate two related approaches to prediction of GS cloud properties from SGS variability: the critical saturation ratio method, and the critical relative humidity method. We find significant correlation between GS cloud water and SGS supersaturation when the critical saturation ratio is set at 100%, as well as a notable relationship between GS cloud water and the width of the SGS total water distribution. Critical relative humidity also compares well with GS cloud water. However, we also find that the assumptions of a static critical saturation ratio of 100% to be unrealistic, as well as a fixed SGS distribution width. Empirical calculations from the ARISE data show a large sensitivity of these SGS variables to GS relative humidity, and so a SGS parameterization allowing them to vary according to GS thermodynamic properties may result in more realistic GS cloud water values.

J. Brant Dodson↗

A Path to Improving Simulated Properties of Low Clouds over the Beaufort Sea using Airborne In Situ Observations of Subgrid-Scale Variability

Arctic low clouds influence the evolution of the Arctic system through their effects on radiative fluxes, boundary layer mixing, stability, turbulence, humidity, and precipitation. Unfortunately, atmospheric models and retrospective analysis (reanalysis) products struggle to accurately simulate the occurrence and properties of low clouds in the Arctic. One of the main reasons for this problem are the possible unrealistic assumptions that models/reanalyses make about the subgrid-scale (SGS) variability of meteorological properties, as well as the relationship between SGS variability and grid-scale (GS) cloud properties. We utilize cloud and thermodynamic data of low level (primarily) liquid clouds collected from two aircraft campaigns conducted over the Beaufort Sea to better understand and characterize this problem. Examining data from the September 2014 Arctic Radiation-IceBridge Sea and Ice Experiment (ARISE) airborne campaign and the 1998 First International Satellite Cloud Climatology Project (ISCCP) Regional Experiment (FIRE)–Arctic Cloud Experiment (ACE) reveals that GS cloud water variability is closely related with SGS distribution of total water (i.e. water vapor + cloud water). We examine the influence of three assumed SGS parameters on the estimation of GS cloud water: the width of the SGS PDF of total water (ΔQT), the shape of the SGS PDF, and the critical saturation ratio (Scrit), which partitions the SGS PDF into water that is available or not for conversion to cloud water. Both Scrit and ΔQT influence predicted GS cloud water strongly. ARISE and FIRE-ACE disagree somewhat in the details, but both campaigns support the possibility that more realistic representations of these two parameters may lead to more realistic GS cloud water. In particular, both parameters are sensitive to GS relative humidity, and so we examine the effects of allowing these quantities to vary as an empirically-derived linear function of GS relative humidity. In contrast, the shape of the PDF has little effect on the predicted cloud water – a surprising result that merits additional investigation.

J Brant Dodson↗

An Analysis of Coordinated Observations from NOAA's Ronald Brown Ship and G-IV Aircraft in a Landfalling Atmospheric River over the North Pacific During CalWater-2015

To gain a more complete observational understanding of atmospheric rivers (ARs) over the data-sparse open ocean, we used a diverse suite of mobile observing platforms deployed on NOAAs Ronald H. Brown (RHB) research vessel and G-IV research aircraft during the CalWater-2015 field campaign to describe the structure and evolution of a long-lived AR modulated by six frontal waves over the northeastern Pacific on 20-25 January 2015. Satellite observations and reanalysis diagnostics provided synoptic-scale context, illustrating the warm, moist southwesterly airstream within the quasi-stationary AR situated between an upper-level trough and ridge. The AR remained offshore of the U.S. West Coast but made landfall across British Columbia where heavy precipitation fell. Forty-seven rawinsondes launched from the RHB provided a comprehensive thermodynamic and kinematic depiction of the AR, including uniquely documenting an upward intrusion of strong water-vapor transport in the low-level moist southwesterly flow during the passage of frontal waves 2 through 6. A collocated 1290-MHz wind profiler showed an abrupt frontal transition from southwesterly to northerly flow below 1 km MSL coinciding with the tail-end of AR conditions. Shipborne radar and disdrometer observations in the AR uniquely captured key microphysical characteristics of shallow warm rain, convection, and deep mixed-phase precipitation. Novel observations of sea-surface fluxes in a midlatitude AR documented persistent ocean-surface evaporation and sensible-heat transfer into the ocean. The G-IV aircraft flew directly over the ship, with dropsonde and radar spatial analyses complementing the temporal depictions of the AR from the RHB. The AR characteristics varied, depending on the location of the cross section relative to the frontal waves.

atmospheric rivers↗

Evaluation of NASA's MERRA Precipitation Product in Reproducing the Observed Trend and Distribution of Extreme Precipitation Events in the United States

This study evaluates the performance of NASA's Modern-Era Retrospective Analysis for Research and Applications (MERRA) precipitation product in reproducing the trend and distribution of extreme precipitation events. Utilizing the extreme value theory, time-invariant and time-variant extreme value distributions are developed to model the trends and changes in the patterns of extreme precipitation events over the contiguous United States during 1979-2010. The Climate Prediction Center (CPC) U.S.Unified gridded observation data are used as the observational dataset. The CPC analysis shows that the eastern and western parts of the United States are experiencing positive and negative trends in annual maxima, respectively. The continental-scale patterns of change found in MERRA seem to reasonably mirror the observed patterns of change found in CPC. This is not previously expected, given the difficulty in constraining precipitation in reanalysis products. MERRA tends to overestimate the frequency at which the 99th percentile of precipitation is exceeded because this threshold tends to be lower in MERRA, making it easier to be exceeded. This feature is dominant during the summer months. MERRA tends to reproduce spatial patterns of the scale and location parameters of the generalized extreme value and generalized Pareto distributions. However, MERRA underestimates these parameters, particularly over the Gulf Coast states, leading to lower magnitudes in extreme precipitation events. Two issues in MERRA are identified: 1) MERRA shows a spurious negative trend in Nebraska and Kansas, which is most likely related to the changes in the satellite observing system over time that has apparently affected the water cycle in the central United States, and 2) the patterns of positive trend over the Gulf Coast states and along the East Coast seem to be correlated with the tropical cyclones in these regions. The analysis of the trends in the seasonal precipitation extremes indicates that the hurricane and winter seasons are contributing the most to these trend patterns in the southeastern United States. In addition, the increasing annual trend simulated by MERRA in the Gulf Coast region is due to an incorrect trend in winter precipitation extremes.

MERRA↗

Planetary Boundary Layer Height Estimates From ICESat-2 and CATS Backscatter Measurements

The lowest layer of the atmosphere in which all human activity occurs is called the Planetary Boundary Layer (PBL). All physical interactions with the surface, such as heat and moisture transport, pollution dispersion and transport happen in this relatively shallow layer. The ability to understand and model the complex interactions that occur in the PBL is very important to air quality, weather prediction and climate modeling. A fundamental and physically important property of the PBL is its thickness or height. This work presents two methods to obtain global PBL height using satellite lidar data from the Ice, Cloud and land Elevation Satellite-2 (ICESat-2) and the Cloud-Aerosol Transport System (CATS). The first method is a straightforward backscatter threshold technique and the second is a machine learning approach known as a Convolutional Neural Network. The PBL height retrievals from the two methods are compared with each other and with PBL height from the NASA GEOS MERRA-2 reanalysis. The lidar-retrieved PBL heights have a high degree of spatial correlation with the model heights but are generally higher over ocean (∼400 m) and over northern hemisphere high latitude regions (∼1,000 m). Over mid-latitude and tropical land areas, the satellite estimated PBL heights agree well with model mid-day estimates. This work demonstrates the feasibility of using satellite lidar backscatter measurements to obtain global PBL height estimates, as well as determining seasonal and regional variability of PBL height.

Stephen P Palm↗

Magmatic Processes That Produced Lunar Fire Fountains

Reanalysis of the Apollo 15 A, B, and C green glass beads from slide 15426,72, led to the discovery of patchy, highly vesicular glass rims adhering to beads. These rims are high in S and Ni and low in MgO, but otherwise compositionally similar to the green glasses. We find that these rims represent a unique melt composition that places constraints on lunar magmatic processes. Combining the compositional data for the vesicular glass rims with new minor element data for the green glass beads leads us to hypothesize that the vesicular rim fluid originated near the green glass source, at about 2.2 GPa. These findings support the theory of a heterogeneous lunar mantle, and suggest that sulfur did not drive the eruption from depth, but that degassing volatiles into the vacuum did drive the final fire fountain eruption.

Elkins-Tanton, Linda T.↗

Antarctic Sea Ice Thickness and Snow-to-Ice Conversion from Atmospheric Reanalysis and Passive Microwave Snow Depth

Passive microwave snow depth, ice concentration, and ice motion estimates are combined with snowfall from the European Centre for Medium Range Weather Forecasting (ECMWF) reanalysis (ERA-40) from 1979-200 1 to estimate the prevalence of snow-to-ice conversion (snow-ice formation) on level sea ice in the Antarctic for April-October. Snow ice is ubiquitous in all regions throughout the growth season. Calculated snow- ice thicknesses fall within the range of estimates from ice core analysis for most regions. However, uncertainties in both this analysis and in situ data limit the usefulness of snow depth and snow-ice production to evaluate the accuracy of ERA-40 snowfall. The East Antarctic is an exception, where calculated snow-ice production exceeds observed ice thickness over wide areas, suggesting that ERA-40 precipitation is too high there. Snow-ice thickness variability is strongly controlled not just by snow accumulation rates, but also by ice divergence. Surprisingly, snow-ice production is largely independent of snow depth, indicating that the latter may be a poor indicator of total snow accumulation. Using the presence of snow-ice formation as a proxy indicator for near-zero freeboard, we examine the possibility of estimating level ice thickness from satellite snow depths. A best estimate for the mean level ice thickness in September is 53 cm, comparing well with 51 cm from ship-based observations. The error is estimated to be 10-20 cm, which is similar to the observed interannual and regional variability. Nevertheless, this is comparable to expected errors for ice thickness determined by satellite altimeters. Improvement in satellite snow depth retrievals would benefit both of these methods.

Markus, Thorsten↗

The Role of DYNAMO in Situ Observations in Improving NASA Ceres-like Daily Surface and Atmospheric Radiative Flux Estimates

The daily surface and atmospheric radiative fluxes from NASA Clouds and the Earths RadiantEnergy System (CERES) Synoptic 1 degree (SYN1deg) Ed3A are among the most widely used data to studycloud-radiative feedback. The CERES SYN1deg data are based on Fu-Liou radiative transfer computations thatuse specific humidity (Q) and air temperature (T) from NASA Global Modeling and Assimilation Office (GMAO)reanalyses as inputs and are therefore subject to the quality of those fields. This study uses in situ Q and Tobservations collected during the Dynamics of the Madden-Julian Oscillation (DYNAMO) field campaign toaugment the input stream used in the NASA GMAO reanalysis and assess the impact on the CERES dailysurface and atmospheric longwave estimates. The results show that the assimilation of DYNAMOobservations considerably improves the vertical profiles of analyzed Q and T over and near DYNAMO stationsby moistening and warming the lower troposphere and upper troposphere and drying and cooling themid-upper troposphere. As a result of these changes in Q and T, the computed CERES daily surface downwardlongwave flux increases by about 5 W m(exp -2), due mainly to the warming and moistening in the lowertroposphere; the computed daily top-of-atmosphere (TOA) outgoing longwave radiation increases by2-3 W m(exp -2) during dry periods only. Correspondingly, the estimated local atmospheric longwave radiativecooling enhances by about 5 W m(exp -2) (7-8 W m(exp -2)) during wet (dry) periods. These changes reduce the bias inthe CERES SYN1deg-like daily longwave estimates at both the TOA and surface and represent animprovement over the DYNAMO region.

DYNAMO↗

A reanalysis of rocket measurements of the ultraviolet dayglow

Rocket measurements of O I 989, 1304, 1356 A and N2 Lyman-Birge-Hopfield band emission in the midlatitude dayglow reported by Gentieu et al. (1979) and Eastes et al. (1985) are reexamined. MSIS-83 atomic oxygen densities, the 1304 and 1356 A excitation cross sections of Zipf and Erdman (1985), and SMM solar 1304 A irradiance measurements are consistent with the observed O I and N2 emission intensities. Atmosphere Explorer E (AE-E) measurements of the solar EUV irradiance near solar maximum are consistent with the 1980 rocket airglow data, but the solar EUV irradiance required to explain the 1978 airglow data is a factor of 1.5 larger than indicated by AE-E. Enhancement of the 1304 A excitation cross section due to radiative entrapment of cascade-feeding photons is much less than the factor of two predicted by the cascade model of Julienne and Davis (1976).

Link, R.↗

The Relation Between Atmospheric Humidity and Temperature Trends for Stratospheric Water

We analyze the relation between atmospheric temperature and water vapor-a fundamental component of the global climate system-for stratospheric water vapor (SWV). We compare measurements of SWV (and methane where available) over the period 1980-2011 from NOAA balloon-borne frostpoint hygrometer (NOAA-FPH), SAGE II, Halogen Occultation Experiment (HALOE), Microwave Limb Sounder (MLS)/Aura, and Atmospheric Chemistry Experiment Fourier Transform Spectrometer (ACE-FTS) to model predictions based on troposphere-to-stratosphere transport from ERA-Interim, and temperatures from ERA-Interim, Modern Era Retrospective-Analysis (MERRA), Climate Forecast System Reanalysis (CFSR), Radiosonde Atmospheric Temperature Products for Assessing Climate (RATPAC), HadAT2, and RICHv1.5. All model predictions are dry biased. The interannual anomalies of the model predictions show periods of fairly regular oscillations, alternating with more quiescent periods and a few large-amplitude oscillations. They all agree well (correlation coefficients 0.9 and larger) with observations for higherfrequency variations (periods up to 2-3 years). Differences between SWV observations, and temperature data, respectively, render analysis of the model minus observation residual difficult. However, we find fairly well-defined periods of drifts in the residuals. For the 1980s, model predictions differ most, and only the calculation with ERA-Interim temperatures is roughly within observational uncertainties. All model predictions show a drying relative to HALOE in the 1990s, followed by a moistening in the early 2000s. Drifts to NOAA-FPH are similar (but stronger), whereas no drift is present against SAGE II. As a result, the model calculations have a less pronounced drop in SWV in 2000 than HALOE. From the mid-2000s onward, models and observations agree reasonably, and some differences can be traced to problems in the temperature data. These results indicate that both SWV and temperature data may still suffer from artifacts that need to be resolved in order to answer the question whether the large-scale flow and temperature field is sufficient to explain water entering the stratosphere.

Fueglistaler, S.↗

Projections of Declining Surface-Water Availability for the Southwestern United States

16 of the CMIP5 models had all the data needed for this work for at least one simulation that was continuous from 1950 to 2040. Details of the models analyzed here are provided in Table S1. The model data analyzed here are available at http://strega.ldeo.columbia.edu:81/expert/home/.naomi/.AR5/.v2/.historical:rcp85/.mmm16/ a. Assessing the climatology of the models Despite increases in horizontal resolution of many models compared to their CMIP3 counterparts none of these models can adequately resolve the topography of the south west United States, such as the Sierra Nevada and Rocky Mountains and the associated orographic precipitation. This requires that caution be used when interpreting the results presented here. To assess the ability of the models to simulate the current hydroclimate, in Figure S1 we show the observed (from the Global Precipitation Climatology Centre gridded rain gauge data, (1)) monthly climatology of precipitation and the same for all the models and the multimodel mean for the California-Nevada, Colorado headwaters and Texas regions. The GPCC data uses rain gauges only and interpolates to regular grids of which we used the 1◦ by 1◦ one. Details of the data set can be found in (2). While the models apparently overestimate precipitation in California and Nevada the seasonal cycle with wet winters and dry summers is very well represented. It is also possible that the rain gauge observations are biased low by inadequately sampling the higher mountain regions. How ever the models might also be expected to underestimate orographic precipitation due to inadequate horizontal resolution. The 25 models are also too wet in the Colorado headwaters region but correctly represent the quite even distribution though the year. The bimodal distribution of precipitation in Texas, with peaks in May and September, and the absolute amounts, are well modeled but with the September peak too weak. The positive precipitation bias translates into a positive runoff bias for the Colorado headwaters as also shown in Figure S1. Here the observed runoff values are taken from simulations of the Variable Infiltration Capacity (VIC) land surface-hydrology model (3) forced by observed meteorology (5) that were conducted as part of the North American Land Data Assimilation System project phase 2 ( (NLDAS-2), http://www.emc.ncep.noaa.gov/mmb/nldas/. Runoff for California-Nevada is better simulated but there is a positive bias over Texas despite no strong precipitation bias. To check whether regional climate models better simulate P and runoff in these regions we analyzed the historical simulation with the Regional Climate Model version 3 driven by the National Centers for Environmental Prediction-Department of Energy Reanalysis 2 available from the North American Regional Climate Change Assessment Program (http://www.narccap.ucar.edu). This model configuration retained these biases in P and runoff although they were reduced in amplitude. Given these varying biases we plot P and P − E changes in actual values but apply the simplest bias correction possible to the runoff and soil moisture values and show the modeled changes in terms of percentages of the 20th Century model climatologies. A thorough assessment of the simulation of North American climate in CMIP5 models is conducted in Sheffield at al. (North American Climate in CMIP5 Experiments. Part I: Evaluation of 20th Century Continental and Regional Climatology, manuscript submit ted to J. Climate, available at http://www.climate.noaa.gov/index.jsp?pg=./cpo pa/ mapp/cmip5 publications.html). Sheffield et al. analyze the climatology of precipitation, surface air temperature, low level winds, moisture fluxes, runoff etc. and conclude that the main features of the hydrological cycle, including characteristics of the atmospheric moisture balance and its seasonality, are captured in the CMP5 models subject to biases in total precipitation amounts. We chose to use all available models instead of selecting some and rejecting others based on an assessment of model realism. This is in accord with the suggestions of Mote et al. for CMIP3 (4) but future work needs to revisit this matter for the case of the CMIP5 ensemble.

Seager, Richard↗

Global Ocean Evaporation Increases Since 1960 in Climate Reanalyses: How Accurate Are They?

AGCMs w/ Specified SSTs (AMIPs) GEOS-5, ERA-20CM Ensembles Incorporate best historical estimates of SST, sea ice, radiative forcing Atmospheric "weather noise" is inconsistent with specified SST. Instantaneous Sfc fluxes can be wrong sign (e.g. Indian Ocean Monsoon, high latitude oceans). Averaging over ensemble members helps isolate SST-forced signal. Reduced Observational Reanalyses: NOAA 20CR V2C, ERA-20C, JRA-55C Incorporate observed Sfc Press (20CR), Marine Winds (ERA-20C) and rawinsondes (JRA-55C) to recover much of true synoptic or weather w/o shock of new sat obs. Comprehensive Reanalyses (MERRA-2) Full suite of observational constraints- both conventional and remote sensing. But... substantial uncertainties owing to evolving satellite observing system. Multi-source Statistically Blended OAFlux, LargeYeager Blend reanalysis, satellite, and ocean buoy information. While climatological biases are removed, non-physical trends or variations in components remain. Satellite Retrievals GSSTF3, SeaFlux, HOAPS3... Global coverage. Retrieved near sfc wind speed, & humidity used with SST to drive accurate bulk aerodynamic flux estimates. Satellite inter-calibration, spacecraft pointing variations crucial. Short record ( late 1987-present). In situ Measurements ICOADS, IVAD, Res Cruises VOS and buoys offer direct measurements. Sparse data coverage (esp south of 30S. Changes in measurement techniques (e.g. shipboard anemometer height).

Climate↗

Estimating Evapotranspiration Using an Observation Based Terrestrial Water Budget

Evapotranspiration (ET) is difficult to measure at the scales of climate models and climate variability. While satellite retrieval algorithms do exist, their accuracy is limited by the sparseness of in situ observations available for calibration and validation, which themselves may be unrepresentative of 500m and larger scale satellite footprints and grid pixels. Here, we use a combination of satellite and ground-based observations to close the water budgets of seven continental scale river basins (Mackenzie, Fraser, Nelson, Mississippi, Tocantins, Danube, and Ubangi), estimating mean ET as a residual. For any river basin, ET must equal total precipitation minus net runoff minus the change in total terrestrial water storage (TWS), in order for mass to be conserved. We make use of precipitation from two global observation-based products, archived runoff data, and TWS changes from the Gravity Recovery and Climate Experiment satellite mission. We demonstrate that while uncertainty in the water budget-based estimates of monthly ET is often too large for those estimates to be useful, the uncertainty in the mean annual cycle is small enough that it is practical for evaluating other ET products. Here, we evaluate five land surface model simulations, two operational atmospheric analyses, and a recent global reanalysis product based on our results. An important outcome is that the water budget-based ET time series in two tropical river basins, one in Brazil and the other in central Africa, exhibit a weak annual cycle, which may help to resolve debate about the strength of the annual cycle of ET in such regions and how ET is constrained throughout the year. The methods described will be useful for water and energy budget studies, weather and climate model assessments, and satellite-based ET retrieval optimization.

Rodell, Matthew↗

Assimilation of Aerosol Observations in the NASA GEOS Model

In the GEOS near real-time system, as well as in MERRA-2 which is the latest reanalysis produced at NASA's Global Modeling Assimilation Office (GMAO), the assimilation of aerosol observations is performed by means of a so-called analysis splitting method. The prognostic model is based on the GEOS model radiatively coupled to GOCART aerosol module and includes assimilation of bias-corrected Aerosol Optical Depth (AOD) at 550 nm from various space-based remote sensing platforms. In line with the transition of the GEOS meteorological data assimilation system to a hybrid Ensemble-Variational formulation, we are updating the aerosol component of our assimilation system to a variational ensemble type of scheme. In this talk we will examine the impact of replacing the current analysis splitting scheme with this new approach. Starting with the assimilation of satellite based single-channel retrievals; we will discuss the impact of this aerosol data assimilation technique on the 3D aerosol distributions by means of innovation statistics and verification against independent datasets such as the Aerosol Robotic Network (AERONET) and surface PM2.5. We will also present preliminary results related to the introduction of new aerosol data types in GEOS, including multi-spectral AOD retrievals.

AERONET↗

A Comprehensive Machine Learning Study to Classify Precipitation Type over Land from Global Precipitation Measurement Microwave Imager (GPM-GMI) Measurements

Precipitation type is a key parameter used for better retrieval of precipitation characteristics as well as to understand the cloud–convection–precipitation coupling processes. Ice crystals and water droplets inherently exhibit different characteristics in different precipitation regimes (e.g., convection, stratiform), which reflect on satellite remote sensing measurements that help us distinguish them. The Global Precipitation Measurement (GPM) Core Observatory’s microwave imager (GMI) and dual-frequency precipitation radar (DPR) together provide ample information on global precipitation characteristics. As an active sensor, the DPR provides an accurate precipitation type assignment, while passive sensors such as the GMI are traditionally only used for empirical understanding of precipitation regimes. Using collocated precipitation type flags from the DPR as the “truth”, this paper employs machine learning (ML) models to train and test the predictability and accuracy of using passive GMI-only observations together with ancillary information from a reanalysis and GMI surface emissivity retrieval products. Out of six ML models, four simple ones (support vector machine, neural network, random forest, and gradient boosting) and the 1-D convolutional neural network (CNN) model are identified to produce 90–94% prediction accuracy globally for five types of precipitation (convective, stratiform, mixture, no precipitation, and other precipitation), which is much more robust than previous similar effort. One novelty of this work is to introduce data augmentation (subsampling and bootstrapping) to handle extremely unbalanced samples in each category. A careful evaluation of the impact matrices demonstrates that the polarization difference (PD), brightness temperature (Tc) and surface emissivity at high-frequency channels dominate the decision process, which is consistent with the physical understanding of polarized microwave radiative transfer over different surface types, as well as in snow and liquid clouds with different microphysical properties. Furthermore, the view-angle dependency artifact that the DPR’s precipitation flag bears with does not propagate into the conical-viewing GMI retrievals. This work provides a new and promising way for future physics-based ML retrieval algorithm development.

machine learning/artificial intelligence↗