Abstract
Severe storms produce hazardous weather phenomena, such as large hail, damaging winds, and tornadoes. However, relationships between convective parameters and confirmed severe weather occurrences are poorly quantified in south-central Brazil. This study explores severe weather reports and measurements from newly available datasets. Hail, damaging wind, and tornado reports are sourced from the PREVOTS project from June 2018 to December 2021, while measurements of convectively induced wind gusts from 1996 to 2019 are obtained from METAR reports and from Brazilās operational network of automated weather stations. Proximal convective parameters were computed from ERA5 reanalysis for these reports and used to perform a discriminant analysis using mixed-layer CAPE and deep-layer shear (DLS). Compared to other regions, thermodynamic parameters associated with severe weather episodes exhibit lower magnitudes in south-central Brazil. DLS displays better performance in distinguishing different types of hazardous weather, but does not discriminate well between distinct severity levels. To address the sensitivity of the discriminant analysis to distinct environmental regimes and hazard types, five different discriminants are assessed. These include discriminants for any severe storm, severe hail only, severe wind gust only, and all environments but broken into āhighā and ālowā CAPE regimes. The best performance of the discriminant analysis is found for the āhighā CAPE regime, followed by the severe wind regime. All discriminants demonstrate that DLS plays a more important role in conditioning Brazilian severe storm environments than other regions, confirming the need to ensure that parameters and discriminants are tuned to local severe weather conditions.
Abstract
Severe storms produce hazardous weather phenomena, such as large hail, damaging winds, and tornadoes. However, relationships between convective parameters and confirmed severe weather occurrences are poorly quantified in south-central Brazil. This study explores severe weather reports and measurements from newly available datasets. Hail, damaging wind, and tornado reports are sourced from the PREVOTS project from June 2018 to December 2021, while measurements of convectively induced wind gusts from 1996 to 2019 are obtained from METAR reports and from Brazilās operational network of automated weather stations. Proximal convective parameters were computed from ERA5 reanalysis for these reports and used to perform a discriminant analysis using mixed-layer CAPE and deep-layer shear (DLS). Compared to other regions, thermodynamic parameters associated with severe weather episodes exhibit lower magnitudes in south-central Brazil. DLS displays better performance in distinguishing different types of hazardous weather, but does not discriminate well between distinct severity levels. To address the sensitivity of the discriminant analysis to distinct environmental regimes and hazard types, five different discriminants are assessed. These include discriminants for any severe storm, severe hail only, severe wind gust only, and all environments but broken into āhighā and ālowā CAPE regimes. The best performance of the discriminant analysis is found for the āhighā CAPE regime, followed by the severe wind regime. All discriminants demonstrate that DLS plays a more important role in conditioning Brazilian severe storm environments than other regions, confirming the need to ensure that parameters and discriminants are tuned to local severe weather conditions.
Abstract
Radiative feedbacks over interannual timescales can be potentially useful for global warming estimation. However, the diversity of the leadālag relationships in global mean surface temperature (GMST) and net radiation flux at the top of the atmosphere (GMTOA) create uncertainty during the estimation of radiative feedbacks. In this study, key physical processes controlling leadālag relationships were elucidated by categorizing preindustrial control simulations of CMIP6 into three groups based on cross correlation values of GMTOA against GMST at Lag 0 and Lag +1 year. The diversity in the leadālag relationships was primarily caused by the climatological state difference of the atmosphere over the equatorial Pacific, which modulated the strength of convective activity and sensitivity of low-level clouds. Diminished atmospheric stability caused enhanced convective activity, more efficient energy release, and smaller lags. In addition, enhanced stability in the lower atmosphere rendered the low-level clouds more sensitive to sea surface temperature changes and considerably delayed the radiative response. The climatological state difference of the atmosphere resulted from model-inherent atmospheric conditions. These findings suggest that the diversity of leadālag relationships of GMST and GMTOA over interannual timescales could represent the characteristics of general atmospheric circulation models and possible solutions of the actual atmosphere, which could also affect long-term feedback features.
Abstract
Radiative feedbacks over interannual timescales can be potentially useful for global warming estimation. However, the diversity of the leadālag relationships in global mean surface temperature (GMST) and net radiation flux at the top of the atmosphere (GMTOA) create uncertainty during the estimation of radiative feedbacks. In this study, key physical processes controlling leadālag relationships were elucidated by categorizing preindustrial control simulations of CMIP6 into three groups based on cross correlation values of GMTOA against GMST at Lag 0 and Lag +1 year. The diversity in the leadālag relationships was primarily caused by the climatological state difference of the atmosphere over the equatorial Pacific, which modulated the strength of convective activity and sensitivity of low-level clouds. Diminished atmospheric stability caused enhanced convective activity, more efficient energy release, and smaller lags. In addition, enhanced stability in the lower atmosphere rendered the low-level clouds more sensitive to sea surface temperature changes and considerably delayed the radiative response. The climatological state difference of the atmosphere resulted from model-inherent atmospheric conditions. These findings suggest that the diversity of leadālag relationships of GMST and GMTOA over interannual timescales could represent the characteristics of general atmospheric circulation models and possible solutions of the actual atmosphere, which could also affect long-term feedback features.
Abstract
Recent observations and numerical simulations have demonstrated the potential for significant interactions between mesoscale eddies and smaller-scale tidally generated internal waves ā also known as internal tides. Here we develop a simple theoretical model that predicts the one-way upscale transfer of energy from internal tides to mesoscale eddies through a critical level mechanism. We find that ā in the presence of a critical level ā the internal tide energy flux into an eddy is partitioned according to the wave frequency Ī© and local inertial frequency f : a fraction of 1 ā f /Ī© is transferred to the eddy kinetic energy while the remainder is viscously dissipated or supports mixing. These predictions are validated by comparison with a suite of numerical simulations. The simulations further show that the wave-driven energisation of the eddies also accelerates the onset of hydrodynamical instabilities and the break down of the eddies, thereby increasing eddy kinetic energy, but reducing eddy lifetimes. Our estimates suggest that in regions of the ocean with both significant eddy fields and internal tidesāsuch as parts of the Gulf Stream and Antarctic Circumpolar Currentāthe critical level effect could drive a ā¼10% per month increase in the kinetic energy of a typical eddy. Our results provide a basis for parameterising internal tide-eddy interactions in global ocean models where they are currently unrepresented.
Abstract
Recent observations and numerical simulations have demonstrated the potential for significant interactions between mesoscale eddies and smaller-scale tidally generated internal waves ā also known as internal tides. Here we develop a simple theoretical model that predicts the one-way upscale transfer of energy from internal tides to mesoscale eddies through a critical level mechanism. We find that ā in the presence of a critical level ā the internal tide energy flux into an eddy is partitioned according to the wave frequency Ī© and local inertial frequency f : a fraction of 1 ā f /Ī© is transferred to the eddy kinetic energy while the remainder is viscously dissipated or supports mixing. These predictions are validated by comparison with a suite of numerical simulations. The simulations further show that the wave-driven energisation of the eddies also accelerates the onset of hydrodynamical instabilities and the break down of the eddies, thereby increasing eddy kinetic energy, but reducing eddy lifetimes. Our estimates suggest that in regions of the ocean with both significant eddy fields and internal tidesāsuch as parts of the Gulf Stream and Antarctic Circumpolar Currentāthe critical level effect could drive a ā¼10% per month increase in the kinetic energy of a typical eddy. Our results provide a basis for parameterising internal tide-eddy interactions in global ocean models where they are currently unrepresented.
Abstract
Based on the Conditional Nonlinear Optimal Perturbation for boundary condition method and Regional Ocean Modeling System (ROMS), this study investigates the influence of wind stress uncertainty on predicting the short-term state transitions of the Kuroshio Extension (KE). The optimal time-dependent wind stress errors that lead to maximum prediction errors are obtained for two KE stable-to-unstable and two reverse transitions, which exhibit local multi-eddies structures with decreasing magnitude as the end time of prediction approaches. The optimal boundary errors initially induce small oceanic errors through Ekman pumping. Subsequently, these errors grow in magnitude as oceanic internal processes take effect, which exerts significant influences on the short-term prediction of the KE state transition process. Specifically, during stable-to-unstable (unstable-to-stable) transitions, the growing error induces an overestimation (underestimation) of the meridional sea surface height gradient across the KE axis, leading to the predicted KE state being more (less) stable. Furthermore, the dynamics mechanism analysis indicates that barotropic instability is crucial for the error growth in the prediction of both the stable-to-unstable and the reverse transition processes due to the horizontal shear of flow field. But work generated by wind stress error plays a more important role in the prediction of the unstable-to-stable transitions because of the synergistic effect of strong wind stress error and strong oceanic error. Eventually, the sensitive areas have been identified based on the optimal boundary errors. Reducing wind stress errors in sensitive areas can significantly improve prediction skills, offering theoretical guidance for devising observational strategies.
Abstract
Based on the Conditional Nonlinear Optimal Perturbation for boundary condition method and Regional Ocean Modeling System (ROMS), this study investigates the influence of wind stress uncertainty on predicting the short-term state transitions of the Kuroshio Extension (KE). The optimal time-dependent wind stress errors that lead to maximum prediction errors are obtained for two KE stable-to-unstable and two reverse transitions, which exhibit local multi-eddies structures with decreasing magnitude as the end time of prediction approaches. The optimal boundary errors initially induce small oceanic errors through Ekman pumping. Subsequently, these errors grow in magnitude as oceanic internal processes take effect, which exerts significant influences on the short-term prediction of the KE state transition process. Specifically, during stable-to-unstable (unstable-to-stable) transitions, the growing error induces an overestimation (underestimation) of the meridional sea surface height gradient across the KE axis, leading to the predicted KE state being more (less) stable. Furthermore, the dynamics mechanism analysis indicates that barotropic instability is crucial for the error growth in the prediction of both the stable-to-unstable and the reverse transition processes due to the horizontal shear of flow field. But work generated by wind stress error plays a more important role in the prediction of the unstable-to-stable transitions because of the synergistic effect of strong wind stress error and strong oceanic error. Eventually, the sensitive areas have been identified based on the optimal boundary errors. Reducing wind stress errors in sensitive areas can significantly improve prediction skills, offering theoretical guidance for devising observational strategies.
Abstract
This study investigates the combined impacts of the Madden-Julian oscillation (MJO) and extratropical anticyclonic Rossby wave breaking (AWB) on sub-seasonal Atlantic tropical cyclone (TC) activity and their physical connections. Our results show that during MJO phases 2-3 (enhanced Indian Ocean convection) and 6-7 (enhanced tropical Pacific convection), there are significant changes in basin-wide TC activity. The MJO and AWB collaborate to suppress basin-wide TC activity during phases 6-7 but not during phases 2-3. During phases 6-7, when AWB occurs, various TC metrics including hurricanes, accumulated cyclone energy and rapid intensification probability decrease by ~50%-80%. Simultaneously, large-scale environmental variables, like vertical wind shear, precipitable water, and sea surface temperatures become extremely unfavorable for TC formation and intensification, compared to periods characterized by suppressed AWB activity during the same MJO phases. Further investigation reveals that AWB events during phases 6-7 occur in concert with the development of a stronger anticyclone in the lower troposphere, which transports more dry, stable extratropical air equatorward, and drives enhanced tropical SST cooling. As a result, individual AWB events in phases 6-7 can disturb the development of surrounding TCs to a greater extent than their phases 2-3 counterparts. The influence of the MJO on AWB over the western subtropical Atlantic can be attributed to the modulation of the convectively-forced Rossby wave source over the tropical eastern Pacific. A significant number of Rossby waves initiated from this region during phases 5-6 propagate into the subtropical North Atlantic, preceding the occurrence of AWB events in phases 6-7.
Abstract
This study investigates the combined impacts of the Madden-Julian oscillation (MJO) and extratropical anticyclonic Rossby wave breaking (AWB) on sub-seasonal Atlantic tropical cyclone (TC) activity and their physical connections. Our results show that during MJO phases 2-3 (enhanced Indian Ocean convection) and 6-7 (enhanced tropical Pacific convection), there are significant changes in basin-wide TC activity. The MJO and AWB collaborate to suppress basin-wide TC activity during phases 6-7 but not during phases 2-3. During phases 6-7, when AWB occurs, various TC metrics including hurricanes, accumulated cyclone energy and rapid intensification probability decrease by ~50%-80%. Simultaneously, large-scale environmental variables, like vertical wind shear, precipitable water, and sea surface temperatures become extremely unfavorable for TC formation and intensification, compared to periods characterized by suppressed AWB activity during the same MJO phases. Further investigation reveals that AWB events during phases 6-7 occur in concert with the development of a stronger anticyclone in the lower troposphere, which transports more dry, stable extratropical air equatorward, and drives enhanced tropical SST cooling. As a result, individual AWB events in phases 6-7 can disturb the development of surrounding TCs to a greater extent than their phases 2-3 counterparts. The influence of the MJO on AWB over the western subtropical Atlantic can be attributed to the modulation of the convectively-forced Rossby wave source over the tropical eastern Pacific. A significant number of Rossby waves initiated from this region during phases 5-6 propagate into the subtropical North Atlantic, preceding the occurrence of AWB events in phases 6-7.
Abstract
National Weather Service (NWS) forecasters have many roles and responsibilities, including communication with core partners throughout the forecast and warning process to ensure that the information they are providing is relevant, understandable, and actionable. While the NWS communicates to many groups, members of the emergency management community are among the most critical partners. However, little is known about the diverse population of emergency managers (EMs) and how they receive, process, and use forecast information. The Extreme Weather and Emergency Management Survey (WxEM) aims to fill this knowledge gap by (1) building a nationwide panel of EMs and (2) fielding routine surveys that include questions of relevance to NWS operations. The panel was built by creating a database with contact information from more than 4,000 EMs across the country. An enrollment survey was sent to the list, and over 700 EMs agreed to participate in the project. Following enrollment, WxEM panelists receive surveys that address how EMs use NWS forecast information three to four times a year. These surveys cover a variety of subjects, with the goal of working with other researchers to develop surveys that address their research needs. By collaborating with other research groups to design short, focused surveys, the WxEM project will reduce the research burden on EMs and, at the same time, increase the quality and comparability of research data in the weather enterprise. The results will be shared with the NWS and the research community, and all data gathered from these surveys will be publicly available.
Abstract
National Weather Service (NWS) forecasters have many roles and responsibilities, including communication with core partners throughout the forecast and warning process to ensure that the information they are providing is relevant, understandable, and actionable. While the NWS communicates to many groups, members of the emergency management community are among the most critical partners. However, little is known about the diverse population of emergency managers (EMs) and how they receive, process, and use forecast information. The Extreme Weather and Emergency Management Survey (WxEM) aims to fill this knowledge gap by (1) building a nationwide panel of EMs and (2) fielding routine surveys that include questions of relevance to NWS operations. The panel was built by creating a database with contact information from more than 4,000 EMs across the country. An enrollment survey was sent to the list, and over 700 EMs agreed to participate in the project. Following enrollment, WxEM panelists receive surveys that address how EMs use NWS forecast information three to four times a year. These surveys cover a variety of subjects, with the goal of working with other researchers to develop surveys that address their research needs. By collaborating with other research groups to design short, focused surveys, the WxEM project will reduce the research burden on EMs and, at the same time, increase the quality and comparability of research data in the weather enterprise. The results will be shared with the NWS and the research community, and all data gathered from these surveys will be publicly available.
Abstract
Because flash drought is a relatively new phenomenon in drought research, defining the concept is critical for scientists and decision-makers. Having detrimental impacts on many sectors, it is important to have a consistent definition and understanding of flash drought, between experts and stakeholders, to provide early warning to the community. This study focuses on onset and progression of conditions and demonstrates the difference in flash drought identification for 15 events across six quantitative definitions of flash drought that use the U.S. Drought Monitor (USDM). Five flash drought events have been studied in the literature while 10 additional events have been perceived as flash drought by stakeholders. The results show that two of six definitions consistently capture the earliest onset of flash drought and include a large percent area in the identification. A qualitative analysis of management challenges and needs determined by stakeholders was completed using survey data. The results found that managing impacts and better communication and education were the top challenges and more data and enhanced and efficient communication as the top needs to better monitor, manage, and respond to flash droughts. The results demonstrate the need for assessing the characteristics of the definitions to enhance communication and monitoring strategies for large and small-scale flash droughts.
Significance Statement
The purpose of this study is to better understand how different numerical flash drought definitions characterize multiple flash drought events and how these definitions are useful in addressing the needs and challenges of stakeholders. This is important because definitions may capture different areas in flash droughts, which can impact how end users identify a flash drought. Further, this study uses events identified by the literature and by people familiar with drought monitoring. From these findings, definitions that capture flash drought earliest would help address the challenge of rapid onset and the need of quicker data. Further, definitions by sector would be beneficial to address the scale of impacts. This study identifies the importance of definitions for early warning systems.
Abstract
Because flash drought is a relatively new phenomenon in drought research, defining the concept is critical for scientists and decision-makers. Having detrimental impacts on many sectors, it is important to have a consistent definition and understanding of flash drought, between experts and stakeholders, to provide early warning to the community. This study focuses on onset and progression of conditions and demonstrates the difference in flash drought identification for 15 events across six quantitative definitions of flash drought that use the U.S. Drought Monitor (USDM). Five flash drought events have been studied in the literature while 10 additional events have been perceived as flash drought by stakeholders. The results show that two of six definitions consistently capture the earliest onset of flash drought and include a large percent area in the identification. A qualitative analysis of management challenges and needs determined by stakeholders was completed using survey data. The results found that managing impacts and better communication and education were the top challenges and more data and enhanced and efficient communication as the top needs to better monitor, manage, and respond to flash droughts. The results demonstrate the need for assessing the characteristics of the definitions to enhance communication and monitoring strategies for large and small-scale flash droughts.
Significance Statement
The purpose of this study is to better understand how different numerical flash drought definitions characterize multiple flash drought events and how these definitions are useful in addressing the needs and challenges of stakeholders. This is important because definitions may capture different areas in flash droughts, which can impact how end users identify a flash drought. Further, this study uses events identified by the literature and by people familiar with drought monitoring. From these findings, definitions that capture flash drought earliest would help address the challenge of rapid onset and the need of quicker data. Further, definitions by sector would be beneficial to address the scale of impacts. This study identifies the importance of definitions for early warning systems.
Abstract
Ensemble data assimilation (DA) is an efficient approach to reduce snow simulation errors by combining observation and land surface modeling. However, there is a small spread between ensemble members of simulated snowpack, which typically occurs for a long time with 100% snow cover fraction (SCF) or snow-free conditions. Here, we apply a hybrid DA method, in which direct insertion (DI) is a supplement of the ensemble square root filter (EnSRF), to assimilate the spaceborne SCF into a land surface model, driven by China Meteorological Administration Land Data Assimilation System high-resolution climate forcings over northern China during the snow season in 2021/22. Compared to the open-loop experiment (without SCF assimilation), the root-mean-square error (RMSE) of SCF is reduced by 6% through the original EnSRF and is even lower (by 14%) in the combined DI and EnSRF (EnSRFDI) experiment. The results reveal the ability of both EnSRF and EnSRFDI to improve the SCF estimation over regions where the snow cover is low, while only EnSRFDI is able to efficiently reduce the RMSE over areas with high SCF. Moreover, the SCF assimilation is also observed to improve the snow depth and soil temperature simulations, with the KlingāGupta efficiency (KGE) increasing at 60% and 56%ā70% stations, respectively, particularly under conditions with near-freezing temperature, in which reliable simulations are typically challenging. Our results demonstrate that the EnSRFDI hybrid method can be applied for the assimilation of spaceborne observational snow cover to improve land surface simulations and snow-related operational products.
Significance Statement
Due to the small spread between the seasonal snowpack of ensemble simulations, ensemble snow cover fraction (SCF) data assimilation (DA) proves to be ineffective. Therefore, we apply a hybrid method that combines the direct insertion (DI) and ensemble square root filter (EnSRF) to assimilate the spaceborne SCF into a land surface model (LSM) driven by high-resolution climate forcings. Our results reveal the applicability of the EnSRFDI to further improve snow cover simulations over regions with high SCF. Furthermore, the DA experiments were validated through a large number of in situ observations from the China Meteorological Administration. The uncertainties of snow depth and soil temperature simulations are also slightly reduced by the SCF DAs, particularly over regions with a poor LSM performance.
Abstract
Ensemble data assimilation (DA) is an efficient approach to reduce snow simulation errors by combining observation and land surface modeling. However, there is a small spread between ensemble members of simulated snowpack, which typically occurs for a long time with 100% snow cover fraction (SCF) or snow-free conditions. Here, we apply a hybrid DA method, in which direct insertion (DI) is a supplement of the ensemble square root filter (EnSRF), to assimilate the spaceborne SCF into a land surface model, driven by China Meteorological Administration Land Data Assimilation System high-resolution climate forcings over northern China during the snow season in 2021/22. Compared to the open-loop experiment (without SCF assimilation), the root-mean-square error (RMSE) of SCF is reduced by 6% through the original EnSRF and is even lower (by 14%) in the combined DI and EnSRF (EnSRFDI) experiment. The results reveal the ability of both EnSRF and EnSRFDI to improve the SCF estimation over regions where the snow cover is low, while only EnSRFDI is able to efficiently reduce the RMSE over areas with high SCF. Moreover, the SCF assimilation is also observed to improve the snow depth and soil temperature simulations, with the KlingāGupta efficiency (KGE) increasing at 60% and 56%ā70% stations, respectively, particularly under conditions with near-freezing temperature, in which reliable simulations are typically challenging. Our results demonstrate that the EnSRFDI hybrid method can be applied for the assimilation of spaceborne observational snow cover to improve land surface simulations and snow-related operational products.
Significance Statement
Due to the small spread between the seasonal snowpack of ensemble simulations, ensemble snow cover fraction (SCF) data assimilation (DA) proves to be ineffective. Therefore, we apply a hybrid method that combines the direct insertion (DI) and ensemble square root filter (EnSRF) to assimilate the spaceborne SCF into a land surface model (LSM) driven by high-resolution climate forcings. Our results reveal the applicability of the EnSRFDI to further improve snow cover simulations over regions with high SCF. Furthermore, the DA experiments were validated through a large number of in situ observations from the China Meteorological Administration. The uncertainties of snow depth and soil temperature simulations are also slightly reduced by the SCF DAs, particularly over regions with a poor LSM performance.
Abstract
The NASA Soil Moisture Active Passive (SMAP) mission Level-4 Soil Moisture (L4_SM) product provides global, 9-km resolution, 3-hourly surface and root-zone soil moisture from April 2015 to the present with a mean latency of 2.5 days from the time of observation. The L4_SM algorithm assimilates SMAP L-band (1.4 GHz) brightness temperature (Tb) observations into the NASA Catchment land surface model as the model is driven with observation-based precipitation. This paper describes and evaluates the use of satellite- and gauge-based precipitation from the NASA Integrated Multi-satellitE Retrievals for the Global Precipitation Measurement (IMERG) products in the L4_SM algorithm beginning with L4_SM Version 6. Specifically, IMERG is used in two ways: (i) The L4_SM precipitation reference climatology is primarily based on IMERG-Final (Version 06B) data, replacing the Global Precipitation Climatology Project Version 2.2 data used in previous L4_SM versions, and (ii) the precipitation forcing outside of North America and the high latitudes is corrected to match the daily totals from IMERG, replacing the gauge-only, daily product or uncorrected weather analysis precipitation used there in earlier L4_SM versions. The use of IMERG precipitation improves the anomaly time series correlation coefficient of L4_SM surface soil moisture (versus independent satellite estimates) by 0.03 in the global average and by up to ā¼0.3 in parts of South America, Africa, Australia, and East Asia, where the quality of the gauge-only precipitation product used in earlier L4_SM versions is poor. The improvements also reduce the time series standard deviation of the Tb observation-minus-forecast residuals from 5.5 K in L4_SM Version 5 to 5.1 K in Version 6.
Significance Statement
Soil moisture links the land surface water, energy, and carbon cycles. NASA Soil Moisture Active Passive (SMAP) satellite observations and observation-based precipitation data are merged into a numerical model of land surface water and energy balance processes to generate the global, 9-km resolution, 3-hourly Level-4 Soil Moisture (L4_SM) data product. The product is available with ā¼2.5-day latency to support Earth science research and applications, such as flood prediction and drought monitoring. We show that a recent L4_SM algorithm update using satellite- and gauge-based precipitation inputs from the NASA Integrated Multi-satellitE Retrievals for the Global Precipitation Measurement (IMERG) products improves the temporal variations in the estimated soil moisture, particularly in otherwise poorly instrumented regions in South America, Africa, Australia, and East Asia.
Abstract
The NASA Soil Moisture Active Passive (SMAP) mission Level-4 Soil Moisture (L4_SM) product provides global, 9-km resolution, 3-hourly surface and root-zone soil moisture from April 2015 to the present with a mean latency of 2.5 days from the time of observation. The L4_SM algorithm assimilates SMAP L-band (1.4 GHz) brightness temperature (Tb) observations into the NASA Catchment land surface model as the model is driven with observation-based precipitation. This paper describes and evaluates the use of satellite- and gauge-based precipitation from the NASA Integrated Multi-satellitE Retrievals for the Global Precipitation Measurement (IMERG) products in the L4_SM algorithm beginning with L4_SM Version 6. Specifically, IMERG is used in two ways: (i) The L4_SM precipitation reference climatology is primarily based on IMERG-Final (Version 06B) data, replacing the Global Precipitation Climatology Project Version 2.2 data used in previous L4_SM versions, and (ii) the precipitation forcing outside of North America and the high latitudes is corrected to match the daily totals from IMERG, replacing the gauge-only, daily product or uncorrected weather analysis precipitation used there in earlier L4_SM versions. The use of IMERG precipitation improves the anomaly time series correlation coefficient of L4_SM surface soil moisture (versus independent satellite estimates) by 0.03 in the global average and by up to ā¼0.3 in parts of South America, Africa, Australia, and East Asia, where the quality of the gauge-only precipitation product used in earlier L4_SM versions is poor. The improvements also reduce the time series standard deviation of the Tb observation-minus-forecast residuals from 5.5 K in L4_SM Version 5 to 5.1 K in Version 6.
Significance Statement
Soil moisture links the land surface water, energy, and carbon cycles. NASA Soil Moisture Active Passive (SMAP) satellite observations and observation-based precipitation data are merged into a numerical model of land surface water and energy balance processes to generate the global, 9-km resolution, 3-hourly Level-4 Soil Moisture (L4_SM) data product. The product is available with ā¼2.5-day latency to support Earth science research and applications, such as flood prediction and drought monitoring. We show that a recent L4_SM algorithm update using satellite- and gauge-based precipitation inputs from the NASA Integrated Multi-satellitE Retrievals for the Global Precipitation Measurement (IMERG) products improves the temporal variations in the estimated soil moisture, particularly in otherwise poorly instrumented regions in South America, Africa, Australia, and East Asia.
Abstract
There is growing use of machine learning algorithms to replicate sub-grid parametrisation schemes in global climate models. Parametrisations rely on approximations, thus there is potential for machine learning to aid improvements. In this study, a neural network is used to mimic the behaviour of the non-orographic gravity wave scheme used in the Met Office climate model, important for stratospheric climate and variability. The neural network is found to require only two of the six inputs used by the parametrisation scheme, suggesting the potential for greater efficiency in this scheme. Use of a one-dimensional mechanistic model is advocated, allowing neural network hyperparameters to be chosen based on emergent features of the coupled system with minimal computational cost, and providing a test bed prior to coupling to a climate model. A climate model simulation, using the neural network in place of the existing parametrisation scheme, is found to accurately generate a quasi-biennial oscillation of the tropical stratospheric winds, and correctly simulate the non-orographic gravity wave variability associated with the El NiƱo Southern Oscillation and stratospheric polar vortex variability. These internal sources of variability are essential for providing seasonal forecast skill, and the gravity wave forcing associated with them is reproduced without explicit training for these patterns.
Abstract
There is growing use of machine learning algorithms to replicate sub-grid parametrisation schemes in global climate models. Parametrisations rely on approximations, thus there is potential for machine learning to aid improvements. In this study, a neural network is used to mimic the behaviour of the non-orographic gravity wave scheme used in the Met Office climate model, important for stratospheric climate and variability. The neural network is found to require only two of the six inputs used by the parametrisation scheme, suggesting the potential for greater efficiency in this scheme. Use of a one-dimensional mechanistic model is advocated, allowing neural network hyperparameters to be chosen based on emergent features of the coupled system with minimal computational cost, and providing a test bed prior to coupling to a climate model. A climate model simulation, using the neural network in place of the existing parametrisation scheme, is found to accurately generate a quasi-biennial oscillation of the tropical stratospheric winds, and correctly simulate the non-orographic gravity wave variability associated with the El NiƱo Southern Oscillation and stratospheric polar vortex variability. These internal sources of variability are essential for providing seasonal forecast skill, and the gravity wave forcing associated with them is reproduced without explicit training for these patterns.
