Uptake and Storage of Anthropogenic CO2 in the Pacific Ocean Estimated Using Two Modeling Approaches

A basin-wide ocean general circulation model(OGCM) of the Pacific Ocean is employed to estimate the uptake and storage of anthropogenic CO 2 using two different simulation approaches.The simulation(named BIO) makes use of a carbon model with biological processes and full thermodynamic equations to calculate surface water partial pressure of CO 2,whereas the other simulation(named PTB) makes use of a perturbation approach to calculate surface water partial pressure of anthropogenic CO 2.The results from the two simulations agree well with the estimates based on observation data in most important aspects of the vertical distribution as well as the total inventory of anthropogenic carbon.The storage of anthropogenic carbon from BIO is closer to the observation-based estimate than that from PTB.The Revelle factor in 1994 obtained in BIO is generally larger than that obtained in PTB in the whole Pacific,except for the subtropical South Pacific.This,to large extent,leads to the difference in the surface anthropogenic CO 2 concentration between the two runs.The relative difference in the annual uptake between the two runs is almost constant during the integration processes after 1850.This is probably not caused by dissolved inorganic carbon(DIC),but rather by a factor independent of time.In both runs,the rate of change in anthropogenic CO 2 fluxes with time is consistent with the rate of change in the growth rate of atmospheric partial pressure of CO 2.     

Nonstationary time series prediction by incorporating external forces

Almost all climate time series have some degree of nonstationarity due to external forces of the observed system. Therefore, these external forces should be taken into account when reconstructing the climate dy- namics. This paper presents a novel technique in predicting nonstationary time series. The main difference of this new technique from some previous methods is that it incorporates the driving forces in the pre- diction model. To appraise its effectiveness, three prediction experiments were carried out using the data generated from some known classical dynamical models and a climate model with multiple external forces. Experimental results indicate that this technique is able to improve the prediction skill effectively.     

GENERAL FORMS OF DYNAMIC EQUATIONS FOR ATMOSPHERE IN NUMERICAL MODELS WITH TOPOGRAPHY

The atmospheric dynamic equations have been transformed from the z-coordinate system into a generalized vertical coordinate system by using a sc-called DDD transformation method. Then the generalized system is assumed being pressure, sigma or incorporated pressure-sigma coordinate system and corresponding equations are obtained with the second-order accuracy. It is pointed out that the usual equations are only of the first-order accuracy when their space-differential terms are approximated by central finite differences. Therefore the usual forms of the equations may result in quite large errors on steep slopes of mountains included in a model.     

EXPERIMENTS OF A REDUCED GRID IN LASG/IAPWORLD OCEAN GENERAL CIRCULATION MODELS (OGCMs)

Due to the decrease in grid size associated with the convergence of meridians toward the poles in spherical coordinates, the time steps in many global climate models with finite-difference method are restricted to be unpleasantly small. To overcome the problem, a reduced grid is introduced to LASG/IAP world ocean general circulation models. The reduced grid is implemented successfully in the coarser resolutions version model L30T63 at first. Then, it is carried out in the improved version model LICOM with finer resolutions. In the experiment with model L30T63, under time step unchanged though, execution time per single model run is shortened significantly owing to the decrease of grid number and filtering execution in high latitudes. Results from additional experiments with L30T63 show that the time step of integration can be quadrupled at most in reduced grid with refinement ratio 3. In the experiment with model LICOM and with the model’s original time step unchanged, the model covered area is extended to the whole globe from its original case with the grid point of North Pole considered as an isolated island and the results of experiment are shown to be acceptable.     

Dependence of the agcm climatology on the method of prescribing surface boundary conditions and its climatological implication

