一个全球耦合模式的ENSO后报试验

A group of seasonal hindcast experiments are conducted using a coupled model known as the Flexible Global Ocean-Atmosphere-Land System Model-gamil1.11 (FGOALS-g1.11) developed at the State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics (LASG). Two steps are included in our El Niño-Southern Oscillation (ENSO) hindcast experiments. The first step is to integrate the coupled GCM with the Sea Surface Temperature (SST) strongly nudged towards the observation from 1971 to 2006. The second step is to remove the SST nudging term. We carried out a one-year hindcast by adopting the initial values from SST nudging experiments from the first step on January 1st, April 1st, July 1st, and October 1st from 1982 to 2005. In the SST nudging experiment, the model can reproduce the observed equatorial thermocline anomalies and zonal wind stress anomalies in the Pacific, which demonstrates that the SST nudging approach can provide realistic atmospheric and oceanic initial conditions for seasonal prediction experiments. The model also demonstrates a high Anomaly Correlation Coefficient (ACC) score for SST in most of the tropical Pacific, Atlantic Ocean, and some Indian Ocean regions with a 3-month lead. Compared with the persistence ACC score, this model shows much higher ACC scores for the Nino3.4 index for a 9-month lead.     

一个全球耦合模式的ENSO后报试验

A group of seasonal hindcast experiments are conducted using a coupled model known as the Flexible Global Ocean-Atmosphere-Land System Modelgamil1.11 (FGOALS-g1.11) developed at the State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics (LASG).Two steps are included in our ElNi o-Southern Oscillation (ENSO) hindcast experiments.The first step is to integrate the coupled GCM with the Sea Surface Temperature (SST) strongly nudged towards the observation from 1971 to 2006.The second step is to remove the SST nudging term.The authors carried out a one-year hindcast by adopting the initial values from SST nudging experiments from the first step on January 1st,April 1st,July 1st,and October 1st from 1982 to 2005.In the SST nudging experiment,the model can reproduce the observed equatorial thermocline anomalies and zonal wind stress anomalies in the Pacific,which demonstrates that the SST nudging approach can provide realistic atmospheric and oceanic initial conditions for seasonal prediction experiments.The model also demonstrates a high Anomaly Correlation Coefficient (ACC) score for SST in most of the tropical Pacific,Atlantic Ocean,and some Indian Ocean regions with a 3-month lead.Compared with the persistence ACC score,this model shows much higher ACC scores for the Ni o-3.4 index for a 9-month lead.     

IMPROVEMENT OF OCEAN DATA ASSIMILATION SYSTEM AND CLIMATE PREDICTION BY ASSIMILATING ARGO DATA

The Argo(Array for Real-time Geostrophic Oceanography) data from 1998 to 2003 were used in the Beijing Climate Center-Global Ocean Data Assimilation System(BCC-GODAS). The results show that the utilization of Argo global ocean data in BCC-GODAS brings about remarkable improvements in assimilation effects. The assimilated sea surface temperature(SST) of BCC-GODAS can well represent the climatological states of observational data. Comparison experiments based on a global coupled atmosphere-ocean general circulation model(AOCGM) were conducted for exploring the roles of ocean data assimilation system with or without Argo data in improving the climate predictability of rainfall in boreal summer. Firstly, the global ocean data assimilation system BCC-GODAS was used to obtain ocean assimilation data under the conditions with or without Argo data. Then, the global coupled atmosphere-ocean general circulation model(AOCGM) was utilized to do hindcast experiments with the two sets of the assimilation data as initial oceanic fields. The simulated results demonstrate that the seasonal predictability of rainfall in boreal summer, particularly in China, increases greatly when initial oceanic conditions with Argo data are utilized. The distribution of summer rainfall in China hindcast by the AOGCM under the condition when Argo data are used is more in accordance with observation than that when no Agro data are used. The area of positive correlation between hindcast and observation enlarges and the hindcast skill of rainfall over China in summer improves significantly when Argo data are used.     

中国业务动力季节预报的进展

