Estimating the components of the sensible heat budget of a tall forest canopy in complex terrain

Ultrasonic wind measurements, sonic temperature and air temperature data at two heights in the advection experiment MORE II were used to establish a complete budget of sensible heat including vertical advection, horizontal advection and horizontal turbulent flux divergence. MORE II took place at the long-term Carbo-Europe IP site in Tharandt, Germany. During the growing period of 2003 three additional towers were established to measure all relevant parameters for an estimation of advective fluxes, primarily of CO₂. Additionally, in relation to other advection experiments, a calculation of the horizontal turbulent flux divergence is proposed and the relation of this flux to atmospheric stability and friction velocity is discussed. In order to obtain a complete budget, different scaling heights for horizontal advection and horizontal turbulent flux divergence are tested. It is shown that neglecting advective fluxes may lead to incorrect results. If advective fluxes are taken into account, the sensible heat budget based upon vertical turbulent flux and storage change only, is reduced by approximately 30%. Additional consideration of horizontal turbulent flux divergence would in turn add 5–10% to this sum (i.e., the sum of vertical turbulent flux plus storage change plus horizontal and vertical advection). In comparison with available energy horizontal advection is important at night whilst horizontal turbulent flux divergence is rather insignificant. Obviously, advective fluxes typically improve poor nighttime energy budget closure and might change ecosystem respiration fluxes considerably.     

A study on the variability of atmospheric and oceanic processes over the Arabian Sea during contrasting monsoons

Summary The atmospheric and oceanic conditions associated with the southwest monsoon during the contrasting monsoon years of 2002 and 2003 over the Arabian Sea have been analyzed in the present study. Early onset of southwesterlies and reduced net heat gain due to low solar radiation were responsible for low sea-surface temperatures (SSTs) over the Arabian Sea during 2002 pre-monsoon (particularly in May). Conversely, light winds and an increased net heat gain set up the pre-monsoon warming in 2003. The development and intensification of deep convection over a large area of the Arabian Sea prior to the onset of the monsoon was observed during 2003, but was absent in 2002. Weak cross equatorial flow and a weak low level jet over the Arabian Sea reduced moisture transport towards the Indian subcontinent in July 2002. This scenario helped to contribute to a prolonged break in monsoon conditions during July. However, no such break in conditions occurred during July 2003. In 2002, the summer monsoon cooling of the Arabian Sea occurred well before July, whereas in 2003 cooling occurred during July. Estimates of wind driven Ekman (horizontal) and vertical transports showed maximum values in the month of June (July) in 2002 (2003). These estimates clearly show the importance of horizontal and vertical advection in the summer cooling of the Arabian Sea. During the southwest monsoon period, the Arabian Sea was warmer in 2003 than in 2002. Late onset of the southwesterlies in June, late cooling of the Arabian Sea in July, and downwelling Rossby wave propagation were responsible for the warm SSTs in 2003. Weak wind stress curl in July dampened the westward propagating sea surface height anomaly signals (Rossby waves) before they reached the western Arabian Sea in 2002, whereas, in 2003 strong wind stress curl enhanced Rossby wave propagation. During the summer monsoon period, subsurface temperatures in the south central Arabian Sea were warmer in 2003 than in 2002, particularly in July and August. Strong Ekman convergence, solar penetration, and downwelling (downward velocities) are responsible for the enhanced subsurface warming in 2003.     

The scaling of flow in vegetative structures

Sensible heat flux data obtained from the U.S. Environmental Protection Agency's Regional Air Pollution Study in St. Louis, Missouri are presented and discussed. Large spatial variations exist in heat flux on both a land-use scale and the urban scale. Arguments based upon empirical data and sampling theory show that estimates of heat flux representative of an upwind fetch l _x require a minimum observation height proportional to z ^3/4. The influence of advection on the magnitude of the heat flux is also explored for both the urban and sub-urban or land-use scales. The data clearly indicate that advection can strongly modulate and even dominate the vertical heat flux above surfaces in areas which maintain large horizontal temperature gradients. The advection contribution is positive for cold air advection and negative for warm air advection, and may result from either the urban heat island or land-use mesoscale features. The depth of advective influence is directly proportional to the horizontal scale of the phenomenon and inversely proportional to horizontal temperature gradient.     

东中国海黑潮海洋锋的季节变化特征及其成因