By using IAP 9L AGCM, two sets of long-term climatological integration have been per-formed with the two different interpolation procedures for generating the daily surface boundary conditions. One interpolation procedure is the so-called “traditional” scheme, for which the daily surface boundary conditions are obtained by linearly interpolating between the observed monthly mean values, however the observed monthly means cannot be preserved after interpolation. The other one is the “new” scheme, for which the daily surface boundary conditions are obtained by linearly interpolating between the "artificial" monthly mean values which are based on, but are dif-ferent from the observed ones, after interpolating with this new scheme, not only the observed monthly mean values are preserved, the time series of the new generated daily values is also more consistent with the observation. Comparison of the model results shows that the differences of the globally or zonally averaged fields between these two integrations are quite small, and this is due to the compensating effect between the different regions. However, the differences of the two patterns (the global or regional geographical distributions), are quite significant, for example, the magni-tude of the difference in the JJA mean rainfall between these two integrations can exceed 2 mm/day over Asian monsoon regions, and the difference in DJF mean surface air temperature can also exceed 2oC over this region. The fact that the model climatology depends quite strongly on the method of prescribing the daily surface boundary conditions suggests that in order to validate the climate model or to predict the short-term climate anomalies, either the " new* interpolation scheme or the high frequency surface boundary conditions (e.g., daily or weekly data instead of the monthly data) should be introduced. Meanwhile, as for the coupled model, the daily coupling scheme between the different component cli?mate models (e.g., atmospheric and oceanic general circulation models) is preferred in order to partly eliminate the “climate drift” problem which may appear during the course of direct coupling.     

27.3-day and 13.6-day atmospheric tide and lunar forcing on atmospheric circulation

An analysis of time variations of the earth‘s length of day (LOD) versus atmospheric geopotential height fields and lunar phase is presented. A strong correlation is found between LOD and geopotential height from which a close relationship is inferred and found between atmospheric circulation and the lunar cycle around the earth. It is found that there is a 27.3-day and 13.6-day east-west oscillation in the atmospheric circulation following the lunar phase change. The lunar revolution around the earth strongly influences the atmospheric circulation. During each lunar cycle around the earth there is, on average, an alternating change of 6.8-day-decrease, 6.8-day-increase, 6.8-day-deerease and 6.8-day-increase in atmospheric zonal wind, atmospheric angular momentum and LOD. The dominant factor producing such an oscillation in atmospheric circulation is the periodic change of lunar declination during the lunar revolution around the earth. The 27.3- day and 13.6-day atmospheric oscillatory phenomenon is akin to a strong atmospheric tide, which is different from the weak atmospheric tides, diurnal and semidiurnal, previously documented in the literature. Also it is different from the tides in the ocean in accordance with their frequency and date of occurrence. Estimation shows that the 27.3-day lunar forcing produces a 1-2 m s^-1 change in atmospheric zonal wind. Therefore, it should be considered in models of atmospheric circulation and short and middle term weather forecasting. The physical mechanism and dynamic processes in lunar forcing on atmospheric circulation are discussed.     

An Impact Study of the Vertical Coordinate on a Non-hydrostatic Mesoscale High-Resolution Model

With the high-speed development of high-powered computer techniques, it is possible that a high-resolution and multi-scale unified numerical model is applied to the operational weather prediction. Some techniques about mesoscale non-hydrostatic numerical weather prediction are addressed, and the impact of the vertical coordinate system is one of them. Based on WRF (Weather Research and Forecasting) model, the influence of vertical coordinates on the non-hydrostatic mesoscale high-resolution model is compared. The results show that the error of various coordinates in lower levels is almost same when we use the geometry height (z) and the pressure (p) to set up a terrain-following coordinate; but the error of the height terrain-following coordinate in higher levels is smaller than that of the pressure terrain-following coordinate. The higher the resolution is, the bigger the error will be. The results of the high-resolution simulation exhibited that the trend of the difference in the two coordinates existed. In addition, the correlative coefficient and standard error are also analysed by the comparison between the forecast fields and the corresponding analysis fields.     

Comparison and Analysis of Two Different Calculation Schemes of SLEVE- Hybrid Coordinate in GRAPES_MESO Model