利用动力模式开展季节到年际的短期气候预测 ,是目前国际上气候预测的发展方向。自 1996年以来 ,经过 8a多的研制和发展 ,国家气候中心已建立起第 1代动力气候模式预测业务系统 ,其中包括 1个全球大气 海洋耦合模式 (CGCM )、1个高分辨率东亚区域气候模式 (RegCM_NCC)和 5个简化的ENSO预测模式 (SAOMS) ,可用于季节—年际时间尺度的全球气候预测 ;全球海气耦合模式与区域气候模式嵌套 ,可以提供高分辨率的东亚区域气候模式制做季节预测。CGCM对 1982～ 2 0 0 0年夏季的历史回报试验表明 ,该模式对热带太平洋海表面温度和东亚区域的季节预测具有较好的预测能力。RegCM NCC的 5a模拟基本上能再现东亚地区主要雨带的季节进展。利用嵌套的区域气候模式RegCM NCC对 1991～ 2 0 0 0年的夏季回报表明 ,在预报主要季节雨带方面有一定技巧。 2 0 0 1～ 2 0 0 3年 ,CGCM和RegCM NCC的实时季节预报与观测相比基本合理。特别是 ,模式成功地预报了 2 0 0 3年梅雨季节长江和黄河之间比常年偏多的降水。SAOMS模式系统的回报试验表明 ,该系统对热带太平洋海表面温度距平有一定的预报能力 ,模式超前 6～ 12个月的回报与观测的相关系数明显高于持续预报。 1997～ 2 0 0 3年 ,SAOMS多模式集合实时预报与观测的相关系数达到     

基于弱耦合资料同化的冬季北大西洋涛动年际变率预测

北大西洋涛动作为冬季北大西洋地区大气环流的主模态之一,其年际变率对全球许多地区气候变率具有重要影响,但目前其预测技巧并不高。采用降维投影四维变分同化方法,在耦合模式中建立了基于全球大气资料的弱耦合资料同化系统,直接同化月平均再分析资料,并进行了年代际后报试验。结果表明,通过耦合资料同化的手段,可以显著提升耦合模式对冬季北大西洋涛动年际变率及其相关的欧洲北部、美国东部、欧亚大陆北部的冬季近地面温度年际变率的后报效果,相关系数均超过0.05显著水平t检验。该后报效果的改进主要与在耦合同化过程中通过耦合模式中自由发展的海-气相互作用将大气的观测信息储存在耦合模式的海洋分量中,改进了冬季北大西洋地区海表温度“三极”型分布的时空变率及其时间序列的后报效果有关。该研究强调了耦合模式初始状态的准确度对提升冬季北大西洋涛动年际变率的后报技巧具有重要作用。     

Global seasonal climate predictability in a two tiered forecast system: part I: boreal summer and fall seasons

This paper shows demonstrable improvement in the global seasonal climate predictability of boreal summer (at zero lead) and fall (at one season lead) seasonal mean precipitation and surface temperature from a two-tiered seasonal hindcast forced with forecasted SST relative to two other contemporary operational coupled ocean–atmosphere climate models. The results from an extensive set of seasonal hindcasts are analyzed to come to this conclusion. This improvement is attributed to: (1) The multi-model bias corrected SST used to force the atmospheric model. (2) The global atmospheric model which is run at a relatively high resolution of 50 km grid resolution compared to the two other coupled ocean–atmosphere models. (3) The physics of the atmospheric model, especially that related to the convective parameterization scheme. The results of the seasonal hindcast are analyzed for both deterministic and probabilistic skill. The probabilistic skill analysis shows that significant forecast skill can be harvested from these seasonal hindcasts relative to the deterministic skill analysis. The paper concludes that the coupled ocean–atmosphere seasonal hindcasts have reached a reasonable fidelity to exploit their SST anomaly forecasts to force such relatively higher resolution two tier prediction experiments to glean further boreal summer and fall seasonal prediction skill.     

Ensemble hindcasts of ENSO events over the past 120 years using a large number of ensembles

Based on an intermediate coupled model (ICM), a probabilistic ensemble prediction system (EPS) has been developed. The ensemble Kalman filter (EnKF) data assimilation approach is used for generating the initial ensemble conditions, and a linear, first-order Markov-Chain SST anomaly error model is embedded into the EPS to provide model-error perturbations. In this study, we perform ENSO retrospective forecasts over the 120 year period 1886–2005 using the EPS with 100 ensemble members and with initial conditi...     

