Assimilating the LAI Data to the VEGAS Model Using the Local Ensemble Transform Kalman Filter: An Observing System Simulation Experiment

Information on the spatial and temporal pat- terns of surface carbon flux is crucial to understanding of source/sink mechanisms and projection of future atmospheric CO2 concentrations and climate. This study presents the construction and implementation of a terrestrial carbon cycle data assimilation system based on a dynamic vegetation and terrestrial carbon model Vegetation-Global-Atmosphere-Soil （VEGAS） with an advanced assimilation algorithm, the local ensemble transform Kalman filter （LETKF, hereafter LETKF-VEGAS）. An observing system simulation experiment （OSSE） framework was designed to evaluate the reliability of this system, and numerical experiments conducted by the OSSE using leaf area index （LAI） observations suggest that the LETKF -VEGAS can improve the estimations of leaf carbon pool and LAI significantly, with reduced root mean square errors and increased correlation coefficients with true values, as compared to a control run without assimilation. Furthermore, the LETKF-VEGAS has the potential to provide more accurate estimations of the net primary productivity （NPP） and carbon flux to atmosphere （CFta）.     

一次典型高压型海雾过程中海上大气波导的数值模拟

2009年4月9—12日黄海海域发生了一次受高压系统影响的海雾过程。利用卫星观测与探空数据、WRF模式(Weather Research and Forecasting Model)对此次海雾过程及相伴的大气波导进行了观测分析与数值模拟。海雾与波导发展可分为3个阶段:(1)大气波导先于海雾存在于黄海海面;受高压下沉影响,黄海上空存在逆温层和较强的湿度梯度,表现为较强的贴海表面波导和非贴海表面波导。(2)海雾始于高压西部,并随高压系统逐渐东移减弱,向黄海北部扩展;辐射冷却虽然使雾顶附近逆温增强,但海雾的机械湍流使其顶部湿度梯度减小,雾顶附近对应弱悬空波导或波导消失。(3)高压系统影响使干空气下沉到雾区导致黄海海雾消散;雾顶附近逆温仍存在,同时湿度梯度增大,黄海上空逐渐变为非贴海表面波导。本研究结果表明:高压系统不仅极易为波导的发生提供有利条件,而且有利于海雾的生成,在海雾演变过程中主要是雾顶水汽梯度的变化导致了波导类型及强度的变化。     

Amodified discrete element model for sea ice dynamics

Considering the discontinuous characteristics of sea ice on various scales,a modified discrete element model（DEM） for sea ice dynamics is developed based on the granular material rheology.In this modified DEM,a soft sea ice particle element is introduced as a self-adjustive particle size function.Each ice particle can be treated as an assembly of ice floes,with its concentration and thickness changing to variable sizes under the conservation of mass.In this model,the contact forces among ice particles are calculated using a viscous-elastic-plastic model,while the maximum shear forces are described with the Mohr-Coulomb friction law.With this modified DEM,the ice flow dynamics is simulated under the drags of wind and current in a channel of various widths.The thicknesses,concentrations and velocities of ice particles are obtained,and then reasonable dynamic process is analyzed.The sea ice dynamic process is also simulated in a vortex wind field.Taking the influence of thermodynamics into account,this modified DEM will be improved in the future work.     

基于卫星高度计的北极海冰厚度变化研究

A modified algorithm taking into account the first year(FY) and multiyear(MY) ice densities is used to derive a sea ice thickness from freeboard measurements acquired by satellite altimetry ICESat(2003–2008). Estimates agree with various independent in situ measurements within 0.21 m. Both the fall and winter campaigns see a dramatic extent retreat of thicker MY ice that survives at least one summer melting season. There were strong seasonal and interannual variabilities with regard to the mean thickness. Seasonal increases of 0.53 m for FY the ice and 0.29 m for the MY ice between the autumn and the winter ICESat campaigns, roughly 4–5 month separation, were found. Interannually, the significant MY ice thickness declines over the consecutive four ICESat winter campaigns(2005–2008) leads to a pronounced thickness drop of 0.8 m in MY sea ice zones. No clear trend was identified from the averaged thickness of thinner, FY ice that emerges in autumn and winter and melts in summer. Uncertainty estimates for our calculated thickness, caused by the standard deviations of multiple input parameters including freeboard, ice density, snow density, snow depth, show large errors more than 0.5 m in thicker MY ice zones and relatively small standard deviations under 0.5 m elsewhere. Moreover, a sensitivity analysis is implemented to determine the separate impact on the thickness estimate in the dependence of an individual input variable as mentioned above. The results show systematic bias of the estimated ice thickness appears to be mainly caused by the variations of freeboard as well as the ice density whereas the snow density and depth brings about relatively insignificant errors.     

基于GPS反射信号技术的渤海海冰实验

传统的海冰检测手段都存在着受天气影响或实时检测成本太高的缺点,为了更实时有效的检测海冰,减轻海冰灾害对于我国的影响,并达到防灾减灾的目的,本文研究了GPS反射信号(GPS-R)技术检测海冰的可能性。通过对2013年2月下旬在中国天津渤海湾的岸基架设GPS信号接收机所接收的观测数据进行分析,使用不同特性的反射海域(几乎全为海水或是海冰)的数据进行比对分析,其结果显示极化比值(反射左旋信号与直射右旋信号的比值)能够明显反演海冰密集度。本文首次将GPS-R技术应用于国内自主设计的渤海海冰检测实验中。     

