Development and Validation of an Evaporation Duct Model.Part Ⅰ:Model Establishment and Sensitivity Experiments

Based on the Coupled Ocean-Atmospheric Response Experiment（COARE）bulk algorithm and the Naval Postgraduate School（NPS）model,a universal evaporation duct（UED）model that can flexibly accommodate the latest improvements in component（such as stability function,velocity roughness,and scalar roughness）schemes for different stratification and wind conditions,is proposed in this paper.With the UED model,the sensitivity of the model-derived evaporation duct height（EDH）to stability function（Ψ）,ocean wave effect under moderate to high wind speeds,and scalar roughness length parameterization,is investigated,and relative contributions of these factors are compared.The results show that the stability function is a key factor influencing the simulated EDH values.Under unstable conditions,the EDH values from stability functions of Fairall et al.（1996）and Hu and Zhang（1992）are generally higher than those from others;while under stable conditions,unreasonably high EDHs can be avoided by use of the stability functions of Hu and Zhang（1992）and Grachev et al.（2007）.Under moderate to high wind speeds,the increase in velocity roughness length z₀ due to consideration of the true ocean wave effect acts to reduce modeled EDH values;this trend is more pronounced under stable conditions.Although the scalar roughness length parameterization has a minor effect on the model-derived EDH,a positive correlation is found between the scalar roughness length z_0qand the model-derived EDH.     

Regional Nested Model Simulations of Present Day and 2 × CO2 Climate over the Central Plains of the U.S.

A nested regional climate model is used to generate a scenario of climate change over the MINK region (Missouri, Iowa, Nebraska, Kansas) due to doubling of carbon dioxide concentration (2 × CO₂) for use in agricultural impact assessment studies. Five-year long present day (control) and 2 × CO₂ simulations are completed at a horizontal grid point spacing of 50 km. Monthly and seasonal precipitation and surface air temperature over the MINK region are reproduced well by the model in the control run, except for an underestimation of both variables during the spring months. The performance of the nested model in the control run is greatly improved compared to a similar experiment performed with a previous version of the nested modeling system by Giorgi et al. (1994). The nested model generally improves the simulation of spatial precipitation patterns compared to the driving general circulation model (GCM), especially during the summer. Seasonal surface warming of 4 to 6 K and seasonal precipitation increases of 6 to 24% are simulated in 2 × CO₂ conditions. The control run temperature biases are smaller than the simulated changes in all seasons, while the precipitation biases are of the same order of magnitude as the simulated changes. Although the large scale patterns of change in the driving GCM and nested RegCM model are similar, significant differences between the models, and substantial spatial variability, occur within the MINK region.     

A Global Spectral Model and Test of Its Performance

A brief introduction is given of a global spectral model, its dynamical framework and diabatic physical processes involved. A number of real forecasts are shown to illustrate the forecasting capability of the model for various weather processes. It can even manage to predict some of those typical weather processes in summer which used to be difficult to forecasters.     

Optimal Parameter and Uncertainty Estimation of a Land Surface Model： Sensitivity to Parameter Ranges and Model Complexities

Most previous land-surface model calibration studies have defined global ranges for their parameters to search for optimal parameter sets. Little work has been conducted to study the impacts of realistic versus global ranges as well as model complexities on the calibration and uncertainty estimates. The primary purpose of this paper is to investigate these impacts by employing Bayesian Stochastic Inversion (BSI) to the Chameleon Surface Model (CHASM). The CHASM was designed to explore the general aspects of land-surface energy balance representation within a common modeling framework that can be run from a simple energy balance formulation to a complex mosaic type structure. The BSI is an uncertainty estimation technique based on Bayes theorem, importance sampling, and very fast simulated annealing.The model forcing data and surface flux data were collected at seven sites representing a wide range of climate and vegetation conditions. For each site, four experiments were performed with simple and complex CHASM formulations as well as realistic and global parameter ranges. Twenty eight experiments were conducted and 50 000 parameter sets were used for each run. The results show that the use of global and realistic ranges gives similar simulations for both modes for most sites, but the global ranges tend to produce some unreasonable optimal parameter values. Comparison of simple and complex modes shows that the simple mode has more parameters with unreasonable optimal values. Use of parameter ranges and model complexities have significant impacts on frequency distribution of parameters, marginal posterior probability density functions, and estimates of uncertainty of simulated sensible and latent heat fluxes.Comparison between model complexity and parameter ranges shows that the former has more significant impacts on parameter and uncertainty estimations.     