Impact of different ocean reanalyses on decadal climate prediction

We study the impact of three ocean state estimates (GECCO, SODA, [ECMWF]-ORA-S3) on decadal predictability in one particular forecast system, the Earth system model from the Max Planck Institute for Meteorology in Hamburg. The forecast procedure follows two steps. First, anomalies of temperature and salinity of the observational estimates are assimilated into our coupled model. Second, the assimilation runs are then used to initialize 10-year-long hindcasts/forecasts starting from each year between 1960 and 2001. The impact of the individual ocean state estimates is evaluated both by the extent to which climate variations from the ocean state estimates are adopted by the forecast system (‘fidelity’) and by the prediction skill of the corresponding hindcast experiments. The evaluation focuses on North Atlantic (NA) sea surface temperature (SST), upper-level (0–700?m) NA ocean heat content (OHC) and the Atlantic meridional overturning circulation (MOC). Regarding fidelity, correlations between observations and the assimilation runs are generally high for NA SST and NA OHC, except for NA OHC in the GECCO assimilation. MOC variations experience strong modifications when GECCO and SODA are assimilated, much less so when assimilating ORA-S3. Regarding prediction skill, when initializing with ORA-S3 and SODA, correlations with observations are high for NA OHC and moderate for NA SST. Correlations in case of GECCO, on the other hand, are high for NA SST and moderate for NA OHC. Relatively high MOC correlations between hindcasts and respective assimilation run appear in the first five years in GECCO in the tropics and subtropics and in ORA-S3 north of 50N. Correlations are largely reduced when the MOC signals are detrended. The trends in the assimilation runs are to some extent artifacts of the assimilation procedure. Hence, our potential predictabilities of the MOC are optimistic estimates of the upper limits of predictability. However, the ORA-S3 reanalysis gives the best results for our forecast system as measured by both overall fidelity of the assimilation procedure and predictions of upper-level OHC in the North Atlantic.     

Impact of Subsurface Entrainment on ENSO Prediction: 1997-98 El Nino

Twenty-one-year hindcasts of sea surface temperature (SST) anomalies in the tropical Pacific were performed to validate the influence of ocean subsurface entrainment on SST prediction.A new hybrid coupled model was used that considered the entrainment of subsurface temperature anomalies into the sea surface.The results showed that predictions were improved significantly in the new coupled model.The predictive correlation skill increased by about 0.2 at a lead time of 9 months,and the root-mean-square (RMS) errors were decreased by nearly 0.2°C in general.A detailed analysis of the 1997-98 El Nio hindcast showed that the new model was able to predict the onset,peak (both time and amplitude),and decay of the 1997-98 strong El Nio event up to a lead time of one year,factors that are not represented well by many other forecast systems.This implies,in terms of prediction,that subsurface anomalies and their impact on the SST are one of the controlling factors in ENSO cycles.Improving the presentation of such effects in models would increase the forecast skill.     

BCC_CSM1.1对10年尺度全球及区域温度的预测研究

近期10～30年时间尺度的年代际预测是第五次耦合模式国际比较计划(CMIP5)重要内容之一。按照CMIP5试验要求, 国家气候中心利用气候系统模式BCC_CSM1.1完成并提交了年代际试验结果。本文评估了该模式年代际试验对10年尺度全球及区域地表温度的预测能力, 并通过与20世纪历史气候模拟试验的对比分析, 研究模式模拟对海洋初始观测状态的依赖程度。分析结果表明:(1)在有、无海洋初始化条件下, 模式均能模拟出1960～2005年间全球10年平均实测地表温度的变暖趋势, 但在有海洋初始化条件下, 可以明显减小BCC_CSM1.1模式模拟的全球升温趋势, 使得年代际试验比历史试验的结果更接近观测值。这一特点在观测资料相对丰富的南北纬50°以内地区更为显著。(2)在年代际试验预测前期, 通过Nudging方法, 利用SODA再分析海洋温度资料对模式进行初始化, 经过前期8～12月的协调后, 模式预测的第1年南北纬50°范围海洋、陆面的平均地表气温接近于观测值(CRUTEM3, HadSST2)。由于模式初值SODA再分析SST资料与HadSST2观测值存在明显的全球大洋系统暖偏差以及模式本身系统偏差的影响, 年代际试验模拟的地表气温在2～7年之内, 从观测SST状态逐渐恢复到模式系统本身状态。在同组Decadal试验中, 陆面和海洋恢复调整的时间长度几乎一致。(3) 从10年平均气候异常在区域尺度上的预报技巧来看, 有、无海洋初始同化对预测结果影响不大, 高预测技巧区主要分布在南半球印度洋中高纬度、热带西太平洋以及热带大西洋区域。(4)SST变化与下垫面热通量密切相关, 在热带和副热带海洋区域, 长波辐射和感热通量是影响10年时间尺度SST变化较大的物理量, 在中高纬度海洋, 洋面温度变化主要受潜热通量的影响相对较大。     

