主流海表风场资料在舟山群岛近海的性能评估

再分析风场资料已广泛应用于我国舟山群岛海域可再生能源评估、海洋灾害预防决策以及港口运维和船舶运输等涉海发展领域,然而不同业务机构所提供的再分析数据在舟山近海的性能表现不一,严重阻碍了此类数据的有效应用。基于2018年全年单点浮标观测资料,综合评价了舟山群岛近海面(10 m)风场的长期变化趋势,并利用误差分析和风玫瑰图等统计工具对6种主流海表风场再分析资料,包括:ECMWF第五代全球大气再分析数据(the 5th generation ECMWF atmospheric reanalysis,ERA5)、NECP第二版全球高分辨率再分析数据(climate forecast system version 2,CFSv2)、美国宇航局物理海洋学分布存档中心的多卫星融合资料(cross-calibrated multi-platform,CCMP)、日本55年再分析数据(Japanese 55-year reanalysis,JRA-55)、第二版现代研究与应用回顾性分析数据(modern-era retrospective analysis for research and applications version 2,MERRA-2)和ECMWF哥白尼大气监测服务再分析数据(the Copernicus Atmosphere Monitoring Service,CAMS)在时间变化特征上进行了对比与评估。研究表明:在综合性能方面,ERA5对风场的再现能力最优,其次为JRA-55;在要素可信度方面,ERA5对风速的再现情况相对较优,而CFSv2的风向再现情况较好;风场产品在不同季节的模拟能力有所差异;不同风场产品在不同风速区间的重构能力也有所不同;在全年风向分布方面,各再分析资料都存在显著的东向偏差。研究结果为不同应用场景下风场资料的选取提供评估依据,并为进一步开发适用于舟山群岛近海的高精度长周期风场数据产品奠定基础。     

中国沿海HY-2A 校正微波辐射计水汽含量数据精度检验方法

为评价“海洋二号”卫星（HaiYang-2A, HY-2A）校正微波辐射计（Calibration Microwave Radiometer, CMR）近海水汽产品精度,以中国沿海全球导航卫星系统（Global Navigation Satellite System, GNSS）业务观测站数据和欧洲中期天气预报中心（European Centre for Medium-Range Weather Forecasts,ECMWF)发布的第五代大气再分析资料（ECMWF Reanalysis 5, ERA5）作为验证数据。首先对选取的GNSS业务观测站数据和CMR水汽含量数据进行时空匹配,两者的观测时间一致、空间范围取为100 km;然后利用精密单点定位方法反演GNSS业务观测站上空的大气可降水量（Precipitable Water Vapor, PWV）,同时对1 h分辨率的ERA5再分析资料内插计算,得到CMR水汽数据点处的ERA5 PWV;最后以GNSS PWV和ERA5 PWV为参考,分析2015年CMR水汽产品精度和偏差时空分布。结果表明,CMR水汽含量和GNSS PWV、ERA5 PWV之间的相关系数r均高于0.96,平均均方根误差分别为3.17 mm和1.58 mm,具有较高的精度;CMR水汽含量相对于GNSS PWV和ERA5 PWV的偏差不随季节变化而变化,但CMR水汽含量数据精度随纬度的增加而有所提高。     

ERA-Interim和ERA5再分析数据在江苏区域的适用性对比研究

基于江苏省73个国家级自动气象观测站和3个高空探测站的观测数据,对比研究了ERA-Interim和ERA5两套再分析数据在江苏区域的适用性。结果表明:无论是地面气象要素还是高空气象要素,ERA5再分析数据的适用性均优于ERA-Interim再分析数据,特别是ERA5再分析数据中的2 m气温及2 m相对湿度相对于ERA-Interim再分析数据体现出系统性的改进,相关系数的提升超过0.4。两套数据的小时降水量与观测之间均存在相对较大的误差,在实际使用中需谨慎。ERA5再分析数据与观测数据在一致性及误差方面存在较为显著的时空变化特征,该数据在夏季的适用性低于其他季节,在空间上夏半年体现了明显的南北分布特征,与之对应,ERA-Interim再分析数据没有显著的时空分布特征。采用新一代的ERA5再分析数据对江苏区域的天气/气候进行研究/预报具有潜在的优势。     