Development and Test of Hydrostatic Extraction Scheme in Spectral Model

The introduction of “hydrostatic extraction” scheme, or “standard stratification approximation”, into spectral model gained some advantages compared with commonly used schemes. However, computational instability may oc-cur for high vertical resolution versions if the stratification parameter C0 taken as a constant. In this paper, the pos-sible cause leading to the instability is discussed and an improved scheme presented where C0 is generalized to be a function of both height and latitudes. Hence the reference atmosphere gets closer to the real atmosphere and the tem-perature deviation field to be expanded becomes smoother everywhere. Test by real case forecasts shows good computational stability of the new scheme and better prediction performance than usual schemes of spectral model.     

Application of Radar Reflectivity Factor in Initializing Cloud-Resolving Mesoscale Model. Part Ⅰ: Retrieval of Microphysical Elements and Vertical Velocity

Assuming that cloud reaches static state in the warm microphysical processes, water vapor mixing ratio(qv), cloud water mixing ratio (qc), and vertical velocity (w) can be calculated from rain water mixing ratio (qr)- Through relation of Z-qr, qr can be retrieved by radar reflectivity factor (Z). Retrieval results indicate that the distributions of mixing ratios of vapor, cloud, rain, and vertical velocity are consistent with radar images, and the three-dimensional spatial structure of the convective cloud is presented. Treating q,v saturated at the echo area, the retrieved qr is about 0.1 g kg-1, qc is always less than 0.3 g kg-1, w is usually below 0.5 m s-1, and rain droplet terminal velocity (vr) is around 5.0 m s-1 in the place where radar reflectivity factor is about 25 dBz; in the place where echo is 45 dBz, the retrieved qr and qc are always about 3.0 g kg-1, w is greater than 5.0 m s-1, and vr is around 7.0 m s-1. In the vertical, the maximum updraft velocity is greater than 3.0 m s-1 at the height of around 5.0 kin, the maximum cloud water content is about 3.0 g kg-1 above 5 km and the maximum rain water content is about 3.0 g kg-1 below 6 kin. Due to the assumption that the cloud is in static state, there will be some errors in the retrieved variables within the clouds which axe rapidly growing or dying-out, and in such cases, more sophisticated radar data control technique will help to improve the retrieval results.     

Application of Radar Reflectivity Factor in Initializing Cloud-Resolving Mesoscale Model. Part Ⅰ： Retrieval of Microphysical Elements and Vertical Velocity

Assuming that cloud reaches static state in the warm microphysical processes, water vapor mixing ratio （qv）, cloud water mixing ratio （qc）, and vertical velocity （w） can be calculated from rain water mixing ratio （qr）. Through relation of Z-qr, qr can be retrieved by radar reflectivity factor （Z）. Retrieval results indicate that the distributions of mixing ratios of vapor, cloud, rain, and vertical velocity are consistent with radar images, and the three-dimensional spatial structure of the convective cloud is presented. Treating qv saturated at the echo area, the retrieved qr is about 0.1 g kg^-1, qc is always less than 0.3 g kg^-1, w is usually below 0.5 m s^-1, and rain droplet terminal velocity （vr） is around 5.0 m s^-1 in the place where radar reflectivity factor is about 25 dBz; in the place where echo is 45 dBz, the retrieved qr and qc are always about 3.0 g kg^-1, w is greater than 5.0 m s^-1, and vr is around 7.0 m s^-1. In the vertical, the maximum updraft velocity is greater than 3.0 m s^-1 at the height of around 5.0 km, the maximum cloud water content is about 3.0 g kg^-1 above 5 km and the maximum rain water content is about 3.0 g kg^-1 below 6 km. Due to the assumption that the cloud is in static state, there will be some errors in the retrieved variables within the clouds which are rapidly growing or dying-out, and in such cases, more sophisticated radar data control technique will help to improve the retrieval results.     

Improvements of a Dynamic Global Vegetation Model and Simulations of Carbon and Water at an Upland-Oak Forest

