北京地区城市环境对云和降水影响的个例数值模拟研究

本文利用中尺度可分辨云模式(WRF,Weather Research and Forecasting)对2011年8月14日北京地区一次强降水过程进行了数值模拟和敏感性试验,研究了城市环境(包括城市气溶胶和城市地表)对北京地区云和降水的影响.研究结果表明:城市气溶胶污染增强和城市地表使得北京地区(城区和周边)降水量减少,对城区平均累积降水量的影响作用分别为-38.92%和-3.4%.降水系统向北京主城区移动过程中,城市气溶胶在城市环境影响降水过程中的作用逐渐减小为85.13%,城市地表的作用增加到14.87%.城市污染气溶胶增强,使得非降水性粒子增多,而降水性粒子减少,这不利于对流的发展增强,使得水汽的垂直输送减弱,导致区域降水量减小.城市地表对强对流的发展也表现为减弱作用,这使得水汽的垂直输送减小,区域降水量减小.     

中尺度模式积云参数化方案的改进及其在暴雨模拟中的应用

本文针对中国暴雨发生发展天气特征,改进和发展了一种适合于描述东亚暴雨的中尺度积云参数化方案.首先基于近年来(1990—2010)江淮流域汛期降水合成分析的基础上,诊断出组织化对流降水环境的动力参数;其次利用该动力参数作为动力控制条件,改进了Kain-Fritsch Eta中尺度积云对流参数化方案;最后利用改进的中尺度积云参数化方案对梅雨期暴雨、华南前汛期暴雨过程进行了数值模拟,结果表明:改进后的中尺度积云参数化方案对上述两次暴雨过程的落区及强度的模拟,均有明显改进.     

北京精细下垫面信息引入对暴雨模拟的影响

首先根据2000年环北京实际的精细下垫面布局资料（500m分辨率）,按美国USGS陆面资料分类标准（25类）对其提供的全球30 s经纬分辨率（≈1 km）下垫面分类资料进行了更新设计.进而针对一个北京夏季暴雨过程,利用10:3.3km双向双重嵌套的MM5V3.6-Noah LSM陆气耦合模式进行24h数值对比试验,研究了北京精细下垫面信息引入对暴雨的影响.分析表明：新设计的陆面资料更真实地反映了环北京区域的下垫面结构,尤其针对北京城区面积迅增特征;同时还修正了原资料将亚洲中纬度区域落叶阔叶林下垫面类型归属为热带(或亚热带)稀疏大草原类型的问题.其在数值天气模式中的引入会对短期暴雨过程的发生发展产生重要影响.对此次暴雨主要降水中心的模拟,12h差值分布范围远达30km以上,中心值相对差异可达30%.研究发现在城市下垫面和大气相互间存在一个重要的相互影响机制,即由于城区面积的扩大会导致自然植被减少,进而会减少地表蒸发及相应局地大气水分供应、加深边界层高度并增强大气水汽混合,这不利于降水的发生发展.     

基于WRF模式的青藏高原斜坡和平台加热影响亚洲夏季风的模拟研究

青藏高原大地形的热力强迫作用对亚洲夏季风的形成和发展具有重要的影响.本文利用较高分辨率的WRF区域模式,探讨了高原不同区域(斜坡和平台)的地形加热分别对南亚夏季风和东亚夏季风的影响.结果表明:高原南部喜马拉雅山脉的斜坡地形加热对其周围局地的环流形势和降水影响十分明显,是南亚夏季风北支分量形成和维持的主导因子,也是斜坡上气流爬坡和降水发生的必要条件.斜坡加热对东亚夏季风也有明显的增强作用,它不仅加强了中国东部低空西南季风环流,还会造成北部南下的异常干冷空气的响应.斜坡上的地形加热作用也是对流层高层暖中心位置维持在斜坡上空的一个重要原因.而高原平台加热对季风环流和降水的影响虽然没有喜马拉雅山脉斜坡加热那么显著,但是对南亚夏季风的影响范围更广,对经向哈得来环流影响更明显,能够调控高原以外更远处热带洋面上的西南季风环流.通过比较高原不同区域地形加热条件下的多种季风指数,进一步表明了高原地形加热对南亚和东亚夏季风均有增强作用,但是高原不同区域的地形加热对两类夏季风子系统又会产生不一样的影响.     

地表拖曳对斜压波地面锋生的影响