海冰对北极海冰边缘区大洋光学观测的影响评估

Diffuse attenuation coefficient(DAC) of sea water is an important parameter in ocean thermodynamics and biology, reflecting the absorption capability of sea water in different layers. In the Arctic Ocean, however, sea ice affects the radiance/irradiance measurements of upper ocean, which results in obvious errors in the DAC calculation. To better understand the impacts of sea ice on the ocean optics observations, a series of in situ experiments were carried out in the summer of 2009 in the southern Beaufort Sea. Observational results show that the profiles of spectral diffuse attenuation coefficients of seawater near ice cover within upper surface of 50 m were not contaminated by the sea ice with a solar zenith angle of 55°, relative azimuth angle of 110°≤φ≤115° and horizontal distance between the sensors and ice edge of greater than 25 m. Based on geometric optics theory, the impact of ice cover could be avoided by adjusting the relative solar azimuth angle in a particular distance between the instrument and ice. Under an overcast sky, ice cover being 25 m away from sensors did not affect the profiles of spectral DACs within the upper 50 m either. Moreover, reliable spectral DACs of seawater could be obtained with sensors completely covered by sea ice.     

长江口北支异常强盐水入侵观测与分析

长江口北支由于径流分流比很小,盐水入侵较强,特别是枯季大潮期盐水甚至倒灌进入南支,影响上海市和江苏省的水源地水质。为了进一步研究北支盐水入侵的规律及影响因素,2014年1月1-9日在北支进行了大小潮同步水文观测。本次观测到了北支异常强盐水入侵:小潮期发生强盐水入侵,且强度大于大潮期。小潮期强盐水入侵导致中下游河段被高浓度盐水控制,盐度从下游B01到中游B02几乎没有变化,且盐度的涨落潮变化几乎消失。分析认为,小潮观测期间长江入河口流量较小、河口潮差不是太小,加上强偏北风,三者的叠加是导致强盐水入侵的主要原因。     

船舶海平面气压观测资料质量控制

随着数值预报模式不断发展,各类观测资料在数值模式中发挥着越来越重要的作用,船舶观测资料是海上3类观测系统 (卫星、飞机、船舶) 资料之一。为了保障船舶观测资料质量并有效地应用于数值模式,该文依据船舶观测资料的时空分布特征,以及2011年1月和7月T639分析场与观测资料对比分析结果,建立了一套船舶海平面气压资料质量控制方案,包括要素极值范围检查、缺测和冗余资料剔除、背景场一致性检查、测站黑名单建立等,并将该方案应用于2011年2月、6月、8月的观测资料。结果表明:船舶观测资料在时间上不连续且空间分布不均匀,会影响时间一致性检查和空间一致性检查质量控制效果;船舶海平面气压观测资料在所有观测要素中资料量最大,但其缺测和冗余资料量约占50%;黑名单资料的质量控制方案能够有效识别和剔除黑名单资料,且有利于对各测站的检查与维修工作;由于五大湖和大奴湖地区地形高度的影响,在背景场一致性检查过程中需对这些区域的背景场资料进行订正。     

探空观测黑名单检查在变分同化系统中的应用

针对探空观测资料使用造成的某些区域GRAPES分析场存在虚假的高、低压系统问题,该文通过对比全球探空资料的位势高度观测与NCEP分析场,统计站点中观测质量较差的时次出现频数,确定探空位势高度观测黑名单。研究表明:500 hPa在印度地区、北大西洋和南极洲附近的探空位势高度观测与NCEP分析场的均方根误差在30 gpm以上的站点较多,且位势高度观测不可靠观测比率为20%以上的站点主要集中这些区域,以上观测站均列入黑名单。文中在GRAPES全球三维变分分析场的质量控制中加入探空位势高度观测黑名单检查,通过6 h分析预报循环试验表明:探空位势高度观测黑名单检查能有效提高分析场质量,GRAPES位势高度分析场在南极洲附近和印度地区有所改善。     

雷电定位系统与人工观测雷暴日数统计比较

为了利用雷电定位系统 (lightning location system,LLS) 资料统计人工观测雷暴日数,采用湖北省2007—2012年LLS监测资料,选取25个气象站为圆心,统计其不同监测半径 (r) 圆区域内LLS监测的雷电日数,并与人工观测雷暴日数进行比较。结果表明:r≤7 km时,LLS监测平均年雷电日数小于人工观测平均年雷暴日数;r≥8 km时, LLS监测平均年雷电日数大于人工观测平均年雷暴日数;r=22 km圆区域内年平均雷电日数可替代最大年雷暴日数。根据r=7 km,r=8 km圆区域内LLS监测的年雷电日数、年平均地闪密度资料,分别采用直接替代法、地闪密度法和该文提出的二元法计算年雷暴日数,结果显示:二元法效果最好。二元法计算的2007—2012年25个站平均年雷暴日数与人工观测相等,平均差异为7.4%;二元法计算的2013年年雷暴日数与人工观测相差0.8 d,平均差异为12.3%。