利用高分辨率海洋和大气再分析资料研究了东中国海黑潮海洋锋的季节变化特征及其成因。研究表明,东海黑潮海洋锋存在明显的季节变化,从冬季到次年春季逐渐增强,并在春季达到最强,初夏以后强度逐渐减弱,9—10月达到全年最弱。通过诊断混合层热流量方程发现,东海黑潮区一年四季均表现为暖的温度平流,有利于海洋锋的形成和维持,该暖平流也存在季节变化并在春季达到最大,对海洋锋在春季达到最强起了重要作用。海气界面净热通量在秋冬季对海洋锋的形成有促进作用,有利于海洋锋增强,而在春夏季则起抑制作用,促进海洋锋消亡。温度垂直输送全年对海洋锋起一定程度的抑制作用。总之,在海温水平平流和海表净热通量的共同作用下导致海洋锋春季达到最强,而夏秋季海表净热通量和温度垂直输送作用的共同作用致使海洋锋减弱并最后消失。因此,海洋的动力和热力共同作用导致了东海黑潮海洋锋的季节变化,其中海温水平平流和海表净热通量对海洋锋的季节变化起主要作用,而温度垂直输送项对海洋锋的发展起抑制作用,但影响相对较小。     

Divergence of turbulent fluxes in the surface layer: case of a coastal city

This study quantifies the processes that take place in the layer between the mean building height and the measurement level of an energy balance micrometeorological tower located in the dense old core of a coastal European city. The contributions of storage, vertical advection, horizontal advection and radiative divergence for heat are evaluated with the available measurements and with a three-dimensional, high-resolution meteorological simulation that had been evaluated against observations. The study focused on a summer period characterized by sea-breeze flows that affect the city. In this specific configuration, it appears that the horizontal advection is the dominant term. During the afternoon when the sea breeze is well established, correction of the sensible heat flux with horizontal heat advection increases the measured sensible heat flux up to 100 W m⁻². For latent heat flux, the horizontal moisture advection converted to equivalent latent heat flux suggests a decrease of 50 W m⁻². The simulation reproduces well the temporal evolution and magnitude of these terms.     

On the upper oceanic heat budget in the south china sea:Annual cycle

l.Intr0ducti0nInrecentyCarstheSouthChinaSea(SCS)hasbecomeoneofthemostimP0rtantregionsinthelocalair-seainteractionresearchbecauseofitssPecialgeographicPOsitionandsemi-encloseddeepbasincharacteristics.TheimP0rtanceoftheSCSisembodiedmainlyinwhichisoneofthekeyheatandmoisturesourcesofatmosphericcirculationineasternAsia.TheonsetandmaintenanceoftheSCSmonsoonarecloselyconnectedwiththelargeheattransportfromtheSCStoair(Yan,l997).Thestudyonspatial-temporalvariationofair-seaheatexchangeintheSCS…     

冬季北大西洋涛动对热带印度洋海表热通量的影响

利用1979—2017年TropFlux海气热通量资料、ERA5再分析资料及HadISST资料,分析了冬季北大西洋涛动（North Atlantic Oscillation,NAO）与同期热带印度洋海气热通量的关系。结果表明,NAO指数与热带印度洋海气净热通量整体上呈负相关,意味着NAO为正位相时,海洋向大气输送热量,其显著区域主要位于热带西印度洋（50°～70°E,10°S～10°N）。净热通量的变化主要依赖于潜热通量和短波辐射的变化;潜热通量和短波辐射在NAO正（负）位相事件期间的贡献率分别为72.96％和61.48％（71.72％和57.06％）。NAO可通过Rossby波列影响印度洋地区局地大气环流,进而影响海气热通量;当NAO为正位相时,波列沿中低纬路径传播至印度洋地区,在阿拉伯海北部对流层高层触发异常反气旋环流。该异常反气旋性环流加强了阿拉伯高压,使得北印度洋偏北风及越赤道气流加强。伴随风速的加强,海面蒸发增强,同时加强的越赤道气流导致热带辐合带强度偏强,深对流加强引起对流层水汽和云量增多,进而引起海表下行短波辐射减少。     

暴雨中尺度环境场特征及积云对流的反馈作用

