MASNUM-WAM海浪模式集合Kalman滤波同化研究-I.风场扰动对海浪模拟影响

模式集合样本的代表性和观测信息的可靠性是制约数据同化效果的重要因素,而前者对海浪模式同化的影响尤为显著。由于海浪模式对初始场的敏感性较弱,来自大气的风输入源函数是海浪的重要能量输入,如何合理地对风输入进行扰动,构造海浪的集合模式运行,是实现和改进海浪模式集合Kalman滤波同化的关键问题。为了实现海浪模式集合运行,本文提出了风场的三种集合扰动方案,分别为：纯随机数、随机场和时间滞后的风场扰动方法。本研究利用2014年1月ECMWF全球风场,基于这三种风场扰动方法开展了集合海浪模式的集合运行实验,并统计分析了海浪特征要素（有效波高）和二维波数谱对风场扰动的响应。结果表明,随机场集合扰动方案所构造的风场集合效果最佳,所得海浪模拟结果的集合样本发散度适中,能够较为合理地反映背景误差的统计特征,可用于进一步的集合Kalman滤波海浪数据同化实验。     

Seismic attenuation estimation using a complete ensemble empirical mode decomposition-based method

The frequency attenuation gradient method can provide important information for hydrocarbon detection. In this paper, a method using Complete Ensemble Empirical Mode Decomposition (CEEMD), Hilbert transform and the least-squares curve-fitting is proposed for seismic attenuation estimation as an effective frequency attenuation gradient estimation approach. We first use CEEMD to obtain the different Intrinsic Mode Functions (IMFs), which have a narrow band and can enhance the physical meaning of instantaneous attributes trace by trace. The time-frequency spectrum, which is computed using a Hilbert transform of each IMF, is represented as a spectrum with a single-peak that has narrow side lobes, which is conducive to frequency attenuation gradient estimation. Second, for each time sample, the frequency-amplitude spectrum of each IMF trace is extracted from the time-frequency spectrum to conduct the attenuation gradient computation. Then, the logarithm operation is performed for each IMF trace. Due to the very narrow bands of some IMFs in some seismic traces, a variable frequency window is adopted along the IMF trace according to the local data characteristics. Finally, the attenuation gradient for each IMF in a seismic trace can be computed using least-squares fitting. A different IMF reflects a seismic trace with a different spatiotemporal scale and can highlight different geologic and stratigraphic information. The correlation weighted average operation is used to highlight some useful details in seismic trace and obtains the attenuation gradient for each seismic trace. Field data examples demonstrate our method and its effectiveness. The proposed method can stably estimate the frequency attenuation gradient.     

概率密度分布法在“一带一路”地区极端月气温评估和预估中的应用

“一带一路”地区人口众多,气候类型复杂,亟待加强区域气候变化风险的认识。文中将该区分成10个区域,基于第五次耦合模式比较计划（CMIP5）中的31个全球模式模拟结果,应用概率密度分布（PDF）方法评估历史阶段（1986—2005年）各模式模拟暖月和冷月气温的能力,挑选并建立较优模式集合,用以预估21世纪中叶（2041—2060年）和21世纪末（2081—2100年）的极端月气温。结果表明,模式对观测中冷月气温距平PDF的模拟水平整体较暖月高。与多模式平均以及中位值相比,较优模式集合方法更适于极端暖/冷月气温的评估。在中等排放RCP4.5情景下,与低纬度地区相比,较优模式模拟中高纬地区未来极端暖/冷月气温的增温幅度的不确定性范围较大。21世纪中叶和21世纪末较优模式模拟的极端暖月气温在地中海增幅整体最大,东南亚增幅整体最小。对较优模式集合预估的极端冷月气温而言,无论是21世纪中叶还是世纪末,北欧增幅整体最大,东南亚增幅整体最小。     

基于CMIP6模式优化集合平均预估21世纪全球陆地生态系统总初级生产力变化