基于ERA5再分析数据的中国邻近海域极端波高特征分析

极端波浪对沿海地区基础设施有着深远的影响,了解它们的变化规律是进行海岸带风险分析和灾害预防的基础。文章基于欧洲中期天气预报中心(European Centre for Medium-Range Weather Forecasting,ECMWF)第五代再分析数据集(ECMWF reanalysis v5,ERA5),对中国邻近海域1979~2018年间极端波高展开时空特征分析,并统计了40 a厄尔尼诺系数,利用广义极值分布(generalized extreme value,GEV)函数探究了近40 a厄尔尼诺现象对我国海域极端波高的影响,结果显示:统计1979~2018年整个研究区域前2%极端波浪年均值在6~10 m间浮动,且整体趋势递增,在四季趋势变化中,春夏极端波高增长趋势较秋冬高,且波动明显,在年际极端波高变化趋势中有较大波动时大多伴随着厄尔尼诺或是拉尼娜现象的发生,结合厄尔尼诺现象对GEV分布中位置参数的影响分布图和极端波高年、季节际趋势变化分布图,厄尔尼诺现象影响大的地区极端波高大多呈现增长趋势,表明厄尔尼诺现象对极端波高有较高影响。     

HRCLDAS与ERA5-Land气温日统计数据在山东地区的对比评估

基于山东省2021年3月—2022年2月1 519个气象观测站2 m气温观测数据,对中国气象局高分辨率陆面数据同化系统(High Resolution China Meteorological Administration Land Data Assimilation System,HRCLDAS)和欧洲中期天气预报中心第五代陆面再分析数据集(ERA5-Land)逐小时2 m气温分析的日统计数据(平均气温、最高气温、最低气温)进行对比评估。结果显示:(1)HRCLDAS/ERA5-Land日统计平均气温、最高气温、最低气温的均方根误差分别为0.1/1.2 ℃、0.6/1.9 ℃、0.4/1.7 ℃,表明HRCLDAS具有更高的精度,且在不同地理区域、不同海拔高度的表现均优于ERA5-Land,大部地区的偏差(-0.5~0.5 ℃)远低于ERA5-Land(-2.0~2.0 ℃)。(2)两套数据对高温及寒潮过程的监测能力对比评估表明,HRCLDAS能够捕捉到大部分的高温以及寒潮过程,其与观测的高温日数及寒潮日数空间分布较为相似,但对影响范围存在一定的低估;ERA5-Land则只能监测到部分高温及寒潮过程,并对高温日数与寒潮日数存在严重的低估。     

菲律宾海盆深层温盐结构与环流变化特征研究

基于ECCO2 (Estimating the Circulation and Climate of the Ocean)、 GLORYS12V1 (Global Ocean Reanalysis and Simulations)、ORA-S5 (Ocean ReAnalysis)三种海洋再分析数据,对比研究了菲律宾海盆深层温盐及环流的季节和年际变化特征。结果表明:三种数据显示的海盆深层温盐季节变化特征基本一致,在3 000～4 000 m水深区域,海水呈春夏两季高温低盐而秋冬季低温高盐特性,4 000 m以下海水温盐季节变化很小;沿西边界,温度与内部有明显差异且季节变化幅度相对较大。沿西边界的输运季节变化特征表现为10月至次年4月输运向南,5—9月输运向北,并且在8月份达最大值;表明存在沿西边界的流动,即菲律宾海盆与南端西卡罗林海盆(West Caroline Basin)之间存在季节性水体交换。海盆深层海水温盐年际变化也十分显著,但不同数据显示的变化特征存在较大差异。EOF和相关分析显示,三种再分析数据的深层位温与ENSO均存在一定相关性,ECCO2的深层位温变化与ENSO的相关性最强。由于长期观测数据较少,再分析数据的结果难以验证,因此目前对年际变化特征的研究仍具有很大的不确定性。     