本文利用一套具有中尺度分辨率的观测资料,对发生在1983年6月下旬的一次长江流域梅雨暴雨过程进行了诊断研究,并比较了该过程的对流降水活跃期和非活跃期的中尺度环境场特征。结果表明,两个时期的动量场存在明显差异,而水热场差异不大,Q₁和Q₂以及涡度收支也存在明显差异。在降水活跃期（暴雨集中期）,积云对流对能量和涡度的垂直输送有着重要作用;在水热收支中,起主要作用的是垂直输送项,潜热加热基本上为抬升冷却所平衡,水汽的垂直输送是积云对流的主要水汽源;在涡度收支中,低层散度项和扭转项制造正涡度,并通过积云对流向上输送,无论在低层或高层积累的正涡度都被平流非线性过程所耗损。      

Aspects of CO2 advection measurements

Observations of vegetation–atmosphere exchange of carbon dioxide (CO₂) by the eddy covariance (EC) technique are limited by difficult conditions such as nighttime and heterogeneous terrain. Thus, advective flux components are included into the net ecosystem exchange (NEE) budget. However, advection measurements are experimentally challenging and do not always help to solve the night flux problem of the EC technique. This study investigates alternative methods for the observation of horizontal advection, in particular horizontal concentration gradients, as well as different approaches to coordinate rotation and vertical advection. Continuous high-frequency measurements of the horizontal CO₂ concentration field are employed and compared to the often used discontinuous sequential sampling. Significant differences were found in the case of 30-min mean concentration values between the conventional discontinuous sampling approach and the complete observation of the time series by continuous sampling. Estimates of vertical advection rely on accurate estimates of vertical wind velocity ( $\emph{w}$ ). Therefore, different approaches to the planar fit coordinate rotation have been investigated. Sector-wise rotation was able to eliminate directional dependencies of mean $\emph{w}$ . Furthermore, the effect of the data set length used for rotation (window length) was investigated and was found to have significant impact on estimates of vertical advection, with larger window lengths yielding about 50% larger vertical advection. A sequential planar fit with controlled window length is proposed to give reproducible results. The different approaches to the measurement and calculation of horizontal and vertical advection presented are applied to data obtained during the exchange processes in mountainous region experiment at the FLUXNET site Waldstein–Weidenbrunnen (DE-Bay). Estimates of NEE including advection are compared to NEE from turbulent and storage flux alone without advection. NEE including vertical advection with sector-wise planar fit rotation and controlled window length and including horizontal advection from continuous gradient measurements, which were comprehensively bias corrected by a new approach, did compare well with the expected night flux error, with meteorological drivers of the fluxes and with soil chamber measurements. Unrealistically large and noisy values of horizontal advection from the conventional discontinuous sampling approach, which lead to unrealistic values of NEE, could be eliminated by the alternative approaches presented. We therefore suggest the further testing of those approaches at other sites in order to improve the accuracy of advection measurements and, subsequently, estimates of NEE.     

Impacts of Sea Surface Temperature on the Interannual Variability of Winter Haze Days in Guangdong Province

The impact of sea surface temperature (SST) on winter haze in Guangdong province (WHDGD) was analyzed on the interannual scale. It was pointed out that the northern Indian Ocean and the northwest Pacific SST play a leading role in the variation of WHDGD. Cold (warm) SST anomalies over the northern Indian Ocean and the Northwest Pacific stimulate the eastward propagation of cold (warm) Kelvin waves through the Gill forced response, causing Ekman convergence (divergence) in the western Pacific, inducing abnormal cyclonic (anticyclonic) circulation. It excites the positive (negative) Western Pacific teleconnection pattern (WP), which results in the temperature and the precipitation decrease (increase) in Guangdong and forms the meteorological variables conditions that are conducive (not conducive) to the formation of haze. ENSO has an asymmetric influence on WHDGD. In El Ni?o (La Ni?a) winters, there are strong (weak) coordinated variations between the northern Indian Ocean, the northwest Pacific, and the eastern Pacific, which stimulate the negative (positive) phase of WP teleconnection. In El Ni?o winters, the enhanced moisture is attributed to the joint effects of the horizontal advection from the surrounding ocean, vertical advection from the moisture convergence, and the increased atmospheric apparent moisture sink (Q2) from soil evaporation. The weakening of the atmospheric apparent heat source (Q1) in the upper layer is not conducive to the formation of inversion stratification. In contrast, in La Ni?a winters, the reduced moisture is attributed to the reduced upward water vapor transport and Q2 loss. Due to the Q1 increase in the upper layer, the temperature inversion forms and suppresses the diffusion of haze.     