利用国际耦合模式比较计划第六阶段(CMIP6)中18个地球系统模式总初级生产力(GPP)模拟数据,基于传统的多模式集合平均(MME)和可靠集合平均方法(REA),在4个未来情景(SSP1-2.6,SSP2-4.5,SSP3-7.0和SSP5-8.5)下预估了21世纪全球陆地生态系统GPP的变化量,并分析了GPP变化的驱动因子。研究结果表明：在4个未来情景下,基于REA方法预估的全球陆地生态系统年GPP在未来时期(2068—2100年)比历史时期(1982—2014年)分别增长了(14.85±3.32)、(28.43±4.97)、(37.66±7.61)和(45.89±9.21)Pg C,其增量大小和不确定性都明显低于MME方法。在4个情景下,大气CO₂浓度增长对GPP变化的贡献最大,基于REA方法计算的贡献占比分别为140%、137%、115%和75%;除SSP5-8.5(24%)外,其他情景下升温均导致全球陆地生态系统GPP降低(-42%、-37%、-16%),部分抵消了CO₂施肥效应的正面贡献。温度的影响存在纬度差异：升温在低纬度地区对GPP有负向贡献,在中高纬度地区为正向贡献。降水和辐射变化对GPP变化的贡献相对较小。     

Comparison of different land-surface perturbation methods in short-range ensemble forecasts

为了比较不同陆面扰动方法对短时集合预报的影响,本研究设计了陆面模式扰动实验(LSMPE),初始土壤湿度扰动实验(ISMPE),陆面-大气耦合系数扰动实验(LCCPE)以及大气扰动对照实验(GEFSPE).结果表明,在三组陆面扰动实验中,LSMPE能代表最大的不确定性且误差最小;ISMPE的离散度要比LCCPE稍大,但是耦合系数的减弱能有效降低区域平均降水量.然而,所有陆面扰动实验产生的离散度均小于GEFSPE,这表明未来有必要将陆而扰动与大气扰动相结合,以使得模式能代表更多的不确定性.     

A CRPS-Based Spatial Technique for the Verification of Ensemble Precipitation Forecasts

Traditional precipitation skill scores are affected by the well-known“double penalty”problem caused by the slight spatial or temporal mismatches between forecasts and observations. The fuzzy (neighborhood) method has been proposed for deterministic simulations and shown some ability to solve this problem. The increasing resolution of ensemble forecasts of precipitation means that they now have similar problems as deterministic forecasts. We developed an ensemble precipitation verification skill score, i.e., the Spatial Continuous Ranked Probability Score (SCRPS), and used it to extend spatial verification from deterministic into ensemble forecasts. The SCRPS is a spatial technique based on the Continuous Ranked Probability Score (CRPS) and the fuzzy method. A fast binomial random variation generator was used to obtain random indexes based on the climatological mean observed frequency, which were then used in the reference score to calculate the skill score of the SCRPS. The verification results obtained using daily forecast products from the ECMWF ensemble forecasts and quantitative precipitation estimation products from the OPERA datasets during June-August 2018 shows that the spatial score is not affected by the number of ensemble forecast members and that a consistent assessment can be obtained. The score can reflect the performance of ensemble forecasts in modeling precipitation and thus can be widely used.     

CMIP5多模式集合对江苏省气候变化模拟评估及情景预估

利用中国气象局所属的2 400余个台站观测资料制作的分辨率为0.25°×0.25°数据集中的气温、降水量资料评估了CMIP5中17个模式对于1961—2004年江苏省气温和降水量空间分布特征的模拟能力,筛选出了5个对江苏省气候特征模拟较好的模式。之后基于5个优选模式集合平均的结果预估了3种典型浓度路径(Representative Concentration Pathways,RCPs)下江苏省2006—2100年的气温和降水量变化趋势。结果表明:(1)全球耦合气候模式对江苏省的气温和降水量空间分布特征具有一定的模拟能力,并且模式集合平均的气温和降水量与观测资料的空间相关系数分别为0.85和0.93;(2)在低浓度路径(RCP2.6)、中浓度路径(RCP4.5)和高浓度路径(RCP8.5)3种温室气体排放情景下,江苏省2006—2100年的地表温度均呈现明显的增温趋势,并且苏北的增温幅度要高于苏南;(3)3种温室气体排放情景下,江苏省未来百年降水量均呈现出北方增多南方减少的趋势;(4)未来百年江苏省降水量随气温变化的趋势并不稳定,RCP2.6和RCP4.5情景下降水量随气温的升高而增加,而RCP8.5情景下降水量随气温的增加而减少。     

