Treatment of discontinuities inside Earth models:Effects on computed coseismic deformations

In this paper,we study how coseismic deformations calculated in 1066 Earth models are affected by how the models treat Earth discontinuities.From the results of applying models 1066A(continuous)and 1066B(discontinuous),we find that the difference in Love numbers of strike-slip and horizontal tensile sources are bigger than dip-slip and vertical tensile sources.Taken collectively,discontinuities have major effects on Green’s functions of four independent sources.For the near-field coseismic deformations of the 2013 Okhotsk earthquake(Mw 8.3),the overall differences between theoretical calculations in vertical displacement,geoid,and gravity changes caused by discontinuities are 10.52 percent,9.07 percent and 6.19 percent,with RMS errors of 0.624 mm,0.029 mm,and 0.063μGal,respectively.The difference in far-field displacements is small,compared with GPS data,and we can neglect this effect.For the shallow earthquake,2011 Tohoku-Oki earthquake(Mw 9.0),the differences in near-field displacements are 0.030 m(N-S),0.093 m(E-W),and 0.025 m(up-down)in our study area with the ARIA slip model,which gives results closer to GPS data than those from the USGS model.The difference in vertical displacements and gravity changes on the Earth’s surface caused by discontinuities are larger than 10 percent.The difference in the theoretical gravity changes at spatially fixed points truncated to degrees 60,as required by GRACE data,is 0.0016μGal and the discrepancy is 11 percent,with the theoretical spatial gravity changes from 1066B closer to observations than from 1066A.The results show that an Earth model with discontinuities in the medium has a large effect on the calculated coseismic deformations.

Spatial Inhomogeneity of Atmospheric CO_(2) Concentration and Its Uncertainty in CMIP6 Earth System Models

This paper provides a systematic evaluation of the ability of 12 Earth System Models(ESMs)participating in the Coupled Model Intercomparison Project Phase 6(CMIP6)to simulate the spatial inhomogeneity of the atmospheric carbon dioxide(CO_(2))concentration.The multi-model ensemble mean(MME)can reasonably simulate the increasing trend of CO_(2) concentration from 1850 to 2014,compared with the observation data from the Scripps CO_(2) Program and CMIP6 prescribed data,and improves upon the CMIP5 MME CO_(2) concentration(which is overestimated after 1950).The growth rate of CO_(2) concentration in the northern hemisphere(NH)is higher than that in the southern hemisphere(SH),with the highest growth rate in the mid-latitudes of the NH.The MME can also reasonably simulate the seasonal amplitude of CO_(2) concentration,which is larger in the NH than in the SH and grows in amplitude after the 1950s(especially in the NH).Although the results of the MME are reasonable,there is a large spread among ESMs,and the difference between the ESMs increases with time.The MME results show that regions with relatively large CO_(2) concentrations(such as northern Russia,eastern China,Southeast Asia,the eastern United States,northern South America,and southern Africa)have greater seasonal variability and also exhibit a larger inter-model spread.Compared with CMIP5,the CMIP6 MME simulates an average spatial distribution of CO_(2) concentration that is much closer to the site observations,but the CMIP6-inter-model spread is larger.The inter-model differences of the annual means and seasonal cycles of atmospheric CO_(2) concentration are both attributed to the differences in natural sources and sinks of CO_(2) between the simulations.

Global air-sea CO_(2) exchange fl ux since 1980s: results from CMIP6 Earth System Models

The ocean could profoundly modulate the ever-increasing atmospheric CO_(2) by air-sea CO_(2) exchange process,which is also able to cause signifi cant changes of physical and biogeochemical properties in return.In this study,we assessed the long-term average and spatial-temporal variability of global air-sea CO_(2) exchange fl ux(F CO_(2))since 1980s basing on the results of 18 Coupled Model Intercomparison Project Phase 6(CMIP6)Earth System Models(ESMs).Our fi ndings indicate that the CMIP6 ESMs simulated global CO_(2) sink in recent three decades ranges from 1.80 to 2.24 Pg C/a,which is coincidence with the results of cotemporaneous observations.What’s more,the CMIP6 ESMs consistently show that the global oceanic CO_(2) sink has gradually intensifi ed since 1980s as well as the observations.This study confi rms the simulated F CO_(2) could reach agreements with the observations in the aspect of primary climatological characteristics,however,the simulation skills of CIMP6 ESMs in diverse open-sea biomes are unevenness.None of the 18 CMIP6 ESMs could reproduce the observed F CO_(2) increasement in the central-eastern tropical Pacifi c and the midlatitude Southern Ocean.Defi ciencies of some CMIP6 ESMs in reproducing the atmospheric pressure systems of the Southern Hemisphere and the El Niño-Southern Oscillation(ENSO)mode of the tropical Pacifi c are probably the major causes.