利用三维非静力中尺度数值模式MM5模拟了干湿大气条件中纬度典型斜压波及其锋面系统的生成与演变过程, 重点讨论地表拖曳对干、湿大气中地面锋结构、锋生过程的影响作用. 研究结果表明, 在干大气中, 地表拖曳力对地面锋锋生具有双向作用, 一方面是锋消作用, 主要体现在地表拖曳力减慢地面锋锋生、地面斜压波系统发展; 另一方面, 地表拖曳力导致强的非地转流形成, 从而延长了冷锋锋生过程维持时间, 有利于冷锋强度增大. 同时地表拖曳力可以造成边界层内锋面近乎垂直于地面, 导致锋前垂直运动增强, 这些结果进一步推广了谈哲敏和伍荣生的理论结果. 在湿大气中, 地表拖曳过程对锋面雨带分布有重要的影响作用, 地表拖曳力可减缓对流上 升, 从而导致地表能量的耗散减缓. 当大气低层湿度较小时, 对流不是很强, 地表拖曳力可减缓地表水汽、能量的迅速耗散, 且在锋后边界层中产生摩擦辐合上升区, 这些上升区可逐渐东移到冷锋前, 补偿了锋前上升带的强度, 有利于冷锋降水的维持. 当大气低层湿度场很强时, 对流发展比较旺盛, 此时地表拖曳对低层水汽与能量的束缚作用相对较弱, 相应地表拖曳对锋面及其降水系统影响较小.     

积云并合扩展层化型积层混合云的数值模拟分析

积层混合云是我国的主要降水云系,也是人工影响天气的主要作业对象,从云降水物理的角度来研究云系的形成和发展维系具有重要的意义.基于积层混合云的重要性,本文从个例研究入手,利用中尺度数值模式WRF,模拟2005年5月17~18日发生在我国西南山区（主要以贵州省为主）的积层混合云降水过程.发现这次过程是由对流云并合扩大层化形成的.云系形成以后,云系附近会不断有对流云生成,并在移动过程中并合进入云系,补充云系发展维系所需的含水量和能量,促使云系不断维持.在积层混合云系的内部,对流云和层状云区不断地发生作用.对流云给周围的层状云不断输送含水量和能量,支持着层状云的发展.云系内部两种云相互作用的结果体现在：对流云内的上升气流速度逐渐渐小,层状云的上升气流速度不断维持,总上升气流面积区扩大.对流云的降水量不断减小,而层状云的降水不断维持,带来了大面积持续时间很长的降水.     

碧利斯(0604)暴雨增幅的云微物理成因

2006年"碧利斯"台风登陆中国大陆后,在湖南、广东和江西三省交界附近地区造成明显暴雨增幅,造成十分严重的灾害,影响巨大.本文利用高分辨率数值模拟资料,从微观云物理过程角度出发,研究了"碧利斯"暴雨增幅发生前和增幅强降水发生时段云微物理特征的差异,探讨了登陆台风暴雨增幅云微物理方面的可能成因,结果指出:暴雨增幅前后,强降水区云微物理特征存在明显不同,与降水强度的明显增强相伴,云中各种水凝物含量也明显增加,其中云冰、雪和霰等固态水凝物的增加尤为显著,冰相过程对地面降水的贡献明显加大,降水云系发展旺盛、高大;云微物理转化率的对比分析发现,暴雨增幅时段,由水汽凝结过程显著增强所带来的云水的增加,主要通过两个途径作用于暴雨增幅:一是通过云中雨水对云水的碰并收集,促进雨水含量显著增加,进而增强地面降水;二是通过云中雪粒子对云水的碰并造成雪粒子含量增加,增加的雪粒子又被云中霰粒子碰并收集造成霰含量增长,进而由霰粒子融化为雨水,并最终作用于地面降水的增幅.文中最后通过分析总结给出了"碧利斯"暴雨增幅云微物理成因示意图.     

气溶胶污染对夏季降水精细数值预报的影响研究

