• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Impacts of Increased SST Resolution on the North Pacific Storm Track in ERA-Interim

    2020-10-15 10:09:36ChaoZHANGHailongLIUJinboXIEChongyinLIandPengfeiLIN
    Advances in Atmospheric Sciences 2020年11期

    Chao ZHANG, Hailong LIU, Jinbo XIE, Chongyin LI,, and Pengfei LIN,3

    1College of Meteorology and Oceanography, National University of Defense Technology, Nanjing 211101, China

    2State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics,Institute of Atmospheric Physics, Chinese Academy of Sciences, Beijing 100029, China

    3College of Earth and Planetary Sciences, University of Chinese Academy of Sciences, Beijing 100049, China

    ABSTRACT This study examines the artificial influence of increasing the SST resolution on the storm track over the North Pacific in ERA-Interim. Along with the mesoscale oceanic eddies and fronts resolved during the high-resolution-SST period, the low-level storm track strengthens northward, reaching more than 30% of the maximum values in the low-resolution-SST period after removing the influence of ENSO. The mesoscale structure firstly imprints on the marine atmospheric boundary layer, which then leads to changes in turbulent heat flux and near-surface convergence, forcing a secondary circulation into the free atmosphere, strengthening the vertical eddy heat, momentum and specific humidity fluxes, and contributing to the enhancement of the storm track. Results from a high-resolution atmospheric model further indicate the changes in the storm track due to the mesoscale SST and their relationship.

    Key words: storm track, mesoscale SST, air?sea interaction, ERA-Interim, CAM4

    1. Introduction

    The European Centre for Medium-Range Weather Forecasts (ECMWF) interim reanalysis (ERA-Interim) dataset is typically applied to detect climate change (Chang and Yau,2016), assess interannual or decadal variability (O’Reilly and Czaja, 2015; Wills and Thompson, 2018), or evaluate against model results (Piazza et al., 2016; Ma et al., 2017;Lee et al., 2018). However, the prescribed sea surface temperature (SST) resolution for the ERA-Interim analysis system has been increased twice since 1979, resulting in a significant impact on the marine atmospheric boundary layer(MABL) (Masunaga et al., 2015, 2016; Parfitt et al., 2017).Masunaga et al. (2015) noted that mesoscale features, such as surface wind convergence and cloudiness, can only be seen during the period with high-resolution SST. Parfitt et al. (2017) showed the sensitivity to the SST resolution of the frontal air?sea interaction in the Gulf Stream region.The authors found that the occurrence of atmospheric cold fronts and precipitation increased up to 30% during the high-resolution period, along with strengthening of the surface storm track.

    The increased spatial resolution of the prescribed SST translates to better resolved mesoscale features of the ocean,especially in some eddy-rich regions such as the Kuroshio Extension area or Gulf Stream area. Different from the large-scale air?sea interaction, the atmosphere is forced by the ocean at the meso scale (Bishop et al., 2017). Recent studies have further suggested the important role of oceanic mesoscale eddies (Ma et al., 2015b, 2017) and fronts (Nakamura et al., 2008; Small et al., 2008; Kuwano-Yoshida and Minobe, 2017) on the storm track. Comparing twin experiments in a regional atmospheric model, one with realistic SST and another with spatial smoothing SST, Ma et al.(2017) showed that removing the mesoscale SST not only led to a reduction in the local storm track intensity, but also forced a significant southward shift in the storm track over the eastern North Pacific. Similar results were also found over the North Atlantic (Piazza et al., 2016). In addition,Woollings et al. (2010) showed that the storm track intensity in the model was overestimated along the coast of North America when the SST resolution was low, which was later confirmed by Small et al. (2014). A recent study by Small et al. (2019b) also pointed out that the storm track location improved in the coupled model with an eddy-resolving ocean model.

    According to the aforementioned studies, the mesoscale SST signals resolved by high-resolution SST could play an important role in the intensity and location of the storm track. Therefore, the increases in SST resolution in ERA-Interim may be able to change the storm track characteristics and introduce artificial signals. Meanwhile, the storm track, which is known to serve as a bridge of air?sea interaction at midlatitudes, could greatly influence weather and climate processes. Such an artificial effect would be very likely to introduce uncertainties and mislead us into reaching inaccurate or wrong conclusions. For example, if we ignore this effect, the long-term trend may be overestimated or underestimated, especially considering global warming.

    In the present study, we attempt to identify and quantify the signals of the storm track due to the increase in the prescribed SST resolution in ERA-Interim. How does its pattern change due to the increased SST resolution? The results will help us to clarify the uncertainties of the ERAInterim dataset in terms of both the spatial pattern and temporal variability of the storm track in the North Pacific, as well as benefit studies of climate change.

    2. Data and methods

    ERA-Interim is available on a 0.75° grid for the period January 1979 to December 2018 (Dee et al., 2011). We mainly use the daily fields, such as SST, sea level pressure(SLP), horizontal wind, specific humidity, and air temperature. In addition, the surface heat fluxes are taken from the forecasting system. The prescribed SST resolution in ERAInterim has increased twice, from 1.0° to 0.5° on 1 January 2002, and from 0.5° to 0.05° on 1 February 2009. The SST dataset with a resolution of 0.5° was from NCEP Real-Time Global SST, and then the Operational Sea Surface Temperature and Sea-Ice Analysis (OSTIA) from February 2009(Dee et al., 2011). Despite the oversampling of the OSTIA product, all of these SST datasets were interpolated to the same model grid, which was the same as the atmospheric model. The latter has remained at the same resolution from 1979 to the present day (ECMWF, 2019). It is also important to note that, except for the SST datasets, the atmospheric physical processes do not change. Here, the period of January 1979 to December 2001 is referred to as the lowresolution (LR) period, while the period of February 2009 to December 2018 is the high-resolution (HR) period.Although the prescribed horizontal resolution of ERAInterim is 0.75°, Masunaga et al. (2015) suggested that mesoscale features were well represented during the ERAInterim HR period, and were overall consistent with high-resolution satellite observations (see their Fig. 1).

    To separate the mesoscale signals, a 5° × 5° spatial boxcar filter was employed (Small et al., 2019a; Zhang et al.,2019). First, we calculated the area average in the 5° × 5°box for the original field (X) at every grid point. Then, the smoothed field was obtained, referred to as . Next, the mesoscale anomalies (X*) were defined as the original field minus the smoothed one. This procedure is denoted as a high-pass boxcar filter and can be expressed as X*= X ?, where < > represents the low-pass spatial filter. We examined our boxcar filter carefully and compared it with other methods, such as the Loess spatial filter used in Ma et al. (2017). The results demonstrated that the mesoscale signals with wavelength less than 500 km could be well separated by this boxcar filter (not shown).

    The storm track in the present study is represented by the meridional eddy heat flux (v′T′; Nakamura et al., 2002;Taguchi et al., 2009; Yao et al., 2018), the synoptic-scale variance in meridional wind (; Hoskins and Hodges, 2002),and eddy kinetic energy at 850 hPa. The prime here denotes the synoptic-scale eddies, which were derived by a 2?8-day Lanczos bandpass filter (Duchon, 1979). The focus of this paper is on midlatitude air?sea interaction, and so the influence of ENSO on the storm track was reduced by the removal of the linear regression on the Ni?o3.4 index (Qiu et al., 2014; Small et al., 2019b). This process can be expressed as follows:

    where X ( t) i s the atmospheric variable,is the ENSO signal, andstands for the part not influenced by ENSO.N, α, and b are the Ni?o3.4 series, regression coefficient,and residual, respectively. Czaja and Frankignoul (2002)showed that the covariance between the atmosphere and SST anomalies reached a maximum after two months in the observation, indicating that the atmospheric response may take a few months to fully develop (Ferreira and Frankignoul,2005). Following previous studies (Frankignoul et al., 2011;Gan and Wu, 2013; Révelard et al., 2016), two-month lagged regression was applied (τ = 2). Figure 1 shows an example of the storm track before and after linearly removing the ENSO impacts.

    To identify the response to the increase in SST resolution in ERA-Interim, we directly computed the difference in the climatological winter mean variables between the HR and LR periods (HR minus LR). To avoid the influences of interannual variabilities in the LR period, the differences were further computed following the method of Parfitt et al.(2017), in which the winter means of the randomly selected eight years during the LR period were also used, denoted as“LRi”, where “i” is the identifier of the different samples.Here, 1000 random samples during the LR period were selected to compute the differences. The average (AVE) and standard deviation (STD) of these differences can be used to measure the effects of the SST resolution and the interannual variability, respectively.

    Fig. 1. (a) Winter-mean storm track measured by the meridional eddy heat flux (v′T′; units: m s?1 K) in the LR period. (b) As in (a) but for the storm track with the influence of ENSO linearly removed.

    Two simulations using the high-resolution Community Atmosphere Model, version 4 (CAM4; Neale et al., 2013)were conducted to further confirm the impact of mesoscale SST on the storm track: one with the eddy-resolving prescribed SST, and the other with the smoothed SST, referred to as the control run (CTRL) and mesoscale-SST-filtered run (MSFR), respectively. In MSFR, the mesoscale SST was filtered out by the low-pass spatial boxcar filter with scales of 5° × 5°. CAM4 was configured on the finitevolume dynamical core with an approximately 25-km horizontal resolution and 26 vertical levels. A five-year integration with the “present-day” greenhouse gas conditions and the monthly c

    l

    imatological SST boundary condition (climatological mean around 2000) was first conducted as a spin-up run. The results of the final day of the spin-up integration were then used as the initial condition of CTRL and MSFR.Meanwhile, the SST boundary condition was taken from a high-resolution coupled model output with a horizontal resolution of 0.1° × 0.1° (Lin et al., 2019). As the horizontal resolution of the daily SST is higher than that of CAM4, it was then interpolated to the CAM4 grid based on the areal conservative remapping method. Next, the remapped daily SST was used to force the atmospheric model. The two simulations were run for six model-years.

    3. Results

    3.1. Changes in SST and its impact on the storm track

    Fig. 2. (a) Differences (colors; unit: °C) in the wintertime climatology of the SST between the HR and LR periods(HR minus LR) and the climatological SST (CI = 3°C, CI means contour interval.) during the HR period. (b)Average (colors; unit: °C) and STD (CI = 0.05°C) of the differences between the HR and LRi periods. (c)Climatological mesoscale SST (unit: °C) during the HR period. (d) As in (c) but for the LR period. The red box defines the KOCR.

    Figure 2a shows the differences in the climatological winter mean SST between the HR and LR periods, which shows significant positive and negative anomalies across the Kuroshio and Oyashio confluence region (KOCR;33°?45°N, 142°?175°E), where the oceanic mesoscale eddies are active (Chelton et al., 2011). The maximum and minimum of the differences over KOCR are 1.96°C and?1.93°C, respectively, while the positive and negative anomalies over the eastern North Pacific between 180°E and 130°W can only reach values of about 0.91°C and ?0.66°C.

    The average of the differences between the HR and each LRi period is shown in Fig. 2b, which is similar to that in Fig. 1a, with the pattern correlation coefficient reaching 0.99 for the whole region. In addition, the maximum and minimum of the mean values (shaded) over KOCR are 2.01°C and ?1.93°C, which are also very close to those in Fig. 2a.This result indicates that the pattern of SST differences is barely affected by the period we selected during LR, and the effects of interannual variability in the LR period are negligible.

    The contours in Fig. 2b show the STD of the differences between the HR and LRi periods. The large values occur in KOCR, indicating a relatively strong natural variability there. Moreover, the value of the maximum STD is only 0.32°C, which is less than the values of the differences in Fig. 2a (maximum of 1.96°C). This result confirms the robustness of the differences against interannual variability. The results of more random samples also show a similar pattern(not shown).

    To separate the mesoscale SST, we applied the highpass 5° × 5° spatial boxcar filter mentioned above to the SST field in both the LR and HR periods, the results of which are shown in Figs. 2c and d, respectively. The highpass-filtered SST anomalies are confined within KOCR,while anomalies of SST over the eastern North Pacific disappear. The pattern of SST differences in KOCR in Fig. 2a is similar to the mesoscale SST during the HR period (Fig. 2c),with the pattern correlation coefficient reaching 0.71. The maximum and minimum of mesoscale SST during the HR period are 1.90°C and ?1.83°C, respectively—only slightly less than the SST differences between the HR and LR period. However, there is almost no mesoscale signal in the high-pass-filtered SST during the LR period (Fig. 2d), with the positive and negative values reaching 0.58°C and?0.67°C in KOCR, respectively. Therefore, we speculate that the SST differences between the HR and LR period over KOCR arose from the increase in the SST horizontal resolution and are mainly dominated by mesoscale SST, while the anomalies in the eastern North Pacific are at a large scale and probably induced by the climate variabilities between the HR and LR periods, such as the Pacific Decadal Oscillation and global warming.

    Compared to that of the LR period, the pattern of the storm track during the HR period, measured by the meridional eddy heat flux () at 850 hPa, shows a northward strengthening. The differences of the storm track show an enhancement over KOCR and extend northeast to the Aleutian Islands (Fig. 3a), while it decreases over the northwestern North Pacific as well as the Gulf of Alaska. The maximum and minimum of the differences in the meridional eddy heat flux are 1.16 m s?1K and ?1.12 m s?1K, respectively, while the maximum of the climatological mean storm track (ENSO influence linearly removed) is 2.95 m s?1K,nearly 30%?40% of the maximum climatology. In addition,the maximum and minimum of the AVE values (shaded)over the North Pacific are 1.15 m s?1K and ?1.12 m s?1K,which are very close to those in Fig. 3a. However, the value of the maximum STD is 0.93 m s?1K, which is comparable to the AVE maximum. This result suggests a subtle influence of the increased SST resolution on the low-level storm track.

    3.2. Local response

    Masunaga et al. (2015) suggested that the mesoscale SST anomalies could imprint on the MABL. Figure 4 shows the mesoscale SST and surface heat flux differences between the HR and LR periods. The turbulent heat flux anomalies (Fig. 4c), defined as the sum of latent heat flux and sensible heat flux, are in phase with the mesoscale SST anomalies, with the pattern correlation coefficient reaching 0.89 over KOCR. Indeed, the changes of turbulent heat flux are dominated by the latent heat flux, accounting for more than 50% over KOCR (Fig. 4). This result suggests that more (less) heat fluxes out of the ocean over the warm(cold) mesoscale SST anomalies. Small et al. (2019a)revealed that, in midlatitude ocean frontal zones, the variability of surface heat fluxes is driven by small-scale SST on monthly time scales. The ocean forcing of heat flux further leads to the changes in the MABL thermal structure, represented by the in-phase variation of the meridional SST gradient and surface temperature gradient, with the pattern correla-tion coefficient reaching 0.88 (Fig. 4d). Changes can then be found in near-surface convergence/divergence, as described below.

    Two mechanisms are usually applied to explain the local atmospheric response to frontal-scale or mesoscale SST. One is the “pressure adjustment mechanism” (PAM;Lindzen and Nigam, 1987), which suggests in-phase spatial coherence among the sign-reversed SST Laplacian, the Laplacian of SLP, and near-surface wind convergence. Minobe et al. (2008) showed that the wind speed convergence is proportional to the SLP Laplacian as follows:

    The other mechanism is vertical momentum mixing(VMM; Wallace et al., 1989), which attributes the changes of surface wind speed to atmospheric instability. The strength can be measured by the linear relationship between wind stress divergence and downwind SST gradients (Ma et al., 2015a; Piazza et al., 2016). Figure 5c shows the differences of downwind SST gradients and wind stress divergence. Here, the downwind SST gradient is calculated by, where V is the 10-m wind vector and k the wind direction. As shown, the shape of the downwind SST gradient is very similar to the wind stress divergence, with a correlation coefficient of 0.79, which demonstrates the VMM. In fact, the relative contributions of the two mechanisms has long been debated (Small et al., 2008).Using a high-resolution atmospheric general circulation model, Koseki and Watanabe (2010) suggested that these two mechanisms make almost equal contributions in the Kuroshio Extension in January. In contrast, the results reported by Chen et al. (2017) showed VMM to be dominant, occupying 60% in the winter Kuroshio Extension region. On the other hand, considering the different atmospheric background wind direction, Schneider and Qiu (2015) and Bai et al. (2019) indicated that the response to SST in the MABL is quietly different, suggesting that the relative importance of PAM and VMM should be studied under different wind conditions.

    The imprints on the MABL suggest that the improved prescribed SST resolution does exert a significant influence on the near-surface atmosphere. But how does the impact penetrate into the free atmosphere and force the storm track in the troposphere? Figure 6a plots the cross-section of vertical motion and convergence at pressure levels along 42°N,and Fig. 6c shows the profiles of the SLP Laplacian and mesoscale SST at the same latitude. It is clear that an upward(downward) motion is induced over the convergence (divergence), forcing a secondary circulation above the oceanic eddies. Moreover, the vertical velocity anomalies are not confined within the MABL and could penetrate into the free atmosphere (above 700 hPa). Also, Fig. 6a implies that the VMM can be operative since the convergence and divergence patches straddle the peaks of high-pass-filtered SST.

    Fig. 5. Differences (colors) between the HR and LR periods for (a) the Laplacian of SST (units: 10?7 K m?2), (b) nearsurface convergence (units: 10?7 s?1), and (c) the downwind SST gradient [units: 10?6°C (100 km)?1]. The contours represent the Laplacian of SLP (CI = 8 × 10?7 Pa m?2) in (a, b)and wind stress divergence (CI = 0.1 N m?3) in (c). For clarity,the zero contour is omitted in all plots.

    As a consequence, the vertical eddy fluxes strengthen over KOCR (Fig. 6b). Cai et al. (2007) suggested that the vertical eddy heat flux is closely related with the baroclinic conversion (BC) between the eddy available potential energy and EKE, which is expressed as follows:C1=(P0/P)Cv/Cp(R/g)

    whereωis vertical velocity and . R and Cp(Cv) represent the gas constant for dry air and the specific heat capacity of dry air at constant pressure (volume),respectively. Figure 6d shows the difference in BC at 850 hPa between the HR and LR periods. The positive anomaly east of Japan further indicates the strengthening of the storm track during the HR period. Besides, the enhancement of vertical eddy specific humidity flux over KOCR, which is associated with latent heat fluxes, can also significantly wetten the atmosphere. As shown in Ma et al. (2017), the moist atmosphere could reduce the static stability and strengthen the baroclinicity in the atmosphere, ultimately leading to the enhanced storm track.

    3.3. Response of the storm track in model results

    Due to the method for calculating the differences in the storm track, changes may be induced by both the influence of the increase in SST resolution and/or the long-term climate variability between the HR and LR periods, most likely on the decadal time scale. As shown in Fig. 2, increasing the SST resolution results in the differences dominating at the meso scale over KOCR. To investigate the role of mesoscale SST in the storm track and to exclude the effects of the climate variability on the decadal time scale, we conducted two experiments using a high-resolution model, CAM4.One was forced by the eddy-resolving prescribed SST and the other by the smoothed SST, denoted as CTRL and MSFR, respectively. The potential differences between the two simulations suggest the influence of mesoscale SST without the long-term climate variability. However, these mesoscale SSTs may be induced by oceanic eddies and fronts. To determine the components of these mesoscale SSTs, the meridional SST gradient was then calculated to represent the oceanic fronts. The climatological mean SST gradient magnitudes for the two simulations were averaged within KOCR to provide a single representative value. The magnitude of the oceanic front decreased only 4% in MSFR compared to that in CTRL. We suspect that the mesoscale structure of the differences in the SST in our case may be mainly induced by the oceanic eddies, rather than the SST front. The influence of SST fronts is beyond the scope of this paper, but is worthy of further investigation in the future.

    Therefore, the response of the storm track to mesoscale SST was examined. Similarly, the influence of ENSO on the storm track was removed. Here, the ENSO signal is defined by the first principal component of monthly SST anomalies in the tropical Pacific between 12.5°N and 12.5°S. Figure 7e shows the differences of winter mean SST in CTRL and MSFR. Across KOCR, the mesoscale SST anomalies dominate, which agrees with the situation during the HR period in ERA-Interim (Fig. 2c). The differences of meridional eddy heat flux between CTRL and MSFR show that the storm track is significantly enhanced across KOCR(Fig. 7f). Moreover, the positive anomalies extend northeast from KOCR to the Gulf of Alaska, while the negative anomalies occur over the eastern and northwestern North Pacific, which shows a similar pattern to the differences in Figs. 3a and c. These results further indicate that the differences in the storm track between the HR and LR periods can be induced by the increase in the prescribed SST resolution in ERA-Interim. However, it should be noted that the specific location and value of the response are not closely consistent with the observation. This may be explained by the fact that the climatology of the storm track in the model is not identical to that in the observation. Also, the shorter period of model results may be another contributory factor for the inconsistency.

    To quantify the magnitude of the changes in storm track, a metric following Small et al. (2014) is introduced,(unit: %), which is expressed as the percentage of the maximum absolute difference between the HR and LR periods to the maximum value of the storm track during the LR period. The S and (ˉ) stand for the variables representing the storm track and the temporal mean during winter. The M is calculated over the domain of the North Pacific (30°?60°N, 140°?240°E). The M value of the meridional eddy heat flux and the variance of the synoptic-scale meridional wind is 39% and 32%,respectively. Then, we calculated the value of M for CAM4,which was roughly 42% of the maximum value in MSFR.

    4. Conclusion and discussion

    This study investigated the artificial changes of storm track induced by the increase in the SST resolution in the ERA-Interim dataset over the North Pacific. Along with the signals of resolved mesoscale oceanic eddies and fronts, the storm track strengthens northward across KOCR during the HR period. The positive anomalies of meridional eddy heat flux or meridional momentum flux extend from east of Japan to the northeastern North Pacific. The differences of storm track between the HR and LR periods reach more than 30%?40% of the maximum values during the LR period after removing the influence of ENSO. In addition, a set of experiments using a high-resolution atmospheric model were used to further examine whether the changes of the North Pacific storm track can be due to the resolved mesoscale SST. The results also indicated that the changes in SST horizontal resolution possibly affect the long-term variab-ility amplitude of the storm track. Therefore, the artificial influence induced by the changes in SST resolution should be considered when analyzing the long-term variability using the ERA-Interim dataset.

    Fig. 7. (a) Winter mean SST (unit: °C) in CTRL. (c) Winter mean SST (units: °C) in MSFR. (e) Differences (colors)of SST between CTRL and MSFR, with the winter mean of SST in MSFR overlaid. (b, d, f) As in (a, c, e) but for the storm track represented by the transient eddy heat flux at 850 hPa (units: m s?1 K). Statistically significant differences at the 95% confidence level according to the bootstrapping test are stippled.

    The mechanism by which the mesoscale SST impacts the storm track can be summarized as follows. The mesoscale imprints were first found in the turbulent heat flux with more (less) heat fluxed out of the ocean over the warm(cold) SST. That will further change the thermal structure in the MABL, and lead to anomalous convergence/divergence at the surface, forcing a secondary circulation above the mesoscale SST anomalies over KOCR. The PAM and VMM are both operative to the response of near-surface convergence.However, our results from reanalysis data could not fully separate the respective contributions of the PAM and VMM.The impact can penetrate into the free atmosphere, enhance the vertical eddy heat, momentum and specific humidity fluxes, and finally strengthen the storm track intensity in the free atmosphere both locally and remotely.

    Note also there is a large-scale positive anomaly of SST in the subtropical North Pacific. Wang et al. (2017) showed that when there is a North Pacific Gyre Oscillation—like positive SST anomaly in the subtropical North Pacific, the intensity of the subarctic frontal zone (SAFZ) gets stronger. This stronger SAFZ tends to lead a northward shift of the baroclinic zone and westerly jet, which also indicates the northward strengthening of the storm track. Recently, studies have also revealed that an intensified subtropical frontal zone (STFZ) would lead to stronger transient eddy activity(Wang et al., 2017, 2019; Chen et al., 2019; Wen et al.,2020). The impact of the STFZ on the intensity of the storm track seems comparable with that to the SAFZ (Yao et al.,2016; Wang et al., 2019). Hence, the relative contribution of the STFZ to the storm track should also be studied in the future. Indeed, based on our model experiments, we have confirmed that only the mesoscale SST anomalies can lead to a similar pattern of changes in the storm track, which indicates that the influences of large-scale SST anomalies may not be as important as the mesoscale SST in this case.Besides, the numerical experiments of Sun et al. (2018) also illustrated that oceanic stochastic forcing due to small-scale SST variability is crucial to the atmospheric circulation above, which can organize well the atmospheric transient eddies and improve the simulated storm track.

    In addition, the differences shown in the model experiments might be due to internal climate variability. Furthermore, the atmosphere is highly chaotic and variable in the midlatitudes, and our limited samples may not be sufficient to capture the internal climate variability. However, previous studies have highlighted oceanic mesoscale structures other than internal climate variability that result in changes of the storm track. Ma et al. (2017) conducted twin ensemble simulations using a regional model, each of which contained 10 members. Their results also showed a similar result to ours. They found a significant weakening of the storm track across KOCR after removing the mesoscale SST, indicating a robust influence of mesoscale SST on the atmosphere rather than the climate variability (see their Fig.5a). In addition, more ensemble simulations using an atmospheric model coupled with a slab-ocean model indicated the ocean mesoscale feedback to the enhancement of the storm track (Jia et al., 2019). Small et al. (2019b) also demonstrated the important role of ocean mesoscale structures on the extratropical storm track using a global coupled model.Based on these valuable works, we focused mainly on the impact of mesoscale SST. Nevertheless, we are planning to investigate the effect of model internal variability in our future work.

    Overall, our analysis and experiments, to some extent,show evidence and an indication that the changes in storm track can be due to the resolved mesoscale SST. However,more experiments are needed to prove that part of the changes in the period mean pattern and the climate variability in the ERA-Interim dataset between LR and HR can be attributed to the resolved ocean mesoscale eddies.

    Acknowledgements. This study was supported by National Key R&D Program for Developing Basic Sciences(2018YFA0605703, 2016YFC1401401), the National Natural Science Foundation of China (Grant Nos. 41490642, 41776030,41806034, 4160501) and the research project of the National University of Defense Technology (ZK20-45 and ZK17-02-010). HL and PL acknowledge the technical support from the National Key Scientific and Technological Infrastructure project “Earth System Science Numerical Simulator Facility” (EarthLab). The ERAInterim data were obtained from the ECMWF at http://apps.ecmwf.int/datasets/data/interim-full-daily/.

    女人精品久久久久毛片| 午夜福利,免费看| 日韩欧美精品免费久久| 人妻系列 视频| 一二三四在线观看免费中文在| 国产精品久久久久久久久免| 国产精品久久久久久久久免| 少妇精品久久久久久久| 午夜福利视频在线观看免费| 国产精品不卡视频一区二区| 亚洲在久久综合| 最近中文字幕2019免费版| 国产黄频视频在线观看| 亚洲欧美精品综合一区二区三区 | 男女国产视频网站| 啦啦啦在线观看免费高清www| 自拍欧美九色日韩亚洲蝌蚪91| 色哟哟·www| 自拍欧美九色日韩亚洲蝌蚪91| 男女高潮啪啪啪动态图| 自拍欧美九色日韩亚洲蝌蚪91| 国产精品欧美亚洲77777| 9色porny在线观看| 久久午夜福利片| 一二三四在线观看免费中文在| 亚洲精品中文字幕在线视频| 亚洲美女搞黄在线观看| 青春草视频在线免费观看| 国产又色又爽无遮挡免| 国产精品不卡视频一区二区| 成人亚洲欧美一区二区av| 亚洲精品自拍成人| 久久久久精品久久久久真实原创| 一本色道久久久久久精品综合| 午夜久久久在线观看| 综合色丁香网| 久久精品国产综合久久久| 亚洲精品第二区| 亚洲国产成人一精品久久久| 免费少妇av软件| 尾随美女入室| 国产精品av久久久久免费| 18在线观看网站| 日产精品乱码卡一卡2卡三| 久久久亚洲精品成人影院| 天堂俺去俺来也www色官网| 热re99久久国产66热| 亚洲综合精品二区| 在现免费观看毛片| 天天躁日日躁夜夜躁夜夜| 欧美人与性动交α欧美软件| 91精品国产国语对白视频| 97人妻天天添夜夜摸| 久久精品夜色国产| 黄色 视频免费看| 亚洲国产精品国产精品| 男人添女人高潮全过程视频| 人妻一区二区av| 夫妻性生交免费视频一级片| 母亲3免费完整高清在线观看 | 亚洲国产av新网站| www.精华液| 两个人看的免费小视频| 黄色怎么调成土黄色| 国产高清不卡午夜福利| tube8黄色片| 精品酒店卫生间| 亚洲久久久国产精品| 久久女婷五月综合色啪小说| 人妻一区二区av| av.在线天堂| 汤姆久久久久久久影院中文字幕| 极品人妻少妇av视频| 制服人妻中文乱码| 天天躁日日躁夜夜躁夜夜| 国产乱来视频区| av线在线观看网站| 国产精品国产av在线观看| 国产一区二区激情短视频 | 亚洲精品,欧美精品| 亚洲一级一片aⅴ在线观看| 久久综合国产亚洲精品| 国产精品久久久久久精品古装| 十八禁网站网址无遮挡| 视频区图区小说| 好男人视频免费观看在线| 亚洲av在线观看美女高潮| 欧美日韩综合久久久久久| av免费在线看不卡| 欧美人与性动交α欧美软件| 精品视频人人做人人爽| 亚洲av中文av极速乱| 成人影院久久| 999精品在线视频| 看免费成人av毛片| 免费在线观看黄色视频的| videossex国产| 国产精品麻豆人妻色哟哟久久| 亚洲欧美精品自产自拍| 欧美中文综合在线视频| 亚洲精品成人av观看孕妇| 亚洲,欧美,日韩| 日本色播在线视频| 黄色毛片三级朝国网站| 亚洲av在线观看美女高潮| 免费观看av网站的网址| 欧美国产精品一级二级三级| 亚洲精品第二区| 在线精品无人区一区二区三| 最黄视频免费看| 热re99久久国产66热| 午夜福利视频在线观看免费| 天天操日日干夜夜撸| 黑人猛操日本美女一级片| 一级黄片播放器| 国产精品二区激情视频| 精品国产乱码久久久久久小说| 好男人视频免费观看在线| 国产精品久久久av美女十八| 欧美精品人与动牲交sv欧美| 在线观看国产h片| 999精品在线视频| 精品视频人人做人人爽| 纯流量卡能插随身wifi吗| av国产精品久久久久影院| 午夜老司机福利剧场| 纯流量卡能插随身wifi吗| 波多野结衣av一区二区av| 久久av网站| 亚洲在久久综合| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | av有码第一页| 超色免费av| 有码 亚洲区| 亚洲综合精品二区| 伊人久久大香线蕉亚洲五| 免费在线观看完整版高清| 美女大奶头黄色视频| 女人精品久久久久毛片| 亚洲欧美一区二区三区国产| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 丝袜脚勾引网站| 1024香蕉在线观看| av福利片在线| 欧美国产精品va在线观看不卡| 蜜桃在线观看..| 波多野结衣一区麻豆| 午夜老司机福利剧场| 久久这里有精品视频免费| 只有这里有精品99| 黄色毛片三级朝国网站| 精品一区二区三卡| av在线观看视频网站免费| 一二三四中文在线观看免费高清| 欧美老熟妇乱子伦牲交| 天天躁夜夜躁狠狠久久av| 久久热在线av| 亚洲伊人色综图| 亚洲精品aⅴ在线观看| 久久午夜综合久久蜜桃| 99久久人妻综合| 亚洲av日韩在线播放| 久久ye,这里只有精品| 亚洲成人一二三区av| 成人亚洲欧美一区二区av| 中文精品一卡2卡3卡4更新| 777久久人妻少妇嫩草av网站| 国产探花极品一区二区| 久久久久久伊人网av| 欧美 亚洲 国产 日韩一| 另类亚洲欧美激情| 观看av在线不卡| 又黄又粗又硬又大视频| 国产精品国产av在线观看| 久久精品国产自在天天线| 久久国内精品自在自线图片| 欧美黄色片欧美黄色片| 性少妇av在线| 水蜜桃什么品种好| 成年人免费黄色播放视频| 成人国语在线视频| 国产乱来视频区| 久久久久久久久久人人人人人人| 国产亚洲欧美精品永久| 国产片特级美女逼逼视频| 一区二区三区激情视频| 精品一区二区免费观看| 在线观看美女被高潮喷水网站| 晚上一个人看的免费电影| 欧美日韩成人在线一区二区| 午夜久久久在线观看| 一本大道久久a久久精品| 久久精品国产a三级三级三级| 午夜福利乱码中文字幕| 男的添女的下面高潮视频| 人人妻人人澡人人看| 黄色怎么调成土黄色| 免费黄网站久久成人精品| 久久久久精品人妻al黑| 极品人妻少妇av视频| av国产精品久久久久影院| 最近手机中文字幕大全| 亚洲国产成人一精品久久久| 精品人妻熟女毛片av久久网站| 免费高清在线观看视频在线观看| 国产视频首页在线观看| 亚洲欧美日韩另类电影网站| 一区二区日韩欧美中文字幕| 亚洲国产精品一区二区三区在线| 女人精品久久久久毛片| 咕卡用的链子| 国产在线一区二区三区精| 日韩熟女老妇一区二区性免费视频| 少妇被粗大猛烈的视频| 久久久久久久精品精品| 中文字幕亚洲精品专区| 亚洲经典国产精华液单| av免费在线看不卡| 欧美国产精品va在线观看不卡| 国产午夜精品一二区理论片| av在线观看视频网站免费| xxxhd国产人妻xxx| 飞空精品影院首页| 老汉色av国产亚洲站长工具| 电影成人av| 一区福利在线观看| 香蕉丝袜av| 久久精品国产亚洲av高清一级| 肉色欧美久久久久久久蜜桃| 国产高清不卡午夜福利| 国产 一区精品| 午夜福利网站1000一区二区三区| 欧美老熟妇乱子伦牲交| 国产一区二区在线观看av| 不卡av一区二区三区| 久热这里只有精品99| 99re6热这里在线精品视频| 老鸭窝网址在线观看| 久久综合国产亚洲精品| 伊人久久大香线蕉亚洲五| 免费高清在线观看日韩| 日日撸夜夜添| 免费观看性生交大片5| 亚洲国产精品国产精品| 母亲3免费完整高清在线观看 | 中文字幕人妻熟女乱码| 卡戴珊不雅视频在线播放| 女的被弄到高潮叫床怎么办| 日韩一区二区视频免费看| 免费人妻精品一区二区三区视频| av国产久精品久网站免费入址| 亚洲精品国产av蜜桃| 亚洲中文av在线| 亚洲美女视频黄频| 一本大道久久a久久精品| 亚洲视频免费观看视频| 两性夫妻黄色片| 人人妻人人添人人爽欧美一区卜| 国产精品麻豆人妻色哟哟久久| 国产高清不卡午夜福利| 一级毛片 在线播放| 国产精品久久久久成人av| 国产成人精品一,二区| 在线看a的网站| 色哟哟·www| 精品少妇一区二区三区视频日本电影 | 在线观看免费日韩欧美大片| 我的亚洲天堂| www.自偷自拍.com| 精品国产露脸久久av麻豆| 国产精品欧美亚洲77777| 亚洲一区中文字幕在线| 99久久人妻综合| 天堂中文最新版在线下载| 日韩av免费高清视频| 丰满饥渴人妻一区二区三| 五月开心婷婷网| 人人妻人人澡人人看| 国语对白做爰xxxⅹ性视频网站| 熟妇人妻不卡中文字幕| 亚洲一码二码三码区别大吗| xxx大片免费视频| 国产精品免费视频内射| 十八禁网站网址无遮挡| 最新中文字幕久久久久| 人妻人人澡人人爽人人| 国产一区有黄有色的免费视频| 婷婷色麻豆天堂久久| 久久97久久精品| 日韩精品免费视频一区二区三区| 亚洲精品,欧美精品| 人妻人人澡人人爽人人| 亚洲人成网站在线观看播放| 哪个播放器可以免费观看大片| 嫩草影院入口| 久久精品aⅴ一区二区三区四区 | 久久影院123| 婷婷成人精品国产| 国产成人精品在线电影| 美女国产视频在线观看| 男人爽女人下面视频在线观看| 尾随美女入室| 日本wwww免费看| 99久国产av精品国产电影| 欧美在线黄色| 男女下面插进去视频免费观看| 成年女人毛片免费观看观看9 | 日本欧美国产在线视频| 一区福利在线观看| 青草久久国产| 99九九在线精品视频| 国产免费又黄又爽又色| 婷婷色麻豆天堂久久| 亚洲精品av麻豆狂野| 久久久久网色| 80岁老熟妇乱子伦牲交| 26uuu在线亚洲综合色| av不卡在线播放| 99re6热这里在线精品视频| 日韩伦理黄色片| 一区二区三区乱码不卡18| 国产不卡av网站在线观看| 日韩伦理黄色片| 成年人午夜在线观看视频| 精品少妇黑人巨大在线播放| 99久久中文字幕三级久久日本| 秋霞在线观看毛片| 久久久久久久久久久久大奶| 日本欧美国产在线视频| 亚洲第一区二区三区不卡| 王馨瑶露胸无遮挡在线观看| 免费大片黄手机在线观看| 亚洲国产精品一区三区| 久久精品aⅴ一区二区三区四区 | 青春草亚洲视频在线观看| 国产女主播在线喷水免费视频网站| 久久 成人 亚洲| 日韩电影二区| 国产黄频视频在线观看| 亚洲欧美一区二区三区久久| 久久这里有精品视频免费| 国产一区二区三区av在线| 男男h啪啪无遮挡| 国产色婷婷99| 亚洲av中文av极速乱| 亚洲成人手机| 午夜影院在线不卡| 成人黄色视频免费在线看| 午夜日韩欧美国产| 久久这里有精品视频免费| 一区在线观看完整版| 啦啦啦中文免费视频观看日本| 亚洲伊人色综图| 制服丝袜香蕉在线| 黑人欧美特级aaaaaa片| 日日啪夜夜爽| 亚洲人成电影观看| 久久韩国三级中文字幕| 久久精品aⅴ一区二区三区四区 | av天堂久久9| 国产亚洲av片在线观看秒播厂| 汤姆久久久久久久影院中文字幕| 王馨瑶露胸无遮挡在线观看| 亚洲人成网站在线观看播放| 久久这里只有精品19| 人妻少妇偷人精品九色| 成人免费观看视频高清| 亚洲内射少妇av| 啦啦啦啦在线视频资源| 午夜福利网站1000一区二区三区| 亚洲国产色片| 麻豆av在线久日| 丝袜脚勾引网站| 一级片免费观看大全| 丝袜脚勾引网站| 电影成人av| 久久久久国产精品人妻一区二区| 久久久亚洲精品成人影院| 水蜜桃什么品种好| 最新的欧美精品一区二区| 一边摸一边做爽爽视频免费| 黄片无遮挡物在线观看| 婷婷色av中文字幕| 最黄视频免费看| 久久久久久人妻| 最新中文字幕久久久久| 高清在线视频一区二区三区| 十分钟在线观看高清视频www| 免费黄频网站在线观看国产| 啦啦啦在线观看免费高清www| 国产在线一区二区三区精| 色婷婷av一区二区三区视频| 精品亚洲成a人片在线观看| 日韩大片免费观看网站| 欧美精品高潮呻吟av久久| 少妇人妻久久综合中文| 国产野战对白在线观看| 最近最新中文字幕免费大全7| 在线观看人妻少妇| 久久久国产一区二区| 午夜免费观看性视频| 伦理电影免费视频| 欧美激情 高清一区二区三区| 久久人人爽av亚洲精品天堂| 国产一区二区激情短视频 | 日韩一本色道免费dvd| 久久人人爽人人片av| 国产黄频视频在线观看| 三上悠亚av全集在线观看| 日韩伦理黄色片| 精品亚洲乱码少妇综合久久| 亚洲国产欧美在线一区| 日韩一区二区视频免费看| 精品一区在线观看国产| 亚洲成av片中文字幕在线观看 | 蜜桃在线观看..| 国产成人av激情在线播放| 捣出白浆h1v1| 9191精品国产免费久久| 久久人妻熟女aⅴ| 下体分泌物呈黄色| 国产精品偷伦视频观看了| 欧美 日韩 精品 国产| 日本wwww免费看| 天天躁狠狠躁夜夜躁狠狠躁| 男女边吃奶边做爰视频| 人人妻人人爽人人添夜夜欢视频| 午夜福利视频精品| 你懂的网址亚洲精品在线观看| 夜夜骑夜夜射夜夜干| 国产福利在线免费观看视频| 久久久久国产精品人妻一区二区| 日本色播在线视频| 韩国精品一区二区三区| 亚洲国产欧美在线一区| 1024香蕉在线观看| 一二三四中文在线观看免费高清| 肉色欧美久久久久久久蜜桃| 久久精品夜色国产| 69精品国产乱码久久久| 国产乱人偷精品视频| 国产男人的电影天堂91| 欧美老熟妇乱子伦牲交| 国产免费视频播放在线视频| 天堂中文最新版在线下载| 一二三四在线观看免费中文在| 国产精品成人在线| 国产在线一区二区三区精| 大码成人一级视频| 久久亚洲国产成人精品v| 中文字幕av电影在线播放| 老司机影院毛片| 成年人午夜在线观看视频| 最新的欧美精品一区二区| 欧美 日韩 精品 国产| 七月丁香在线播放| 国产av一区二区精品久久| 久久人人爽av亚洲精品天堂| 女人久久www免费人成看片| 一区二区三区四区激情视频| 你懂的网址亚洲精品在线观看| 欧美日韩精品网址| 国产熟女午夜一区二区三区| 亚洲国产精品国产精品| 97在线视频观看| 国产精品二区激情视频| 精品卡一卡二卡四卡免费| 高清黄色对白视频在线免费看| 人妻系列 视频| 日韩精品有码人妻一区| 男人添女人高潮全过程视频| 9色porny在线观看| 成人亚洲欧美一区二区av| 久久久久精品性色| 两性夫妻黄色片| 伊人久久国产一区二区| 国产亚洲午夜精品一区二区久久| 欧美日本中文国产一区发布| 亚洲四区av| 国产极品天堂在线| 熟妇人妻不卡中文字幕| 日韩av免费高清视频| 中文字幕人妻丝袜一区二区 | 少妇精品久久久久久久| 欧美国产精品一级二级三级| 一区福利在线观看| 黄色怎么调成土黄色| 国产精品av久久久久免费| 一二三四在线观看免费中文在| 观看av在线不卡| 欧美av亚洲av综合av国产av | 蜜桃国产av成人99| 一本色道久久久久久精品综合| 久久久欧美国产精品| 91在线精品国自产拍蜜月| 免费观看在线日韩| 99热全是精品| 精品久久久久久电影网| av视频免费观看在线观看| 亚洲,一卡二卡三卡| 久久精品国产亚洲av高清一级| 久久久精品区二区三区| 色吧在线观看| 免费在线观看完整版高清| 亚洲国产最新在线播放| 男女边摸边吃奶| 久久人人97超碰香蕉20202| 久久久国产一区二区| 久久精品国产亚洲av涩爱| 又黄又粗又硬又大视频| 两性夫妻黄色片| 老熟女久久久| 国产片内射在线| 亚洲经典国产精华液单| 国产精品免费视频内射| 久久久精品国产亚洲av高清涩受| 久久婷婷青草| 在线观看免费高清a一片| 男女下面插进去视频免费观看| 国产一级毛片在线| 亚洲综合色网址| 欧美人与善性xxx| 国产在视频线精品| 欧美精品一区二区免费开放| 久久久久久久亚洲中文字幕| 黑人欧美特级aaaaaa片| av.在线天堂| 另类精品久久| 亚洲成国产人片在线观看| 欧美黄色片欧美黄色片| 一本大道久久a久久精品| 久久99蜜桃精品久久| 精品人妻一区二区三区麻豆| 成人漫画全彩无遮挡| 国产精品av久久久久免费| 青春草国产在线视频| 精品一区二区免费观看| 国产精品不卡视频一区二区| 亚洲精品国产av成人精品| 街头女战士在线观看网站| 免费黄频网站在线观看国产| 国产成人av激情在线播放| 制服丝袜香蕉在线| 久久久久久人人人人人| 一级片'在线观看视频| 黄色配什么色好看| 国产片特级美女逼逼视频| 最黄视频免费看| 久久这里有精品视频免费| 热re99久久国产66热| 亚洲av中文av极速乱| 国产日韩欧美亚洲二区| 黄片播放在线免费| 狠狠婷婷综合久久久久久88av| 天堂俺去俺来也www色官网| 少妇的丰满在线观看| 不卡视频在线观看欧美| 亚洲欧美中文字幕日韩二区| 又大又黄又爽视频免费| 午夜免费男女啪啪视频观看| 亚洲综合精品二区| 街头女战士在线观看网站| 免费少妇av软件| 一二三四中文在线观看免费高清| 亚洲激情五月婷婷啪啪| 亚洲精品一区蜜桃| 黄片播放在线免费| 国产成人精品无人区| 亚洲精品国产av蜜桃| 精品一区在线观看国产| 香蕉丝袜av| 亚洲国产欧美日韩在线播放| 色哟哟·www| 亚洲成人av在线免费| 啦啦啦啦在线视频资源| av网站在线播放免费| 久久精品人人爽人人爽视色| 国产日韩欧美视频二区| 精品久久蜜臀av无| 女人久久www免费人成看片| 精品少妇一区二区三区视频日本电影 | 国产极品天堂在线| 国产精品麻豆人妻色哟哟久久| 一区二区三区四区激情视频| 亚洲精品自拍成人| 久久久a久久爽久久v久久| 亚洲 欧美一区二区三区| 纵有疾风起免费观看全集完整版| 亚洲图色成人| 成人漫画全彩无遮挡| 高清黄色对白视频在线免费看| 九九爱精品视频在线观看| 成人18禁高潮啪啪吃奶动态图| 国产成人精品久久久久久| 999久久久国产精品视频| 中文字幕最新亚洲高清| 黑丝袜美女国产一区| 老汉色∧v一级毛片| 国产精品 欧美亚洲| 久久国产精品大桥未久av| 韩国高清视频一区二区三区| 欧美日韩亚洲国产一区二区在线观看 | 97人妻天天添夜夜摸| 亚洲综合色惰| 欧美日本中文国产一区发布| 在线免费观看不下载黄p国产| 各种免费的搞黄视频| 国产精品.久久久| 性色av一级| 午夜久久久在线观看| 精品国产乱码久久久久久男人| 18禁动态无遮挡网站| av天堂久久9| 国产女主播在线喷水免费视频网站| 国产乱人偷精品视频|