潜热加热纳近方法在地面降水资料同化中的应用

基于国家气象中心GRAPES_Meso高分辨率区域模式,针对中尺度数值预报模式中预报雨带形成滞后问题,研究了潜热加热纳近方法在地面降水资料同化中的应用,以期提高短时数值天气预报的水平。2013年6月20日—7月20日的初步试验结果表明:通过调整模式潜热加热廓线,可以改进初始场中温、湿、风等要素的合理分布,增加降水区的对流不稳定性;潜热加热纳近方法可以缩短模式的调整适应 (spin-up) 时间,改进短时降水预报的落区和强度,提高3 h,6 h,12 h的降水预报TS,ETS评分;与传统的冷潜热加热纳近的试验结果相比,改进的暖潜热加热纳近试验对降水落区和强度的预报更接近观测,但强降水中心范围略大。     

快速更新循环同化预报系统的汛期试验与分析

为了支持短时临近预报,利用新一代数值预报技术,结合高时空分辨的新一代探测所提供的观测信息,开发了基于快速更新循环同化的数值预报系统。基于GRAPES预报模式及其三维变分,开发了逐时循环同化,包括雷达、卫星、飞机、常规地面和探空等观测资料的同化模块,并采用nudging技术,引入雷达回波,订正模式的初始云水、雨水等信息,实现快速更新预报。多种测试和汛期连续试验表明,系统运行稳定可靠。逐时同化场合理,与实况基本一致。通过1个月的滚动预报综合分析和个例预报分析均表明,预报稳定有效,与观测分布基本一致,初步具备开展短时临近预报的能力。     

EXPERIMENTS ON ASSIMILATION OF INITIAL VALUES IN NUMERICAL PREDICTION OF A WARM-SECTOR PRECIPITATION IN SOUTH CHINA

In order to understand the impact of initial conditions upon prediction accuracy of short-term forecast and nowcast of precipitation in South China, four experiments i.e. a control, an assimilation of conventional sounding and surface data, testing with nudging rainwater data and the assimilation of radar-derived radial wind, are respectively conducted to simulate a case of warm-sector heavy rainfall that occurred over South China, by using the GRAPES_MESO model. The results show that (1) assimilating conventional surface and sounding observations helps improve the 24-h rainfall forecast in both the area and order of magnitude; (2) nudging rainwater contributes to a significant improvement of nowcast, and (3) the assimilation of radar-derived radial winds distinctly improves the 24-h rainfall forecast in both the area and order of magnitude. These results serve as significant technical reference for the study on short-term forecast and nowcast of precipitation over South China in the future.     

Evaluation of ocean data assimilation in CAS-ESM-C: Constraining the SST field

A weakly coupled assimilation system, in which SST observations are assimilated into a coupled climate model(CASESM-C) through an ensemble optimal interpolation scheme, was established. This system is a useful tool for historical climate simulation, showing substantial advantages, including maintaining the atmospheric feedback, and keeping the oceanic fields from drifting far away from the observation, among others. During the coupled model integration, the bias of both surface and subsurface oceanic fields in the analysis can be reduced compared to unassimilated fields. Based on 30 model years of output from the system, the climatology and interannual variability of the climate system were evaluated. The results showed that the system can reasonably reproduce the climatological global precipitation and SLP, but it still suffers from the double ITCZ problem. Besides, the ENSO footprint, which is revealed by ENSO-related surface air temperature, geopotential height and precipitation during El Ni ?no evolution, is basically reproduced by the system. The system can also simulate the observed SST–rainfall relationships well on both interannual and intraseasonal timescales in the western North Pacific region, in which atmospheric feedback is crucial for climate simulation.     