本文基于北京地区快速更新循环同化和预报系统(BJ-RUC),通过对北京2014年汛期的32次降水过程进行数值模拟和敏感性试验,研究了气溶胶污染引起的云滴数增加对北京夏季降水精细数值预报的影响.研究结果表明:BJ-RUC系统中默认的云滴数浓度偏低,采用根据环境模式预测结果计算的云滴数浓度后(BJ-RUC-Nc算例)预报效果更接近实况;气溶胶增多可以增加或减少降水.具体表现在:(1)2014年汛期北京地区降水主要集中在东北区域(平谷、密云、顺义)、城区次之,数值模式能较好的反映出降水的整体分布和趋势;(2)BJ-RUCNc对于小雨、中雨、暴雨的预报优于BJ-RUC,更接近观测;(3)气溶胶浓度增加,当水汽供应充足(不足)时,会促进(抑制)降水性粒子形成,降水效率提高(降低),降水增加(减少).     

长江流域梅雨锋暴雨过程的中尺度结构个例分析

利用外场试验资料, 用双多普勒雷达技术和径向速度场分析方法, 研究了2002年7月22~23日发生在长江流域一次暴雨过程的中尺度结构动力特征和演变过程. 结果表明, 在西南-东北取向1000 km长的暴雨雨带中, 存在有许多尺度在20~50 km大小的β或γ中尺度强回波带或回波团, 在长江中游, 混合性强降水雨带在长200 km以上的低空切变线上形成; 在切变线南侧的低空西南急流和北侧的偏东气流共同作用下形成上升气流, 对流云得到发展, 切变线低空风场的扰动、中尺度切变和β中尺度辐合是造成对流发展的原因; 新回波常常在老回波右后侧生成, 并移向西南气流区, 从而得到充足的水汽, 这种回波发展旺盛, 持续时间长. β中尺度对流系统常常伴有尺度更小的中尺度涡旋和中尺度辐合等γ中尺度结构, 这些γ中尺度结构在强对流的发展过程中也起了很重要的作用.     

风垂直切变对中尺度地形对流降水影响的研究

针对长江中下游中尺度地形特点以及暴雨过程发生发展期间风垂直切变的主要观测特征,设计了一系列中尺度地形的三维理想数值试验,分析了干大气地形流和重力波特征,探讨了条件不稳定湿大气地形对流降水的模态分布,在此基础上研究了圆形、直线风垂直切变和切变厚度对中尺度地形对流降水强度和模态分布的影响.结果发现:在 F_r≈1的干大气条件下,气流遇到地形后分支、绕流和爬升现象同时存在,地形激发的重力波在水平和垂直方向上传播,其在迎风坡、背风坡、地形上游和下游的振幅不同,并组织出不同强度的垂直上升运动.在F_r > 1的条件不稳定湿大气下,地形对流降水主要存在三种模态,即迎风坡和背风坡准静止对流降水以及地形下游移动性对流降水,地形对流降水的形成与重力波在低层组织的上升运动密切相关.风垂直切变对地形对流降水的强度和模态分布有重要作用,其中圆形风垂直切变(风随高度旋转)不仅影响地形下游对流降水系统的移动方向,而且影响迎风坡和背风坡山脚处对流降水中心的分布和强度;直线风垂直切变(风随高度无旋转)主要影响地形对流降水的移动速度和强度.风随高度自下而上顺(逆)时针旋转,地形对流系统向下游传播时向右(左)偏移.风垂直切变主要通过影响地形重力波的结构和传播以及对流系统的形成、移动方向和速度,来影响地形对流降水的模态分布,其中对流层中低层的风垂直切变对地形对流降水强度和模态分布有重要影响.     

有限水域近水面气流演化特征试验研究

气流作为湖泊、湿地等有限水域的主要驱动力之一,其演化特征是研究水-气间能量与物质传递的基础,决定了水域水环境与水生态格局,具有重要研究意义.采用室内风洞水槽,通过设置不同试验风速,研究近水面气流与风速和吹程的响应关系与演化特征,并将光滑壁面条件下的气流特征作为试验对照组,分析时间域内风速垂向分布形式、摩阻风速、特征粗糙...     

Airside velocity measurements over the wind-sheared water surface using particle image velocimetry