The calculation scheme of the smoothed-level and hybrid (SLEVE-hybrid for short) coordinates in numerical forecasting model is not limited to one. It is divided into the semi-analytical scheme and the finite differential scheme in terms of the various differential methods of the coordinate deformation variables. Comparing the dynamic equation and the long-time batch simulation results of the two schemes, the present study draws the following conclusions. The first- order finite difference accuracy of the coordinate deformation variables in the finite differential scheme is theoretically lower than that in the semi-analytical scheme. The larger the vertical gradient of the layer thickness is, the larger the relative errors of the finite differential scheme are. The long-time batch simulation test in the GRAPES model dynamic core demonstrates that the bias of the temperature and the geopotential height in the semi-analytical scheme is smaller under the default layering, while the simulation difference of the two schemes is greatly reduced when the layering is more uniform.     

The Development of a Nonhydrostatic Global Spectral Model

With the development of numerical weather prediction technology,the traditional global hydrostatic models used in many countries of the world for operational weather forecasting and numerical simulations of general circulation have become more and more unfit for high-impact weather prediction.To address this,it is important to invest in the development of global nonhydrostatic models.Few existing nonhydrostatic global models use consistently the grid finite difference scheme for the primitive equations of dynamical cores,which can subsequently degrade the accuracy of the calculations.A new nonhydrostatic global spectral model,which utilizes the Eulerian spectral method,is developed here from NCAR Community Atmosphere Model 3.0(CAM3.0).Using Janjic's hydrostatic/nonhydrostatic method,a global nonhydrostatic spectral method for the primitive equations has been formulated and developed.In order to retain the integrity of the nonhydrostatic equations,the atmospheric curvature correction and eccentricity correction are considered. In this paper,the Held-Suarez idealized test and an idealized baroclinic wave test are first carried out,which shows that the nonhydrostatic global spectral model has similar climate states to the results of many other global models for long-term idealized integration,as well as better simulation ability for short-term idealized integration.Then,a real case experiment is conducted using the new dynamical core with the full physical parameterizations of subgrid-scale physical processes.The 10-day numerical integration indicates a decrease in systematic error and a better simulation of zonal wind,temperature,and 500-hPa height.     

Properties of High-Order Finite Difference Schemes and Idealized Numerical Testing

Construction of high-order difference schemes based on Taylor series expansion has long been a hot topic in computational mathematics, while its application in comprehensive weather models is still very rare. Here, the properties of high-order finite difference schemes are studied based on idealized numerical testing, for the purpose of their application in the Global/Regional Assimilation and Prediction System(GRAPES) model. It is found that the pros and cons due to grid staggering choices diminish with higher-order schemes based on linearized analysis of the one-dimensional gravity wave equation. The improvement of higher-order difference schemes is still obvious for the mesh with smooth varied grid distance. The results of discontinuous square wave testing also exhibits the superiority of high-order schemes. For a model grid with severe non-uniformity and non-orthogonality, the advantage of high-order difference schemes is inapparent, as shown by the results of two-dimensional idealized advection tests under a terrain-following coordinate. In addition, the increase in computational expense caused by high-order schemes can be avoided by the precondition technique used in the GRAPES model. In general, a high-order finite difference scheme is a preferable choice for the tropical regional GRAPES model with a quasi-uniform and quasi-orthogonal grid mesh.     

A Numerical Study on Effects of Land-Surface Heterogeneity from ‘Combined Approach’ on Atmospheric Process Part II: Coupling-Model Simulations

Two land surface schemes, one the standard Biosphere / Atmosphere Transfer Scheme Version le (B0Z) and the other B1Z based on B0Z and heterogeneously-treated by ‘combined approach’ , were coupled to the meso-scale model MM4, respectively. Through the calculations of equations from the companion paper, parameters representing land surface heterogeneity and suitable for the coupling models were found out. Three cases were simulated for heavy rainfalls during 36 hours, and the sensitivity of short-term weath-er modeling to the land surface heterogeneity was tested. Through the analysis of the simulations of the three heavy rainfalls, it was demonstrated that B1Z, compared with B0Z, could more realistically reflect the features of the land surface heterogeneity, therefore could more realistically reproduce the circulation and precipitation amount in the heavy rainfall processes of the three cases. This shows that even short-term weather is sensitive to the land surface heterogeneity, which is more obvious with time passing, and whose influence is more pronounced in the lower layer and gradually extends to the middle and upper layer. Through the analysis of these simulations with B1Z, it is suggested that the bulk effect of smaller-scale fluxes (i.e., the momentum, water vapor and sensible heat fluxes) near the significantly-heterogeneous Land surface is to change the larger-scale (i.e., meso-scale) circulation, and then to influence the development of the low-level jets and precipitation. And also, the complexity of the land-atmosphere interaction was shown in these simulations.     