华南秋季大尺度大气水汽汇时空演变特征

用1958—2004年NCEP/NCAR再分析资料和中国160站月降水量资料分析了华南秋季大尺度大气水汽汇的时空变化特征。结果表明,华南中西部、东部地区是华南秋季水汽汇的2个主要变异中心区。华南中西部地区秋季水汽汇与该地区降水一样,以年际尺度变化为主;而该地区蒸发量的年代际变化比年际变化还稍显著。华南东部地区秋季降水、水汽汇和蒸发都存在明显的年际和年代际变化特征。如果我国南方上空出现向东北(向西南)的水汽通量距平,则会导致华南上空的水汽汇偏强(偏弱)。     

台风条件下ERA5再分析数据对中国近海适用性评估

2019年第9号台风“利奇马”是自1949年以来在浙江登陆的第三强台风,对中国近海沿岸地区造成了巨大的经济损失以及人员伤亡。为了研究ERA5再分析数据在渤黄东海的适用性,本文利用统计方法,针对台风“利奇马”期间位于渤海（QF104）、黄海（QF108）和东海（QF5003）的三套浮标观测数据对ECMWF最新推出的ERA5再分析数据进行评估分析。评估结果表明：1） ERA5再分析风速、风向、有效波高数据与浮标观测数据的匹配度较好,平均周期与海面温度数据次之;2）ERA5再分析平均波向数据相对平滑,而实测数据波动性大,两者之间的相关性有待进一步提高。     

ERA-Interim和ERA5再分析风资料在中国近海的适用性对比研究

全球再分析海面风资料在波浪模拟和风能资源评估等研究中发挥着重要作用，但风场资料种类繁多，且准确性在不同海域差异较大，使用时需要进行适用性分析。本文基于欧洲中期天气预报中心的ERA5和ERA-Interim再分析风场，利用多个站点的实测数据，分析了其在中国近海的适用性，并将再分析风场输入FVCOM-SWAVE波浪模型，对比了它们在常风天和台风天对波浪模拟的效果。结果表明：（1）常风天条件下ERA5和ERA-Interim资料在中国近海表现相似，风速较实测值略偏大，均能基本反映海表面风场变化和平均风速分布，吻合度指标在各站点均超过0.9;(2) ERA5对台风的模拟显著优于ERA-Interim，能较好模拟台风风速结构，对不同台风模拟精度差异大，整体上会低估台风风速；（3）风场质量是造成波浪模拟误差的主要原因之一，ERA5和ERA-Interim均能较好地模拟常海况下的波浪变化情况，而在台风浪的模拟中ERA5更优，“双台风”现象对风速和波浪的模拟准确度影响大。     

山东省冰雹频数月际变化的时空特征与成因探究

基于观测数据,统计分析了1999—2017年山东省冰雹频数月际变化的时空特征;基于大气再分析数据,通过对大尺度环流客观分型与物理量场的合成分析,探究了月际变化的成因。统计结果显示,山东省冰雹频数存在显著的月际变化,逐月冰雹落区具有明显的空间差异。成因探究结果表明:1)约88%的冰雹发生于低槽槽前,或槽后西北气流的环流形势下,此两种形势发生频率之和的月际变化与冰雹频数的月际变化呈现一致的起伏特征; 2)水汽与热力条件影响了冰雹发生频数,对4—6月和9—10月而言,充足的水汽与强不稳定对应了较高的冰雹频数,其中6月冰雹频数最高归因于水汽含量与K指数均达到最大,且距平正异常最小; 3) 4—9月冰雹落区的空间分布与低层水汽输送及辐合存在较好的对应关系,而10月对应关系不明显。     

An assessment of Arctic cloud water paths in atmospheric reanalyses

