Search NASA⌕ Search

SEARCH · Search NASA

Results for “large ensembles”

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 469 records · Page 26

A Taxonomy-Based Approach to Shed Light on the Babel of Mathematical Models for Rice Simulation

For most biophysical domains, differences in model structures are seldom quantified. Here, we used a taxonomy-based approach to characterise thirteen rice models. Classification keys and binary attributes for each key were identified, and models were categorised into five clusters using a binary similarity measure and the unweighted pair-group method with arithmetic mean. Principal component analysis was performed on model outputs at four sites. Results indicated that (i) differences in structure often resulted in similar predictions and (ii) similar structures can lead to large differences in model outputs. User subjectivity during calibration may have hidden expected relationships between model structure and behaviour. This explanation, if confirmed, highlights the need for shared protocols to reduce the degrees of freedom during calibration, and to limit, in turn, the risk that user subjectivity influences model performance.

model parameterisation↗

Structure and Parameter Uncertainty in Centennial Projections of Forest Community Structure and Carbon Cycling

Secondary forest regrowth shapes community succession and biogeochemistry for decades, including in the Upper Great Lakes region. Vegetation models encapsulate our understanding of forest function, and whether models can reproduce multi‐decadal succession patterns is an indication of our ability to predict forest responses to future change. We test the ability of a vegetation model to simulate C cycling and community composition during 100 years of forest regrowth following stand‐replacing disturbance, asking (a) Which processes and parameters are most important to accurately model Upper Midwest forest succession? (b) What is the relative importance of model structure versus parameter values to these predictions? We ran ensembles of the Ecosystem Demography model v2.2 with different representations of processes important to competition for light. We compared the magnitude of structural and parameter uncertainty and assessed which sub‐model–parameter combinations best reproduced observed C fluxes and community composition. On average, our simulations underestimated observed net primary productivity (NPP) and leaf area index (LAI) after 100 years and predicted complete dominance by a single plant functional type (PFT). Out of 4,000 simulations, only nine fell within the observed range of both NPP and LAI, but these predicted unrealistically complete dominance by either early hardwood or pine PFTs. A different set of seven simulations were ecologically plausible but under‐predicted observed NPP and LAI. Parameter uncertainty was large; NPP and LAI ranged from ~0% to >200% of their mean value, and any PFT could become dominant. The two parameters that contributed most to uncertainty in predicted NPP were plant–soil water conductance and growth respiration, both unobservable empirical coefficients. We conclude that (a) parameter uncertainty is more important than structural uncertainty, at least for ED‐2.2 in Upper Midwest forests and (b) simulating both productivity and plant community composition accurately without physically unrealistic parameters remains challenging for demographic vegetation models.

canopy radiative transfer↗

Numerical Simulations of TRMM LBA, TOGA, COARE, GATE, ARM and PRESTORM Convective Systems: Sensitivity tests on Microphysical Processes

The 3D Goddard Cumulus Ensemble (GCE) model was utilized to examine the behavior and response of simulated deep tropical cloud systems that occurred over the west Pacific warm pool region, the Atlantic ocean and the central United States. The periods chosen for simulation were convectively active periods during TOGA-COARE (February 22 1993, December 11-17, 1992; December 19-28, February 9-13, 1993), GATE (September 4, 1974), LBA (January 26 and February 23, 1998), ARM (1997 IOP) and PRESTORM (June 11, 1985). We will examine differences in the microphysics for both warm rain and ice processes (evaporation /sublimation and condensation/ deposition), Q1 (Temperature), Q2 (Water vapor) and Q3 (momentum both U and V) budgets for these three convective events from different large-scale environments. The contribution of stratiform precipitation and its relationship to the vertical shear of the large-scale horizontal wind will also be examined. New improvements to the GCE model (i.e., microphysics: 4ICE two moments and 3ICE one moment; advection schemes) as well as their sensitivity to the model results will be discussed. Preliminary results indicated that various microphysical schemes could have a major impact on stratiform formation as well as the size of convective systems. However, they do not change the major characteristics of the convective systems, such as: arc shape, strong rotational circulation on both ends of system, heavy precipitation along the leading edge of systems.