Mixed layer thermodynamics of the Southern South China Sea

Seasonal and inter-annual variability of the mixed layer temperature in the Southern South China Sea (SSCS) is investigated using a regional ocean circulation model simulation. The mixed layer depth (MLD) over the SSCS exhibits a strong seasonal signal with deeper MLDs during the northeast and southwest monsoons. The main factor that drives the mixed layer temperature variation in the SSCS is the air-sea heat fluxes, with vertical ocean processes acting as a relatively weak negative feedback. In general, the budget analysis demonstrates a net balance between the vertical ocean processes and surface heat flux during the pre-monsoon and southwest monsoon. Northeast monsoon period is noted by an offsetting of surface heat flux, horizontal and vertical ocean processes. The first dominant mode of mixed layer temperature inter-annual variability in the SSCS shows significant correlation (0.34) with the El Nino phenomenon in the Pacific Ocean and is best correlated (0.67) with a lag of 5 months.     

The Influence of Advection on the Short Term CO 2 -Budget in and Above a Forest Canopy

An experimental micrometeorological set-up was established at the CARBOEURO-FLUX site in Tharandt, Germany, to measure all relevant variables for the calculation of the vertical and horizontal advective fluxes of carbon dioxide. The set-up includes two auxiliary towers to measure horizontal and vertical CO₂ and H₂O gradients through the canopy, and to make ultrasonic wind measurements in the trunk space. In combination with the long-term flux tower an approximately even-sided prism with a typical side-length of 50 m was established. It is shown that under stable (nighttime) conditions the mean advective fluxes have magnitudes on the same order as the daily eddy covariance (EC) flux, which implies that they play a significant, but not yet fully understood, role in the carbon budget equation. The two advective fluxes are opposite and seem to cancel each other at night (at least for these measurements). During the day, vertical advection tends to zero, while horizontal advection is still present implying a flow of CO₂ out of the control volume. From our measurements, a mean daily gain of 2.2 gC m^–2 d^–1 for the horizontal advection and a mean daily loss of 2.5 gC m^–2d^–1 for the vertical advection is calculated for a period of 20 days. However the large scatter of the advective fluxes has to be further investigated. It is not clear yet whether the large variability is natural or due to measurement errors and conceptual deficiencies of the experiment. Similar results are found in the few comparable studies.     

Measurements of the turbulent heat flux over the sea

The turbulent heat flux was measured with two instruments simultaneously over the Baltic Sea by means of the eddy-correlation method. In one observational period, a small but noticeable divergence in heat flux was found, which may be explained by the advection of colder air. The parameterization of heat flux by the bulk method leads to a value for C _Hof 1 × 10^–3.     

Thermal effect of the sea breeze on the structure of the boundary layer and the heat budget over land

The daytime boundary-layer heating process and the air-land heat budget were investigated over the coastal sea-breeze region by means of observations over the Sendai plain in Japan during the summer. In this area, the onset of the sea breeze begins at the coast around 0900 LST, intruding about 35 km inland by late afternoon. The cold sea breeze creates a temperature difference of over 10°C between the coastal and inland areas in the afternoon. On the other hand, warm air advection due to the combination of the counter-sea breeze and land-to-sea synoptic wind occurs in the layer above the cold sea breeze in the coastal region. Owing to this local warm air advection, there is no significant difference in the daytime heating rate over the entire atmospheric boundary layer between the coastal and inland areas. The sensible heat flux from the land surface gradually decreases as distance from the coastline increases, being mainly attributed to the cold sea breeze. The daytime mean cold air advection due to the sea breeze is estimated asQ _adv ^local =–29 W m^–2 averaged over the sea breeze region (035 km from the coastline). This value is 17% of the surface sensible heat fluxH over the same region. The results of a two-dimensional numerical model show that the value ofQ _adv ^local /H is strongly affected by the upper-level synoptic wind direction. The absolute value ofQ _adv ^local /H becomes smaller when the synoptic wind has the opposite direction of the sea breeze. This condition occurred during the observations used in the present study.     

The Influence of Hilly Terrain on Canopy-Atmosphere Carbon Dioxide Exchange