The interest in the development and improvement of dynamic global vegetation models （DGVMs）, which have the potential to simulate fluxes of carbon, water and nitrogen, along with changes in the vegetation dynamics, within an integrated system, has been increasing. In this paper, some numerical schemes and a higher resolution soil texture dataset were employed to improve the Sheffield Dynamic Global Vegetation Model （SDGVM）. Using eddy covariance-based measurements, we then tested the standard version of the SDGVM and the modified version of the SDGVM. Detailed observations of daily carbon and water fluxes made at the upland oak forest on the Walker Branch Watershed in Tennessee, USA offered a unique opportunity for these comparisons. The results revealed that the modified version of the SDGVM did a reasonable job of simulating the carbon and water flux and the variation of soil water content （SWC）. However, at the end of the growing season, it failed to simulate the effect of the limitations on the soil respiration dynamics and as a result underestimated this respiration. It was also noted that the modified version overestimated the increase in the SWC following summer rainfall, which was attributed to an inadequate representation of the ground water and thermal cycle.     

Implementation of a Surface Runoff Model with Horton and Dunne Mechanisms into the Regional Climate Model RegCM_NCC

A surface runoff parameterization scheme that dynamically represents both Horton and Dunne runoff generation mechanisms within a model grid cell together with a consideration of the subgrid-scaie soil heterogeneity, is implemented into the National Climate Center regional climate model （RegCM_NCC）. The effects of the modified surface runoff scheme on RegCMANCC performance are tested with an abnormal heavy rainfall process which occurred in summer 1998. Simulated results show that the model with the original surface runoff scheme （noted as CTL） basically captures the spatial pattern of precipitation, circulation and land surface variables, but generally overestimates rainfall compared to observations. The model with the new surface runoff scheme （noted as NRM） reasonably reproduces the distribution pattern of various variables and effectively diminishes the excessive precipitation in the CTL. The processes involved in the improvement of NRM-simulated rainfall may be as follows： with the new surface runoff scheme, simulated surface runoff is larger, soil moisture and evaporation （latent heat flux） are decreased, the available water into the atmosphere is decreased; correspondingly, the atmosphere is drier and rainfall is decreased through various processes. Therefore, the implementation of the new runoff scheme into the RegCMANCC has a significant effect on results at not only the land surface, but also the overlying atmosphere.     

Sensitivity of the Grid-point Atmospheric Model of IAP LASG （GAMILI.I.0） Climate Simulations to Cloud Droplet Effective Radius and Liquid Water Path

This paper documents a study to examine the sensitivity to cloud droplet effective radius and liquid water path and the alleviation the energy imbalance at the top of the atmosphere and at the surface in the latest version of the Grid-point Atmospheric Model of the State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics （LASG）, Institute of Atmospheric Physics （IAP） （GAMIL1.1.0）. Considerable negative biases in all flux components, and thus an energy imbalance, are found in GAMIL1.1.0. In order to alleviate the energy imbalance, two modifications, namely an increase in cloud droplet effective radius and a decrease in cloud liquid water path, have been made to the cloud properties used in GAMIL. With the increased cloud droplet effective radius, the single scattering albedo of clouds is reduced, and thus the reflection of solar radiation into space by clouds is reduced and the net solar radiation flux at the top of the atmosphere is increased. With the reduced cloud optical depth, the net surface shortwave radiation flux is increased, causing a net warming over the land surface. This results in an increase in both sensible and latent heat fluxes over the land regions, which is largely balanced by the increased terrestrial radiation fluxes. Consequently, the energy balance at the top of atmosphere and at the surface is achieved with energy flux components consistent with available satellite observations.     

Evaluation of Haney-Type Surface ThermalBoundary Conditions Using a CoupledAtmosphere and Ocean Model

A coupled atmosphere-ocean model developed at the Institute for Space Studies at NASA Goddard Space Flight Center (Russell et al., 1995) was used to verify the validity of Haney-type surface thermal boundary condition, which linearly connects net downward surface heat flux Q to air / sea temperature difference △T by a relaxation coefficient k. The model was initiated from the National Centers for Environmental Prediction (NCEP) atmospheric observations for 1 December 1977, and from the National Ocean Data Center (NODC) global climatological mean December temperature and salinity fields at 1° ×1° resolution. The time step is 7.5 minutes. We integrated the model for 450 days and obtained a complete model-generated global data set of daily mean downward net surface flux Q, surface air temperature TA,and sea surface temperature To. Then, we calculated the cross-correlation coefficients (CCC) between Q and △T. The ensemble mean CCC fields show (a) no correlation between Q and △T in the equatorial regions, and (b) evident correlation (CCC≥ 0.7) between Q and △T in the middle and high latitudes.Additionally, we did the variance analysis and found that when k= 120 W m-2K-1, the two standard deviations, σQ and σk△T, are quite close in the middle and high latitudes. These results agree quite well with a previous research (Chu et al., 1998) on analyzing the NCEP re-analyzed surface data, except that a smaller value of k (80 W m-2K-1) was found in the previous study.     

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 Hybrid Coupled Ocean–Atmosphere Model and Its Simulation of ENSO and Atmospheric Responses