Tao, W.-K.↗

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↗

Analysis of Buzz in a Supersonic Inlet

A dual-stream, low-boom supersonic inlet designed for use on a small, Mach 1.6 aircraft was tested experimentally in the 8- by 6-Foot Supersonic Wind Tunnel (SWT) at the NASA Glenn Research Center (GRC). The tests showed that the inlet had good recovery and stable operation over large mass flow range. The inlet went into buzz at mass flows well below that needed for engine operation, and the experiments generated a wealth of data during buzz. High frequency response pressure measurements and high-speed schlieren videos were recorded for many buzz events. The objective of the present work was to use computational fluid dynamics (CFD) to predict some of the experimental data taken during buzz, compare those predictions to the experimental data, and to use both datasets to explain the physics of the buzz cycle. The calculations were done with the Wind-US CFD code using a second-order time-accurate differencing scheme and the SST turbulence model. Computed Mach number contours were compared with schlieren images, and ensemble-averaged unsteady pressures were compared to data. The results showed that the buzz cycle consisted partly of spike buzz, an unsteady oscillation of the main shock at the spike tip while the inlet pressure dropped, and partly of choked flow while the inlet repressurized. Most of the results could be explained by theory proposed by Dailey in 1954, but did not support commonly used acoustic resonance explanations.

Chima, Rodrick V.↗

Subseasonal-to-Seasonal Hindcast Skill Assessment of Ridging Events Related to Drought Over the Western United States

Persistent atmospheric ridging events centered near the western United States are associated with widespread precipitation deficits and meteorological drought. Due to the relatively low skill of- dynamical models in forecasting precipitation on subseasonal‐to‐seasonal (S2S) time scales across this region, forecasts of ridging are explored in this study as a potential bridge for early warning drought prediction. To assess skill, we evaluate deterministic and probabilistic S2S hindcasts (out to 6 week lead time) of ridging events and geopotential height anomalies over the western United States in five ensemble hindcast systems. Prediction skill for ridging events is shown to be highly variable across models. For some models, longer‐time‐averaged patterns of geopotential height anomalies across the first 6 weeks can be skillfully simulated, when evaluated probabilistically against climatology. The most skillful models show modest skill in forecasting above normal ridging occurrences at lead times of Weeks 3–4 and Weeks 5–6,with some sensitivity to the specific ridge location and method for determining skill. Using the European Centre for Medium‐Range Weather Forecasts (ECMWF) model as a case study, longer lead time forecast busts are shown to often occur under extended periods of highly active La Nina‐like tropical convection. Model errors at simulating these tropical convection features may have a disproportionately large impact and degrade downstream forecasts over the western United States. Despite the documented model shortcomings, our results highlight an opportunity to target skillful ridging forecasts from dynamical models for improving early warning drought forecasting at S2S lead times.

S2S↗

SMOS Soil Moisture Data Assimilation in the NASA Land Information System: Impact on LSM Initialization and NWP Forecasts

Land surface models are important components of numerical weather prediction (NWP) models, partitioning incoming energy into latent and sensitive heat fluxes that affect boundary layer growth and destabilization. During warm-season months, diurnal heating and convective initiation depend strongly on evapotranspiration and available boundary layer moisture, which are substantially affected by soil moisture content. Therefore, to properly simulate warm-season processes in NWP models, an accurate initialization of the land surface state is important for accurately depicting the exchange of heat and moisture between the surface and boundary layer. In this study, soil moisture retrievals from the Soil Moisture and Ocean Salinity (SMOS) satellite radiometer are assimilated into the Noah Land Surface Model via an Ensemble Kalman Filter embedded within the NASA Land Information System (LIS) software framework. The output from LIS-Noah is subsequently used to initialize runs of the Weather Research and Forecasting (WRF) NWP model. The impact of assimilating SMOS retrievals is assessed by initializing the WRF model with LIS-Noah output obtained with and without SMOS data assimilation. The southeastern United States is used as the domain for a preliminary case study. During the summer months, there is extensive irrigation in the lower Mississippi Valley for rice and other crops. The irrigation is not represented in the meteorological forcing used to drive the LIS-Noah integration, but the irrigated areas show up clearly in the SMOS soil moisture retrievals, resulting in a case with a large difference in initial soil moisture conditions. The impact of SMOS data assimilation on both Noah soil moisture fields and on short-term (0-48 hour) WRF weather forecasts will be presented.