The role of Arctic clouds in the recent rapid Arctic warming has attracted much attention. However, Arctic cloud water paths(CWPs) from reanalysis datasets have not been well evaluated. This study evaluated the CWPs as well as LWPs(cloud liquid water paths) and IWPs(cloud ice water paths) from five reanalysis datasets(MERRA-2,MERRA, ERA-Interim, JRA-55, and ERA5) against the COSP(Cloud Feedback Model Intercomparison Project Observations Simulator Package) output for MODIS from the MERRA-2 CSP(COSP satellite simulator) collection(defined as M2 Modis in short). Averaged over 1980–2015 and over the Arctic region(north of 60°N), the mean CWPs of these five datasets range from 49.5 g/m~2(MERRA) to 82.7 g/m~2(ERA-Interim), much smaller than that from M2 Modis(140.0 g/m~2). However, the spatial distributions of CWPs, show similar patterns among these reanalyses, with relatively small values over Greenland and large values over the North Atlantic. Consistent with M2 Modis, these reanalyses show larger LWPs than IWPs, except for ERA-Interim. However, MERRA-2 and MERRA underestimate the ratio of IWPs to CWPs over the entire Arctic, while ERA-Interim and JRA-55 overestimate this ratio. ERA5 shows the best performance in terms of the ratio of IWPs to CWPs. All datasets exhibit larger CWPs and LWPs in summer than in winter. For M2 Modis, IWPs hold seasonal variation similar with LWPs over the land but opposite over the ocean. Following the Arctic warming, the trends in LWPs and IWPs during 1980～2015 show that LWPs increase and IWPs decrease across all datasets, although not statistically significant. Correlation analysis suggests that all datasets have similar interannual variability. The study further found that the inclusion of re-evaporation processes increases the humidity in the atmosphere over the land and that a more realistic liquid/ice phase can be obtained by independently treating the liquid and ice water contents.     

Interannual variation of atmospheric mass and the southern oscillation

Two reanalysis datasets, one generated by the Goddard Laboratory for Atmospheres for 1982–1993 and the other generated by the National Centers for Environmental Prediction for 1982–1995, are used to examine the relationship between the Southern Oscillation (SO) and the interannual variation of atmospheric mass. Both reanalyses show that atmospheric mass increases (decreases) during the positive (negative) SO phase. Atmospheric mass consists of dry air and moisture. Since dry mass is conserved, the interannual variation of atmospheric mass results from the variation of water vapor pressure. Thus, global atmospheric hydrological processes are analyzed to illustrate how the SO affects the interannual variation of atmospheric mass. During the positive (negative) SO phase, water vapor is converged (diverged) toward (out of) the central-eastern tropical Pacific [where sea surface temperatures (SSTs) are higher (lower) than normal] to maintain (suppress) cumulus convection in that area. An anomalous east-west Walker circulation straddling the Dateline is driven by the anomalous cumulus convection in this region to create positive (negative) surface pressure anomalies over the western tropical Pacific-Indian Ocean, which result in an increase (decrease) in atmospheric mass.     

“21世纪海上丝绸之路”风能资源时空变化评估

利用1979—2021 年的 ERA5 再分析资料,采用经验正交函数分解法、Mann-Kendall 趋势检验法等统计方法,对“21 世纪海上丝绸之路”相关海区的海表风场与风能密度的空间分布特征、季节变化特征以及长期变化趋势进行分析。结果表明：(1)研究海域风能密度在不同季节表现出很大的空间差异,夏季的阿拉伯海和孟加拉湾,冬季的中国南海,以及全年的热带南印度洋风能资源都极为丰富。(2)研究时段内,中国南海北部及附近海域、阿拉伯海西部、孟加拉湾西部以及热带西北印度洋风能密度等级整体较高。(3)研究海域的风能密度以年变化特征为主,其中中国南海风能密度的季节变幅最大且在春、秋两季表现出明显的转换特征。(4)在研究海区中,结合水深条件与风能密度时空变化特征的评估结果,可以重点关注台湾海峡、吕宋海峡、中南半岛东南沿海、阿拉伯海西部近岸海域及热带西北印度洋近岸大陆架海域风能资源的开发利用,加强其他海域风能资源的储备。此研究可为“21 世纪海上丝绸之路”风能资源的中长期开发规划提供依据。     