利用ARGO资料改进ENSO和我国夏季降水气候预测

全球海洋ARGO资料的获取为气候预测的研究提供了前所未有的海洋资料。该文首先利用ARGO资料改进次表层参数化方案后的Zebiak-Cane海洋模式, 与统计大气模式耦合, 建立了热带动力海洋-统计大气耦合模式。通过比较应用和未应用ARGO资料改进的海洋模式, 进行了耦合模式的长期回报试验。结果表明:ARGO资料的应用极大地改善了耦合模式对热带太平洋海面温度异常的预测能力, 提前3个月和6个月的回报结果都有很大的改进, 基本上回报出了Ni?o3.4区海面温度异常的演变特征, 对厄尔尼诺和拉尼娜都能够给出较准确的回报, 回报结果与观测之间的相关性在整个热带太平洋区域明显提高。该文还利用国家气候中心 (NCC) 全球海气耦合模式, 对我国夏季降水进行了个例和多年季度回报模拟试验, 探讨了包含和不包含ARGO观测资料的同化资料作为初始场对我国夏季降水预测的影响, 表明采用带有ARGO观测资料的海洋同化初始场, 回报的我国夏季降水分布形式与观测更一致, 回报结果与观测之间的正相关区域变大, 对我国夏季降水的回报水平比采用没有ARGO观测资料的海洋同化初始场时有明显提高。     

Dynamic-analogue correction of the decadal change of East Asian summer precipitation in the late 1990s

This paper systematically evaluates the deviations that appear in the hindcasts of the East Asian summer precipitation (EASP) decadal change in the late 1990s in two global coupled models (BCC_CGCM and BCC_CSM). The possible causes for the deviations between the model hindcasts and observations are analyzed. The results show that the hindcasts of EASP by BCC_CGCM and BCC_CSM deviate from observations, with the anomaly correlation coefficient (ACC) being -0.01 and -0.09 for the two models, respectively. The SST anomalies in North and West Pacific and the SST index values predicted by the two models also deviate from the observations, indicating that inconsistent SST fields may be the key factor leading to the deviation in the prediction of the EASP decadal shift. Thus, a dynamic-analogue scheme is proposed to correct the precipitation hindcasts by using SSTs, where SST and EASP are highly correlated, to select historical analogue cases. Cross validations show that the average ACC of the temporal-latitude distribution of the EASP between the corrected hindcasts and observations is 0.18 for BCC_CGCM and 0.02 for BCC_CSM; both are much higher than the uncorrected hindcasts. Applying the dynamic-analogue correction scheme in both models successfully improves prediction of the EASP decadal change in the late 1990s.     

Effects of Spectral Nudging on Oceanic States in a Coarse-Resolution Model

A 1-degree global model is used to investigate the skill of spectral nudging at coarse resolution by performing two numerical experiments, one with spectral nudging and the other without. In the spectral nudging experiment, the model temperature and salinity are nudged to an observed climatological monthly-mean field. The study compares the model mean state, as well as the interannual and decadal variability of oceanic quantities with observations, (e.g., sea surface height (SSH) and sea surface temperature (SST)). Spectral nudging is found to be effective in constraining model drift from the observed mean state of temperature and salinity in the global ocean, which has been reported in previous studies. The present study further shows that spectral nudging significantly improves the model skill of topostrophy (a measure of currents flowing along the topography) in water depth below 2000?m with no clear improvement elsewhere. Despite its known ability to damp oceanic variability at various time scales, spectral nudging can still represent the interannual and decadal variability of SSH and SST well, to a degree comparable to the other experiment.     

Impact of assimilation of conventional and satellite meteorological observations on the numerical simulation of a Bay of Bengal Tropical Cyclone of November 2008 near Tamilnadu using WRF model