Coupling of the Calculated Freezing and Thawing Front Parameterization in the Earth System Model CAS-ESM

The soil freezing and thawing process affects soil physical properties,such as heat conductivity,heat capacity,and hydraulic conductivity in frozen ground regions,and further affects the processes of soil energy,hydrology,and carbon and nitrogen cycles.In this study,the calculation of freezing and thawing front parameterization was implemented into the earth system model of the Chinese Academy of Sciences(CAS-ESM)and its land component,the Common Land Model(CoLM),to investigate the dynamic change of freezing and thawing fronts and their effects.Our results showed that the developed models could reproduce the soil freezing and thawing process and the dynamic change of freezing and thawing fronts.The regionally averaged value of active layer thickness in the permafrost regions was 1.92 m,and the regionally averaged trend value was 0.35 cm yr–1.The regionally averaged value of maximum freezing depth in the seasonally frozen ground regions was 2.15 m,and the regionally averaged trend value was–0.48 cm yr–1.The active layer thickness increased while the maximum freezing depth decreased year by year.These results contribute to a better understanding of the freezing and thawing cycle process.

Simulations of dissolved oxygen concentration in CMIP5 Earth system models

The climatologies of dissolved oxygen concentration in the ocean simulated by nine Earth system models(ESMs) from the historical emission driven experiment of CMIP5(Phase 5 of the Climate Model Intercomparison Project) are quantitatively evaluated by comparing the simulated oxygen to the WOA09 observation based on common statistical metrics. At the sea surface, distribution of dissolved oxygen is well simulated by all nine ESMs due to well-simulated sea surface temperature(SST), with both globally-averaged error and root mean square error(RMSE) close to zero, and both correlation coefficients and normalized standard deviation close to 1. However, the model performance differs from each other at the intermediate depth and deep ocean where important water masses exist. At the depth of 500 to 1 000 m where the oxygen minimum zones(OMZs) exist, all ESMs show a maximum of globally-averaged error and RMSE, and a minimum of the spatial correlation coefficient. In the ocean interior, the reason for model biases is complicated, and both the meridional overturning circulation(MOC) and the particulate organic carbon flux contribute to the biases of dissolved oxygen distribution. Analysis results show the physical bias contributes more. Simulation bias of important water masses such as North Atlantic Deep Water(NADW), Antarctic Bottom Water(AABW) and North Pacific Intermediate Water(NPIW) indicated by distributions of MOCs greatly affects the distributions of oxygen in north Atlantic, Southern Ocean and north Pacific, respectively.Although the model simulations of oxygen differ greatly from each other in the ocean interior, the multi-model mean shows a better agreement with the observation.

Spatial Variation in CO_(2) Concentration Improves the Simulated Surface Air Temperature Increase in the Northern Hemisphere

The increasing concentration of atmospheric CO_(2) since the Industrial Revolution has affected surface air temperature.However,the impact of the spatial distribution of atmospheric CO_(2) concentration on surface air temperature biases remains highly unclear.By incorporating the spatial distribution of satellite-derived atmospheric CO_(2) concentration in the Beijing Normal University Earth System Model,this study investigated the increase in surface air temperature since the Industrial Revolution in the Northern Hemisphere(NH) under historical conditions from 1976-2005.In comparison with the increase in surface temperature simulated using a uniform distribution of CO_(2),simulation with a nonuniform distribution of CO_(2)produced better agreement with the Climatic Research Unit(CRU) data in the NH under the historical condition relative to the baseline over the period 1901-30.Hemispheric June-July-August(JJA) surface air temperature increased by 1.28℃ ±0.29℃ in simulations with a uniform distribution of CO_(2),by 1.00℃±0.24℃ in simulations with a non-uniform distribution of CO_(2),and by 0.24℃ in the CRU data.The decrease in downward shortwave radiation in the non-uniform CO_(2) simulation was primarily attributable to reduced warming in Eurasia,combined with feedbacks resulting from increased leaf area index(LAI) and latent heat fluxes.These effects were more pronounced in the non-uniform CO_(2)simulation compared to the uniform CO_(2) simulation.Results indicate that consideration of the spatial distribution of CO_(2)concentration can reduce the overestimated increase in surface air temperature simulated by Earth system models.