建立η坐标系下的假不压模式及其数值试验

滤除声波的大气运动方程中不包含声波,基于滤除声波方程建立的数值模式可以用较大的时间步长进行数值积分。Durran在1989年提出了一种新的滤除声波的方法,命名为“假不可压”方程,该方程考虑了温度扰动引起的密度变化,忽略了气压扰动引起的密度变化。本文根据Durran提出的假不可压理论,推导出了一组地形追随坐标下的通量形式的假不可压方程。该方程在形式上与WRF（Weather Research and Forecasting）模式中ARW（Advanced Research WRF）动力框架的控制方程非常接近。我们进一步将推导出的假不可压控制方程改写到了WRF模式中,建立了基于WRF模式框架的假不可压模式。用构建的假不可压模式和WRF模式做了两组对比试验:湿热泡对流试验和重力流试验。比较两种模式的模拟结果,可以看出假不可压模式的模拟结果与WRF模式的模拟结果非常接近,说明在WRF模式框架下建立的假不可压模式是合理可信的。     

Stabilization of thermohaline circulation by wind-driven and vertical diffusive salt transport

 The stability of the thermohaline circulation is investigated using an ocean general circulation model coupled to a simple atmospheric model. The atmospheric model is so developed that it represents the wind stress and the freshwater flux more realistically than existing energy balance models. The coupled model can reproduce the realistic deep ocean circulation without any flux adjustment. Effects of the wind stress and the vertical diffusion on the thermohaline circulation are studied by conducting various experiments with the coupled model. The Ekman upwelling between 60^∘N and 90^∘N brings up salt to the sea surface, while the compensation flow of the Ekman transport and the wind-driven gyre circulation between 30^∘N and 60^∘N carry salt horizontally to the high latitudes. By carrying out experiments where the wind stress is completely or partly removed, it is demonstrated that either of the vertical or the horizontal salt transport prevents the halocline formation at high latitudes and maintains the thermohaline circulation. For an experiment in which the vertical diffusivity is enhanced at high latitudes, it is shown that the vertical diffusion at high latitudes also prevents the halocline formation and stabilizes the thermohaline circulation. It is also shown that the value of the vertical diffusivity at high latitude affects the existence of the multiple equilibria of the thermohaline circulation. Received: 26 April 2000 / Accepted: 10 January 2001     

Global high resolution versus Limited Area Model climate change projections over Europe: quantifying confidence level from PRUDENCE results

Four high resolution atmospheric general circulation models (GCMs) have been integrated with the standard forcings of the PRUDENCE experiment: IPCC-SRES A2 radiative forcing and Hadley Centre sea surface temperature and sea-ice extent. The response over Europe, calculated as the difference between the 2071–2100 and the 1961–1990 means is compared with the same diagnostic obtained with nine Regional Climate Models (RCM) all driven by the Hadley Centre atmospheric GCM. The seasonal mean response for 2m temperature and precipitation is investigated. For temperature, GCMs and RCMs behave similarly, except that GCMs exhibit a larger spread. However, during summer, the spread of the RCMs—in particular in terms of precipitation—is larger than that of the GCMs. This indicates that the European summer climate is strongly controlled by parameterized physics and/or high-resolution processes. The temperature response is larger than the systematic error. The situation is different for precipitation. The model bias is twice as large as the climate response. The confidence in PRUDENCE results comes from the fact that the models have a similar response to the IPCC-SRES A2 forcing, whereas their systematic errors are more spread. In addition, GCM precipitation response is slightly but significantly different from that of the RCMs.     