The objective of this study is to examine the impact of assimilation of conventional and satellite data on the prediction of a severe cyclonic storm that formed in the Bay of Bengal during November 2008 with the four-dimensional data assimilation (FDDA) technique. The Weather Research and Forecasting (WRF ARW) model was used to study the structure, evolution, and intensification of the storm. Five sets of numerical simulations were performed using the WRF. In the first one, called Control run, the National Centers for Environmental Prediction (NCEP) Final Analysis (FNL) was used for the initial and boundary conditions. In the remaining experiments available observations were used to obtain an improved analysis and FDDA grid nudging was performed for a pre-forecast period of 24 h. The second simulation (FDDAALL) was performed with all the data of the Quick Scatterometer (QSCAT), Special Sensor Microwave Imager (SSM/I) winds, conventional surface, and upper air meteorological observations. QSCAT wind alone was used in the third simulation (FDDAQSCAT), the SSM/I wind alone in the fourth (FDDASSMI) and the conventional observations alone in the fifth (FDDAAWS). The FDDAALL with assimilation of all observations, produced sea level pressure pattern closely agreeing with the analysis. Examination of various parameters indicated that the Control run over predicted the intensity of the storm with large error in its track and landfall position. The assimilation experiment with QSCAT winds performed marginally better than the one with SSM/I winds due to better representation of surface wind vectors. The FDDAALL outperformed all the simulations for the intensity, movement, and rainfall associated with the storm. Results suggest that the combination of land-based surface, upper air observations along with satellite winds for assimilation produced better prediction than the assimilation with individual data sets.     

The use of successive dynamic initialization by nudging to simulate cyclogenesis during GALE IOP 1

Summary Synoptic/diagnostic case studies have increasingly come to rely on numerical simulations started from some initial state after which the model generated fields receive no further information from observed data. The purpose of this paper is to demonstrate the use of a dynamic data assimilation technique based on nudging to create a dynamically consistent high-resolution four-dimensional data set that can be used for synoptic diagnostic studies. The nudging technique is applied in the Goddard Mesoscale Atmospheric Simulation System (GMASS) using the 3-h radiosonde data collected during GALE IOP 1. A unique aspect of this application is nudging toward data analyses for which the areal coverage shifts with time. One of the two nudging simulations assimilates surface pressure in addition to the temperature, mixing ratio, and wind components. The nudging values are determined by linear interpolation between 3-h observation times. Assuming a linear variation of the assimilated value in time leads to estimates of the nudging coefficients which take into account the accuracy of the observations.Both nudging simulations are more accurate in terms of S hand root-mean-square error (RMSE) scores than a control sunulation without successive initialization. The nudging simulation with surface pressure is more accurate than the nudging simulation without surface pressure assimilation for this case. The simulation with surface pressure nudging captures the surface cyclogenesis and the associated strong rise-fall couplet in the 500 hPa height field. It also exhibits the strongest ageostrophic flow and exit region vertical circulation associated with a jet streak on the western side of the intensifying upper-level trough.The data sets made possible by the dynamic assimilation/ simulation cycles are dynamically consistent, have high spatial and temporal resolution and are ideally suited for diagnostic studies. Examples presented include the evolution of the ageostrophic flow associated with the exit region of an upper-level jet propagating toward the base of an intensifying trough with increasing cyclonic curvature of the flow. The nudging simulation with surface pressure provides the resolution and accuracy required to depict the rapid transformation (within a 12-h period) of the exit region ageostrophic flow from predominantly cross contour to along contour as the jet streak approaches the base of the trough.With 19 Figures     

TRMM卫星的闪电观测资料在中尺度数值模式中的Nudging同化应用研究

在以往研究的基础上,本文提出一种利用TRMM卫星闪电观测资料改善中尺度数值模式水物质初始场的Nudging同化技术,这是把物理初始化方法与Nudging同化方法相结合的一种技术.根据TRMM卫星闪电观测水平分布,对水汽和云凝结物进行Nudging同化,对积云对流参数化方案进行调整,改变模式初始水汽和云凝结物场的含量,以...