A global heat and freshwater forcing dataset for ocean models

A global dataset based on the ECMWF Re-Analyses (ERA) is presented that can be used as surface boundary conditions for ocean models with sea-ice components. The definition of these conditions is based on bulk formulae. To study the mean ocean circulation, a mean annual cycle on a daily basis was constructed from ERA for all relevant parameters including wind stress. Continental runoff is considered by using information about the catchment areas of the rivers and about the main drainage basins. The bulk formulae were extended by using sea ice concentration.To estimate meridional heat transports (MHT) and to avoid any drift in ocean model simulations, the heat and fresh water budgets have been closed by applying an inverse procedure to fine-tune the fluxes towards observed transports. To improve the MHTs on the Southern Hemisphere the winds and the short wave radiation at southern higher latitudes should be corrected. Furthermore, tests were performed concerning short wave radiation which was increased in the tropics and decreased in the subsidence zones.The heat and fresh water fluxes are assessed by using a scheme of Macdonald and Wunsch based on hydrographic sections. The net heat fluxes of ERA and of the forcing dataset are consistent with the heat flux divergences and convergences estimated by this scheme except for parts of the South Atlantic and the Indian Ocean sector of the Southern Ocean where none of these datasets is consistent with these estimates. In the subtropical South Indian Ocean the forcing dataset is consistent with these estimates while ERA are not. The flux components of ERA and the forcing dataset were compared to several observational datasets (SRB, SOC, HOAPS, GPCP, and CMAP). For each component, at least one of these datasets (especially HOAPS) supports the effects of the inverse procedure and the bulk formulae almost globally with some regional exceptions: short wave radiation in the tropical oceans and the subtropical North Atlantic, latent heat flux at higher latitudes, and precipitation in the northern North Atlantic.Comparisons to the NCAR/NCEP Re-Analyses (NRA) (versions 1 and 2) and the ECHAM model in place of ERA lead to similar results. In the North Atlantic the net heat fluxes of the model based datasets approach the hydrographic estimate with increasing resolution. Applied to any ocean/sea-ice model and compared to ERA, the forcing dataset would induce only a relative small net sea-surface buoyancy loss.A comparison of the forcing dataset to measurements made using one buoy deployed in the western Pacific warm pool and five buoys deployed in the subduction region of the Northeast Atlantic shows that at the site of the first buoy the net heat fluxes of the forcing dataset are in poorer agreement than those of ERA. At the sites of two subduction buoys both datasets show the same level of agreement within the error bars specified. At the sites of the three remaining subduction buoys the forcing dataset shows a marginal improvement on ERA.     

Water-vapor content in the atmosphere over European Russia during the summer 2010 fires

We study the water vapor (WV) content over European Russia (ER) during the period of forest and peatbog fires in July–August 2010 using total column water vapor observations from MODIS instruments (both Aqua and Terra platforms) as well as aerological data and NCEP/NCAR reanalysis. It is found that the spatial distribution of total column water vapor (TCWV) over ER in this period was anomalous, with the WV excess in the north of the territory and its deficit in the south of ER. The relationship between WV variations, atmospheric dynamics and the fire situation is analyzed. Along with the processes of the WV advection and evaporation we evaluate the contribution of pyrogenic emission of WV in spatial-temporal evolution of WV over ER during wildfires. The changes of water vapor at different heights in the troposphere and stratosphere are investigated. The results of a comparative analysis of WV contents during the periods of summertime atmospheric blockings in 1972 and 2010 are also presented. The near-infrared total-column precipitable water MODIS products (L3) are validated by upper-air radiosonde data.     

4个代表站点的珊瑚δ18O数据对热带印度洋气候变率的反演能力分析*