Effect of tropospheric temperature change on the zonal mean circulation and SH winter extratropical cyclones

This study aims to understand the mechanisms which cause an overall reduction of SH extratropical cyclone activity with a slight increase in the high latitudes in a warmer climate simulated in general circulation models (GCMs) with increasing CO₂. For this purpose, we conducted idealized model experiments by forcing warm temperature anomalies to the areas where climate change models exhibit local maximum warming—the tropics in the upper troposphere and the polar regions in the lower troposphere—simultaneously and separately. The Melbourne University atmospheric GCM (R21) coupled with prescribed SST was utilized for the experiments. Our results demonstrate that the reduction of SH extratropical cyclone frequency and depth in the midlatitudes but the slight increase in the high latitudes suggested in climate change models result essentially from the tropical upper tropospheric warming. With this tropical warming, the enhanced static stability which decreases baroclinicity in the low and midlatitudes turns out to be a major contributor to the decrease of cyclone activity equatorward of 45°S whereas the increased meridional temperature gradient in the high latitudes seems an important mechanism for the increase of cyclone activity over 50°–60°S.     

Design and Numerical Simulation of an Arctic Ocean Circulation and Thermodynamic Sea－Ice Model

In this paper, the first version of a new Arctic Ocean circulation and thermodynamic sea-ice model is presented by the authors based on the framework of a twenty-layer World Oceanic general circulation model developed by Zhang et al. in 1994, The model’s domain covers the Arctic Ocean and Greenland-Norwegian Seas with the horizon-tal resolution of 200 km × 200 km on a stereographic projection plane. In vertical, the model uses the Eta-coordinate (Sigma modified to have quasi-horizontal coordinate surfaces) and has ten unevenly-spaced layers to cover the deep-est water column of 3000 m. Two 150-year integrations of coupling the ocean circulation model with the sea-ice model have been performed with seasonally cyclic surface boundary conditions. The only difference between the two experiments is in the model’s geography. Some preliminary analyses of the experimental results have been done fo-cused on the following aspects: (1) surface layer temperature, salinity and current; (2) the "Atlantic Layer"; (3) sea-ice cover and its seasonal variation. In comparison with the available observational data, these results are accept-able with reasonable accuracy.     

A study on combining global and regional climate model results for generating climate scenarios of temperature and precipitation for the Netherlands

Climate scenarios for the Netherlands are constructed by combining information from global and regional climate models employing a simplified, conceptual framework of three sources (levels) of uncertainty impacting on predictions of the local climate. In this framework, the first level of uncertainty is determined by the global radiation balance, resulting in a range of the projected changes in the global mean temperature. On the regional (1,000–5,000 km) scale, the response of the atmospheric circulation determines the second important level of uncertainty. The third level of uncertainty, acting mainly on a local scale of 10 (and less) to 1,000 km, is related to the small-scale processes, like for example those acting in atmospheric convection, clouds and atmospheric meso-scale circulations—processes that play an important role in extreme events which are highly relevant for society. Global climate models (GCMs) are the main tools to quantify the first two levels of uncertainty, while high resolution regional climate models (RCMs) are more suitable to quantify the third level. Along these lines, results of an ensemble of RCMs, driven by only two GCM boundaries and therefore spanning only a rather narrow range in future climate predictions, are rescaled to obtain a broader uncertainty range. The rescaling is done by first disentangling the climate change response in the RCM simulations into a part related to the circulation, and a residual part which is related to the global temperature rise. Second, these responses are rescaled using the range of the predictions of global temperature change and circulation change from five GCMs. These GCMs have been selected on their ability to simulate the present-day circulation, in particular over Europe. For the seasonal means, the rescaled RCM results obey the range in the GCM ensemble using a high and low emission scenario. Thus, the rescaled RCM results are consistent with the GCM results for the means, while adding information on the small scales and the extremes. The method can be interpreted as a combined statistical–dynamical downscaling approach, with the statistical relations based on regional model output.