A Brief Introduction to BNU-HESM1.0 and Its Earth Surface Temperature Simulations

Integrated assessment models and coupled earth system models both have their limitations in understanding the interactions between human activity and the physical earth system. In this paper, a new human--earth system model, BNU- HESM1.0, constructed by combining the economic and climate damage components of the Dynamic Integrated Model of Climate Change and Economy to the BNU-ESM model, is introduced. The ability of BNU-HESM1.0 in simulating the global CO2 concentration and surface temperature is also evaluated. We find that, compared to observation, BNU-HESM1.0 underestimates the global CO2 concentration and its rising trend during 1965-2005, due to the uncertainty in the economic components. However, the surface temperature simulated by BNU-HESM1.0 is much closer to observation, resulting from the overestimates of surface temperature by the original BNU-ESM model. The uncertainty of BNU-ESM falls within the range of present earth system uncertainty, so it is the economic and climate damage component of BNU-HESM 1.0 that needs to be improved through further study. However, the main purpose of this paper is to introduce a new approach to investigate the complex relationship between human activity and the earth system. It is hoped that it will inspire further ideas that prove valuable in guiding human activities appropriate for a sustainable future climate.

Global Climate Internal Variability in a 2000-year Control Simulation with Community Earth System Model(CESM)

Using the low-resolution （T31, equivalent to 3.75°× 3.75°） version of the Community Earth System Model （CESM） from the National Center for Atmospheric Research （NCAR）, a global climate simulation was carried out with fixed external forcing factors （1850 Common Era. （C.E.） conditions） for the past 2000 years. Based on the simulated results, spatio-temporal structures of surface air temperature, precipitation and internal variability, such as the E1 Nifio-Southem Oscillation （ENSO）, the Atlantic Multi-decadal Oscilla- tion （AMO）, the Pacific Decadal Oscillation （PDO）, and the North Atlantic Oscillation （NAO）, were compared with reanalysis datasets to evaluate the model performance. The results are as follows： 1） CESM showed a good performance in the long-term simulation and no significant climate drift over the past 2000 years; 2） climatological patterns of global and regional climate changes simulated by the CESM were reasonable compared with the reanalysis datasets; and 3） the CESM simulated internal natural variability of the climate system performs very well. The model not only reproduced the periodicity of ENSO, AMO and PDO events but also the 3-8 years vari- ability of the ENSO. The spatial distribution of the CESM-simulated NAO was also similar to the observed. However, because of weaker total irradiation and greenhouse gas concentration forcing in the simulation than the present, the model performances had some differences from the observations. Generally, the CESM showed a good performance in simulating the global climate and internal natu- ral variability of the climate system. This paves the way for other forced climate simulations for the past 2000 years by using the CESM.

Major Modes of Short-Term Climate Variability in the Newly Developed NUIST Earth System Model(NESM)

A coupled earth system model（ESM） has been developed at the Nanjing University of Information Science and Technology（NUIST） by using version 5.3 of the European Centre Hamburg Model（ECHAM）, version 3.4 of the Nucleus for European Modelling of the Ocean（NEMO）, and version 4.1 of the Los Alamos sea ice model（CICE）. The model is referred to as NUIST ESM1（NESM1）. Comprehensive and quantitative metrics are used to assess the model＇s major modes of climate variability most relevant to subseasonal-to-interannual climate prediction. The model＇s assessment is placed in a multi-model framework. The model yields a realistic annual mean and annual cycle of equatorial SST, and a reasonably realistic precipitation climatology, but has difficulty in capturing the spring–fall asymmetry and monsoon precipitation domains. The ENSO mode is reproduced well with respect to its spatial structure, power spectrum, phase locking to the annual cycle, and spatial structures of the central Pacific（CP）-ENSO and eastern Pacific（EP）-ENSO; however, the equatorial SST variability,biennial component of ENSO, and the amplitude of CP-ENSO are overestimated. The model captures realistic intraseasonal variability patterns, the vertical-zonal structures of the first two leading predictable modes of Madden–Julian Oscillation（MJO）, and its eastward propagation; but the simulated MJO speed is significantly slower than observed. Compared with the T42 version, the high resolution version（T159） demonstrates improved simulation with respect to the climatology, interannual variance, monsoon–ENSO lead–lag correlation, spatial structures of the leading mode of the Asian–Australian monsoon rainfall variability, and the eastward propagation of the MJO.