An experimental investigation of the airflow structure in the near surface region over the wind-sheared air–water interface is reported. The two-dimensional velocity fields in a plane perpendicular to the water surface were measured using particle image velocimetry (PIV) technique over a wind speed range from 1.5 to 4.4 m s⁻¹. The results show a reduction in the mean velocity magnitudes and the tangential stresses when gravity waves appear on the surface. An enhanced vorticity layer was observed immediately above the water surface that extended to a height of approximately 2 cm. The vorticity was enhanced by an order of magnitude, and the energy dissipation rate was enhanced by a factor of 7 in this layer at all wind speeds. The vertical profiles of Reynolds stress, energy production, and dissipation indicate the contribution of surface waves in the enhanced transfer of momentum and energy between the two fluids. The results in this study show that the flow dynamics in a layer immediately adjacent to the water surface whose thickness is of the order of the significant wave height is significantly different from that at greater heights. Thus, it is concluded that the quantitative investigation of the flow in the immediate vicinity of the interface is vital for an improved understanding of the heat, mass, and momentum exchange between air and water. The present study demonstrates that PIV is an effective technique to accurately measure the velocity fields in this region.     

Numerical Analysis of Groundwater Ridging Processes Considering Water‐Air Flow in a Hillslope

In this study, a water‐air two‐phase flow model was employed to investigate the formation, extension, and dissipation of groundwater ridging induced by recharge events in a hypothetical hillslope‐riparian zone, considering interactions between the liquid and gas phases in soil voids. The simulation results show that, after a rain begins, the groundwater table near the stream is elevated instantaneously and significantly, thereby generating a pressure gradient driving water toward both the stream (the discharge of groundwater to the stream) and upslope (the extension of groundwater ridging into upslope). Meanwhile, the airflow upslope triggered by the advancing wetting front moves downward gradually. Therefore, the extension of groundwater ridging into upslope and the downward airflow interact within a certain region. After the rain stops, groundwater ridging near the stream declines quickly while the airflow in the lower part of upslope is still moving into the hillslope. Thus, the airflow upslope mitigates the dissipation of groundwater ridging. Additionally, the development of groundwater ridging under different conditions, including rain intensity, intrinsic permeability, capillary fringe height, and initial groundwater table, was analyzed. Changes in intrinsic permeability affect the magnitude of groundwater ridging near the stream, as well as the downward speed of airflow, thereby generating highly complex responses. The capillary fringe is not a controlling factor but an influence factor on the formation of groundwater ridging, which is mainly related to the antecedent moisture. It was demonstrated that groundwater ridging also occurs where an unsaturated zone occurs above the capillary fringe with a subsurface lateral flow.     

The cloud-microphysical cause of torrential rainfall amplification associated with Bilis (0604)

After its landfall in China’s mainland in 2006, Typhoon Bilis brought about torrential rainfall amplification at the edge of Guangdong, Jiangxi, and Hunan provinces, causing severe disasters. From a cloud-microphysical perspective, we discuss the differences of cloud-microphysical processes before and during the precipitation amplification and possible causes of the rainfall amplification by using high-resolution simulation data. The results show that the cloud-microphysical characteristics during the above two periods are significantly different. With the distinct increase in the rainfall intensity, the cloud hydrometeor contents increase markedly, especially those of the ice-phase hydrometeors including ice, snow and graupel, contributing more to the surface rainfall. The clouds develop highly and vigorously. Comparisons of conversion rates of the cloud hydrometeors between the above two periods show that the distinct increases in the cloud water content caused by the distinct enhancement of the water vapor condensation rate contribute to the surface rainfall mainly in two ways. First, the rain water content increases significantly by accretion of cloud water by rain water, which thus contributes to the surface rainfall. Second, the accretion of cloud water by snow increases significantly the content of snow, which is then converted to graupel by accretion of snow by graupel. And then the graupel melts into rain water, which subsequently contributes to the surface rainfall amplification. In summary, a flow chart is given to clarify the cloud-microphysical cause of the torrential rainfall amplification associated with Bilis.     

Sensitivity of model parameterizations for simulated latent heat flux at the snow surface for complex mountain sites

