{"title":"Geofluid object workbench (GeoFLOW) for atmospheric dynamics in the approach to exascale: Spectral element formulation and CPU performance","authors":"D. Rosenberg, B. Flynt, M. Govett, I. Jankov","doi":"10.1175/mwr-d-22-0250.1","DOIUrl":"https://doi.org/10.1175/mwr-d-22-0250.1","url":null,"abstract":"\u0000A new software framework using a well-established high-order spectral element discretization is presented for solving the compressible Navier–Stokes equations for purposes of research in atmospheric dynamics in bounded and unbounded limited-area domains, with a view toward capturing spatiotemporal intermittency that may be particularly challenging to attain using low order schemes. A review of the discretization is provided, emphasizing properties such as the matrix product formalism and other design considerations that will facilitate its effective use on emerging exascale platforms, and a new geometry-independent, element boundary exchange method is described to maintain continuity. A variety of test problems are presented that demonstrate accuracy of the implementation primarily in wave-dominated or transitional flow regimes; conservation properties are also demonstrated. A strong scaling CPU study in a three-dimensional domain without using threading shows an average parallel efficiency of ≳ 99% up to 2×104 MPI tasks that is not affected negatively by expansion polynomial order. On-node performance is also examined and reveals that, while the primary numerical operations achieve their theoretical arithmetic intensity, the application performance is largely limited by available memory bandwidth.","PeriodicalId":18824,"journal":{"name":"Monthly Weather Review","volume":" ","pages":""},"PeriodicalIF":3.2,"publicationDate":"2023-07-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"48799795","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"A Quantile-Conserving Ensemble Filter Framework. Part II: Regression of Observation Increments in a Probit and Probability Integral Transformed Space","authors":"Jeffrey L. Anderson","doi":"10.1175/mwr-d-23-0065.1","DOIUrl":"https://doi.org/10.1175/mwr-d-23-0065.1","url":null,"abstract":"\u0000Traditional ensemble Kalman filter data assimilation methods make implicit assumptions of Gaussianity and linearity that are strongly violated by many important Earth system applications. For instance, bounded quantities like the amount of a tracer and sea ice fractional coverage cannot be accurately represented by a Gaussian which is unbounded by definition. Nonlinear relations between observations and model state variables abound. Examples include the relation between a remotely sensed radiance and the column of atmospheric temperatures, or the relation between cloud amount and water vapor quantity. Part 1 of this paper described a very general data assimilation framework for computing observation increments for non-Gaussian prior distributions and likelihoods. These methods can respect bounds and other non-Gaussian aspects of observed variables. However, these benefits can be lost when observation increments are used to update state variables using the linear regression that is part of standard ensemble Kalman filter algorithms. Here, regression of observation increments is performed in a space where variables are transformed by the probit and probability integral transforms, a specific type of Gaussian anamorphosis. This method can enforce appropriate bounds for all quantities and deal much more effectively with nonlinear relations between observations and state variables. Important enhancements like localization and inflation can be performed in the transformed space. Results are provided for idealized bivariate distributions and for cycling assimilation in a low-order dynamical system. Implications for improved data assimilation across Earth system applications are discussed.","PeriodicalId":18824,"journal":{"name":"Monthly Weather Review","volume":" ","pages":""},"PeriodicalIF":3.2,"publicationDate":"2023-07-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"47876325","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Relationships Between 10 Years of Radar-Observed Supercell Characteristics and Hail Potential","authors":"C. Homeyer, E. M. Murillo, M. Kumjian","doi":"10.1175/mwr-d-23-0019.1","DOIUrl":"https://doi.org/10.1175/mwr-d-23-0019.1","url":null,"abstract":"\u0000Supercell storms are commonly responsible for severe hail, which is the costliest severe storm hazard in the United States and elsewhere. Radar observations of such storms are common and have been leveraged to estimate hail size and severe hail occurrence. However, many established relationships between radar-observed storm characteristics and severe hail occurrence have been found using data for few storms and in isolation from other radar metrics. This study leverages a 10-year record of polarimetric Doppler radar observations in the United States to evaluate and compare radar observations of thousands of severe hail-producing supercells based on their maximum hail size. In agreement with prior studies, it is found that increasing hail size relates to increasing volume of high (≥50 dBZ) radar reflectivity, increasing mid-altitude mesocyclone rotation (azimuthal shear), increasing storm-top divergence, and decreased differential reflectivity and co-polar correlation coefficient at low levels (mostly below the environmental 0°C level). New insights include increasing vertical alignment of the storm mesocyclone with increasing hail size and a Doppler velocity spectrum width minimum aloft near storm center that increases in area with increasing hail size and is argued to indicate increasing updraft width. To complement the extensive radar analysis, near-storm environments from reanalyses are compared and indicate that the greatest environmental differences exist in the middle troposphere (within the hail growth region), especially the wind speed perpendicular to storm motion. Recommendations are given for future improvements to radar-based hail-size estimation.","PeriodicalId":18824,"journal":{"name":"Monthly Weather Review","volume":" ","pages":""},"PeriodicalIF":3.2,"publicationDate":"2023-07-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"46784160","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
W. Yanase, Udai Shimada, N. Kitabatake, Eigo Tochimoto
{"title":"Tropical Transition of Tropical Storm Kirogi (2012) over the Western North Pacific: Synoptic analysis and meso-scale simulation","authors":"W. Yanase, Udai Shimada, N. Kitabatake, Eigo Tochimoto","doi":"10.1175/mwr-d-22-0190.1","DOIUrl":"https://doi.org/10.1175/mwr-d-22-0190.1","url":null,"abstract":"\u0000Tropical transition (TT) is a cyclogenesis process in which a baroclinic disturbance is transformed into a tropical cyclone. Many studies have analyzed TT events over the North Atlantic. This study assesses TT processes from a possible subtropical cyclone to Tropical Storm Kirogi at relatively high latitude over the western North Pacific in an environment of enhanced baroclinicity in August 2012. Analyses based on satellite observation, the JRA-55 reanalysis, and a simulation with 2.5 km horizontal grid spacing demonstrate three stages during the TT: the baroclinic, intermediate, and convective stages. Over the baroclinic stage, Kirogi had an asymmetric comma-shaped cloud pattern with convection in the northern and eastern parts of the cyclone. This convection is attributed to quasi-geostrophic forcing and frontogenesis associated with advection of warm and moist air. Vorticity locally generated by this convection was advected to the cyclone center by cyclone-relative northerly flow. Kirogi also had a shallow warm-core structure due to the interaction with an upper-level cold trough extending from the mid-latitudes. In the intermediate stage, the warm and moist air in the lower troposphere and the cold trough in the upper troposphere wrapped around Kirogi. In the convective stage, Kirogi attained characteristics of a typical tropical cyclone with convection concentrated near the cyclone center and a deep warm-core structure. These results demonstrate that baroclinic processes can directly trigger formation of a tropical storm at relatively high latitudes over the western North Pacific in a similar manner to that over the North Atlantic.","PeriodicalId":18824,"journal":{"name":"Monthly Weather Review","volume":" ","pages":""},"PeriodicalIF":3.2,"publicationDate":"2023-07-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"49177734","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Benjamin W. Green, E. Sinsky, Shan Sun, V. Tallapragada, G. Grell
{"title":"Sensitivities of Subseasonal Unified Forecast System Simulations to Changes in Parameterizations of Convection, Cloud Microphysics, and Planetary Boundary Layer","authors":"Benjamin W. Green, E. Sinsky, Shan Sun, V. Tallapragada, G. Grell","doi":"10.1175/mwr-d-22-0338.1","DOIUrl":"https://doi.org/10.1175/mwr-d-22-0338.1","url":null,"abstract":"\u0000NOAA has been developing a fully-coupled Earth system model under the Unified Forecast System framework which will be responsible for global (ensemble) predictions at lead times of 0-35 days. The development has involved several prototype runs consisting of bimonthly initializations over a 7-year period for a total of 168 cases.\u0000This study leverages these existing (baseline) prototypes to isolate the impact of substituting (one-at-a-time) parameterizations for convection, microphysics, and planetary boundary layer on 35-day forecasts. Through these physics sensitivity experiments, it is found that no particular configuration of the subseasonal-length coupled model is uniformly better or worse, based on several metrics including mean-state biases and skill scores for the Madden-Julian Oscillation, precipitation, and 2-m temperature. Importantly, the spatial patterns of many “first-order” biases (e.g., impact of convection on precipitation) are remarkably similar between the end of the first week and weeks 3-4, indicating that some subseasonal biases may be mitigated through tuning at shorter timescales. This result, while shown for the first time in the context of subseasonal prediction with different physics schemes, is consistent with findings in climate models that some mean-state biases evident in multi-year averages can manifest in only a few days. An additional convective parameterization test using a different baseline shows that attempting to generalize results between or within modeling systems may be misguided. The limitations of generalizing results when testing physics schemes are most acute in modeling systems that undergo rapid, intense development from myriad contributors – as is the case in (quasi) operational environments.","PeriodicalId":18824,"journal":{"name":"Monthly Weather Review","volume":" ","pages":""},"PeriodicalIF":3.2,"publicationDate":"2023-06-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"44424108","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Kota Endo, A. Monahan, J. Bessac, H. Christensen, N. Weitzel
{"title":"Robustness of the stochastic parameterization of sub-grid scale wind variability in sea-surface fluxes","authors":"Kota Endo, A. Monahan, J. Bessac, H. Christensen, N. Weitzel","doi":"10.1175/mwr-d-22-0319.1","DOIUrl":"https://doi.org/10.1175/mwr-d-22-0319.1","url":null,"abstract":"\u0000High-resolution numerical models have been used to develop statistical models of the enhancement of sea surface fluxes resulting from spatial variability of sea-surface wind. In particular, studies have shown that the flux enhancement is not a deterministic function of the resolved state. Previous studies focused on single geographical areas or used a single high-resolution numerical model. This study extends the development of such statistical models by considering six different high-resolution models, four different geographical regions, and three different ten-day periods, allowing for a systematic investigation of the robustness of both the deterministic and stochastic parts of the data-driven parameterization. Results indicate that the deterministic part, based on regressing the unresolved normalized flux onto resolved scale normalized flux and precipitation, is broadly robust across different models, regions, and time periods. The statistical features of the stochastic part of the model (spatial and temporal autocorrelation and parameters of a Gaussian process fit to the regression residual) are also found to be robust and not strongly sensitive to the underlying model, modelled geographical region, or time period studied. Best-fit Gaussian process parameters display robust spatial heterogeneity across models, indicating potential for improvements to the statistical model. These results illustrate the potential for the development of a generic, explicitly stochastic parameterization of sea-surface flux enhancements dependent on wind variability.","PeriodicalId":18824,"journal":{"name":"Monthly Weather Review","volume":" ","pages":""},"PeriodicalIF":3.2,"publicationDate":"2023-06-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"42921001","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Matthew S. Van Den Broeke, Matthew B. Wilson, Cynthia A. Van Den Broeke, Devon J. Healey, Michael J. Wood, Raychel E. Nelson
{"title":"Polarimetric Radar Observations of a Long-lived Supercell and Associated Tornadoes on 10–11 December 2021","authors":"Matthew S. Van Den Broeke, Matthew B. Wilson, Cynthia A. Van Den Broeke, Devon J. Healey, Michael J. Wood, Raychel E. Nelson","doi":"10.1175/mwr-d-22-0330.1","DOIUrl":"https://doi.org/10.1175/mwr-d-22-0330.1","url":null,"abstract":"\u0000We present environmental and polarimetric radar observations of a long-lived December supercell which tracked approximately 750 km from Arkansas to northern Kentucky. The storm was associated with two long-track EF4 tornadoes, one of which was among the longest-tracked tornadoes recorded in the United States. The supercell’s life cycle is documented from 2000 UTC on 10 December 2021 – 0700 UTC on 11 December 2021, using data from five operational polarimetric weather radars. After convection initiation in central Arkansas, it took nearly 4 hours for a supercell to develop. Afterward, the storm’s ZDR column and arc became anomalously large leading up to genesis of the first EF4 tornado. During this time, the storm’s environment had moderate convective available potential energy (CAPE) and strong deep-layer shear. A cell interaction at about 0200 UTC disrupted the supercell updraft, weakening the ZDR arc and column and initiating the largest radar-implied hailfall event observed with this storm. The remnant circulation associated with the first EF4 tornado did not fully dissipate, and it appeared to merge with the low-level mesocyclone on the nose of a rear flank downdraft surge likely initiated by the hailfall. It is hypothesized that this merger was important to the intensification of the storm’s second EF4 tornado, which lasted nearly 3 hours and traveled approximately 267 km. During the second EF4 tornado the storm experienced decreasing CAPE and increasing storm relative helicity. Increasing interactions with other cells eventually weakened the storm, and its original updraft was obscured before the storm’s remnants dissipated in northern Kentucky.","PeriodicalId":18824,"journal":{"name":"Monthly Weather Review","volume":" ","pages":""},"PeriodicalIF":3.2,"publicationDate":"2023-06-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"49300595","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"The influence of convection initiation strength on subsequent simulated supercell evolution","authors":"Matthew D. Flournoy, E. Rasmussen","doi":"10.1175/mwr-d-22-0069.1","DOIUrl":"https://doi.org/10.1175/mwr-d-22-0069.1","url":null,"abstract":"\u0000Recent studies have shown how very small differences in the background environment of a supercell can yield different outcomes, particularly in terms of tornado production. In this study, we use a novel convection initiation (CI) technique to simulate six supercells with a focus on their early development. Each experiment is identical except for the strength of thermal forcing for the initial convection initiation. Each experiment yields a mature supercell, but differences in storm-scale characteristics like updraft speed, cold pool temperature deficit, and vertical vorticity development abound. Of these, the time when the mid-level updraft strengthens is most strongly related to initiation strength, with stronger thermal forcing favoring quicker updraft development. The same is true for the low-level updraft, with the additional relationship that stronger thermal forcing also tends to yield stronger low-level updrafts for around the first 2 hrs of the simulations. The experiments with faster updraft development tend to be associated with more rapid surface vortex intensification; however, cold pool evolution differs between simulations with weaker vs. stronger thermal forcing. Stronger thermal forcing also yields deviant, rightward storm motion earlier in the supercell’s life cycle that remains more consistent for the duration of the simulation. These results highlight the range of supercellular outcomes that are possible across a background environment due to differences in storm-scale initiation strength. They are also of potential importance for predicting the paths and tornado potential of supercells in real time.","PeriodicalId":18824,"journal":{"name":"Monthly Weather Review","volume":" ","pages":""},"PeriodicalIF":3.2,"publicationDate":"2023-06-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"45815107","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Chuntao Liu, Laufey Jörgensdóttir, Paul Walter, G. Morris, J. Flynn, P. Kucera
{"title":"On the detection of icing conditions at altitude in conjunction with mesoscale convective complexes using balloon sondes","authors":"Chuntao Liu, Laufey Jörgensdóttir, Paul Walter, G. Morris, J. Flynn, P. Kucera","doi":"10.1175/mwr-d-23-0062.1","DOIUrl":"https://doi.org/10.1175/mwr-d-23-0062.1","url":null,"abstract":"\u0000Balloon-borne radiosondes are launched twice daily at coordinated times worldwide to assist with weather forecasting. Data collection from each flight is usually terminated when the balloon bursts at an altitude above 20 km. This paper highlights cases where the balloon's turnaround occurs at lower altitudes and is associated with ice formation on the balloon, a weather condition of interest to aviation safety. Four examples of such cases are shown, where the balloon oscillates between 3-6 km altitude before rising to high altitudes and bursting. This oscillation is due to the accumulation and melting of ice on the balloon, causing the pattern to repeat multiple times. An analysis of National Weather Service radiosonde data over a five-year period and a global data set from the National Centers for Environmental Information from 1980 to 2020 identified that 0.18% of soundings worldwide satisfied these criteria. This indicates that weather conditions important to aviation safety are not rare in the worldwide database. We recommend that soundings that show descent at altitudes lower than typically expected continue to be tracked, particularly given that these up-down oscillating soundings can provide valuable information for weather forecasting on days with significant precipitation and icing conditions that might lead to aviation safety concerns.","PeriodicalId":18824,"journal":{"name":"Monthly Weather Review","volume":" ","pages":""},"PeriodicalIF":3.2,"publicationDate":"2023-06-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"42479527","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Observed air–sea turbulent heat flux anomalies during the onset of the South China Sea summer monsoon in 2021","authors":"Xiangzhou Song, Xinyue Wang, Wenbo Cai, Xuehan Xie","doi":"10.1175/mwr-d-22-0314.1","DOIUrl":"https://doi.org/10.1175/mwr-d-22-0314.1","url":null,"abstract":"\u0000This study presents observational findings of air–sea turbulent heat flux anomalies during the onset of the South China Sea summer monsoon (SCSSM) in 2021 and explains the mechanism for high-resolution heat flux variations. Turbulent heat flux discrepancies are not uniform throughout the basin but indicate a significant regional disparity in the South China Sea (SCS), which also experiences evident year-to-year variability. Based on buoy- and cruise-based air–sea measurements, high-temporal-resolution (less than hourly) anomalies in the latent heat flux during the SCSSM burst are unexpectedly determined by sea-air humidity differences instead of wind effects under near-neutral and mixed marine atmospheric boundary layer (MABL) stability conditions. However, latent heat anomalies are mainly induced by wind speed under changing MABL conditions. The sensible heat flux is much weaker, with its anomalies dominated by sea-air temperature differences regardless of the boundary layer condition. The observational results are used to examine the discrepancies in turbulent heat fluxes and associated air–sea variables in reanalysis products. The comparisons indicate that latent and sensible heat fluxes in the reanalysis are overestimated by approximately 55 Wm−2 and 3 Wm−2, respectively. These overestimations are mainly induced by higher estimates of sea-air humidity/temperature differences. The relative humidity is underestimated by approximately 4.2% in the two high-resolution reanalysis products. The higher SST (near-surface specific humidity) and lower air temperature (specific air humidity) eventually lead to higher estimates of sea-air humidity/temperature differences (1.75 g·kg−1/0.25 °C), which are the dominant factors controlling the variations in the air–sea turbulent heat fluxes.","PeriodicalId":18824,"journal":{"name":"Monthly Weather Review","volume":" ","pages":""},"PeriodicalIF":3.2,"publicationDate":"2023-06-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"48691799","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}