Evaluation of boundary layer cloud parameterizations in the ECHAM5 general circulation model using CALIPSO and CloudSat satellite data
Explore the source record for details and available documents.
SEARCH · Search NASA
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.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Retrievals of the isotopic composition of water vapor from the Aura Tropospheric Emission Spectrometer (TES) have unique value in constraining moist processes in climate models. Accurate comparison between simulated and retrieved values requires that model profiles that would be poorly retrieved are excluded, and that an instrument operator be applied to the remaining profiles. Typically, this is done by sampling model output at satellite measurement points and using the quality flags and averaging kernels from individual retrievals at specific places and times. This approach is not reliable when the model meteorological conditions influencing retrieval sensitivity are different from those observed by the instrument at short time scales, which will be the case for free-running climate simulations. In this study, we describe an alternative, categorical approach to applying the instrument operator, implemented within the NASA GISS ModelE general circulation model. Retrieval quality and averaging kernel structure are predicted empirically from model conditions, rather than obtained from collocated satellite observations. This approach can be used for arbitrary model configurations, and requires no agreement between satellite-retrieved and model meteorology at short time scales. To test this approach, nudged simUlations were conducted using both the retrieval-based and categorical operators. Cloud cover, surface temperature and free-tropospheric moisture content were the most important predictors of retrieval quality and averaging kernel structure. There was good agreement between the D fields after applying the retrieval-based and more detailed categorical operators, with increases of up to 30 over the ocean and decreases of up to 40 over land relative to the raw model fields. The categorical operator performed better over the ocean than over land, and requires further refinement for use outside of the tropics. After applying the TES operator, ModelE had D biases of 8 over ocean and 34 over land compared to TES D, which were less than the biases using raw model D fields.
Climate models have traditionally only represented heat and water fluxes within relatively shallow soil layers, but there is increasing interest in the possible role of heat and water exchanges with the deeper subsurface. Here, we integrate an idealized 50m deep aquifer into the land surface module of the GISS ModelE general circulation model to test the influence of aquifer-soil moisture and heat exchanges on climate variables. We evaluate the impact on the modeled climate of aquifer-soil heat and water fluxes separately, as well as in combination. The addition of the aquifer to ModelE has limited impact on annual-mean climate, with little change in global mean land temperature, precipitation, or evaporation. The seasonal amplitude of deep soil temperature is strongly damped by the soil-aquifer heat flux. This not only improves the model representation of permafrost area but propagates to the surface, resulting in an increase in the seasonal amplitude of surface air temperature of >1K in the Arctic. The soil-aquifer water and heat fluxes both slightly decrease interannual variability in soil moisture and in landsurface temperature, and decrease the soil moisture memory of the land surface on seasonal to annual timescales. The results of this experiment suggest that deepening the modeled land surface, compared to modeling only a shallower soil column with a no-flux bottom boundary condition, has limited impact on mean climate but does affect seasonality and interannual persistence.
We formulate a new second-order closure turbulence model by employing a recent closure for the pressure-temperature correlation at the equation level. As a result, we obtain new heat flux equations that avoid the long-standing issue of a finite critical Richardson number. The new, structurally simpler model improves on the Mellor-Yamada 1982 and Galperin et al. 1988 models; key feature includes enhanced mixing under stable conditions facilitating agreement with observational, experimental and high-resolution numerical data sets. The model predicts a planetary boundary layer height deeper than predicted by models with low critical Richardson numbers, as demonstrated in single column model runs of the GISS ModelE general circulation model.
The suitability of applying the Mellor and Yamada (1974, 1982) Level 2.5 second-order turbulence closure model to general circulation models is investigated by examining not only the scheme's simulation of fully (or nearly fully) developed turbulence, but also its simulation of rapidly growing or strongly decaying turbulence. The behavior of the model is presented over its entire domain of definition, with special consideration given to the pathologies of the model. The model is then modified for the case of growing turbulence to rectify some of its physical shortcomings for that case, and to remove the pathologies that prohibit its use in a general circulation model. The performance of the modified Level 2.5 model is compared to the performance of various other modified versions through the numerical simulation for a growing convective PBL. The results show that the modified Level 2.5 model is a viable candidate for the prediction of turbulence and the simulation of the PBL in general circulation models.
Atmospheric general circulation model (AGCM) cloud parameterizations generally include an assumption about the subgrid-scale probability distribution function (PDF) of total water and its vertical profile. In the present study, the Atmospheric Infrared Sounder (AIRS) monthly-mean cloud amount and relative humidity fields are used to compute a proxy for the second moment of an AGCM total water PDF called the RH01 diagnostic, which is the AIRS mean relative humidity for cloud fractions of 0.1 or less. The dependence of the second moment on horizontal grid resolution is analyzed using results from a high-resolution global model simulation.The AIRS-derived RH01 diagnostic is generally larger near the surface than aloft, indicating a narrower PDF near the surface, and varies with the type of underlying surface. High-resolution model results show that the vertical structure of profiles of the AGCM PDF second moment is unchanged as the grid resolution changes from 200 to 100 to 50 km, and that the second-moment profiles shift toward higher values with decreasing grid spacing.Several Goddard Earth Observing System, version 5 (GEOS-5), AGCM simulations were performed with several choices for the profile of the PDF second moment. The resulting cloud and relative humidity fields were shown to be quite sensitive to the prescribed profile, and the use of a profile based on the AIRS-derived proxy results in improvements relative to observational estimates. The AIRS-guided total water PDF profiles, including their dependence on underlying surface type and on horizontal resolution, have been implemented in the version of the GEOS-5 AGCM used for publicly released simulations.
Mixing coefficients due to clear-air turbulence are estimated from turbulence observations from aircraft, and from large-scale dissipation estimates from the large-scale energy budgets. Maximum coefficients occur near middle-latitude jet streams, and eddy viscosity there is of order of 10 sq m/sec; eddy conductivity is estimated to be about ten times smaller. These coefficients are introduced into the 12-layer general circulation model of the National Center of Atmospheric Research. They produce an apparently significant, though small reduction in maximum speed of the jet, and a reduction in eddy energy. Further, the stratospheric polar-night jet is produced at about the correct location with about the correct intensity.
The atmospheric general circulation models which are being used as components of climate models rely on their boundary layer parameterizations to produce realistic simulations of the surface turbulent fluxes of sensible heat. moisture. and momentum: of the boundary-layer depth over which these fluxes converge: of boundary layer cloudiness: and of the interactions of the boundary layer with the deep convective clouds that grow upwards from it. Two current atmospheric general circulation models are used as examples to show how these requirements are being addressed: these are version 3 of the Community Climate Model. which has been developed at the U.S. National Center for Atmospheric Research. and the Colorado State University atmospheric general circulation model. The formulations and results of both models are discussed. Finally, areas for future research are suggested.
A version of an atmospheric general circulation model was vectorized to run on a CDC STAR 100. The numerical model was coded and run in two different vector languages, CDC and LRLTRAN. A factor of 10 speed improvement over an IBM 360/95 was realized. Efficient use of the STAR machine required some redesigning of algorithms and logic. This precludes the application of vectorizing compilers on the original scalar code to achieve the same results. Vector languages permit a more natural and efficient formulation for such numerical codes.
General Circulation Model (GCM) studies of the atmospheric response to change boundary conditions are discussed. Results are reported on an extensive series of numerical studies based on the National Center for Atmospheric Research (NCAR) Community Climate Model (CCM) general circulation model. In these studies the authors determined the response to systematic changes in atmospheric CO2 ranging from 100 to 1000 ppm; to changes in the prescribed sea surface temperature (SST) in the Gulf of Mexico, such as occurred during the deglaciation phase of the last ice age; to changes in soil moisture over North America; and to changes in sea ice extent in the Southern Hemisphere. Study results show that the response of surface temperature and other variables is nearly logarithmic, with lower levels of CO2 implying greater sensitivity of the atmospheric state to changes in CO2. It was found that the surface temperature of the Gulf of Mexico exerts considerable control over the storm track and behavior of storm systems over the North Atlantic through its influence on evaporation and the source of latent heat. It was found that reductions in soil moisture can play a significant role in amplifying and maintaining North American drought, particularly when a negative soil moisture anomaly prevails late in the spring.
Earth system models (ESMs) and general circulation models (GCMs) are heavily used to provide inputs to sectoral impact and multisector dynamic models, which include representations of energy, water, land, economics, and their interactions. Therefore, representing the full range of model uncertainty, scenario uncertainty, and interannual variability that ensembles of these models capture is critical to the exploration of the future co-evolution of the integrated human–Earth system. The pre-eminent source of these ensembles has been the Coupled Model Intercomparison Project (CMIP). With more modeling centers participating in each new CMIP phase, the size of the model archive is rapidly increasing, which can be intractable for impact modelers to effectively utilize due to computational constraints and the challenges of analyzing large datasets. In this work, we present a method to select a subset of the latest phase, CMIP6, featuring models for use as inputs to a sectoral impact or multisector dynamics models, while prioritizing preservation of the range of model uncertainty, scenario uncertainty, and interannual variability in the full CMIP6 ensemble results. This method is intended to help impact modelers select climate information from the CMIP archive efficiently for use in downstream models that require global coverage of climate information. This is particularly critical for large-ensemble experiments of multisector dynamic models that may be varying additional features beyond climate inputs in a factorial design, thus putting constraints on the number of climate simulations that can be used. We focus on temperature and precipitation outputs of CMIP6 models, as these are two of the most used variables among impact models, and many other key input variables for impacts are at least correlated with one or both of temperature and precipitation (e.g., relative humidity). Besides preserving the multi-model ensemble variance characteristics, we prioritize selecting CMIP6 models in the subset that preserve the very likely distribution of equilibrium climate sensitivity values as assessed by the latest Intergovernmental Panel on Climate Change (IPCC) report. This approach could be applied to other output variables of climate models and, possibly when combined with emulators, offers a flexible framework for designing more efficient experiments on human-relevant climate impacts. It can also provide greater insight into the properties of existing CMIP6 models.
Martian electron density profiles provided by the Mars Global Surveyor (MGS) Radio Science (RS) experiment over the 95-200 km altitude range indicate what the height of the electron peak and the longitudinal structure of the peak height are sensitive indicators of the physical state of the Mars lower and upper atmospheres. The present analysis is carried out on five sets of occultation profiles, all at high solar zenith angles (SZA). Variations spanning 2 Martian years are investigated near aphelion conditions at high northern latitudes (64.7 - 77.6 N) making use of four of these data sets. A mean ionospheric peak height of 133.5 - 135 km is obtained near SZA = 78 - 82 deg.; a corresponding mean peak density of 7.3 - 8.5 x l0(exp 4)/ qu cm is also measured during solar moderate conditions at Mars. Strong wave number 2 - 3 oscillations in peak heights are consistently observed as a function of longitude over the 2 Martian years. These observed ionospheric features are remarkably similar during aphelion conditions 1 Martian year apart. This year-to-year repeatability in the thermosphere-ionosphere structure is consistent with that observed in multiyear aphelion temperature data of the Mars lower atmosphere. Coupled Mars general circulation model (MGCM) and Mars thermospheric general circulation model (MTGCM) codes are run for Mars aphelion conditions, yielding mean and longitude variable ionospheric peak heights that reasonably match RS observations. A tidal decomposition of MTGCM thermospheric densities shows that observed ionospheric wave number 3 features are linked to a non-migrating tidal mode with semidiurnal period (sigma = 2) and zonal wave number 1 (s = -1) characteristics. The height of this photochemically determined ionospheric peak should be monitored regularly.
We present a description of the ModelE2 version of the Goddard Institute for Space Studies (GISS) General Circulation Model (GCM) and the configurations used in the simulations performed for the Coupled Model Intercomparison Project Phase 5 (CMIP5). We use six variations related to the treatment of the atmospheric composition, the calculation of aerosol indirect effects, and ocean model component. Specifically, we test the difference between atmospheric models that have noninteractive composition, where radiatively important aerosols and ozone are prescribed from precomputed decadal averages, and interactive versions where atmospheric chemistry and aerosols are calculated given decadally varying emissions. The impact of the first aerosol indirect effect on clouds is either specified using a simple tuning, or parameterized using a cloud microphysics scheme. We also use two dynamic ocean components: the Russell and HYbrid Coordinate Ocean Model (HYCOM) which differ significantly in their basic formulations and grid. Results are presented for the climatological means over the satellite era (1980-2004) taken from transient simulations starting from the preindustrial (1850) driven by estimates of appropriate forcings over the 20th Century. Differences in base climate and variability related to the choice of ocean model are large, indicating an important structural uncertainty. The impact of interactive atmospheric composition on the climatology is relatively small except in regions such as the lower stratosphere, where ozone plays an important role, and the tropics, where aerosol changes affect the hydrological cycle and cloud cover. While key improvements over previous versions of the model are evident, these are not uniform across all metrics.
A general circulation model is presented which uses quadratically conservative, fourth order horizontal space differences on an unstaggered grid and second order vertical space differences with a forward-backward or a smooth leap frog time scheme to solve the primitive equations of motion. The dynamic equations for motion, finite difference equations, a discussion of the structure and flow chart of the program code, a program listing, and three relevent papers are given.
A general circulation model for numerical weather prediction, analysis, and forecast impact experiments, stratospheric analysis, and extended range weather forcasting was developed.