Topography influences many aspects of forest-atmosphere carbon exchange; yet only a small number of studies have considered the role of topography on the structure of turbulence within and above vegetation and its effect on canopy photosynthesis and the measurement of net ecosystem exchange of CO₂ (N_ee) using flux towers. Here, we focus on the interplay between radiative transfer, flow dynamics for neutral stratification, and ecophysiological controls on CO₂ sources and sinks within a canopy on a gentle cosine hill. We examine how topography alters the forest-atmosphere CO₂ exchange rate when compared to uniform flat terrain using a newly developed first-order closure model that explicitly accounts for the flow dynamics, radiative transfer, and nonlinear eco physiological processes within a plant canopy. We show that variation in radiation and airflow due to topography causes only a minor departure in horizontally averaged and vertically integrated photosynthesis from their flat terrain values. However, topography perturbs the airflow and concentration fields in and above plant canopies, leading to significant horizontal and vertical advection of CO₂. Advection terms in the conservation equation may be neglected in flow over homogeneous, flat terrain, and then N_ee = F_c, the vertical turbulent flux of CO₂. Model results suggest that vertical and horizontal advection terms are generally of opposite sign and of the same order as the biological sources and sinks. We show that, close to the hilltop, F_c departs by a factor of three compared to its flat terrain counterpart and that the horizontally averaged F_c-at canopy top differs by more than 20% compared to the flat-terrain case.     

Response of the tropical Indian Ocean SST to decay phase of La Niña and associated processes

Delayed impact of El Niño on Tropical Indian Ocean (TIO) Sea Surface Temperature (SST) variations and associated physical mechanisms are well documented by several studies. However, TIO SST evolution during the decay phase of La Niña and related processes are not adequately addressed before. Strong cooling associated with La Niña decay over the TIO could influence climate over the Indian Oceanic rim including Indian summer monsoon circulation and remotely northwest Pacific circulation. Thus understanding the TIO basin-wide cooling and related physical mechanisms during decaying La Niña years is important. Composite analyses revealed that negative SST anomalies allied to La Niña gradually dissipate from its mature phase (winter) till subsequent summer in central and eastern Pacific. In contrast, magnitude of negative SST anomalies in TIO, induced by La Niña, starts increasing from winter and attains their peak values in early summer. It is found that variations in heat flux play an important role in SST cooling over the central and eastern equatorial Indian Ocean, Bay of Bengal and part of Arabian Sea from late winter to early summer during the decay phase of La Niña. Ocean dynamical processes are mainly responsible for the evolution of southern TIO SST cooling. Strong signals of westward propagating upwelling Rossby waves between 10°S to 20°S are noted throughout (the decaying phase of La Niña) spring and summer. Anomalous cyclonic wind stress curl to the south of the equator is responsible for triggering upwelling Rossby waves over the southeastern TIO. Further, upwelling Rossby waves are also apparent in the Arabian Sea from spring to summer and partly contributing to the SST cooling. Heat budget analysis reveals that negative SST/MLT (mixed layer temperature) anomalies over the Arabian Sea are mostly controlled by heat flux from winter to spring and vertical advection plays an important role during early summer. Vertical and horizontal advection terms primarily contribute to the SST cooling anomalies over southern TIO and the Bay of Bengal cooling is primarily dominated by heat flux. Further we have discussed influence of TIO cooling on local rainfall variations.     

寒潮过程中风浪对黄海海气热量通量和动量通量影响研究

采用2009—2013年CFSR(Climate Forecast System Reanalysis)大气和海洋再分析资料对黄海海气间热量通量和动量通量的特征进行统计分析,并通过FVCOMSWAVE浪流耦合模式对典型寒潮过程中风浪的影响效果进行模拟研究与对比分析。统计结果显示,通量受海表大风、海气温差及海洋环流等因子影响,秋冬季节强烈,春夏季节相对较弱,在寒潮活跃的冷季该海域的海流处于弱流期,风浪对海面通量的作用明显增强。海温特征也显示冷季的不稳定性显著强于暖季,因此该海域冷季具有更强的海气热量通量。沿岸站点的比较显示,南部吕泗站面向更开阔的东海海域,其平均波高高出北部20%左右。这与沿海南部通量强于北部特征对应。数值模拟显示,在寒潮过程中,海气界面热量通量和动量通量输送比多年月平均状态显著增强,动量通量增大1~5倍,热量通量增大1~6倍。寒潮过程入海冷锋走向、强度、移动方向显著影响海面热量通量和动量通量大值区的分布。偏北路寒潮纬向型冷锋入海,其强度东部大于西部,造成通量大值区形成在黄海东北部,而偏西路寒潮经向型冷锋入海,其强度南部大于北部,造成通量大值区形成在黄海南部。同时偏北路径寒潮强度大于偏西路径,海气动量通量响应较偏西路径强约25%,热量通量强约50%。耦合风浪作用的模拟显示,海气间热量通量和动量通量明显增大,对不同强度风浪,浪高增加1.5倍,动量通量最大值增大约2倍,热量通量增大10~160 W/m2;浪高减弱至0.5倍,动量通量最大值则减弱约40%,热量通量减小10~55 W/m2。冷锋及其驱动的风浪强烈影响区域海气通量时空特征。     