Blankenship, Clay↗

Larger Drought and Flood Hazards and Adverse Impacts on Population and Economic Productivity Under 2.0 than 1.5°C Warming

Climate change may have major influences on surface runoff, which would consequently result in important implications for terrestrial ecosystems and human well‐being. At global scale there is limited understanding of these issues with respect to the warming targets stipulated in the Paris Agreement. Here we use a well‐established hydrological model (Variable Infiltration Capacity [VIC]) forced with a representative ensemble of latest climate projections from four global circulation models (GCMs) to estimate potential future changes in runoff and Terrestrial Ecosystem Water Retention (TEWR), as well as changes in extreme runoff and their impacts on population, and overall gross domestic product (GDP) worldwide. Results suggest that annual runoff generally would have larger increases, while annual TEWR generally would have larger decreases under the 2.0°C warming scenario as opposed to 1.5°C warming scenario. Global mean warming of 2°C versus 1.5°C would lead to more distinct spatial patterns in runoff change, with a general shift of the runoff distribution towards more extreme low runoff in Mexico, western United States, Western Europe, southeastern China, West Siberian Plain and more extreme high runoff in Alaska, northern Canada, and large parts of Asia. More people and GDP would be exposed to extreme low runoff decrease, extreme high runoff increase, extreme low runoff decrease as well as extreme high runoff increase under a higher warming scenario. This study differentiates hydrological impacts between the two warming scenarios and illustrates higher runoff, lower TEWR, larger potential drought and flood hazards and adverse impacts on population and GDP under 2°C than 1.5°C.

Paris Agreement↗

Assessment of Arctic and Antarctic Sea Ice Predictability in CMIP5 Decadal Hindcasts

This paper examines the ability of coupled global climate models to predict decadal variability of Arctic and Antarctic sea ice. We analyze decadal hindcasts/predictions of 11 Coupled Model Intercomparison Project Phase 5 (CMIP5) models. Decadal hindcasts exhibit a large multimodel spread in the simulated sea ice extent, with some models deviating significantly from the observations as the predicted ice extent quickly drifts away from the initial constraint. The anomaly correlation analysis between the decadal hindcast and observed sea ice suggests that in the Arctic, for most models, the areas showing significant predictive skill become broader associated with increasing lead times. This area expansion is largely because nearly all the models are capable of predicting the observed decreasing Arctic sea ice cover. Sea ice extent in the North Pacific has better predictive skill than that in the North Atlantic (particularly at a lead time of 3-7 years), but there is a reemerging predictive skill in the North Atlantic at a lead time of 6-8 years. In contrast to the Arctic, Antarctic sea ice decadal hindcasts do not show broad predictive skill at any timescales, and there is no obvious improvement linking the areal extent of significant predictive skill to lead time increase. This might be because nearly all the models predict a retreating Antarctic sea ice cover, opposite to the observations. For the Arctic, the predictive skill of the multi-model ensemble mean outperforms most models and the persistence prediction at longer timescales, which is not the case for the Antarctic. Overall, for the Arctic, initialized decadal hindcasts show improved predictive skill compared to uninitialized simulations, although this improvement is not present in the Antarctic.

sea ice↗

A New Look into the Effect of Large Drops on Radiative Transfer Process