Effects of Earth's curvature and radial heterogeneity in dislocation studies:Case studies of the 2008 Wenchuan earthquake and the 2004 Sumatra earthquake

Recently, effects of Earth＇s curvature and radial heterogeneity on coseismic deformations are often investi- gated based on the 2004 Sumatra earthquake. However, such effects are strongly related to earthquake types. As a low dip angle event, the 2004 Sumatra earthquake is not a good seismic case for such a topic since the effects for moderate dip angle events are much bigger. In this study, the half-space and spherical dislocation theories are used, respectively, to calculate co- seismic displacements caused by the 2008 Wenchuan earthquake and the 2004 Sumatra earthquake. Effects of Earth＇s curva- ture and stratification are investigated through the discrepancies of results calculated using the two dislocation theories. Re- sults show that the effects of Earth＇s curvature and stratification for the 2008 Wenchuan earthquake are much larger than those for the 2004 Sumatra earthquake. Ignoring the effects will cause errors up to 100%-200% in far field displacements for a moderate dip angle event like the 2008 Wenchuan earthquake. Such great effects are much bigger than those conclusions of previous studies. Besides, comparison with observations verifies that spherical dislocation theories yield better results than half-space ones in far fields.

Regional earth system modeling：review and future directions

The authors review recent advances in the development of coupled Regional Earth System Models （RESMs）,a field that is still in its early stages.To date,coupled regional atmosphere-ocean-sea ice,atmosphere-aerosol and atmosphere-biosphere models have been developed,but they have been applied onlyto limited regional settings.Much more work is thus needed to assess their transferability to a wide range of settings.Future challenges in regional climate modeling are identified,including the development of fully coupled RESMs encompassing not only atmosphere,ocean,cryosphere,biosphere,chemosphere,but also the human component in a fully interactive way.

Quantifying the impacts of fire aerosols on global terrestrial ecosystem productivity with the fully-coupled Earth system model CESM

Fire is a global phenomenon and a major source of aerosols from the terrestrial biosphere to the atmosphere.Most previous studies quantified the effect of fire aerosols on climate and atmospheric circulation,or on the regional and site-scale terrestrial ecosystem productivity.So far,only one work has quantified their global impacts on terrestrial ecosystem productivity based on offline simulations,which,however,did not consider the impacts of aerosol–cloud interactions and aerosol–climate feedbacks.This study quantitatively assesses the influence of fire aerosols on the global annual gross primary productivity(GPP)of terrestrial ecosystems using simulations with the fully coupled global Earth system model CESM1.2.Results show that fire aerosols generally decrease GPP in vegetated areas,with a global total of−1.6 Pg C yr^−1,mainly because fire aerosols cool and dry the land surface and weaken the direct photosynthetically active radiation(PAR).The exception to this is the Amazon region,which is mainly due to a fire-aerosol-induced wetter land surface and increased diffuse PAR.This study emphasizes the importance of the influence of fire aerosols on climate in quantifying global-scale fire aerosols’impacts on terrestrial ecosystem productivity.

Earth System Model FGOALS-s2: Coupling a Dynamic Global Vegetation and Terrestrial Carbon Model with the Physical Climate System Model

