• <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/.

    亚洲性夜色夜夜综合| 一进一出抽搐动态| 亚洲av电影不卡..在线观看| 在线国产一区二区在线| 日韩人妻高清精品专区| 国产视频一区二区在线看| 日韩中字成人| 亚洲最大成人av| 国产精品久久视频播放| 欧美xxxx性猛交bbbb| 国产高清三级在线| 在线天堂最新版资源| avwww免费| 18禁在线播放成人免费| 麻豆av噜噜一区二区三区| 亚洲三级黄色毛片| 插逼视频在线观看| 精品一区二区免费观看| 精品久久久久久久末码| 精品99又大又爽又粗少妇毛片| 亚洲,欧美,日韩| 亚洲精品日韩在线中文字幕 | 国产一区亚洲一区在线观看| 嫩草影院入口| 99热这里只有精品一区| 亚洲成人av在线免费| 看非洲黑人一级黄片| 欧美性猛交╳xxx乱大交人| 国产人妻一区二区三区在| 国产 一区 欧美 日韩| 亚洲av二区三区四区| 18禁黄网站禁片免费观看直播| 狂野欧美激情性xxxx在线观看| 国产欧美日韩一区二区精品| 精品国内亚洲2022精品成人| 男人舔奶头视频| 哪里可以看免费的av片| 最近的中文字幕免费完整| 午夜久久久久精精品| 白带黄色成豆腐渣| 男女做爰动态图高潮gif福利片| 一区二区三区四区激情视频 | 天天一区二区日本电影三级| 久久婷婷人人爽人人干人人爱| av在线播放精品| 三级毛片av免费| 可以在线观看毛片的网站| 夜夜看夜夜爽夜夜摸| 日韩精品有码人妻一区| 亚洲高清免费不卡视频| 级片在线观看| 一个人看视频在线观看www免费| 亚洲精品一区av在线观看| 国产精品三级大全| 精品欧美国产一区二区三| 中文在线观看免费www的网站| 黄色一级大片看看| 亚洲一级一片aⅴ在线观看| 男女视频在线观看网站免费| 最近的中文字幕免费完整| 搡女人真爽免费视频火全软件 | 日韩欧美三级三区| 在线观看美女被高潮喷水网站| 久久午夜福利片| 秋霞在线观看毛片| 精品99又大又爽又粗少妇毛片| 久久久精品大字幕| 18+在线观看网站| 婷婷亚洲欧美| 欧美三级亚洲精品| 国产精品久久久久久亚洲av鲁大| 欧美bdsm另类| 亚洲熟妇中文字幕五十中出| 中国国产av一级| 精品久久久久久久末码| 国产日本99.免费观看| 日本在线视频免费播放| 观看免费一级毛片| 日日摸夜夜添夜夜添av毛片| 三级男女做爰猛烈吃奶摸视频| 亚洲激情五月婷婷啪啪| 日日撸夜夜添| 99在线人妻在线中文字幕| 日韩欧美在线乱码| 亚洲av五月六月丁香网| 成年女人看的毛片在线观看| 日日啪夜夜撸| 亚洲一级一片aⅴ在线观看| 免费观看人在逋| 日韩高清综合在线| 99热6这里只有精品| 亚洲精品乱码久久久v下载方式| 久久久精品大字幕| 看免费成人av毛片| 亚洲精品日韩在线中文字幕 | 亚洲美女黄片视频| 久久久久久久久久黄片| 亚洲精品国产av成人精品 | 精品久久久久久久久av| 麻豆av噜噜一区二区三区| 亚洲欧美日韩东京热| 亚洲国产精品sss在线观看| 69人妻影院| 免费黄网站久久成人精品| 网址你懂的国产日韩在线| 久久久久久伊人网av| 国模一区二区三区四区视频| 小蜜桃在线观看免费完整版高清| 日本与韩国留学比较| 日日干狠狠操夜夜爽| 校园春色视频在线观看| 久久精品91蜜桃| 国产成年人精品一区二区| 麻豆国产av国片精品| 国产伦精品一区二区三区四那| 日韩av不卡免费在线播放| 欧美一区二区亚洲| 一级av片app| 亚洲av二区三区四区| 99热只有精品国产| 在线看三级毛片| 特级一级黄色大片| 久久久久久伊人网av| 99在线视频只有这里精品首页| 看片在线看免费视频| 国产伦精品一区二区三区视频9| 中国美女看黄片| 国产精品一区二区免费欧美| 97碰自拍视频| 99九九线精品视频在线观看视频| 色av中文字幕| 欧美三级亚洲精品| 听说在线观看完整版免费高清| 晚上一个人看的免费电影| 99久久九九国产精品国产免费| 91av网一区二区| 亚洲精品影视一区二区三区av| 欧美国产日韩亚洲一区| 欧美一区二区精品小视频在线| 听说在线观看完整版免费高清| 国产精品久久久久久久久免| 变态另类成人亚洲欧美熟女| 亚洲精品国产成人久久av| aaaaa片日本免费| 久久精品国产亚洲av涩爱 | 久久韩国三级中文字幕| 久久久精品大字幕| 日本免费一区二区三区高清不卡| 免费观看人在逋| 免费黄网站久久成人精品| 99视频精品全部免费 在线| 久久精品影院6| 免费在线观看成人毛片| 三级经典国产精品| 丰满乱子伦码专区| 欧美国产日韩亚洲一区| 蜜臀久久99精品久久宅男| 亚洲av第一区精品v没综合| av黄色大香蕉| 97碰自拍视频| 俺也久久电影网| 少妇的逼水好多| 亚洲精品成人久久久久久| 欧美成人一区二区免费高清观看| 波多野结衣巨乳人妻| 欧美又色又爽又黄视频| 精品久久久久久久末码| 久久久久免费精品人妻一区二区| 国产精品免费一区二区三区在线| 国产日本99.免费观看| 国产精品国产高清国产av| 国产精品人妻久久久影院| 美女 人体艺术 gogo| 久久精品国产亚洲av涩爱 | 精品久久久久久久末码| 赤兔流量卡办理| 婷婷色综合大香蕉| 国产真实乱freesex| 久久亚洲国产成人精品v| 亚洲欧美成人综合另类久久久 | 成人午夜高清在线视频| 啦啦啦观看免费观看视频高清| 国产乱人视频| 亚洲精品影视一区二区三区av| 一级毛片aaaaaa免费看小| 人人妻人人澡欧美一区二区| 亚洲婷婷狠狠爱综合网| 成人永久免费在线观看视频| 亚洲内射少妇av| 蜜桃亚洲精品一区二区三区| 成人av一区二区三区在线看| 嫩草影院入口| 成人精品一区二区免费| 日韩欧美一区二区三区在线观看| 性欧美人与动物交配| 插逼视频在线观看| 少妇裸体淫交视频免费看高清| 久久精品国产鲁丝片午夜精品| 欧美一区二区精品小视频在线| av视频在线观看入口| 亚洲人成网站高清观看| 国产不卡一卡二| 精品一区二区三区视频在线| 国内精品一区二区在线观看| 我要看日韩黄色一级片| 欧美一级a爱片免费观看看| 久久国产乱子免费精品| 毛片女人毛片| 99久久九九国产精品国产免费| 有码 亚洲区| 99热这里只有精品一区| 变态另类成人亚洲欧美熟女| 国产高清激情床上av| 免费看美女性在线毛片视频| 狂野欧美白嫩少妇大欣赏| 成人综合一区亚洲| 日本与韩国留学比较| 亚洲精品国产av成人精品 | 天天躁日日操中文字幕| 97超碰精品成人国产| 久久久久久伊人网av| 国产视频内射| 成人毛片a级毛片在线播放| 日日干狠狠操夜夜爽| 一级毛片aaaaaa免费看小| 一a级毛片在线观看| 亚洲国产精品成人久久小说 | 日韩制服骚丝袜av| 久久久国产成人精品二区| 国产一级毛片七仙女欲春2| 搡老妇女老女人老熟妇| 免费观看精品视频网站| 久久中文看片网| 最新在线观看一区二区三区| 99久久久亚洲精品蜜臀av| 两个人视频免费观看高清| 亚洲三级黄色毛片| 国产精品免费一区二区三区在线| 亚洲国产精品久久男人天堂| 亚洲成人精品中文字幕电影| 久99久视频精品免费| 亚洲精品一卡2卡三卡4卡5卡| 激情 狠狠 欧美| 中文在线观看免费www的网站| 最好的美女福利视频网| 国产淫片久久久久久久久| 黄色一级大片看看| 校园人妻丝袜中文字幕| 成人综合一区亚洲| 免费不卡的大黄色大毛片视频在线观看 | 少妇的逼水好多| 亚洲综合色惰| 内地一区二区视频在线| 日本与韩国留学比较| 精品一区二区三区视频在线| 亚洲第一区二区三区不卡| 69av精品久久久久久| 看片在线看免费视频| 国产精品一及| 女人被狂操c到高潮| 成人毛片a级毛片在线播放| 夜夜爽天天搞| 熟妇人妻久久中文字幕3abv| 亚洲精品国产成人久久av| 丝袜喷水一区| 婷婷六月久久综合丁香| 精品日产1卡2卡| 久久精品人妻少妇| 一本久久中文字幕| 人人妻人人澡人人爽人人夜夜 | 天堂网av新在线| 精品久久久久久久久亚洲| 国产成人aa在线观看| 十八禁国产超污无遮挡网站| www日本黄色视频网| 最近2019中文字幕mv第一页| 国产午夜福利久久久久久| 岛国在线免费视频观看| 成人精品一区二区免费| 精品久久国产蜜桃| 国产在线精品亚洲第一网站| 日韩欧美一区二区三区在线观看| 欧美成人精品欧美一级黄| 神马国产精品三级电影在线观看| 国产片特级美女逼逼视频| 在线观看66精品国产| 亚洲七黄色美女视频| 国产亚洲精品综合一区在线观看| 国产亚洲91精品色在线| 午夜精品一区二区三区免费看| 男人舔女人下体高潮全视频| 免费看av在线观看网站| 亚洲精品一区av在线观看| 午夜爱爱视频在线播放| 欧美区成人在线视频| 欧美激情在线99| 国产精品人妻久久久久久| 国产成年人精品一区二区| 日本在线视频免费播放| 国产精品av视频在线免费观看| 国产精品一二三区在线看| 1000部很黄的大片| 精品一区二区三区视频在线| 国产免费一级a男人的天堂| 天天躁日日操中文字幕| 一个人看的www免费观看视频| 美女xxoo啪啪120秒动态图| 国产伦精品一区二区三区四那| 久久久久久久久久黄片| 久久久久久久亚洲中文字幕| 在线免费观看的www视频| 内地一区二区视频在线| 91av网一区二区| 麻豆精品久久久久久蜜桃| 精品人妻偷拍中文字幕| 国产人妻一区二区三区在| 永久网站在线| 欧美色视频一区免费| 国产亚洲精品久久久久久毛片| 久久韩国三级中文字幕| 亚洲经典国产精华液单| 少妇猛男粗大的猛烈进出视频 | 最新中文字幕久久久久| 99久国产av精品| 亚洲乱码一区二区免费版| 熟女人妻精品中文字幕| 精品人妻偷拍中文字幕| 热99re8久久精品国产| 91久久精品电影网| 国产真实乱freesex| 日本黄大片高清| 俺也久久电影网| 在线观看66精品国产| 丝袜美腿在线中文| 婷婷精品国产亚洲av| 亚洲七黄色美女视频| 欧美激情国产日韩精品一区| 日本a在线网址| 国产伦在线观看视频一区| a级毛片免费高清观看在线播放| 啦啦啦韩国在线观看视频| 精品午夜福利在线看| 男女那种视频在线观看| 精品午夜福利在线看| 又黄又爽又刺激的免费视频.| 菩萨蛮人人尽说江南好唐韦庄 | 一边摸一边抽搐一进一小说| 国产熟女欧美一区二区| 亚洲最大成人手机在线| 12—13女人毛片做爰片一| 国产精品嫩草影院av在线观看| 国产视频一区二区在线看| 亚洲精品久久国产高清桃花| 三级经典国产精品| 在线观看午夜福利视频| 99视频精品全部免费 在线| 欧美xxxx黑人xx丫x性爽| 国产毛片a区久久久久| 一区二区三区四区激情视频 | 日韩精品有码人妻一区| 校园春色视频在线观看| 少妇丰满av| 老司机福利观看| 久久人人精品亚洲av| 少妇人妻精品综合一区二区 | av中文乱码字幕在线| 简卡轻食公司| 亚洲激情五月婷婷啪啪| 色视频www国产| 99国产极品粉嫩在线观看| 久久精品国产亚洲网站| 美女大奶头视频| 国产精品久久久久久久久免| 黄片wwwwww| 如何舔出高潮| 少妇丰满av| 亚洲美女搞黄在线观看 | 成人欧美大片| 丰满乱子伦码专区| 亚洲性久久影院| 黄色视频,在线免费观看| 色5月婷婷丁香| 观看美女的网站| 亚洲欧美日韩东京热| 国产高清有码在线观看视频| 色综合色国产| 中国美女看黄片| 插逼视频在线观看| 麻豆成人午夜福利视频| 老熟妇仑乱视频hdxx| 亚洲国产色片| av天堂在线播放| 一级黄片播放器| 波多野结衣高清无吗| 非洲黑人性xxxx精品又粗又长| 神马国产精品三级电影在线观看| 男人狂女人下面高潮的视频| 久久这里只有精品中国| 亚洲人成网站高清观看| 国产精品日韩av在线免费观看| 亚洲人成网站在线观看播放| 综合色丁香网| 国产精品一区二区三区四区免费观看 | 欧美潮喷喷水| av在线天堂中文字幕| 国产亚洲精品av在线| 乱人视频在线观看| 国产毛片a区久久久久| 六月丁香七月| 亚洲一区二区三区色噜噜| 亚洲欧美日韩高清专用| 国产精品一区二区性色av| 国产美女午夜福利| 亚洲乱码一区二区免费版| 欧美中文日本在线观看视频| 男女做爰动态图高潮gif福利片| 国内久久婷婷六月综合欲色啪| 丰满的人妻完整版| 真实男女啪啪啪动态图| 国产一区二区激情短视频| 97超碰精品成人国产| 校园人妻丝袜中文字幕| 国产高清不卡午夜福利| 亚洲国产色片| 亚洲乱码一区二区免费版| 精品日产1卡2卡| 国产精品久久久久久精品电影| 老师上课跳d突然被开到最大视频| 国产精品,欧美在线| 国产又黄又爽又无遮挡在线| 国产一区二区亚洲精品在线观看| 尾随美女入室| 免费搜索国产男女视频| 高清午夜精品一区二区三区 | 99久久精品一区二区三区| 中文字幕熟女人妻在线| 国产精品不卡视频一区二区| 亚洲七黄色美女视频| av福利片在线观看| 最近手机中文字幕大全| 久久精品人妻少妇| 晚上一个人看的免费电影| 日韩欧美一区二区三区在线观看| 联通29元200g的流量卡| 99九九线精品视频在线观看视频| 三级经典国产精品| 国产片特级美女逼逼视频| 午夜福利在线在线| 蜜桃久久精品国产亚洲av| 久久久久久久久久黄片| 18+在线观看网站| 亚洲七黄色美女视频| 成人亚洲欧美一区二区av| 国产欧美日韩一区二区精品| 村上凉子中文字幕在线| 免费观看在线日韩| 国产片特级美女逼逼视频| 国产av在哪里看| 国产精品伦人一区二区| 久久99热这里只有精品18| 国产色婷婷99| 亚洲av一区综合| h日本视频在线播放| 国产亚洲精品久久久com| 又粗又爽又猛毛片免费看| 男人狂女人下面高潮的视频| 久久精品国产清高在天天线| 波多野结衣巨乳人妻| 日本爱情动作片www.在线观看 | 乱系列少妇在线播放| 日韩欧美精品v在线| 成年av动漫网址| 国产精品免费一区二区三区在线| 日日摸夜夜添夜夜添小说| 欧美性猛交黑人性爽| 特级一级黄色大片| 97热精品久久久久久| 欧美一区二区精品小视频在线| 亚洲一区高清亚洲精品| 一级a爱片免费观看的视频| 久久草成人影院| 97人妻精品一区二区三区麻豆| 日本一二三区视频观看| 亚洲最大成人中文| 免费观看人在逋| 亚洲性久久影院| 久久天躁狠狠躁夜夜2o2o| 日本黄大片高清| 国产精品一二三区在线看| 久久精品国产亚洲av香蕉五月| 啦啦啦韩国在线观看视频| 少妇裸体淫交视频免费看高清| 成人亚洲欧美一区二区av| 亚洲成人久久爱视频| 亚洲精品成人久久久久久| 国产高清三级在线| 偷拍熟女少妇极品色| 一级黄色大片毛片| 午夜激情福利司机影院| 国产精品av视频在线免费观看| 99热精品在线国产| 国产高清有码在线观看视频| 日本黄大片高清| 性欧美人与动物交配| 99久久中文字幕三级久久日本| 看黄色毛片网站| 毛片一级片免费看久久久久| 国产高清视频在线播放一区| 午夜福利成人在线免费观看| 国产成年人精品一区二区| 久久久久免费精品人妻一区二区| 97热精品久久久久久| 91午夜精品亚洲一区二区三区| 中文字幕免费在线视频6| 亚洲自偷自拍三级| 国产探花极品一区二区| 男人和女人高潮做爰伦理| 国产真实伦视频高清在线观看| 丝袜喷水一区| 欧美日韩乱码在线| 国产成人精品久久久久久| 干丝袜人妻中文字幕| 精品乱码久久久久久99久播| 国产伦在线观看视频一区| 俺也久久电影网| 人人妻人人澡欧美一区二区| 久久人妻av系列| 欧美一区二区亚洲| 欧美日韩一区二区视频在线观看视频在线 | 一级毛片我不卡| 人人妻,人人澡人人爽秒播| 亚洲精品乱码久久久v下载方式| 国产一区二区在线av高清观看| 黄色一级大片看看| 国产精品美女特级片免费视频播放器| 婷婷色综合大香蕉| 可以在线观看毛片的网站| 嫩草影院新地址| 国产午夜福利久久久久久| 成人二区视频| 欧美精品国产亚洲| 国产毛片a区久久久久| 成人性生交大片免费视频hd| 黄片wwwwww| av专区在线播放| 性插视频无遮挡在线免费观看| 国产精品综合久久久久久久免费| 99热这里只有精品一区| 亚洲无线在线观看| 99视频精品全部免费 在线| 久久99热6这里只有精品| 97超视频在线观看视频| 内射极品少妇av片p| 国产精华一区二区三区| 久久人人精品亚洲av| 国产成人a区在线观看| 天美传媒精品一区二区| 免费在线观看影片大全网站| 99热这里只有精品一区| 综合色av麻豆| 国产精品久久久久久久电影| 亚洲精品456在线播放app| 久久精品影院6| 白带黄色成豆腐渣| 一级黄片播放器| 美女黄网站色视频| 欧美+日韩+精品| 毛片女人毛片| 看非洲黑人一级黄片| 免费一级毛片在线播放高清视频| 69av精品久久久久久| 97碰自拍视频| 老女人水多毛片| 少妇裸体淫交视频免费看高清| 色哟哟哟哟哟哟| 亚洲人成网站在线播| 淫秽高清视频在线观看| 你懂的网址亚洲精品在线观看 | 韩国av在线不卡| 日韩欧美三级三区| 看免费成人av毛片| 少妇人妻精品综合一区二区 | 乱人视频在线观看| 国产成人精品久久久久久| 欧美另类亚洲清纯唯美| 综合色av麻豆| 精品一区二区三区视频在线观看免费| 女的被弄到高潮叫床怎么办| 成人特级av手机在线观看| 亚洲av二区三区四区| 夜夜夜夜夜久久久久| 国产亚洲91精品色在线| 色综合色国产| 日本撒尿小便嘘嘘汇集6| 国产一区亚洲一区在线观看| 人妻制服诱惑在线中文字幕| 乱系列少妇在线播放| 高清日韩中文字幕在线| 日韩,欧美,国产一区二区三区 | av卡一久久| 在线观看av片永久免费下载| 国产又黄又爽又无遮挡在线| av视频在线观看入口| 国产精品久久久久久亚洲av鲁大| 日本一本二区三区精品| 日本a在线网址| 97碰自拍视频| 久久久精品大字幕| 色播亚洲综合网| 男人舔奶头视频| 亚洲av成人精品一区久久| 岛国在线免费视频观看| 蜜臀久久99精品久久宅男| 青春草视频在线免费观看| 精品久久久久久久末码|