Recent studies indicate that a cloudy atmosphere absorbs more solar radiation than any current 1D or 3D radiation model can predict. The excess absorption is not large, perhaps 10-15 W/sq m or less, but any such systematic bias is of concern since radiative transfer models are assumed to be sufficiently accurate for remote sensing applications and climate modeling. The most natural explanation would be that models do not capture real 3D cloud structure and, as a consequence, their photon path lengths are too short. However, extensive calculations, using increasingly realistic 3D cloud structures, failed to produce photon paths long enough to explain the excess absorption. Other possible explanations have also been unsuccessful so, at this point, conventional models seem to offer no solution to this puzzle. The weakest link in conventional models is the way a size distribution of cloud particles is mathematically handled. Basically, real particles are replaced with a single average particle. This "ensemble assumption" assumes that all particle sizes are well represented in any given elementary volume. But the concentration of larger particles can be so low that this assumption is significantly violated. We show how a different mathematical route, using the concept of a cumulative distribution, avoids the ensemble assumption. The cumulative distribution has jumps, or steps, corresponding to the rarer sizes. These jumps result in an additional term, a kind of Green's function, in the solution of the radiative transfer equation. Solving the cloud radiative transfer equation with the measured particle distributions, described in a cumulative rather than an ensemble fashion, may lead to increased cloud absorption of the magnitude observed.

Marshak, Alexander↗

Impact of a Regional US Drought on Land and Atmospheric Carbon

The impacts of drought on regional land and atmospheric carbon are still poorly understood. Here we quantify the impact of a regional US drought on land carbon fluxes (Gross Primary Production, or GPP, and Net Biosphere Production, or NBP) and atmospheric carbon (CO2) by imposing an idealized 3-month meteorological drought in an ensemble of coupled land-atmosphere climate simulations. The imposed drought, applied to the lower Mississippi River Valley (~500,000 km (exp 2)), leads to a 23% GPP reduction in the drought area in the month immediately following the drought’s termination. The drought also caused GPP reductions in some remote areas through drought-induced impacts on remote meteorology, particularly the areas adjacent to the imposed drought. In the remote areas, the induced precipitation changes are responsible for most of the anomalous land productivity. The impact of the drought-induced meteorological anomalies on GPP is greater than that of the CO2 anomalies by at least an order of magnitude. While their impact on GPP is secondary, the drought-induced atmospheric CO2 anomalies near the land surface can be as large as 3.57ppm. The significant CO2 anomalies cover an area up to three times of that of the imposed drought, suggesting that atmospheric transport needs to be considered in the interpretation of drought-induced CO2 anomalies in the atmosphere. The imposed drought also leads to column-averaged CO2 increases of up to 0.78ppm, which is at the edge of the uncertainty from single soundings of current greenhouse gas (GHG) observing satellites.

Eunjee Lee↗

An Efficient GPU-Accelerated Multi-Source Global Fit Pipeline for LISA Data Analysis

The large-scale analysis task of deciphering gravitational wave signals in the LISA data stream will be difficult, requiring a large amount of computational resources and extensive development of computational methods. Its high dimensionality, multiple model types, and complicated noise profile require a global fit to all parameters and input models simultaneously. In this work, we detail our global fit algorithm, called “Erebor,” designed to accomplish this challenging task. It is capable of analysing current state-of-the-art datasets and then growing into the future as more pieces of the pipeline are completed and added. We describe our pipeline strategy, the algorithmic setup, and the results from our analysis of the LDC2A Sangria dataset, which contains Massive Black Hole Binaries, compact Galactic Binaries, and a parameterized noise spectrum whose parameters are unknown to the user. The Erebor algorithm includes three unique and very useful contributions: GPU acceleration for enhanced computational efficiency; ensemble MCMC sampling with multiple MCMC walkers per temperature for better mixing and parallelized sample creation; and special online updates to reversible-jump (or trans-dimensional) sampling distributions to ensure sampler mixing and accurate initial estimates for detectable sources in the data. We recover posterior distributions for all 15 (6) of the injected MBHBs in the LDC2A training (hidden) dataset. We catalog ∼12000 Galactic Binaries (∼8000 as high confidence detections) for both the training and hidden datasets. All of the sources and their posterior distributions are provided in publicly available catalogs.

LISA global fit↗

Efficient GPU-Accelerated MultiSource Global Fit Pipeline for LISA Data Analysis