Earth System Models （ESMs） are fundamental tools for understanding climate-carbon feedback. An ESM version of the Flexible Global Ocean-Atmosphere-Land System model （FGOALS） was recently developed within the IPCC AR5 Coupled Model Intercomparison Project Phase 5 （CMIP5） modeling framework, and we describe the development of this model through the coupling of a dynamic global vegetation and terrestrial carbon model with FGOALS-s2. The performance of the coupled model is evaluated as follows. The simulated global total terrestrial gross primary production （GPP） is 124.4 PgC yr-I and net pri- mary production （NPP） is 50.9 PgC yr-1. The entire terrestrial carbon pools contain about 2009.9 PgC, comprising 628.2 PgC and 1381.6 PgC in vegetation and soil pools, respectively. Spatially, in the tropics, the seasonal cycle of NPP and net ecosystem production （NEP） exhibits a dipole mode across the equator due to migration of the monsoon rainbelt, while the seasonal cycle is not so significant in Leaf Area Index （LAI）. In the subtropics, especially in the East Asian monsoon region, the seasonal cycle is obvious due to changes in temperature and precipitation from boreal winter to summer. Vegetation productivity in the northern mid-high latitudes is too low, possibly due to low soil moisture there. On the interannual timescale, the terrestrial ecosystem shows a strong response to ENSO. The model- simulated Nifio3.4 index and total terrestrial NEP are both characterized by a broad spectral peak in the range of 2-7 years. Further analysis indicates their correlation coefficient reaches -0.7 when NEP lags the Nifio3.4 index for about 1-2 months.

A GOCE only gravity model GOSG01S and the validation of GOCE related satellite gravity models

We compile the GOCE-only satellite model GOSG01S complete to spherical harmonic degree of 220 using Satellite Gravity Gradiometry （SGG） data and the Satellite-to-Satellite Tracking （SST） observations along the GOCE orbit based on applying a least-squares analysis. The diagonal components （Vxx, Vyy, Vzz） of the gravitational gradient tensor are used to form the system of observation equations with the band-pass ARMA filter. The point-wise acceleration observations （ax, ay, az） along the orbit are used to form the system of observation equations up to the maximum spherical harmonic degree/order 130. The analysis of spectral accuracy characteristics of the newly derived gravitational model GOSG01S and the existing models GOTIM04S, GODIR04S, GOSPW04S and JYY_GOCE02S based on their comparison with the ultrahigh degree model EIGEN-6C2 reveals a significant consistency at the spectral window approximately between 80 and 190 due to the same period SGG data used to compile these models. The GOCE related satellite gravity models GOSG01S, GOTIM05S, GODIR05S, GOTIM04S, GODIR04S, GOSPW04S, JYY_- GOCE02S, EIGEN-6C2 and EGM2008 are also validated by using GPS-leveling data in China and USA. According to the truncation at degree 200, the statistic results show that all GGMs have very similar differences at GPS-leveling points in USA, and all GOCE related gravity models have better performance than EGM2008 in China. This suggests that all these models provide much more information on the gravity field than EGM2008 in areas with low terrestrial gravity coverage. And STDs of height anomaly differences in China for the selected truncation degrees show that GOCE has improved the accuracy of the global models beyond degree 90 and the accuracies of the models improve from 24 cm to 16 cm. STDs of geoid height differences in USA show that GOSG01S model has best consistency comparing with GPSleveling data for the frequency band of the degree between 20 and 160.

A Novel Method for Locally Updating an Earth Model While Geosteering

Earth models are important tools for support of decision making processes for optimal exploitation of subsurface resources. For geosteering and other real-time processes where time is a major constraint, effective model management is decisive for optimal decision support. During drilling, subsurface information is received which should optimally be used to modify the 3D earth model. Today this model is typically not altered during the operation. We discuss the principles of a novel method that enables a populated earth model grid to be locally modified when the topology (connectivity) of the geological structure is locally altered. The method also allows local updates of the grid resolution. The modelled volume is split into closed regions by the structural model. Each region is individually discretized and obtains its own subgrid. Properties are stored in separate functions, e.g. for each layer, and transferred into each subgrid via a mapping. A local update of the geological structure implies that only subgrids in regions that are directly affected by the updated structure must be discarded and rebuilt, and the rest of the populated earth model grid is retained. Our focus is on decision support for optimal well placement while geosteering. The proposed method aims to manage multiple model realizations that are never fixed and always locally updated with the most recent measurements and interpretations in real-time, and where each realization is always kept at an optimal resolution.

STUDIES ON POLYMERIZATION MECHANISM OF CONJUGATED DIENE WITH η~3-ALLYL RARE EARTH MODEL COMPLEX

(η^-C_H_5)_2CeCl_5Mg_2(tmed)_2 was used as a catalyst to polymerize isoprene,Terminal ends of polyisoprene were identified as -CH_2-CH=CH_,which indicates that the insertion reaction of the monomer occurred between η^-allyl and cerium ion and gives a direct evidence for the η^-allyl mechanism of conjugated diene polymerzation with rare earth coordination catalysts.

