Browse
Abstract
High-resolution oceanic precipitation estimates are needed to increase our understanding of and ability to monitor ocean–atmosphere coupled processes. Satellite multisensor precipitation products such as IMERG provide global precipitation estimates at relatively high resolution (0.1°, 30 min), but the resolution at which IMERG precipitation estimates are considered reliable is coarser than the nominal resolution of the product itself. In this study, we examine the ability of the Rainfall Autoregressive Model (RainFARM) statistical downscaling technique to produce ensembles of precipitation fields at relatively high spatial and temporal resolution when applied to spatially and temporally coarsened precipitation fields from IMERG. The downscaled precipitation ensembles are evaluated against in situ oceanic rain-rate observations collected by passive aquatic listeners (PALs) in 11 different ocean domains. We also evaluate IMERG coarsened to the same resolution as the downscaled fields to determine whether the process of coarsening then downscaling improves precipitation estimates more than averaging IMERG to coarser resolution only. Evaluations were performed on individual months, seasons, by ENSO phase, and based on precipitation characteristics. Results were inconsistent, with downscaling improving precipitation estimates in some domains and time periods and producing worse performance in others. While the results imply that the performance of the downscaled precipitation estimates is related to precipitation characteristics, it is still unclear what characteristics or combinations thereof lead to the most improvement or consistent improvement when applying RainFARM to IMERG.
Abstract
High-resolution oceanic precipitation estimates are needed to increase our understanding of and ability to monitor ocean–atmosphere coupled processes. Satellite multisensor precipitation products such as IMERG provide global precipitation estimates at relatively high resolution (0.1°, 30 min), but the resolution at which IMERG precipitation estimates are considered reliable is coarser than the nominal resolution of the product itself. In this study, we examine the ability of the Rainfall Autoregressive Model (RainFARM) statistical downscaling technique to produce ensembles of precipitation fields at relatively high spatial and temporal resolution when applied to spatially and temporally coarsened precipitation fields from IMERG. The downscaled precipitation ensembles are evaluated against in situ oceanic rain-rate observations collected by passive aquatic listeners (PALs) in 11 different ocean domains. We also evaluate IMERG coarsened to the same resolution as the downscaled fields to determine whether the process of coarsening then downscaling improves precipitation estimates more than averaging IMERG to coarser resolution only. Evaluations were performed on individual months, seasons, by ENSO phase, and based on precipitation characteristics. Results were inconsistent, with downscaling improving precipitation estimates in some domains and time periods and producing worse performance in others. While the results imply that the performance of the downscaled precipitation estimates is related to precipitation characteristics, it is still unclear what characteristics or combinations thereof lead to the most improvement or consistent improvement when applying RainFARM to IMERG.
Abstract
The expansion of the boreal forest poleward is a potentially important driver of feedbacks between the land surface and Arctic climate. A growing body of work has highlighted the importance of differences in evaporative resistance between different possible future Arctic land covers, which in turn alters humidity and cloudiness in the boundary layer, for these feedbacks. While thus far this problem has been studied primarily with complex Earth system models, we turn to a locally focused, idealized model capable of diagnosing and testing the sensitivity of first-order processes connecting vegetation, the atmospheric boundary layer, and low clouds in this critical region. This allows us to benchmark the mechanisms and results at the center of predictions from larger-scale simulations. A surface dominated by broadleaf trees, characterized by higher albedo and lower surface evaporative resistance, drives cooling and moistening of the boundary layer relative to a surface of needleleaf trees, characterized by lower albedo and higher surface evaporative resistance. Differences in evaporative resistance between these hypothetical Arctic vegetation covers are of equal importance to changes in albedo for the initial response of the boundary layer to boreal expansion, even with our idealized approach. However, compensation between the elevation of the lifting condensation level (LCL) and more rapid growth of the mixed layer over higher evaporative resistance surfaces can minimize changes in the favorability of shallow clouds over different land cover types under some conditions. We then perform two tests on the sensitivity of this compensating effect, to changes in water availability, represented first by a reduction in boundary layer humidity and then by both a reduction in humidity and soil moisture available to our vegetation surface. Finally, given the importance of this potential LCL–mixed-layer height compensation in our idealized modeling results, we look to determine its relevance in observational data from a field campaign in boreal Finland. These observations do confirm that such a coupling plays an important role in cumulus-topped boundary layers over a needleleaf forest surface. While our results confirm some underlying mechanisms at the center of prior work with Earth system models, they also provide motivation for future work to constrain the impact of boreal forest expansion. This will include both large eddy simulations to examine the impact of processes and feedbacks not resolved by a mixed-layer model, as well as a more systematic evaluation and comparison of relevant observations at the site in Finland and sites from prior boreal field campaigns.
Significance Statement
Clouds and vegetation are both important components of the climate system that interact across a range of scales. These interactions are central to understanding how changes at the land surface feedback on climate. For example, if a forest expands or recedes, diagnosing how that will impact clouds will determine whether you predict warming or cooling temperatures from that shift in the forest area. These predictions are often made with complex Earth system models, but we look to a more idealized representation of the land–atmosphere system to diagnose how shallow clouds should respond to changes in surface properties with different scenarios of boreal forest expansion at a more foundational level. This both grounds our understanding of previous analysis and provides helpful direction for future studies of this relevant and impactful land cover change.
Abstract
The expansion of the boreal forest poleward is a potentially important driver of feedbacks between the land surface and Arctic climate. A growing body of work has highlighted the importance of differences in evaporative resistance between different possible future Arctic land covers, which in turn alters humidity and cloudiness in the boundary layer, for these feedbacks. While thus far this problem has been studied primarily with complex Earth system models, we turn to a locally focused, idealized model capable of diagnosing and testing the sensitivity of first-order processes connecting vegetation, the atmospheric boundary layer, and low clouds in this critical region. This allows us to benchmark the mechanisms and results at the center of predictions from larger-scale simulations. A surface dominated by broadleaf trees, characterized by higher albedo and lower surface evaporative resistance, drives cooling and moistening of the boundary layer relative to a surface of needleleaf trees, characterized by lower albedo and higher surface evaporative resistance. Differences in evaporative resistance between these hypothetical Arctic vegetation covers are of equal importance to changes in albedo for the initial response of the boundary layer to boreal expansion, even with our idealized approach. However, compensation between the elevation of the lifting condensation level (LCL) and more rapid growth of the mixed layer over higher evaporative resistance surfaces can minimize changes in the favorability of shallow clouds over different land cover types under some conditions. We then perform two tests on the sensitivity of this compensating effect, to changes in water availability, represented first by a reduction in boundary layer humidity and then by both a reduction in humidity and soil moisture available to our vegetation surface. Finally, given the importance of this potential LCL–mixed-layer height compensation in our idealized modeling results, we look to determine its relevance in observational data from a field campaign in boreal Finland. These observations do confirm that such a coupling plays an important role in cumulus-topped boundary layers over a needleleaf forest surface. While our results confirm some underlying mechanisms at the center of prior work with Earth system models, they also provide motivation for future work to constrain the impact of boreal forest expansion. This will include both large eddy simulations to examine the impact of processes and feedbacks not resolved by a mixed-layer model, as well as a more systematic evaluation and comparison of relevant observations at the site in Finland and sites from prior boreal field campaigns.
Significance Statement
Clouds and vegetation are both important components of the climate system that interact across a range of scales. These interactions are central to understanding how changes at the land surface feedback on climate. For example, if a forest expands or recedes, diagnosing how that will impact clouds will determine whether you predict warming or cooling temperatures from that shift in the forest area. These predictions are often made with complex Earth system models, but we look to a more idealized representation of the land–atmosphere system to diagnose how shallow clouds should respond to changes in surface properties with different scenarios of boreal forest expansion at a more foundational level. This both grounds our understanding of previous analysis and provides helpful direction for future studies of this relevant and impactful land cover change.
Abstract
Organized systems of deep convective clouds are often associated with high-impact weather and changes in such systems may have implications for climate sensitivity. This has motivated the derivation of many organization indices that attempt to measure the level of deep convective aggregation in models and observations. Here we conduct a comprehensive review of existing methodologies and highlight some of their drawbacks, such as only measuring organization in a relative sense, being biased toward particular spatial scales, or being very sensitive to the details of the calculation algorithm. One widely used metric, I org, uses statistics of nearest-neighbor distances between convective storms to address the first of these concerns, but we show here that it is insensitive to organization beyond the meso-β scale and very contingent on the details of the implementation. We thus introduce a new and complementary metric, L org, based on all-pair convective storm distances, which is also an absolute metric that can discern regular, random, and clustered cloud scenes. It is linearly sensitive to spatial scale in most applications and robust to the implementation methodology. We also derive a discrete form suited to gridded data and provide corrections to account for cyclic boundary conditions and finite, open boundary domains of nonequal aspect ratios. We demonstrate the use of the metric with idealized synthetic configurations, as well as model output and satellite rainfall retrievals in the tropics. We claim that this new metric usefully supplements the existing family of indices that can help to understand convective organization across spatial scales.
Significance Statement
The clustering and organization of convection is associated with high-impact weather and changes could impact climate sensitivity, but no consensus exists on how to best measure organization. Here we suggest a new metric that is robust to the calculation details and can classify scenes as random, clustered, or regular. This new metric can therefore account for spacing of organized convective systems and convective storms on scales spanning tens of kilometers to the entire tropics. We suggest that the new metric L org addresses many shortcomings of existing measures and can act as a useful additional tool to further understanding of convective organization.
Abstract
Organized systems of deep convective clouds are often associated with high-impact weather and changes in such systems may have implications for climate sensitivity. This has motivated the derivation of many organization indices that attempt to measure the level of deep convective aggregation in models and observations. Here we conduct a comprehensive review of existing methodologies and highlight some of their drawbacks, such as only measuring organization in a relative sense, being biased toward particular spatial scales, or being very sensitive to the details of the calculation algorithm. One widely used metric, I org, uses statistics of nearest-neighbor distances between convective storms to address the first of these concerns, but we show here that it is insensitive to organization beyond the meso-β scale and very contingent on the details of the implementation. We thus introduce a new and complementary metric, L org, based on all-pair convective storm distances, which is also an absolute metric that can discern regular, random, and clustered cloud scenes. It is linearly sensitive to spatial scale in most applications and robust to the implementation methodology. We also derive a discrete form suited to gridded data and provide corrections to account for cyclic boundary conditions and finite, open boundary domains of nonequal aspect ratios. We demonstrate the use of the metric with idealized synthetic configurations, as well as model output and satellite rainfall retrievals in the tropics. We claim that this new metric usefully supplements the existing family of indices that can help to understand convective organization across spatial scales.
Significance Statement
The clustering and organization of convection is associated with high-impact weather and changes could impact climate sensitivity, but no consensus exists on how to best measure organization. Here we suggest a new metric that is robust to the calculation details and can classify scenes as random, clustered, or regular. This new metric can therefore account for spacing of organized convective systems and convective storms on scales spanning tens of kilometers to the entire tropics. We suggest that the new metric L org addresses many shortcomings of existing measures and can act as a useful additional tool to further understanding of convective organization.
Abstract
Stratocumulus occur in closed- or open-cell states, which tend to be associated with high or low cloud cover and the absence or presence of precipitation, respectively. Thus, the transition between these states has substantial implications for the role of this cloud type in Earth’s radiation budget. In this study, we analyze transitions between these states using an ensemble of 127 large-eddy simulations, covering a wide range of conditions. Our analysis is focused on the behavior of these clouds in a cloud fraction (fc ) scene albedo (A) phase space, which has been shown in previous studies to be a useful framework for interpreting system behavior. For the transition from closed to open cells, we find that precipitation creates narrower clouds and scavenges cloud droplets for all fc . However, precipitation decreases the cloud depth for fc > 0.8 only, causing a rapid decrease in A. For fc < 0.8, the cloud depth actually increases due to mesoscale organization of the cloud field. As the cloud deepening balances the effects of cloud droplet scavenging in terms of influence on A, changes in A are determined by the decreasing fc only, causing a linear decrease in A for fc < 0.8. For the transition from open to closed cells, we find that longwave radiative cooling drives the cloud development, with cloud widening dominating for fc < 0.5. For fc > 0.5, clouds begin to deepen gradually due to the decreasing efficiency of lateral expansion. The smooth switch between cloud widening and deepening leads to a more gentle change in A compared to the transitions under precipitating conditions.
Significance Statement
By reflecting a substantial fraction of solar shortwave radiation back to space, shallow clouds constitute a major cooling agent in Earth’s radiation budget. To constrain this effect, a profound understanding of cloud cover and cloud albedo is necessary. In this study, we analyze the processes that drive the variability in these cloud properties in stratocumulus clouds, a very common cloud type covering approximately 20% of the globe. For these clouds, we show that changes from low to high or high to low cloud cover are different due to the underlying cloud micro- and macrophysics, elucidating this crucial aspect of aerosol–cloud–climate interactions.
Abstract
Stratocumulus occur in closed- or open-cell states, which tend to be associated with high or low cloud cover and the absence or presence of precipitation, respectively. Thus, the transition between these states has substantial implications for the role of this cloud type in Earth’s radiation budget. In this study, we analyze transitions between these states using an ensemble of 127 large-eddy simulations, covering a wide range of conditions. Our analysis is focused on the behavior of these clouds in a cloud fraction (fc ) scene albedo (A) phase space, which has been shown in previous studies to be a useful framework for interpreting system behavior. For the transition from closed to open cells, we find that precipitation creates narrower clouds and scavenges cloud droplets for all fc . However, precipitation decreases the cloud depth for fc > 0.8 only, causing a rapid decrease in A. For fc < 0.8, the cloud depth actually increases due to mesoscale organization of the cloud field. As the cloud deepening balances the effects of cloud droplet scavenging in terms of influence on A, changes in A are determined by the decreasing fc only, causing a linear decrease in A for fc < 0.8. For the transition from open to closed cells, we find that longwave radiative cooling drives the cloud development, with cloud widening dominating for fc < 0.5. For fc > 0.5, clouds begin to deepen gradually due to the decreasing efficiency of lateral expansion. The smooth switch between cloud widening and deepening leads to a more gentle change in A compared to the transitions under precipitating conditions.
Significance Statement
By reflecting a substantial fraction of solar shortwave radiation back to space, shallow clouds constitute a major cooling agent in Earth’s radiation budget. To constrain this effect, a profound understanding of cloud cover and cloud albedo is necessary. In this study, we analyze the processes that drive the variability in these cloud properties in stratocumulus clouds, a very common cloud type covering approximately 20% of the globe. For these clouds, we show that changes from low to high or high to low cloud cover are different due to the underlying cloud micro- and macrophysics, elucidating this crucial aspect of aerosol–cloud–climate interactions.
Abstract
Annual spring and summer runoff from western Colorado is relied upon by 40 million people, six states, and two countries. Cool season precipitation and snowpack have historically been robust predictors of seasonal runoff in western Colorado. Forecasts made with this information allow water managers to plan for the season ahead. Antecedent hydrological conditions, such as root zone soil moisture and groundwater storage, and weather conditions following peak snowpack also impact seasonal runoff. The roles of such factors were scrutinized in 2020 and 2021: seasonal runoff was much lower than expectations based on snowpack values alone. We investigate the relative importance of meteorological and hydrological conditions occurring before and after the snowpack season in predicting seasonal runoff in western Colorado. This question is critical because the most effective investment strategy for improving forecasts depends on if errors arise before or after the snowpack season. This study is conducted using observations from the Snow Telemetry Network, root zone soil moisture and groundwater data from the Western Land Data Assimilation Systems, and a random forest–based statistical forecasting framework. We find that on average, antecedent root zone soil moisture and groundwater storage values do not add significant skill to seasonal water supply forecasts in western Colorado. In contrast, using precipitation and temperature data after the time of peak snowpack improves water supply forecasts significantly. The 2020 and 2021 runoffs were hampered by dry conditions both before and after the snowpack season. Both antecedent soil moisture and spring/summer precipitation data improved water supply forecast accuracy in these years.
Significance Statement
Seasonal water supply forecasts in western Colorado are highly valuable because spring and summer runoff from this region helps support the water supply of 40 million people. Accurate forecasts improve the management of the region’s water. Heavy investments have been made in improving our ability to monitor antecedent hydrological conditions in western Colorado, such as root zone soil moisture and groundwater. However, results from this study indicate that the largest source of uncertainty in western Colorado runoff forecasts is future weather. Therefore, improved subseasonal-to-seasonal weather forecasts for western Colorado are what is most needed to improve regional water supply forecasts, and the ability to properly manage western Colorado water.
Abstract
Annual spring and summer runoff from western Colorado is relied upon by 40 million people, six states, and two countries. Cool season precipitation and snowpack have historically been robust predictors of seasonal runoff in western Colorado. Forecasts made with this information allow water managers to plan for the season ahead. Antecedent hydrological conditions, such as root zone soil moisture and groundwater storage, and weather conditions following peak snowpack also impact seasonal runoff. The roles of such factors were scrutinized in 2020 and 2021: seasonal runoff was much lower than expectations based on snowpack values alone. We investigate the relative importance of meteorological and hydrological conditions occurring before and after the snowpack season in predicting seasonal runoff in western Colorado. This question is critical because the most effective investment strategy for improving forecasts depends on if errors arise before or after the snowpack season. This study is conducted using observations from the Snow Telemetry Network, root zone soil moisture and groundwater data from the Western Land Data Assimilation Systems, and a random forest–based statistical forecasting framework. We find that on average, antecedent root zone soil moisture and groundwater storage values do not add significant skill to seasonal water supply forecasts in western Colorado. In contrast, using precipitation and temperature data after the time of peak snowpack improves water supply forecasts significantly. The 2020 and 2021 runoffs were hampered by dry conditions both before and after the snowpack season. Both antecedent soil moisture and spring/summer precipitation data improved water supply forecast accuracy in these years.
Significance Statement
Seasonal water supply forecasts in western Colorado are highly valuable because spring and summer runoff from this region helps support the water supply of 40 million people. Accurate forecasts improve the management of the region’s water. Heavy investments have been made in improving our ability to monitor antecedent hydrological conditions in western Colorado, such as root zone soil moisture and groundwater. However, results from this study indicate that the largest source of uncertainty in western Colorado runoff forecasts is future weather. Therefore, improved subseasonal-to-seasonal weather forecasts for western Colorado are what is most needed to improve regional water supply forecasts, and the ability to properly manage western Colorado water.
Abstract
Satellites provide a useful way of estimating rainfall where the availability of in situ data is low but their indirect nature of estimation means there can be substantial biases. Consequently, the assimilation of in situ data is an important step in improving the accuracy of the satellite rainfall analysis. The effectiveness of this step varies with gauge density, and this study investigated the effectiveness of statistical interpolation (SI), also known as optimal interpolation (OI), on a monthly time scale when gauge density is extremely low using Papua New Guinea (PNG) as a study region. The topography of the region presented an additional challenge to the algorithm. An open-source implementation of SI was developed on Python 3 and confirmed to be consistent with an existing implementation, addressing a lack of open-source implementation for this classical algorithm. The effectiveness of the analysis produced by this algorithm was then compared to the pure satellite analysis over PNG from 2001 to 2014. When performance over the entire study domain was considered, the improvement from using SI was close to imperceptible because of the small number of stations available for assimilation and the small radius of influence of each station (imposed by the topography present in the domain). However, there was still value in using OI as performance around each of the stations was noticeably improved, with the error consistently being reduced along with a general increase in the correlation metric. Furthermore, in an operational context, the use of OI provides an important function of ensuring consistency between in situ data and the gridded analysis.
Significance Statement
The blending of satellite and gauge rainfall data through a process known as statistical interpolation (SI) is known to be capable of producing a more accurate dataset that facilitates better estimation of rainfall. However, the performance of this algorithm over a domain such as Papua New Guinea, where gauge density is extremely low, is not often explored. This study reveals that, although an improvement over the entire Papua New Guinea domain was slight, the algorithm is still valuable as there was a consistent improvement around the stations. Additionally, an adaptable and open-source version of the algorithm is provided, allowing users to blend their own satellite and gauge data and create better geospatial datasets for their own purposes.
Abstract
Satellites provide a useful way of estimating rainfall where the availability of in situ data is low but their indirect nature of estimation means there can be substantial biases. Consequently, the assimilation of in situ data is an important step in improving the accuracy of the satellite rainfall analysis. The effectiveness of this step varies with gauge density, and this study investigated the effectiveness of statistical interpolation (SI), also known as optimal interpolation (OI), on a monthly time scale when gauge density is extremely low using Papua New Guinea (PNG) as a study region. The topography of the region presented an additional challenge to the algorithm. An open-source implementation of SI was developed on Python 3 and confirmed to be consistent with an existing implementation, addressing a lack of open-source implementation for this classical algorithm. The effectiveness of the analysis produced by this algorithm was then compared to the pure satellite analysis over PNG from 2001 to 2014. When performance over the entire study domain was considered, the improvement from using SI was close to imperceptible because of the small number of stations available for assimilation and the small radius of influence of each station (imposed by the topography present in the domain). However, there was still value in using OI as performance around each of the stations was noticeably improved, with the error consistently being reduced along with a general increase in the correlation metric. Furthermore, in an operational context, the use of OI provides an important function of ensuring consistency between in situ data and the gridded analysis.
Significance Statement
The blending of satellite and gauge rainfall data through a process known as statistical interpolation (SI) is known to be capable of producing a more accurate dataset that facilitates better estimation of rainfall. However, the performance of this algorithm over a domain such as Papua New Guinea, where gauge density is extremely low, is not often explored. This study reveals that, although an improvement over the entire Papua New Guinea domain was slight, the algorithm is still valuable as there was a consistent improvement around the stations. Additionally, an adaptable and open-source version of the algorithm is provided, allowing users to blend their own satellite and gauge data and create better geospatial datasets for their own purposes.
Abstract
Severe floods and droughts, including their back-to-back occurrences (weather whiplash), have been increasing in frequency and severity around the world. Improved understanding of systematic changes in hydrological extremes is essential for preparation and adaptation. In this study, we identified and quantified extreme wet and dry events globally by applying a clustering algorithm to terrestrial water storage (TWS) data from the Gravity Recovery and Climate Experiment (GRACE) and GRACE Follow-On (FO). The most intense events, ranked using an intensity metric, often reflect impacts of large-scale oceanic oscillations such as El Niño–Southern Oscillation and consequences of climate change. The severity of both wet and dry events, represented by standardized TWS anomalies, increased significantly in most cases, likely associated with intensification of wet and dry weather regimes in a warmer world, and consequently, exhibited strongest correlation with global temperature. In the Dry climate, the number of wet events decreased while the number of dry events increased significantly, suggesting a drying trend that may be attributed to climate variability and possible increases in irrigation and reliance on groundwater. In the Continental climate where temperature has risen faster than global average, dry events increased significantly. Characteristics of extreme events often showed strong correlations with global temperature, especially when averaged over all climates. These results suggest changes in hydrological extremes and underscore the importance of quantifying total water storage changes when studying hydrological extremes. Extending the GRACE/FO record, which spans 2002 to the present, is essential to continuously tracking changes in TWS and hydrological extremes.
Abstract
Severe floods and droughts, including their back-to-back occurrences (weather whiplash), have been increasing in frequency and severity around the world. Improved understanding of systematic changes in hydrological extremes is essential for preparation and adaptation. In this study, we identified and quantified extreme wet and dry events globally by applying a clustering algorithm to terrestrial water storage (TWS) data from the Gravity Recovery and Climate Experiment (GRACE) and GRACE Follow-On (FO). The most intense events, ranked using an intensity metric, often reflect impacts of large-scale oceanic oscillations such as El Niño–Southern Oscillation and consequences of climate change. The severity of both wet and dry events, represented by standardized TWS anomalies, increased significantly in most cases, likely associated with intensification of wet and dry weather regimes in a warmer world, and consequently, exhibited strongest correlation with global temperature. In the Dry climate, the number of wet events decreased while the number of dry events increased significantly, suggesting a drying trend that may be attributed to climate variability and possible increases in irrigation and reliance on groundwater. In the Continental climate where temperature has risen faster than global average, dry events increased significantly. Characteristics of extreme events often showed strong correlations with global temperature, especially when averaged over all climates. These results suggest changes in hydrological extremes and underscore the importance of quantifying total water storage changes when studying hydrological extremes. Extending the GRACE/FO record, which spans 2002 to the present, is essential to continuously tracking changes in TWS and hydrological extremes.
Abstract
North Atlantic sea surface temperature (SST) variability plays a critical role in modulating the climate system. However, characterizing patterns of North Atlantic SST variability and diagnosing the associated mechanisms is challenging because they involve coupled atmosphere–ocean interactions with complex spatiotemporal relationships. Here we address these challenges by applying a time-evolving self-organizing map approach to a long preindustrial coupled control simulation and identify a variety of 10-yr spatiotemporal evolutions of winter SST anomalies, including but not limited to those associated with the North Atlantic Oscillation–Atlantic multidecadal variability (NAO–AMV)-like interactions. To assess mechanisms and atmospheric responses associated with various SST spatiotemporal evolutions, composites of atmospheric and oceanic variables associated with these evolutions are investigated. Results show that transient-eddy activities and atmospheric circulation responses exist in almost all the evolutions that are closely correlated to the details of the SST pattern. In terms of the mechanisms responsible for generating various SST evolutions, composites of ocean heat budget terms demonstrate that contributions to upper-ocean temperature tendency from resolved ocean advection and surface heat fluxes rarely oppose each other over 10-yr periods in the subpolar North Atlantic. We further explore the potential for predictability for some of these 10-yr SST evolutions that start with similar states but end with different states. However, we find that these are associated with abrupt changes in atmospheric variability and are unlikely to be predictable. In summary, this study broadly investigates the atmospheric responses to and the mechanisms governing the North Atlantic SST evolutions over 10-yr periods.
Significance Statement
Climate variability in the North Atlantic Ocean has wide-ranging impacts on global and regional climate. However, the processes involved include interactions between the ocean and atmosphere that vary across both space and time, making it challenging to characterize and predict. Using a novel machine learning approach, this study identifies various time evolutions of North Atlantic sea surface temperature patterns over 10-yr periods. This includes evolutions with similar start states but different trajectories, which have important implications for predictability. Furthermore, we investigate the mechanisms responsible for these evolutions and how different sea surface temperature patterns affect atmospheric circulation through small-scale atmospheric disturbances. These new insights into the complex ocean–atmosphere interactions over time are critical for improving decadal prediction skill.
Abstract
North Atlantic sea surface temperature (SST) variability plays a critical role in modulating the climate system. However, characterizing patterns of North Atlantic SST variability and diagnosing the associated mechanisms is challenging because they involve coupled atmosphere–ocean interactions with complex spatiotemporal relationships. Here we address these challenges by applying a time-evolving self-organizing map approach to a long preindustrial coupled control simulation and identify a variety of 10-yr spatiotemporal evolutions of winter SST anomalies, including but not limited to those associated with the North Atlantic Oscillation–Atlantic multidecadal variability (NAO–AMV)-like interactions. To assess mechanisms and atmospheric responses associated with various SST spatiotemporal evolutions, composites of atmospheric and oceanic variables associated with these evolutions are investigated. Results show that transient-eddy activities and atmospheric circulation responses exist in almost all the evolutions that are closely correlated to the details of the SST pattern. In terms of the mechanisms responsible for generating various SST evolutions, composites of ocean heat budget terms demonstrate that contributions to upper-ocean temperature tendency from resolved ocean advection and surface heat fluxes rarely oppose each other over 10-yr periods in the subpolar North Atlantic. We further explore the potential for predictability for some of these 10-yr SST evolutions that start with similar states but end with different states. However, we find that these are associated with abrupt changes in atmospheric variability and are unlikely to be predictable. In summary, this study broadly investigates the atmospheric responses to and the mechanisms governing the North Atlantic SST evolutions over 10-yr periods.
Significance Statement
Climate variability in the North Atlantic Ocean has wide-ranging impacts on global and regional climate. However, the processes involved include interactions between the ocean and atmosphere that vary across both space and time, making it challenging to characterize and predict. Using a novel machine learning approach, this study identifies various time evolutions of North Atlantic sea surface temperature patterns over 10-yr periods. This includes evolutions with similar start states but different trajectories, which have important implications for predictability. Furthermore, we investigate the mechanisms responsible for these evolutions and how different sea surface temperature patterns affect atmospheric circulation through small-scale atmospheric disturbances. These new insights into the complex ocean–atmosphere interactions over time are critical for improving decadal prediction skill.
Abstract
Previous studies have indicated that boreal winter-to-spring sea surface temperature anomalies (SSTA) over the tropical Atlantic or Indian Ocean can trigger the central-Pacific (CP) type of ENSO in the following winter due to winds over the western Pacific. Here, with the aid of observational data and CMIP5 model simulations, we demonstrate that the ability of the winter-to-spring north tropical Atlantic (NTA) SSTA or Indian Ocean Basin (IOB) mode to initiate CP ENSO events in the following winter may strongly depend on each other. Most warming events of the IOB and NTA, which are followed by CP La Niña events, are concomitant. The synergistic effect of the IOB and NTA SSTA may produce greater CP ENSO events in the subsequent winter via Walker circulation adjustments. The impacts between warming and cooling events of the IOB and NTA SSTA are asymmetric. IOB and NTA warmings appear to contribute to the subsequent CP La Niña development, which is much greater than IOB and NTA cooling contributing to CP El Niño. Overall, a combination of the IOB and NTA SSTA precursors may improve predictions of La Niña events.
Significance Statement
Although boreal winter-to-spring sea surface temperature anomalies over the tropical Atlantic or Indian Ocean can trigger central-Pacific (CP) ENSO in the following winter, it is not yet clear whether the effects of these two basins are independent. The purpose of this study is to better understand the joint effect of these two basins on CP ENSO events. We demonstrate that the ability of the north tropical Atlantic (NTA) SSTA to initiate CP ENSO events in the following winter may strongly depend on the state of the Indian Ocean Basin mode (IOB). The synergistic impact of these two basins may produce stronger CP ENSO events. These results highlight the role of three-ocean interactions in ENSO diversity and prediction.
Abstract
Previous studies have indicated that boreal winter-to-spring sea surface temperature anomalies (SSTA) over the tropical Atlantic or Indian Ocean can trigger the central-Pacific (CP) type of ENSO in the following winter due to winds over the western Pacific. Here, with the aid of observational data and CMIP5 model simulations, we demonstrate that the ability of the winter-to-spring north tropical Atlantic (NTA) SSTA or Indian Ocean Basin (IOB) mode to initiate CP ENSO events in the following winter may strongly depend on each other. Most warming events of the IOB and NTA, which are followed by CP La Niña events, are concomitant. The synergistic effect of the IOB and NTA SSTA may produce greater CP ENSO events in the subsequent winter via Walker circulation adjustments. The impacts between warming and cooling events of the IOB and NTA SSTA are asymmetric. IOB and NTA warmings appear to contribute to the subsequent CP La Niña development, which is much greater than IOB and NTA cooling contributing to CP El Niño. Overall, a combination of the IOB and NTA SSTA precursors may improve predictions of La Niña events.
Significance Statement
Although boreal winter-to-spring sea surface temperature anomalies over the tropical Atlantic or Indian Ocean can trigger central-Pacific (CP) ENSO in the following winter, it is not yet clear whether the effects of these two basins are independent. The purpose of this study is to better understand the joint effect of these two basins on CP ENSO events. We demonstrate that the ability of the north tropical Atlantic (NTA) SSTA to initiate CP ENSO events in the following winter may strongly depend on the state of the Indian Ocean Basin mode (IOB). The synergistic impact of these two basins may produce stronger CP ENSO events. These results highlight the role of three-ocean interactions in ENSO diversity and prediction.
Abstract
This study focuses on the application of two standard inflow turbulence generation methods for growing convective boundary layer (CBL) simulations: the recycle–rescale (R-R) and the digital filter–based (DF) methods, which are used in computational fluid dynamics. The primary objective of this study is to expand the applicability of the R-R method to simulations of thermally driven CBLs. This method is called the extended R-R method. However, in previous studies, the DF method has been extended to generate potential temperature perturbations. This study investigated whether the extended DF method can be applied to simulations of growing thermally driven CBLs. In this study, idealized simulations of growing thermally driven CBLs using the extended R-R and DF methods were performed. The results showed that both extended methods could capture the characteristics of thermally driven CBLs. The extended R-R method reproduced turbulence in thermally driven CBLs better than the extended DF method in the spectrum and histogram of vertical wind speed. However, the height of the thermally driven CBL was underestimated in about 100 m compared with the extended DF method. Sensitivity experiments were conducted on the parameters used in the extended DF and R-R methods. The results showed that underestimation of the length scale in the extended DF method causes a shortage of large-scale turbulence components. The other point suggested by the results of the sensitivity experiments is that the length of the driver region in the extended R-R method should be sufficient to reproduce the spanwise movement of the roll vortices.
Significance Statement
Inflow turbulence generation methods for large-eddy simulation (LES) models are crucial for the better downscaling of meteorological mesoscale models (RANS models) to microscale models (LES models). Various CFD methods have been developed, but few have been applied to simulations of thermally driven convective boundary layers (CBLs). To address this problem, we focused on a method that recycles turbulence [the recycle–rescale (R-R) method] and another method that synthetically generates turbulence [the digital filter–based (DF) method]. This study extends the R-R method to manage turbulence in thermally driven CBLs. In addition, this study investigated the applicability of the DF method to thermally driven CBL simulations. Both extended methods are effective for downscaling experiments and capture the characteristics of thermally driven CBLs.
Abstract
This study focuses on the application of two standard inflow turbulence generation methods for growing convective boundary layer (CBL) simulations: the recycle–rescale (R-R) and the digital filter–based (DF) methods, which are used in computational fluid dynamics. The primary objective of this study is to expand the applicability of the R-R method to simulations of thermally driven CBLs. This method is called the extended R-R method. However, in previous studies, the DF method has been extended to generate potential temperature perturbations. This study investigated whether the extended DF method can be applied to simulations of growing thermally driven CBLs. In this study, idealized simulations of growing thermally driven CBLs using the extended R-R and DF methods were performed. The results showed that both extended methods could capture the characteristics of thermally driven CBLs. The extended R-R method reproduced turbulence in thermally driven CBLs better than the extended DF method in the spectrum and histogram of vertical wind speed. However, the height of the thermally driven CBL was underestimated in about 100 m compared with the extended DF method. Sensitivity experiments were conducted on the parameters used in the extended DF and R-R methods. The results showed that underestimation of the length scale in the extended DF method causes a shortage of large-scale turbulence components. The other point suggested by the results of the sensitivity experiments is that the length of the driver region in the extended R-R method should be sufficient to reproduce the spanwise movement of the roll vortices.
Significance Statement
Inflow turbulence generation methods for large-eddy simulation (LES) models are crucial for the better downscaling of meteorological mesoscale models (RANS models) to microscale models (LES models). Various CFD methods have been developed, but few have been applied to simulations of thermally driven convective boundary layers (CBLs). To address this problem, we focused on a method that recycles turbulence [the recycle–rescale (R-R) method] and another method that synthetically generates turbulence [the digital filter–based (DF) method]. This study extends the R-R method to manage turbulence in thermally driven CBLs. In addition, this study investigated the applicability of the DF method to thermally driven CBL simulations. Both extended methods are effective for downscaling experiments and capture the characteristics of thermally driven CBLs.