The snowcover energy balance is typically dominated by net radiation and sensible and latent heat fluxes. Validation of the two latter components is rare and often difficult to undertake at complex mountain sites. Latent heat flux, the focus of this paper, is the primary coupling mechanism between the snow surface and the atmosphere. It accounts for the critical exchange of mass (sublimation or condensation), along with the associated snowcover energy loss or gain. Measured and modelled latent heat fluxes at a wind‐exposed and wind‐sheltered site were compared to evaluate variability in model parameters. A well‐tested and well‐validated snowcover energy balance model, Snobal, was selected for this comparison because of previously successful applications of the model at these sites and because of the adjustability of the parameters specific to latent heat transfer within the model. Simulated latent heat flux and snow water equivalent (SWE) were not sensitive to different formulations of the stability profile functions associated with heat transfer calculations. The model parameters of snow surface roughness length and active snow layer thickness were used to improve latent heat flux simulations while retaining accuracy in the simulation of the SWE at an exposed and sheltered study site. Optimal parameters for simulated latent heat flux and SWE were found at the exposed site with a shorter roughness length and thicker active layer, and at the sheltered site with a longer roughness length and thinner active layer. These findings were linked to physical characteristics of the study sites and will allow for adoption into other snow models that use similar parameters. Physical characteristics of wind exposure and cover could also be used to distribute critical parameters in a spatially distributed modelling domain and aid in parameter selection for application to other watersheds where detailed information is not available. Copyright © 2012 John Wiley & Sons, Ltd.     

Impact of the surface temperature and vertical shear of zonal wind on the dynamics of a simple two-layer model of the atmosphere

The paper presents and analyzes, from the point of view of smooth dynamic systems theory, a two-layer baroclinic model of the troposphere in geostrophic approximation. The model describes airflow in β-channel within the tropospheric part of the main Hadley circulation cell. It enables to obtain, after application of the Galerkin method, a fairly simple low-parametric dynamic system describing the phenomena of non-linear interactions, bifurcations and blocking in the atmosphere. This enables to take into consideration such basic factors influencing the atmospheric dynamics like the heat exchange within the surface, orography, vertical variability of zonal wind and hydrostatic stability. Impact of zonal thermal variability of the surface and vertical shear of zonal wind in the troposphere on the orographic bifurcation was investigated and the oscillation character in the dynamic system after Hopf bifurcation of the second kind was analyzed. Additionally, the model dynamics was investigated in conditions including momentum forcing in the upper and lower parts of the troposphere and excluding orographic interaction, as well as in the conditions of thermal interaction between the troposphere and the surface for the vertical shear of zonal wind in both tropospheric layers. Impact of the mean zonal wind in the troposphere on the properties of model dynamics was assessed. It was proved that zonally varied surface temperature and layered mean zonal wind in the atmosphere are the parameters that have basic influence on the model dynamics. They cause numerous bifurcations and strongly influence the periods of oscillations of the model variables. They are often Hopf bifurcations of the second kind during which tropospheric states fairly distant from the ones before the bifurcations are generated. This significantly influences the model predictability.     

The cold rain‐on‐snow event of June 2013 in the Canadian Rockies — characteristics and diagnosis

The June 2013 flood in the Canadian Rockies featured rain‐on‐snow (ROS) runoff generation at alpine elevations that contributed to the high streamflows observed during the event. Such a mid‐summer ROS event has not been diagnosed in detail, and a diagnosis may help to understand future high discharge‐producing hydrometeorological events in mountainous cold regions. The alpine hydrology of the flood was simulated using a physically based model created with the modular cold regions hydrological modelling platform. The event was distinctive in that, although at first, relatively warm rain fell onto existing snowdrifts inducing ROS melt; the rainfall turned to snowfall as the air mass cooled and so increased snowcover and snowpacks in alpine regions, which then melted rapidly from ground heat fluxes in the latter part of the event. Melt rates of existing snowpacks were substantially lower during the ROS than during the relatively sunny periods preceding and following the event as a result of low wind speeds, cloud cover and cool temperatures. However, at the basin scale, melt volumes increased during the event as a result of increased snowcover from the fresh snowfall and consequent large ground heat contributions to melt energy, causing snowmelt to enhance rainfall–runoff by one fifth. Flow pathways also shifted during the event from relatively slow sub‐surface flow prior to the flood to an even contribution from sub‐surface and fast overland flow during and immediately after the event. This early summer, high precipitation ROS event was distinctive for the impact of decreased solar irradiance in suppressing melt rates, the contribution of ground heat flux to basin scale snowmelt after precipitation turned to snowfall, the transition from slow sub‐surface to fast overland flow runoff as the sub‐surface storage saturated and streamflow volumes that exceeded precipitation. These distinctions show that summer, mountain ROS events should be considered quite distinct from winter ROS and can be important contributors to catastrophic events. Copyright © 2016 John Wiley & Sons, Ltd.     

Point vortex dynamics for coupled surface/interior QG and propagating heton clusters in models for ocean convection

Abstract