Calculation of Interaction Parameters from Immiscible Phase Diagram of Alkali Metal or Alkali Earth Metal-Halide System by Means of Subregular Solution Model

In this paper, the interaction parameters in the subregular solution model, λ1 and λ2, are regarded as a linear function of temperature, T. Therefore, the molar excess Gibbs energy of A-B binary system may be reexpressed as follows:Gm^E=xAxB[(λ11+λ12T)+(λ21+λ22T)xB]The calculation of the model parameters, λ11, λ12, λ21and λ22, was carried out numerically from the phase diagrams for 11 alkali metal-alkali halide or alkali earth metal-halide systems. In addition, artificial neural network trained by known data has been used to predict the values of these model parameters. The predicted results are in good agreement with the .calculated ones. The applicability of the subregular solution model to the alkali metal-alkali halide or alkali earth metal-halide systems were tested by comparing the available experimental composition along the boundary of miscibility gap with the calculated ones which were obtained by using genetic algorithm. The good agreement between the calculated and experimental results across the entire liquidus is valid evidence in support of the model.

Earth Ecosystem Theory: Ⅲ Structure and Interface

By applying the segmentation and composition principles of ecosystem,the structure and interface of earth ecosystem are fully elaborated,and corresponding model diagrams are established to display its composition and interconversion. The research finds that earth ecosystem includes three sub-ecosystems： terrestrial ecosystem,marine ecosystem and atmospheric ecosystem. The impact sequence of sub-ecosystems on earth ecosystem is atmospheric ecosystem 〉 marine ecosystem 〉 terrestrial ecosystem. Based on the segmentation and composition principles of ecosystem and the application of earth ecosystem structure and interface,we give new explanations to Taoism and Taiji with the eight diagrams,which contributes to theoretical research of natural science.

Construction of Earth's gravitational field model from CHAMP,GRACE and GOCE data

The basic principle of spectral combination method is discussed,and the general expressions of the spectral weight and spectral combination of the united-processing of various types of gravimetric data are shown.What＇s more,based on degree error RMS of potential coefficients,the detailed expressions of spectral combination formulae and the corresponding spectral weights in the Earth＇s gravitational field model（EGM） determination using GOCE ＋ GRACE and CHAMP ＋ GRACE ＋ GOCE are derived.The fundamental situation that ulux-champ2013 s,tongji-GRACE01,go-cons-gcf-2-tim-r5 constructed respectively by CHAMP,GRACE,GOCE data and go-cons-gcf-2-dir-r5 constructed by syncretic processing of GRACE,GOCE and LAGEOS data are explained briefly,the degree error RMS,cumulative geoid height error and cumulative gravity anomaly error of these models are calculated.A syncretic model constructed from CHAMP,GRACE and GOCE data,which is expressed by champ ＋ grace ＋ goce,is obtained based on spectral combination method.Experimentation results show that the precision of CHAMP data model is the lowest in satellite-only models,so it is not needed in the determination of syncretic models.The GRACE data model can improve the GOCE data model in medium-long wavelength,so the overall precision of syncretic model can be improved.Consequently,as many types of gravimetric data as possible should be combined together in the data processing in order to strengthen the quality and reliability with widening scope and improve the precision and spatial resolution of the computational results.

SST effect on the pre-monsoon intraseasonal oscillation over the South China Sea based on atmospheric-coupled GCM comparison

The role of sea surface temperature(SST)variability in the pre-monsoonal(April to July)intraseasonal oscillation(ISO)over the South China Sea(SCS)is investigated using the Community Earth System Model Version 2(CESM2).An Atmospheric Model Intercomparison Project(AMIP)simulation forced by daily sea surface temperatures(SSTs)derived from a parallel coupled general circulation model(CGCM)run was compared with observations and the mother coupled simulation.In the coupled model,the SST warming leads the peak convection about 1/4 period as in observations.The paralell uncoupled model fails to simulate this phase relationship,implying the importance of air-sea coupling in reproducing realistic ISO.Due to the near-quadrature phase relationship between SST and precipitation ISOs during the ISO events,it is difficult to distinguish the active/passive role of SST from observations alone.Significant correlation in intraseasonal precipitation between the daily SST-forced AMIP and mother CGCM runs indicates that SST plays a role in driving the atmospheric ISO.