A new hybrid coupled model(HCM) is presented in this study, which consists of an intermediate tropical Pacific Ocean model and a global atmospheric general circulation model. The ocean component is the intermediate ocean model(IOM)of the intermediate coupled model(ICM) used at the Institute of Oceanology, Chinese Academy of Sciences(IOCAS). The atmospheric component is ECHAM5, the fifth version of the Max Planck Institute for Meteorology atmospheric general circulation model. The HCM integrates its atmospheric and oceanic components by using an anomaly coupling strategy. A100-year simulation has been made with the HCM and its simulation skills are evaluated, including the interannual variability of SST over the tropical Pacific and the ENSO-related responses of the global atmosphere. The model shows irregular occurrence of ENSO events with a spectral range between two and five years. The amplitude and lifetime of ENSO events and the annual phase-locking of SST anomalies are also reproduced realistically. Despite the slightly stronger variance of SST anomalies over the central Pacific than observed in the HCM, the patterns of atmospheric anomalies related to ENSO,such as sea level pressure, temperature and precipitation, are in broad agreement with observations. Therefore, this model can not only simulate the ENSO variability, but also reproduce the global atmospheric variability associated with ENSO, thereby providing a useful modeling tool for ENSO studies. Further model applications of ENSO modulations by ocean–atmosphere processes, and of ENSO-related climate prediction, are also discussed.     

Effect of Counter-Gradient in the Computation of Turbulent Fluxes of Heat and Moisture in a Regional Model

The counter-gradient terms in the computations of turbulent fluxes of heat and moisture have been included in the PBL parameterization of a regional model for monsoon prediction. Results show that inclusion of counter-gradient terms has a marginal impact in the prediction of large scale monsoon circulation and rainfall rates.     

Development of a Model for Water and Heat Exchange Between the Atmosphere and a Water Body

A model for studying the heat and mass exchange between the atmosphere and a water body is developed, in which the phase change process of water freezing in winter and melting in summer and the function of the convective mixing process are taken into consideration. The model uses enthalpy rather than temperature as the predictive variable. It helps to set up governing equations more concisely, to deal with the phase change process more easily, and make the numerical scheme simpler. The model is verified by observed data from Lake Kinneret for a non-frozen lake in summer time, and Lake Lower Two Medicine for a frozen lake in winter time. Reasonably good agreements between the model simulations and observed data indicate that the model can serve as a component for a water body in a land surface model. In order to more efficiently apply the scheme in a climate system model, a sensitivity study of various division schemes with less layers in the vertical direction in the water body is conducted. The results of the study show that the division with around 10 vertical layers could produce a prediction accuracy that is comparable to the fine division with around 40 layers.     

Model Study on the Interannual Variability of Asian Winter Monsoon and Its Influence

ＭｏｄｅｌＳｔｕｄｙｏｎｔｈｅＩｎｔｅｒａｎｎｕａｌＶａｒｉａｂｉｌｉｔｙｏｆＡｓｉａｎＷｉｎｔｅｒＭｏｎｓｏｏｎａｎｄＩｔｓＩｎｆｌｕｅｎｃｅＪｉＬｉｒｅｎ（纪立人），ＳｕｎＳｈｕｑｉｎｇ（孙淑清）ＩｎｓｔｉｔｕｔｅｏｆＡｔｍｏｓｐｈｅｒｉｃＰｈｙ...     

A Coupled Dynamical-Radiational Model of Stratocumulus

A model dealing with interactions between the air and low stratiform clouds is presented based on the mixed-layer model Lilly (1968) pioneered and on Deardorffs three dimensional numerical model results. Its main new aspects lie in 1) consideration of the natures of both the atmosphere and cloud; 2) a new entrainment velocity scheme with few arbitrary assumptions; 3) transition from one-mixed layer to two-mixed layer model; and 4) parameterization of radiation and precipitation calculations.The model results for radiation, moisture, and heat turbulent fluxes turn out to be in good agreement with those calculated or observed by Kawa (1988), Nicholls (1984), and Schmets et al. (1981) in California, the North Sea, and the North Atlantic, respectively.Basically, this paper furnishes the theoretical basis for a model to address questions concerning the time-evolution of thermodynamical profiles both in cloud and out of cloud. The applications of this model wil be in a separate paper.