The large-scale analysis task of deciphering gravitational-wave signals in the LISA data stream will be difficult, requiring a large amount of computational resources and extensive development of computational methods. Its high dimensionality, multiple model types, and complicated noise profile require a global fit to all parameters and input models simultaneously. In this work, we detail our global fit algorithm, called “Erebor,” designed to accomplish this challenging task. It is capable of analyzing current state-of-the-art datasets and then growing into the future as more pieces of the pipeline are completed and added. We describe our pipeline strategy, the algorithmic setup, and the results from our analysis of the LDC2A Sangria dataset, which contains massive black hole binaries, compact galactic binaries, and a parametrized noise spectrum whose parameters are unknown to the user. The Erebor algorithm includes three unique and very useful contributions: GPU acceleration for enhanced computational efficiency; ensemble Markov Chain Monte Carlo (MCMC) sampling with multiple MCMC walkers per temperature for better mixing and parallelized sample creation; and special online updates to reversible-jump (or transdimensional) sampling distributions to ensure sampler mixing and accurate initial estimates for detectable sources in the data.We recover posterior distributions for all 15 (6) of the injected massive black hole binaries (MBHB) in the LDC2A training (hidden) dataset. We catalog ∼12000 galactic binaries (∼8000 as high confidence detections) for both the training and hidden datasets. All of the sources and their posterior distributions are provided in publicly available catalogs.

LISA↗

Lower-Stratospheric Tracer Trends and Variability from Reanalyses

Previous studies have robustly identified a decrease since 1998 in lower stratospheric ozone in the Northern Hemisphere (NH). While this ozone decrease is qualitatively explained as resulting from changes in the large-scale circulation, there is not yet a quantitative mechanistic explanation of these changes. Here we explore the drivers of recent ozone changes using two different configurations of the Goddard Earth Observing System (GEOS) general circulation model. The first configuration of GEOS includes a full chemistry module and is constrained with meteorological fields from the MERRA-2 reanalysis. This configuration (M2GMI) is used to analyze an idealized tracer that covaries closely with ozone on interannual and decadal timescales, revealing that recent ozone decreases in the NH subtropics are associated with a poleward expansion of upwelling in the NH LS, with reduced (enhanced) downwelling over northern subtropics (midlatitudes). The second configuration of GEOS is a free-running version of the GEOS Chemistry-Climate Model (CCM) that is used to perform a ten-member ensemble of free-running simulations. Comparisons of the two configurations reveal that, while the free-running model can produce negative ozone changes in the NH LS, the magnitude of these changes is significantly weaker, relative to both M2GMI and MERRA-2; moreover, these weaker ozone decreases are consistent with weaker simulated changes in the residual circulation. Finally, we examine the GEOS model results in the broader context of the hindcast simulations performed as part of Phase 1 of the Chemistry Climate Modeling Initiative. We show that the majority of the free-running simulations considered here also exhibit weaker long-term residual circulation changes, compared to reanalyses.

Wargan, Kris↗

Tropical Oceanic Precipitation Processes Over Warm Pool: 2D and 3D Cloud Resolving Model Simulations

Rainfall is a key link in the hydrologic cycle as well as the primary heat source for the atmosphere. The vertical distribution of convective latent-heat release modulates the large-scale circulations of the topics. Furthermore, changes in the moisture distribution at middle and upper levels of the troposphere can affect cloud distributions and cloud liquid water and ice contents. How the incoming solar and outgoing longwave radiation respond to these changes in clouds is a major factor in assessing climate change. Present large-scale weather and climate model simulate processes only crudely, reducing confidence in their predictions on both global and regional scales. One of the most promising methods to test physical parameterizations used in General Circulation Models (GCMs) and climate models is to use field observations together with Cloud Resolving Models (CRMs). The CRMs use more sophisticated and physically realistic parameterizations of cloud microphysical processes, and allow for their complex interactions with solar and infrared radiative transfer processes. The CRMs can reasonably well resolve the evolution, structure, and life cycles of individual clouds and clouds systems. The major objective of this paper is to investigate the latent heating, moisture and momentum budgets associated with several convective systems developed during the TOGA COARE IFA - westerly wind burst event (late December, 1992). The tool for this study is the Goddard Cumulus Ensemble (GCE) model which includes a 3-class ice-phase microphysics scheme.