On the mechanism of the seasonal variability of SST in the tropical Indian Ocean

A general form of an equation that "explicitly" diagnoses SST change is derived. All other equations in wide use are its special case. Combining with the data from an ocean general circulation model (MOM2) with an integration of 10 years (1987-1996), the relative importances of various processes that determine seasonal variations of SST in the tropical Indian Ocean are compared mainly for January, April, July and October. The main results are as follows. (1) The net surface heat flux is the most important factor affecting SST over the Arabian Sea, the Bay of Bengal and the region south of the equator in January; in April, its influence covers almost the whole region studied; whereas in July and October, this term shows significance only in the regions south of 10°S and north of the equator, respectively. (2) The horizontal advection dominates in the East African-Arabian coast and the region around the equator in January and July; in October, the region is located south of 10°S. (3) The entrainment is s     

Modulation of El Nino-Southern Oscillation by Freshwater Flux and Salinity Variability in the Tropical Pacific

The El Nin o-Southern Oscillation (ENSO) is modulated by many factors; most previous studies have emphasized the roles of wind stress and heat flux in the tropical Pacific. Freshwater flux (FWF) is another environmental forcing to the ocean; its effect and the related ocean salinity variability in the ENSO region have been of increased interest recently. Currently, accurate quantifications of the FWF roles in the climate remain challenging; the related observations and coupled ocean-atmosphere modeling involve large elements of uncertainty. In this study, we utilized satellite-based data to represent FWF-induced feedback in the tropical Pacific climate system; we then incorporated these data into a hybrid coupled ocean-atmosphere model (HCM) to quantify its effects on ENSO. A new mechanism was revealed by which interannual FWF forcing modulates ENSO in a significant way. As a direct forcing, FWF exerts a significant influence on the ocean through sea surface salinity (SSS) and buoyancy flux (Q B ) in the western-central tropical Pacific. The SSS perturbations directly induced by ENSO-related interannual FWF variability affect the stability and mixing in the upper ocean. At the same time, the ENSO-induced FWF has a compensating effect on heat flux, acting to reduce interannual Q B variability during ENSO cycles. These FWF-induced processes in the ocean tend to modulate the vertical mixing and entrainment in the upper ocean, enhancing cooling during La Nin a and enhancing warming during El Nin o, respectively. The interannual FWF forcing-induced positive feedback acts to enhance ENSO amplitude and lengthen its time scales in the tropical Pacific coupled climate system.     

On the subarctic North Atlantic cooling due to global warming

Various ocean reanalysis data reveal that the subarctic Atlantic sea surface temperature (SST) has been cooling during the twentieth century. A similar cooling pattern is found in the doubling CO₂ experiment obtained from the CMIP3 (coupled model intercomparison project third phase) compared to the pre-industrial experiment. Here, in order to investigate the main driver of this cooling, we perform the heat budget analysis on the subarctic Atlantic upper ocean temperature. The net surface heat flux associated with the increased concentration of greenhouse gases heats the subarctic ocean surface. In the most of models, the longwave radiation, latent heat flux, and sensible heat flux exert a warming effect, and the shortwave radiation exerts a cooling effect. On the other hand, the thermal advection by the meridional current reduces the subarctic upper ocean temperature in all models. This cold advection is attributed to the weakening of the meridional overturning circulation, which is related to the reduction in the ocean surface density. In particular, greater warming of the surface air than of the sea surface results in the reduction of surface evaporation and thereby enhanced freshening of the ocean surface water, while precipitation change was smaller than evaporation change. The thermal advections by both the wind-driven Ekman current and the density-driven geostrophic current contribute to cooling in most of the models, where the heat transport by the geostrophic current tends to be larger than that by the Ekman current.