The Relationship between Deterministic and Ensemble Mean Forecast Errors Revealed by Global and Local Attractor Radii

It has been demonstrated that ensemble mean forecasts, in the context of the sample mean, have higher forecasting skill than deterministic(or single) forecasts. However, few studies have focused on quantifying the relationship between their forecast errors, especially in individual prediction cases. Clarification of the characteristics of deterministic and ensemble mean forecasts from the perspective of attractors of dynamical systems has also rarely been involved. In this paper, two attractor statistics—namely, the global and local attractor radii(GAR and LAR, respectively)—are applied to reveal the relationship between deterministic and ensemble mean forecast errors. The practical forecast experiments are implemented in a perfect model scenario with the Lorenz96 model as the numerical results for verification. The sample mean errors of deterministic and ensemble mean forecasts can be expressed by GAR and LAR, respectively, and their ratio is found to approach2~(1/2) with lead time. Meanwhile, the LAR can provide the expected ratio of the ensemble mean and deterministic forecast errors in individual cases.     

Sub-seasonal to Seasonal Hindcasts of Stratospheric Sudden Warming by BCC_CSM1.1(m):A Comparison with ECMWF

This study focuses on model predictive skill with respect to stratospheric sudden warming(SSW) events by comparing the hindcast results of BCC_CSM1.1(m) with those of the ECMWF's model under the sub-seasonal to seasonal prediction project of the World Weather Research Program and World Climate Research Program. When the hindcasts are initiated less than two weeks before SSW onset, BCC_CSM and ECMWF show comparable predictive skill in terms of the temporal evolution of the stratospheric circumpolar westerlies and polar temperature up to 30 days after SSW onset. However, with earlier hindcast initialization, the predictive skill of BCC_CSM gradually decreases, and the reproduced maximum circulation anomalies in the hindcasts initiated four weeks before SSW onset replicate only 10% of the circulation anomaly intensities in observations. The earliest successful prediction of the breakdown of the stratospheric polar vortex accompanying SSW onset for BCC_CSM(ECMWF) is the hindcast initiated two(three) weeks earlier. The predictive skills of both models during SSW winters are always higher than that during non-SSW winters, in relation to the successfully captured tropospheric precursors and the associated upward propagation of planetary waves by the model initializations. To narrow the gap in SSW predictive skill between BCC_CSM and ECMWF, ensemble forecasts and error corrections are performed with BCC_CSM. The SSW predictive skill in the ensemble hindcasts and the error corrections are improved compared with the previous control forecasts.     

GRAPES全球集合预报系统模式扰动随机动能补偿方案初步探究

为了体现次网格尺度能量升尺度转换过程中存在的不确定性, 文中将随机动能补偿(Stochastic Kinetic Energy Backscatter, SKEB)方案应用于GRAPES(Global/Regional Assimilation and Prediction System)全球集合预报系统(GRAPES-GEPS), 以更好地表征模式误差并且增大集合离散度。使用的SKEB方案基于具有一定时、空相关特征的随机型以及由数值扩散导致的局地动能耗散率来构造随机流函数强迫。并根据流函数与水平风速旋转分量的关系, 将SKEB方案中的流函数强迫转化为适用于GRAPES全球模式的水平风速扰动。结果表明, SKEB方案的使用一方面能够提高GRAPES对大气动能谱的模拟能力; 另一方面能够改善GRAPES-GEPS的集合离散度与集合平均误差的关系, 增加了集合离散度, 并在一定程度上减小了集合平均误差, 尤其是在热带地区这种改进更为显著。而且该方案使得热带地区连续分级概率评分(CRPS评分)显著减小。就降水预报而言, 从Brier评分与相对作用特征面积(AROC, Area under the Relative Operating Characteristics)的结果来看, SKEB方案有助于改善中国地区小雨[0.1 mm, 10 mm)、中雨[10 mm, 25 mm)与大雨[25 mm, 50 mm)量级降水的概率预报技巧, 而对暴雨[50 mm, ∞)量级降水预报技巧影响很小(24 h降水量)。总体上, 模式扰动随机动能补偿方案提高了GRAPES-GEPS的概率预报技巧。