珊瑚中的δ¹⁸O含量可以高精度地记录热带海洋中的气候变率及气候变化信息, 从而弥补器测观测时间长度有限的缺陷, 为反演过去长期的气候状况提供了可能。热带印度洋的气候模态通过海气相互作用影响周边区域甚至全球的气候, 具有重要的研究意义。本文对比分析了1880—1999年间热带印度洋4个站点(肯尼亚、坦桑尼亚、塞舌尔、明打威)的珊瑚δ¹⁸O数据, 研究了其对热带印度洋气候的反演情况。结果显示, 珊瑚δ¹⁸O对海表温度(sea surface temperature, SST)的长期变化趋势反演良好, 珊瑚δ¹⁸O与局地SST在季节循环中保持一致, 且在低温季节对SST的变化更为敏感。此外, 4个站点的珊瑚δ¹⁸O对印度洋海盆模态(Indian Ocean Basin Mode, IOBM)和印度洋偶极子模态(Indian Ocean Dipple Mode, IODM)有一定的表征能力, 并且δ¹⁸O记录的热带印度洋年际变率受年代际变率调制。本文的研究结果表明, 珊瑚δ¹⁸O数据对于了解古气候时期热带印度洋气候变率具有重要的指示意义, 但同时也需要综合考虑不同站点的信息来完整地反映印度洋在历史时期的气候变化情况。     

Seasonal and interannual variability of oceanic carbon cycling in the western and central tropical-subtropical pacific: A physical-biogeochemical modeling study

The oceanic carbon cycle in the tropical-subtropical Pacific is strongly affected by various physical processes with different temporal and spatial scales, yet the mechanisms that regulate air-sea CO₂ flux are not fully understood due to the paucity of both measurement and modeling. Using a 3-D physical-biogeochemical model, we simulate the partial pressure of CO₂ in surface water (pCO_2sea) and air-sea CO₂ flux in the tropical and subtropical regions from 1990 to 2004. The model reproduces well the observed spatial differences in physical and biogeochemical processes, such as: (1) relatively higher sea surface temperature (SST), and lower dissolved inorganic carbon (DIC) and pCO_2sea in the western than in the central tropical-subtropical Pacific, and (2) predominantly seasonal and interannual variations in the subtropical and tropical Pacific, respectively. Our model results suggest a non-negligible contribution of the wind variability to that of the air-sea CO₂ flux in the central tropical Pacific, but the modeled contribution of 7% is much less than that from a previous modeling study (30%; McKinley et al., 2004). While DIC increases in the entire region SST increases in the subtropical and western tropical Pacific but decreases in the central tropical Pacific from 1990 to 2004. As a result, the interannual pCO_2sea variability is different in different regions. The pCO_2sea temporal variation is found to be primarily controlled by SST and DIC, although the role of salinity and total alkalinity, both of which also control pCO_2sea, need to be elucidated by long-term observations and eddy-permitting models for better estimation of the interannual variability of air-sea CO₂ flux.     

潜热通量异常对西北太平洋热带气旋活动影响的机理研究

使用1958—2001年的中国气象局台风年鉴资料和欧洲中心再分析资料,初步揭示了潜热通量异常对西北太平洋(含南海)热带气旋活动的可能影响机理。统计分析发现,西北太平洋(含南海)热带气旋年频数与北太平洋副热带(简称关键区)的潜热通量在过去40余年中均表现为显著的减弱趋势。诊断分析表明,关键区的潜热通量通过低层信风向西的水汽输送—整个西北太平洋热带气旋活动区水汽低层辐合上升而凝结—潜热释放这一链条改变西北太平洋(含南海)大气环境场条件(包括中层湿度场、高低层涡度场和高低层散度场),进而调制热带气旋活动。使用SAMIL模式进行关键区内潜热通量加倍和减半的敏感性试验,进一步证实了关键区潜热通量异常对西北太平洋热带气旋活动调制作用的可能机理,而对南海热带气旋活动影响较小。由此可推断,在过去40余年中,北太平洋副热带中部区潜热通量的减弱趋势可能是造成西北太平洋(含南海)热带气旋频数下降的主要原因之一。