Tao, W.-K.↗

EnKF OSSE Experiments Assessing the Impact of HIRAD Wind Speed and HIWRAP Radial Velocity Data on Analysis of Hurricane Karl (2010)

Previous studies (e.g., Zhang et al. 2009, Weng et al. 2011) have shown that radial velocity data from airborne and ground-based radars can be assimilated into ensemble Kalman filter (EnKF) systems to produce accurate analyses of tropical cyclone vortices, which can reduce forecast intensity error. Recently, wind speed data from SFMR technology has also been assimilated into the same types of systems and has been shown to improve the forecast intensity of mature tropical cyclones. Two instruments that measure these properties were present during the NASA Genesis and Rapid Intensification Processes (GRIP) field experiment in 2010 which sampled Hurricane Karl, and will next be co-located on the same aircraft for the subsequent NASA HS3 experiment. The High Altitude Wind and Rain Profiling Radar (HIWRAP) is a conically scanning Doppler radar mounted upon NASAs Global Hawk unmanned aerial vehicle, and the usefulness of its radial velocity data for assimilation has not been previously examined. Since the radar scans from above with a fairly large fixed elevation angle, it observes a large component of the vertical wind, which could degrade EnKF analyses compared to analyses with data taken from lesser elevation angles. The NASA Hurricane Imaging Radiometer (HIRAD) is a passive microwave radiometer similar to SFMR, and measures emissivity and retrieves hurricane surface wind speeds and rain rates over a much wider swath. Thus, this study examines the impact of assimilating simulated HIWRAP radial velocity data into an EnKF system, simulated HIRAD wind speed, and HIWRAP+HIRAD with the Weather Research and Forecasting (WRF) model and compares the results to no data assimilation and also to the Truth from which the data was simulated for both instruments.

Albers, Cerese↗

Tropical Oceanic Precipitation Processes over Warm Pool: 2D and 3D Cloud Resolving Model Simulations

Rainfall is a key link in the hydrologic cycle as well as the primary heat source for the atmosphere, The vertical distribution of convective latent-heat release modulates the large-scale circulations of the tropics, Furthermore, changes in the moisture distribution at middle and upper levels of the troposphere can affect cloud distributions and cloud liquid water and ice contents. How the incoming solar and outgoing longwave radiation respond to these changes in clouds is a major factor in assessing climate change. Present large-scale weather and climate models simulate cloud processes only crudely, reducing confidence in their predictions on both global and regional scales. One of the most promising methods to test physical parameterizations used in General Circulation Models (GCMS) and climate models is to use field observations together with Cloud Resolving Models (CRMs). The CRMs use more sophisticated and physically realistic parameterizations of cloud microphysical processes, and allow for their complex interactions with solar and infrared radiative transfer processes. The CRMs can reasonably well resolve the evolution, structure, and life cycles of individual clouds and cloud systems, The major objective of this paper is to investigate the latent heating, moisture and momenti,im budgets associated with several convective systems developed during the TOGA COARE IFA - westerly wind burst event (late December, 1992). The tool for this study is the Goddard Cumulus Ensemble (CCE) model which includes a 3-class ice-phase microphysical scheme, The model domain contains 256 x 256 grid points (using 2 km resolution) in the horizontal and 38 grid points (to a depth of 22 km depth) in the vertical, The 2D domain has 1024 grid points. The simulations are performed over a 7 day time period. We will examine (1) the precipitation processes (i.e., condensation/evaporation) and their interaction with warm pool; (2) the heating and moisture budgets in the convective and stratiform regions; (3) the cloud (upward-downward) mass fluxes in convective and stratiform regions; (4) characteristics of clouds (such as cloud size, updraft intensity and cloud lifetime) and the comparison of clouds with Radar observations. Differences and similarities in organization of convection between simulated 2D and 3D cloud systems. Preliminary results indicated that there is major differences between 2D and 3D simulated stratiform rainfall amount and convective updraft and downdraft mass fluxes.

Tao, W.- K.↗

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↗