The dynamic behavior of baroclinic point vortices in two-layer quasi-geostrophic flow provides a compact model for studying the transport of heat in a variety of geophysical flows including recent heton models for open ocean convection as a response to spatially localized intense surface cooling. In such heton models, the exchange of heat with the region external to the compact cooling region reaches a statistical equilibrium through the propagation of tilted heton clusters. Such tilted heton clusters are aggregates of cyclonic vortices in the upper layer and anti-cyclonic vortices in the lower layer which collectively propagate almost as an elementary tilted heton pair even though the individual vortices undergo shifts in their relative locations. One main result in this paper is a mathematical theorem demonstrating the existence of large families of long-lived propagating heton clusters for the two-layer model in a fashion compatible to a remarkable degree with the earlier numerical simulations. Two-layer quasi-geostrophic flow is an idealization of coupled surface/interior quasi-geostrophic flow. The second family of results in this paper involves the systematic development of Hamiltonian point vortex dynamics for coupled surface/interior QG with an emphasis on propagating solutions that transport heat. These are novel vortex systems of mixed species where surface heat particles interact with quasi-geostrophic point vortices. The variety of elementary two-vortex exact solutions that transport heat include two surface heat particles of opposite strength, tilted pairs of a surface heat particle coupled to an interior vortex of opposite strength and two interior tilted vortices of opposite strength at different depths. The propagation speeds of the tilted elementary hetons in the coupled surface/interior QG model are compared and contrasted with those in the simpler two-layer heton models. Finally, mathematical theorems are presented for the existence of large families of propagating long-lived tilted heton clusters for point vortex solutions in coupled surface/interior QG flow.     

Lagrangian circulation and forbidden zone on the West Florida Shelf

This paper presents some recent results of drifters released on the West Florida Shelf during 1996–1997 and compares with the numerical model results of the wind-driven circulation. Using satellite tracked surface drifters during the one year period from February 1996 to February 1997, a drifter free region, called the “forbidden zone”, is found over the southern portion of the West Florida Shelf. This finding is consistent with historical drift bottle data and with a recent numerical model study of the West Florida Shelf circulation response to climatological wind forcing. Direct drifter simulations by numerical model during March 1996 show a good agreement with both the in situ ADCP current observation and drifter observation. Three mechanisms are proposed for the observed Lagrangian features. The primarily dynamic mechanism is the along-shore wind forcing, which induces a coastal jet that tends to leave the coast and the bottom onshore and near surface offshore transports. The second one is the convergent coastal geometry and bottom topography for the southward flow in central shelf near Tampa Bay that enforces the coastal jet and the bottom and near surface transport. The last is a kinematic one, simply due to the short along-shore Lagrangian excursion, driven by the typical synoptic weather systems. Thus near surface shelf waters over the north may not reach the southern coast of the West Florida. Implication is that surface hazard such as oil spill that may occur outside of the southern West Florida shelf may not greatly impact the southern coastal region except Florida Keys. However, the biological and chemical patches over the north that may occur in the water column such as red tides still can easily reach the southern coastal region through the subsurface and bottom waters.     

Air–sea interactions during strong winter extratropical storms

A high-resolution, regional coupled atmosphere–ocean model is used to investigate strong air–sea interactions during a rapidly developing extratropical cyclone (ETC) off the east coast of the USA. In this two-way coupled system, surface momentum and heat fluxes derived from the Weather Research and Forecasting model and sea surface temperature (SST) from the Regional Ocean Modeling System are exchanged via the Model Coupling Toolkit. Comparisons are made between the modeled and observed wind velocity, sea level pressure, 10 m air temperature, and sea surface temperature time series, as well as a comparison between the model and one glider transect. Vertical profiles of modeled air temperature and winds in the marine atmospheric boundary layer and temperature variations in the upper ocean during a 3-day storm period are examined at various cross-shelf transects along the eastern seaboard. It is found that the air–sea interactions near the Gulf Stream are important for generating and sustaining the ETC. In particular, locally enhanced winds over a warm sea (relative to the land temperature) induce large surface heat fluxes which cool the upper ocean by up to 2 °C, mainly during the cold air outbreak period after the storm passage. Detailed heat budget analyses show the ocean-to-atmosphere heat flux dominates the upper ocean heat content variations. Results clearly show that dynamic air–sea interactions affecting momentum and buoyancy flux exchanges in ETCs need to be resolved accurately in a coupled atmosphere–ocean modeling framework.