Search Results
You are looking at 1 - 10 of 37 items for
- Author or Editor: Lin Zhao x
- Refine by Access: All Content x
Abstract
Performance of global climate models (GCMs) is strongly affected by the cumulus parameterization (CP) used. Similar to the approach in GFDL AM4, a double-plume CP, which unifies the deep and shallow convection in one framework, is implemented and tested in the NCAR Community Atmospheric Model version 5 (CAM5). Based on the University of Washington (UW) shallow convection scheme, an additional plume was added to represent the deep convection. The shallow and deep plumes share the same cloud model, but use different triggers, fractional mixing rates, and closures. The scheme was tested in single-column, short-term hindcast, and AMIP simulations. Compared with the default combination of the Zhang–McFarlane scheme and UW scheme in CAM5, the new scheme tends to produce a top-heavy mass flux profile during the active monsoon period in the single-column simulations. The scheme increases the intensity of tropical precipitation, closer to TRMM observations. The new scheme increased subtropical marine boundary layer clouds and high clouds over the deep tropics, both in better agreement with observations. Sensitivity tests indicate that regime-dependent fractional entrainment rates of the deep plume are desired to improve tropical precipitation distribution and upper troposphere temperature. This study suggests that a double-plume approach is a promising way to combine shallow and deep convections in a unified framework.
Abstract
Performance of global climate models (GCMs) is strongly affected by the cumulus parameterization (CP) used. Similar to the approach in GFDL AM4, a double-plume CP, which unifies the deep and shallow convection in one framework, is implemented and tested in the NCAR Community Atmospheric Model version 5 (CAM5). Based on the University of Washington (UW) shallow convection scheme, an additional plume was added to represent the deep convection. The shallow and deep plumes share the same cloud model, but use different triggers, fractional mixing rates, and closures. The scheme was tested in single-column, short-term hindcast, and AMIP simulations. Compared with the default combination of the Zhang–McFarlane scheme and UW scheme in CAM5, the new scheme tends to produce a top-heavy mass flux profile during the active monsoon period in the single-column simulations. The scheme increases the intensity of tropical precipitation, closer to TRMM observations. The new scheme increased subtropical marine boundary layer clouds and high clouds over the deep tropics, both in better agreement with observations. Sensitivity tests indicate that regime-dependent fractional entrainment rates of the deep plume are desired to improve tropical precipitation distribution and upper troposphere temperature. This study suggests that a double-plume approach is a promising way to combine shallow and deep convections in a unified framework.
Abstract
High-resolution global climate models (GCMs) have been increasingly utilized for simulations of the global number and distribution of tropical cyclones (TCs), and how they might change with changing climate. In contrast, there is a lack of published studies on the sensitivity of TC genesis to parameterized processes in these GCMs. The uncertainties in these formulations might be an important source of uncertainty in the future projections of TC statistics.
This study investigates the sensitivity of the global number of TCs in present-day simulations using the Geophysical Fluid Dynamics Laboratory High Resolution Atmospheric Model (GFDL HIRAM) to alterations in physical parameterizations. Two parameters are identified to be important in TC genesis frequency in this model: the horizontal cumulus mixing rate, which controls the entrainment into convective cores within the convection parameterization, and the strength of the damping of the divergent component of the horizontal flow. The simulated global number of TCs exhibits nonintuitive response to incremental changes of both parameters. As the cumulus mixing rate increases, the model produces nonmonotonic response in global TC frequency with an initial sharp increase and then a decrease. However, storm mean intensity rises monotonically with the mixing rate. As the strength of the divergence damping increases, the model produces a continuous increase of global number of TCs and hurricanes with little change in storm mean intensity. Mechanisms for explaining these nonintuitive responses are discussed.
Abstract
High-resolution global climate models (GCMs) have been increasingly utilized for simulations of the global number and distribution of tropical cyclones (TCs), and how they might change with changing climate. In contrast, there is a lack of published studies on the sensitivity of TC genesis to parameterized processes in these GCMs. The uncertainties in these formulations might be an important source of uncertainty in the future projections of TC statistics.
This study investigates the sensitivity of the global number of TCs in present-day simulations using the Geophysical Fluid Dynamics Laboratory High Resolution Atmospheric Model (GFDL HIRAM) to alterations in physical parameterizations. Two parameters are identified to be important in TC genesis frequency in this model: the horizontal cumulus mixing rate, which controls the entrainment into convective cores within the convection parameterization, and the strength of the damping of the divergent component of the horizontal flow. The simulated global number of TCs exhibits nonintuitive response to incremental changes of both parameters. As the cumulus mixing rate increases, the model produces nonmonotonic response in global TC frequency with an initial sharp increase and then a decrease. However, storm mean intensity rises monotonically with the mixing rate. As the strength of the divergence damping increases, the model produces a continuous increase of global number of TCs and hurricanes with little change in storm mean intensity. Mechanisms for explaining these nonintuitive responses are discussed.
Abstract
Using observed and reanalysis data, the pronounced interdecadal variations of Lake Qinghai (LQH) water levels and associated climate factors were diagnosed. From the 1960s to the early 2000s, the water level of LQH in the Tibetan Plateau has experienced a continual decline of 3 m but has since increased considerably. A water budget analysis of the LQH watershed suggested that the water vapor flux divergence 
Abstract
Using observed and reanalysis data, the pronounced interdecadal variations of Lake Qinghai (LQH) water levels and associated climate factors were diagnosed. From the 1960s to the early 2000s, the water level of LQH in the Tibetan Plateau has experienced a continual decline of 3 m but has since increased considerably. A water budget analysis of the LQH watershed suggested that the water vapor flux divergence
Abstract
Using abundant rainfall gauge measurements and Global Precipitation Mission (GPM) data, spatial patterns of rainfall diurnal cycles and their seasonality over high mountain Asia (HMA) were examined. Spatial distributions of rainfall diurnal cycles over the HMA have a prominent seasonality regulated by circulations at different spatiotemporal scales, within which large regional contrasts are embedded. Rainfall diurnal variability is relatively weak in the premonsoon season, with larger amplitude over the western HMA, the southeastern HMA, as well as southern periphery regions, characterized by a dominant late afternoon to morning rainfall preference. The pattern of rainfall spatial distributions is closely related to the midlatitude westerlies. Both the mean rainfall and amplitudes of diurnal cycles become more pronounced with the advance of monsoon season but weaken during postmonsoon. The widespread late afternoon to night pattern over HMA migrating with seasonal atmospheric circulation is consistent with the lifetime of convective systems, which become active from the afternoon due to radiative heating and decay during the night. Stationary terrain-dependent night-to-morning rainfall patterns are visible in those east–west-orientated valleys over HMA and the Qaidam basin throughout the seasons. This salient geographical dependence is associated with local circulation produced by the strong differential thermal conditions over mountains and valleys, which can lift the warm moist air at the mouth of the valley and trigger nocturnal convection.
Significance Statement
The main purpose of this study is to explore how spatial patterns of rainfall diurnal cycles over high mountain Asia vary with the seasons. Our results show that the widespread late afternoon to night rainfall over high mountain Asia migrating with seasonal atmospheric circulation is consistent with the lifetime of convective systems. Stationary terrain-dependent night-to-morning rainfall patterns are visible in those east–west-orientated valleys over high mountain Asia and the Qaidam basin throughout the seasons. These results highlight the importance of large-scale atmospheric circulation and local circulation on precipitation, which is critical for water resources over high mountain Asia.
Abstract
Using abundant rainfall gauge measurements and Global Precipitation Mission (GPM) data, spatial patterns of rainfall diurnal cycles and their seasonality over high mountain Asia (HMA) were examined. Spatial distributions of rainfall diurnal cycles over the HMA have a prominent seasonality regulated by circulations at different spatiotemporal scales, within which large regional contrasts are embedded. Rainfall diurnal variability is relatively weak in the premonsoon season, with larger amplitude over the western HMA, the southeastern HMA, as well as southern periphery regions, characterized by a dominant late afternoon to morning rainfall preference. The pattern of rainfall spatial distributions is closely related to the midlatitude westerlies. Both the mean rainfall and amplitudes of diurnal cycles become more pronounced with the advance of monsoon season but weaken during postmonsoon. The widespread late afternoon to night pattern over HMA migrating with seasonal atmospheric circulation is consistent with the lifetime of convective systems, which become active from the afternoon due to radiative heating and decay during the night. Stationary terrain-dependent night-to-morning rainfall patterns are visible in those east–west-orientated valleys over HMA and the Qaidam basin throughout the seasons. This salient geographical dependence is associated with local circulation produced by the strong differential thermal conditions over mountains and valleys, which can lift the warm moist air at the mouth of the valley and trigger nocturnal convection.
Significance Statement
The main purpose of this study is to explore how spatial patterns of rainfall diurnal cycles over high mountain Asia vary with the seasons. Our results show that the widespread late afternoon to night rainfall over high mountain Asia migrating with seasonal atmospheric circulation is consistent with the lifetime of convective systems. Stationary terrain-dependent night-to-morning rainfall patterns are visible in those east–west-orientated valleys over high mountain Asia and the Qaidam basin throughout the seasons. These results highlight the importance of large-scale atmospheric circulation and local circulation on precipitation, which is critical for water resources over high mountain Asia.
Abstract
Confidence and uncertainty issues of simulations were seldom evaluated in previous studies although the climate models are widely used. This study evaluates the performance of the CMIP6-HighResMIP simulations in presenting long-term variability of tropical cyclone (TC) genesis frequency (TCGF) and track density (TCTD) and quantifies the relative contributions of internal and external forcing to TC activities during the 1950–2014. There is overall poor model performance in simulating long-term changes in TC activities over the Northern Hemisphere, including interdecadal variabilities and long-term linear trends. The simulated long-term changes in TCGF and TCTD over the eastern North Pacific (ENP) in six high-resolution models show opposite characteristics to the observations. Moreover, most models cannot capture the variabilities of TCGF and TCTD over the western part of the western North Pacific (WNP) and northern part of the North Atlantic (NA). However, these models show a high degree of confidence in reproducing the interdecadal variabilities and linear trends of TCGF and TCTD over the eastern part of the WNP and the tropical NA. Quantitative evaluations further show that there are the opposite relative contributions of long-term climate variabilities to TCGF and TCTD changes over the ENP between the observations and the multimodel ensemble mean, followed by large model biases over the western WNP and the northern NA, but relatively consistent contributions over the southern NA and the Caribbean. These results help us cope with contrasting and consistent future TC changes among the model projections.
Significance Statement
While climate models have been widely used to project future changes in tropical cyclone (TC) activity, few studies have examined to what extent we can trust these model projections. We used the CMIP6-HighResMIP simulations to quantify the model biases in presenting TC activity, and evaluate the relative contributions of internal and external forcing to TC activities. In general, the HighResMIP has large discrepancies in representing longer-term climate variability of TC activity. However, the models can capture well TC activity over the eastern part of the western North Pacific and tropical Atlantic, which is attributed to good performance of models in reproducing the relationship between long-term climate variabilities beyond interannual scale and TC activity. These results highlight confidence and uncertainty in future TC changes among the model projections.
Abstract
Confidence and uncertainty issues of simulations were seldom evaluated in previous studies although the climate models are widely used. This study evaluates the performance of the CMIP6-HighResMIP simulations in presenting long-term variability of tropical cyclone (TC) genesis frequency (TCGF) and track density (TCTD) and quantifies the relative contributions of internal and external forcing to TC activities during the 1950–2014. There is overall poor model performance in simulating long-term changes in TC activities over the Northern Hemisphere, including interdecadal variabilities and long-term linear trends. The simulated long-term changes in TCGF and TCTD over the eastern North Pacific (ENP) in six high-resolution models show opposite characteristics to the observations. Moreover, most models cannot capture the variabilities of TCGF and TCTD over the western part of the western North Pacific (WNP) and northern part of the North Atlantic (NA). However, these models show a high degree of confidence in reproducing the interdecadal variabilities and linear trends of TCGF and TCTD over the eastern part of the WNP and the tropical NA. Quantitative evaluations further show that there are the opposite relative contributions of long-term climate variabilities to TCGF and TCTD changes over the ENP between the observations and the multimodel ensemble mean, followed by large model biases over the western WNP and the northern NA, but relatively consistent contributions over the southern NA and the Caribbean. These results help us cope with contrasting and consistent future TC changes among the model projections.
Significance Statement
While climate models have been widely used to project future changes in tropical cyclone (TC) activity, few studies have examined to what extent we can trust these model projections. We used the CMIP6-HighResMIP simulations to quantify the model biases in presenting TC activity, and evaluate the relative contributions of internal and external forcing to TC activities. In general, the HighResMIP has large discrepancies in representing longer-term climate variability of TC activity. However, the models can capture well TC activity over the eastern part of the western North Pacific and tropical Atlantic, which is attributed to good performance of models in reproducing the relationship between long-term climate variabilities beyond interannual scale and TC activity. These results highlight confidence and uncertainty in future TC changes among the model projections.
Abstract
The Indian Ocean witnessed a weak positive Indian Ocean dipole (IOD) event from the boreal summer to autumn in 2015, while an extreme El Niño occurred over the tropical Pacific. This was different from the case in 1997/98, when an extreme El Niño and the strongest IOD took place simultaneously. The analysis here suggests that the unique sea surface temperature anomaly (SSTA) pattern of El Niño in 2015 might have contributed to the weak IOD that year. El Niño in 2015 had a complex SSTA pattern, with positive warming over the central and eastern tropical Pacific. Such a combination of the classic El Niño (also known as cold-tongue El Niño) and the recently identified central Pacific El Niño (also known as El Niño Modoki II) had opposite remote influences on the tropical Indian Ocean. The classic El Niño reduced the strength of the Walker circulation over the tropical Indian Ocean, but this was offset by El Niño Modoki II. This study points out that the IOD can be strongly modulated by combined El Niño types in some circumstances, as in 2015.
Abstract
The Indian Ocean witnessed a weak positive Indian Ocean dipole (IOD) event from the boreal summer to autumn in 2015, while an extreme El Niño occurred over the tropical Pacific. This was different from the case in 1997/98, when an extreme El Niño and the strongest IOD took place simultaneously. The analysis here suggests that the unique sea surface temperature anomaly (SSTA) pattern of El Niño in 2015 might have contributed to the weak IOD that year. El Niño in 2015 had a complex SSTA pattern, with positive warming over the central and eastern tropical Pacific. Such a combination of the classic El Niño (also known as cold-tongue El Niño) and the recently identified central Pacific El Niño (also known as El Niño Modoki II) had opposite remote influences on the tropical Indian Ocean. The classic El Niño reduced the strength of the Walker circulation over the tropical Indian Ocean, but this was offset by El Niño Modoki II. This study points out that the IOD can be strongly modulated by combined El Niño types in some circumstances, as in 2015.
Abstract
Based on a new version of the Geophysical Fluid Dynamics Laboratory (GFDL) coupled model, the Madden–Julian oscillation (MJO) prediction skill in boreal wintertime (November–April) is evaluated by analyzing 11 years (2003–13) of hindcast experiments. The initial conditions are obtained by applying a simple nudging technique toward observations. Using the real-time multivariate MJO (RMM) index as a predictand, it is demonstrated that the MJO prediction skill can reach out to 27 days before the anomaly correlation coefficient (ACC) decreases to 0.5. The MJO forecast skill also shows relatively larger contrasts between target strong and weak cases (32 versus 7 days) than between initially strong and weak cases (29 versus 24 days). Meanwhile, a strong dependence on target phases is found, as opposed to relative skill independence from different initial phases. The MJO prediction skill is also shown to be about 29 days during the Dynamics of the MJO/Cooperative Indian Ocean Experiment on Intraseasonal Variability in Year 2011 (DYNAMO/CINDY) field campaign period. This model’s potential predictability, the upper bound of prediction skill, extends out to 42 days, revealing a considerable unutilized predictability and a great potential for improving current MJO prediction.
Abstract
Based on a new version of the Geophysical Fluid Dynamics Laboratory (GFDL) coupled model, the Madden–Julian oscillation (MJO) prediction skill in boreal wintertime (November–April) is evaluated by analyzing 11 years (2003–13) of hindcast experiments. The initial conditions are obtained by applying a simple nudging technique toward observations. Using the real-time multivariate MJO (RMM) index as a predictand, it is demonstrated that the MJO prediction skill can reach out to 27 days before the anomaly correlation coefficient (ACC) decreases to 0.5. The MJO forecast skill also shows relatively larger contrasts between target strong and weak cases (32 versus 7 days) than between initially strong and weak cases (29 versus 24 days). Meanwhile, a strong dependence on target phases is found, as opposed to relative skill independence from different initial phases. The MJO prediction skill is also shown to be about 29 days during the Dynamics of the MJO/Cooperative Indian Ocean Experiment on Intraseasonal Variability in Year 2011 (DYNAMO/CINDY) field campaign period. This model’s potential predictability, the upper bound of prediction skill, extends out to 42 days, revealing a considerable unutilized predictability and a great potential for improving current MJO prediction.
Abstract
The attenuation-based rainfall estimator is less sensitive to the variability of raindrop size distributions (DSDs) than conventional radar rainfall estimators. For the attenuation-based quantitative precipitation estimation (QPE), the key is to accurately estimate the horizontal specific attenuation A H , which requires a good estimate of the ray-averaged ratio between A H and specific differential phase K DP, also known as the coefficient α. In this study, a variational approach is proposed to optimize the coefficient α for better estimates of A H and rainfall. The performance of the variational approach is illustrated using observations from an S-band operational weather radar with rigorous quality control in south China, by comparing against the α optimization approach using a slope of differential reflectivity Z DR dependence on horizontal reflectivity factor Z H . Similar to the Z DR-slope approach, the variational approach can obtain the optimized α consistent with the DSD properties of precipitation on a sweep-to-sweep basis. The attenuation-based hourly rainfall estimates using the sweep-averaged α values from these two approaches show comparable accuracy when verified against the gauge measurements. One advantage of the variational approach is its feasibility to optimize α for each radar ray, which mitigates the impact of the azimuthal DSD variabilities on rainfall estimation. It is found that, based on the optimized α for radar rays, the hourly rainfall amounts derived from the variational approach are consistent with gauge measurements, showing lower bias (1.0%), higher correlation coefficient (0.92), and lower root-mean-square error (2.35 mm) than the results based on the sweep-averaged α.
Abstract
The attenuation-based rainfall estimator is less sensitive to the variability of raindrop size distributions (DSDs) than conventional radar rainfall estimators. For the attenuation-based quantitative precipitation estimation (QPE), the key is to accurately estimate the horizontal specific attenuation A H , which requires a good estimate of the ray-averaged ratio between A H and specific differential phase K DP, also known as the coefficient α. In this study, a variational approach is proposed to optimize the coefficient α for better estimates of A H and rainfall. The performance of the variational approach is illustrated using observations from an S-band operational weather radar with rigorous quality control in south China, by comparing against the α optimization approach using a slope of differential reflectivity Z DR dependence on horizontal reflectivity factor Z H . Similar to the Z DR-slope approach, the variational approach can obtain the optimized α consistent with the DSD properties of precipitation on a sweep-to-sweep basis. The attenuation-based hourly rainfall estimates using the sweep-averaged α values from these two approaches show comparable accuracy when verified against the gauge measurements. One advantage of the variational approach is its feasibility to optimize α for each radar ray, which mitigates the impact of the azimuthal DSD variabilities on rainfall estimation. It is found that, based on the optimized α for radar rays, the hourly rainfall amounts derived from the variational approach are consistent with gauge measurements, showing lower bias (1.0%), higher correlation coefficient (0.92), and lower root-mean-square error (2.35 mm) than the results based on the sweep-averaged α.
Abstract
Strong tropical cyclones often undergo eyewall replacement cycles that are accompanied by concentric eyewalls and/or rapid intensity changes while the secondary eyewall contracts radially inward and eventually replaces the inner eyewall. To the best of our knowledge, the only documented partial/incomplete tertiary eyewall has been mostly inferred from two-dimensional satellite images or one-dimensional aircraft flight-level measurements that can be regarded as indirect and tangential. This study presents the first high spatial and temporal resolution Doppler radar observations of a tertiary eyewall formation event in Typhoon Usagi (2013) over a 14-h time period before it makes landfall. The primary (tangential) and secondary (radial) circulations of Usagi deduced from the Ground-Based Velocity Track Display (GBVTD) methodology clearly portrayed three distinct axisymmetric maxima of radar reflectivity, tangential wind, vertical velocity, and vertical vorticity. Usagi’s central pressure steadily deepened during the contraction of the secondary and tertiary eyewalls until the tertiary eyewall hit the coast of southeast China, which erminated the intensification of the storm.
Abstract
Strong tropical cyclones often undergo eyewall replacement cycles that are accompanied by concentric eyewalls and/or rapid intensity changes while the secondary eyewall contracts radially inward and eventually replaces the inner eyewall. To the best of our knowledge, the only documented partial/incomplete tertiary eyewall has been mostly inferred from two-dimensional satellite images or one-dimensional aircraft flight-level measurements that can be regarded as indirect and tangential. This study presents the first high spatial and temporal resolution Doppler radar observations of a tertiary eyewall formation event in Typhoon Usagi (2013) over a 14-h time period before it makes landfall. The primary (tangential) and secondary (radial) circulations of Usagi deduced from the Ground-Based Velocity Track Display (GBVTD) methodology clearly portrayed three distinct axisymmetric maxima of radar reflectivity, tangential wind, vertical velocity, and vertical vorticity. Usagi’s central pressure steadily deepened during the contraction of the secondary and tertiary eyewalls until the tertiary eyewall hit the coast of southeast China, which erminated the intensification of the storm.
Abstract
A global atmospheric model with roughly 50-km horizontal grid spacing is used to simulate the interannual variability of tropical cyclones using observed sea surface temperatures (SSTs) as the lower boundary condition. The model’s convective parameterization is based on a closure for shallow convection, with much of the deep convection allowed to occur on resolved scales. Four realizations of the period 1981–2005 are generated. The correlation of yearly Atlantic hurricane counts with observations is greater than 0.8 when the model is averaged over the four realizations, supporting the view that the random part of this annual Atlantic hurricane frequency (the part not predictable given the SSTs) is relatively small (<2 hurricanes per year). Correlations with observations are lower in the east, west, and South Pacific (roughly 0.6, 0.5, and 0.3, respectively) and insignificant in the Indian Ocean. The model trends in Northern Hemisphere basin-wide frequency are consistent with the observed trends in the International Best Track Archive for Climate Stewardship (IBTrACS) database. The model generates an upward trend of hurricane frequency in the Atlantic and downward trends in the east and west Pacific over this time frame. The model produces a negative trend in the Southern Hemisphere that is larger than that in the IBTrACS.
The same model is used to simulate the response to the SST anomalies generated by coupled models in the World Climate Research Program Coupled Model Intercomparison Project 3 (CMIP3) archive, using the late-twenty-first century in the A1B scenario. Results are presented for SST anomalies computed by averaging over 18 CMIP3 models and from individual realizations from 3 models. A modest reduction of global and Southern Hemisphere tropical cyclone frequency is obtained in each case, but the results in individual Northern Hemisphere basins differ among the models. The vertical shear in the Atlantic Main Development Region (MDR) and the difference between the MDR SST and the tropical mean SST are well correlated with the model’s Atlantic storm frequency, both for interannual variability and for the intermodel spread in global warming projections.
Abstract
A global atmospheric model with roughly 50-km horizontal grid spacing is used to simulate the interannual variability of tropical cyclones using observed sea surface temperatures (SSTs) as the lower boundary condition. The model’s convective parameterization is based on a closure for shallow convection, with much of the deep convection allowed to occur on resolved scales. Four realizations of the period 1981–2005 are generated. The correlation of yearly Atlantic hurricane counts with observations is greater than 0.8 when the model is averaged over the four realizations, supporting the view that the random part of this annual Atlantic hurricane frequency (the part not predictable given the SSTs) is relatively small (<2 hurricanes per year). Correlations with observations are lower in the east, west, and South Pacific (roughly 0.6, 0.5, and 0.3, respectively) and insignificant in the Indian Ocean. The model trends in Northern Hemisphere basin-wide frequency are consistent with the observed trends in the International Best Track Archive for Climate Stewardship (IBTrACS) database. The model generates an upward trend of hurricane frequency in the Atlantic and downward trends in the east and west Pacific over this time frame. The model produces a negative trend in the Southern Hemisphere that is larger than that in the IBTrACS.
The same model is used to simulate the response to the SST anomalies generated by coupled models in the World Climate Research Program Coupled Model Intercomparison Project 3 (CMIP3) archive, using the late-twenty-first century in the A1B scenario. Results are presented for SST anomalies computed by averaging over 18 CMIP3 models and from individual realizations from 3 models. A modest reduction of global and Southern Hemisphere tropical cyclone frequency is obtained in each case, but the results in individual Northern Hemisphere basins differ among the models. The vertical shear in the Atlantic Main Development Region (MDR) and the difference between the MDR SST and the tropical mean SST are well correlated with the model’s Atlantic storm frequency, both for interannual variability and for the intermodel spread in global warming projections.
