Nonhydrostatic Modeling System Research Articles

On the Role of the Meridional Jet and Horizontal Potential Vorticity Dipole in the Iowa Derecho of 10 August 2020

Abstract On 10 August 2020, a derecho caused widespread damage across Iowa and Illinois. Des Moines station data show that the arrival of the gust front was characterized by an abrupt shift to northerly flow, exceeding 22 m s−1 for ∼20 min. To test the hypothesis that this northerly jet is associated with a horizontal potential vorticity (PV) dipole in the lower troposphere, we investigated the structure of PV in the University of Wisconsin Nonhydrostatic Modeling System (UWNMS) and of absolute vorticity in High-Resolution Rapid Refresh (HRRR) forecast analyses. This structure is described here for the first time. The negative PV member coincides with the downdraft, while the positive PV member coincides with the updraft, with a northerly jet between. The westerly inflow jet descends anticyclonically in the downdraft, joining with northerly flow from the surface anticyclone. The resulting northerly outflow jet creates the trailing comma-shaped radar echo. The speed of propagation of the derecho is similar to the westerly wind maximum in the 3–5-km layer associated with the approaching synoptic cyclone, which acts as a steering level for resonant amplification. Idealized diagrams and 3D isosurfaces illustrate the commonality of the PV dipole/northerly jet structure. Differences in this structure among the three model states are related to low-level wind shear theory. The PV dipole coincides with the pattern of diabatic stretching tendency, which shifts westward and downward relative to the updraft/downdraft with increasing tilt. The PV dipole can contribute toward dynamical stability in a derecho. Significance Statement The purpose of this work is to investigate the structure of potential vorticity (PV) in the lower troposphere in a derecho. It is found that a northerly outflow jet occurs between an east–west-oriented horizontal PV dipole, which is described here for the first time. The negative PV member coincides with the downdraft and is inertially unstable, while the positive PV member coincides with the updraft. This work contributes toward the theory of resonant structures and longevity. The 3–5-km westerly inflow layer constitutes a steering level, which controls propagation speed despite differences in structure. The degree of westward tilt with height is related to the pattern of forcing by diabatic stretching in producing the PV dipole.

Grid refinement in ICON v2.6.4

Abstract. This article describes the implementation of grid refinement in the ICOsahedral Nonhydrostatic (ICON) modeling system. It basically follows the classical two-way nesting approach known from widely used mesoscale models like MM5 or WRF, but it differs in the way feedback from fine grids to coarser grids is applied. Moreover, the ICON implementation supports vertical nesting in the sense that the upper boundary of a nested domain may be lower than that of its parent domain. Compared to the well-established implementations on quadrilateral grids, new methods had to be developed for interpolating the lateral boundary conditions from the parent domain to the child domain(s). These are based on radial basis functions (RBFs) and partly apply direct reconstruction of the prognostic variables at the required grid points, whereas gradient-based extrapolation from parent to child grid points is used in other cases. The runtime flow control is written such that limited-area domains can be processed identically to nested domains except for the lateral boundary data supply. To demonstrate the functionality and quality of the grid nesting in ICON, idealized tests based on the Jablonowski–Williamson test case (Jablonowski and Williamson, 2006) and the Schär mountain wave test case (Schär et al., 2002) are presented. The results show that the numerical disturbances induced at the nest boundaries are small enough to be negligible for real applications. This is confirmed by experiments closely following the configuration used for operational numerical weather prediction at DWD, which demonstrate that a regional refinement over Europe has a significant positive impact on the forecast quality in the Northern Hemisphere.

Establishing a non-hydrostatic global atmospheric modeling system at 3-km horizontal resolution with aerosol feedbacks on the Sunway supercomputer of China

During the era of global warming and highly urbanized development, extreme and high impact weather as well as air pollution incidents influence everyday life and might even cause the incalculable loss of life and property. Despite the vast development of atmospheric models, there still exist substantial numerical forecast biases objectively. To accurately predict extreme weather, severe air pollution, and abrupt climate change, numerical atmospheric model requires not only to simulate meteorology and atmospheric compositions simultaneously involving many sophisticated physical and chemical processes but also at high spatiotemporal resolution. Global integrated atmospheric simulation at spatial resolutions of a few kilometers remains challenging due to its intensive computational and input/output (I/O) requirement. Through multi-dimension-parallelism structuring, aggressive and finer-grained optimizing, manual vectorizing, and parallelized I/O fragmenting, an integrated Atmospheric Model Across Scales (iAMAS) was established on the new Sunway supercomputer platform to significantly increase the computational efficiency and reduce the I/O cost. The global 3-km atmospheric simulation for meteorology with online integrated aerosol feedbacks with iAMAS was scaled to 39,000,000 processor cores and achieved the speed of 0.82 simulation day per hour (SDPH) with routine I/O, which enabled us to perform 5-day global weather forecast at 3-km horizontal resolution with online natural aerosol impacts. The results demonstrate the promising future that the increasing of spatial resolution to a few kilometers with online integrated aerosol feedbacks may significantly improve the global weather forecast.

On the Formation of Tropopause Folds and Constituent Gradient Enhancement Near Westerly Jets

AbstractThe role of differential advection in creating tropopause folds and strong constituent gradients near midlatitude westerly jets is investigated using the University of Wisconsin Non-hydrostatic Modeling System (UWNMS). Dynamical structures are compared with aircraft observations through a fold and subpolar jet (SPJ) during RF04 of the Stratosphere-Troposphere Analyses of Regional Transport (START08) campaign. The observed distribution of water vapor and ozone during RF04 provides evidence of rapid transport in the SPJ, enhancing constituent gradients above relative to below the intrusion. The creation of a tropopause fold by quasi-isentropic differential advection on the upstream side of the trough is described. This fold was created by a southward jet streak in the SPJ, where upper tropospheric air displaced the tropopause eastward in the 6-10 km layer, thereby overlying stratospheric air in the 3-6 km layer. The subsequent superposition of the subtropical and subpolar jets is also shown to result from quasi-isentropic differential advection.The occurrence of low values of ozone, water vapor, and potential vorticity on the equatorward side of the SPJ can be explained by convective transport of low-ozone air from the boundary layer, dehydration in the updraft, and detrainment of inertially-unstable air in the outflow layer. An example of rapid juxtaposition with stratospheric air in the jet core is shown for RF01. The net effect of upstream convective events is suggested as a fundamental cause of the strong constituent gradients observed in midlatitude jets. Idealized diagrams illustrate the role of differential advection in creating tropopause folds and constituent gradient enhancement.

Numerical Simulation of Tropical Cyclone Mora Using a Regional Coupled Ocean-Atmospheric Model

Tropical storms (TS) are highly sensitive to the behavior of nearby oceans, and the only way to improve their predictability is to include all the possible factors contributing to their genesis and motion. The present study evaluates the response of air–sea coupling on TS Mora, which occurred during the pre-monsoon season of 2017. Mora is simulated using a non-hydrostatic regional atmospheric modeling system coupled to a regional ocean modeling system model. The observations show a low-level Rossby wave in the atmosphere during the storm genesis. There was also evidence of increased sea surface height simultaneously in the form of a downwelling Kelvin wave at the west coast of the Bay of Bengal. The present study emphasizes the relationship between near-surface oceanic and atmospheric circulation patterns. We have tried to elucidate the process with the help of dynamical equations and model solutions. We also analyze the pivotal role of the atmospheric Rossby wave in the formation of the coastal Kelvin wave. The Kelvin wave seems to produce a warm pool required to form a low-pressure system, which further gave rise to the tropical storm Mora. The changes in the temperature and salinity profile within the ocean subsurface layers during storm passage are further used for inferences regarding the vertical structure of the upper ocean due to wind-generated mixing. The model analysis reveals that the ocean and atmospheric counterparts responsible for storm genesis are well recorded in the coupled model when compared to the observations. It was found that air–sea coupling in the model made it possible to realistically capture ocean and atmospheric responses.

On the Structure and Formation of UTLS PV Dipole/Jetlets in Tropical Cyclones by Convective Momentum Surges

AbstractThe structure and origin of mesoscale jets and associated potential vorticity (PV) dipoles in the upper troposphere and lower stratosphere (UTLS) in tropical cyclones (TCs) are investigated. UTLS PV dipole/jetlets, which occurred in Talas (2011), Edouard (2014), and Ita (2014), are simulated with the University of Wisconsin Nonhydrostatic Modeling System (UWNMS). PV dipoles are confined to the UTLS, where the jetlets oppose the ambient anticyclonic flow. They form ~100–250 km from the eye in convective asymmetries and are characterized by surges of air that accelerate in the updraft, overshoot, and extend radially outward. In these cases, the outflow jet merges with the subtropical westerly jet. Analysis of the structure of UTLS PV dipole/jetlets led to a new physical interpretation for their formation, based on the difference in momentum between the updraft and air in the UTLS: the convective momentum transport hypothesis. This view is complementary to the vorticity tilting hypothesis. A jetlet will form whenever an updraft carries horizontal winds to a level with different wind. Schematic diagrams show how to predict jetlet orientation based on horizontal speeds in the updraft and UTLS ambient air. In TCs, horizontal winds in the updraft are cyclonic, so a UTLS jetlet will be cyclonic and oppose the ambient flow. Each jetlet creates an anticyclonic, inertially unstable PV member, which lies radially outward. Estimates of terms in the PV conservation equation support the hypothesis that the dipoles arise from the curl of shear stress. Convective asymmetries associated with PV dipole/jetlets can significantly modify TC evolution by local thermodynamic acceleration.

On the Similarity of Lower-Stratospheric Potential Vorticity Dipoles above Tropical and Midlatitude Deep Convection

Abstract Simulations of the effects of deep convection on the structure of potential vorticity (PV) in the upper troposphere and lower stratosphere (UTLS) have shown that a common signature in the presence of ambient horizontal vorticity is a horizontal PV dipole. Here, the relationship between convection and PV structures in the UTLS in Tropical Cyclone Talas and the extratropical “Super Tuesday” cyclone is investigated with the University of Wisconsin Nonhydrostatic Modeling System (UWNMS). Dipoles of potential temperature in the UTLS are interpreted as an upward deflection of the ambient flow over the updraft (cold), followed by subsidence in its lee (warm), aligned with the wind direction. PV dipoles larger than ±20 PV units (1 PVU = 10−6 K kg−1 m2 s−1) are identified, with typical vertical and horizontal extents of ~3 and ~200 km, respectively, and lifetimes up to 12 h. Confirming the findings of Chagnon and Gray, it is found that horizontal PV dipoles are related to vortex tilting, where horizontally oriented vorticity associated with vertical shear of the ambient wind is bent into a horseshoe shape by the updraft, yielding a PV dipole. This suggests that theta dipoles are perpendicular to PV dipoles and that “low PV lies to the left of the wind shear,” or, in the case of tropical cyclones, “low PV lies radially outward.” Mesoscale jets occur between the dipoles, which oppose the ambient anticyclonic flow. During the extratropical transition of Talas, convective PV anomalies evolved under synoptic-scale deformation into a pair of PV streamers, which modified the midlatitude westerly jet.

On the Relationship between Inertial Instability, Poleward Momentum Surges, and Jet Intensifications near Midlatitude Cyclones

Abstract This study explores the role of inertial instability in poleward momentum surges and “flare ups” of the subpolar jet near midlatitude cyclones. Two cases are simulated with the University of Wisconsin Nonhydrostatic Modeling System to investigate the mechanisms involved in jet accelerations downstream of quasi-stationary “digging troughs.” Deep convection along the cold front leads to regions of inertial instability in the upper troposphere, which are intimately linked to jet accelerations. Terms in the zonal and meridional wind equations following the motion are evaluated for a selected air parcel within the inertially unstable region. A two-stage synoptic evolution is diagnosed, which is a characteristic signature of inertial instability. First, meridional flow accelerates following the motion, because of the subgeostrophic zonal flow and strong northward pressure gradient force (a statement of inertial instability). Second, supergeostrophic poleward flow leads to zonal acceleration and a jet flare-up. Inertial instability thus effectively displaces a westerly jet maximum poleward relative to inertially stable conditions. The structure of the poleward surge involves a distinctive “head” of high angular momentum, with the region of inertial instability enclosing the jet maximum and a core of strongly negative potential vorticity inside the surge. Departures from angular momentum–conserving profiles during meridional displacement are interpreted in terms of the pressure gradient force and degree of inertial stability. Inertial instability reduces the resulting zonal wind profile relative to angular momentum conservation but provides a significant poleward displacement of the resulting zonal wind maximum.

ICON–ART 1.0 – a new online-coupled model system from the global to regional scale

Abstract. We present the first stage of a new online-coupled global to regional-scale modeling framework for the simulation of the spatiotemporal evolution of aerosols and trace gases. The underlying meteorological model is the new nonhydrostatic model system ICON (ICOsahedral Nonhydrostatic) which allows a local grid refinement with two-way interactions between the grids. We develop the extension ART (Aerosol and Reactive Trace gases) with the goal of simulating interactions between trace substances and the state of the atmosphere. Within this paper, we present the basic equations and give an overview of the physical parameterizations as well as numerical methods we use. First applications of the new model system for trace gases, monodisperse particles and polydisperse particles are shown. The simulated distribution of two very short-lived substances (VSLS), bromoform (CHBr3) and dibromomethane (CH2Br2) reflecting the fast upward transport shows a good agreement with observations and previous model studies. Also, the shape of the simulated tropical profiles is well reproduced. As an example for the treatment of monodisperse particles we present the simulated ash plume of the Eyjafjallajökull eruption in April 2010. Here, a novel approach for the source function is applied. The pattern of the simulated distribution of volcanic ash particles shows a good agreement with previous studies. As an example for the treatment of a polydisperse aerosol, where number densities and mass concentrations are accounted for, we simulated the annual emissions of sea salt. We obtain a total emission flux of 26.0 Pg yr−1 and an emission flux of particles with diameter less than 10 μm of 7.36 Pg yr−1.

On the Role of Inertial Instability in Stratosphere–Troposphere Exchange near Midlatitude Cyclones*

Abstract In simulations of midlatitude cyclones with the University of Wisconsin Nonhydrostatic Modeling System (UWNMS), mesoscale regions with large negative absolute vorticity commonly occur in the upper troposphere and lower stratosphere (UTLS), overlying thin layers of air with stratospheric values of ozone and potential vorticity (PV). These locally enhanced stratosphere–troposphere exchange (STE) events are related to upstream convection by tracing negative equivalent potential vorticity (EPV) anomalies along back trajectories. Detailed agreement between the patterns of negative absolute vorticity, PV, and EPV—each indicators of inertial instability in the UTLS—is shown to occur in association with enhanced STE signatures. Results are presented for two midlatitude cyclones in the upper Midwest, where convection develops between the subpolar and subtropical jets. Mesoscale regions of negative EPV air originate upstream in the boundary layer. As they are transported through convection, EPV becomes increasingly negative toward the tropopause. In association with the arrival of each large negative EPV anomaly, a locally enhanced poleward surge of the subpolar jet occurs, characterized by high turbulent kinetic energy and low Richardson number. Isosurfaces of wind speed show that gravity waves emanating from inertially unstable regions connect with and modulate the subpolar and subtropical jets simultaneously. Inertially unstable convective outflow surges can facilitate STE locally by fostering poleward acceleration in the UTLS, with enhanced folding of tropospheric air over stratospheric air underneath the poleward-moving jet.

An improved dynamic core for a non-hydrostatic model system on the Yin-Yang grid

A 3D dynamic core of the non-hydrostatic model GRAPES (Global/Regional Assimilation and Prediction System) is developed on the Yin-Yang grid to address the polar problem and to enhance the computational efficiency. Three-dimensional Coriolis forcing is introduced to the new core, and full representation of the Coriolis forcing makes it straightforward to share code between the Yin and Yang subdomains. Similar to that in the original GRAPES model, a semi-implicit semi-Lagrangian scheme is adopted for temporal integration and advection with additional arrangement for cross-boundary transport. Under a non-centered second-order temporal and spatial discretization, the dry nonhydrostatic frame is summarized as the solution of an elliptical problem. The resulting Helmholtz equation is solved with the Generalized Conjugate Residual solver in cooperation with the classic Schwarz method. Even though the coefficients of the equation are quite different from those in the original model, the computational procedure of the new core is just the same. The bi-cubic Lagrangian interpolation serves to provide Dirichlet-type boundary conditions with data transfer between the subdomains. The dry core is evaluated with several benchmark test cases, and all the tests display reasonable numerical stability and computing performance. Persistency of the balanced flow and development of both the mountain-induced Rossby wave and Rossby-Haurwitz wave confirms the appropriate installation of the 3D Coriolis terms in the semi-implicit semi-Lagrangian dynamic core on the Yin-Yang grid.

Introducing Variable-Step Topography (VST) coordinates within dynamically constrained Nonhydrostatic Modeling System (NMS). Part 1: VST formulation within NMS host model framework

A Variable-Step Topography (VST) surface coordinate system is introduced into a dynamically constrained, scalable, nonhydrostatic atmospheric model for reliable simulations of flows over both smooth and steep terrain without sacrificing dynamical integrity over either type of surface. Backgrounds of both terrain-following and step coordinate model developments are presented before justifying the turn to a VST approach within an appropriately configured host model. In this first part of a two-part sequence of papers, the full formulation of the VST model, prefaced by a description of the framework of its apposite host, i.e., a re-tooled Nonhydrostatic Modeling System (NMS), are presented. [The second part assesses the performance and benefits of the new VST coordinate system in conjunction with seven orthodox obstacle flow problems.] The NMS is a 3-dimensional, nonhydrostatic cloud-mesoscale model, designed for integrations from plume-cloud scales out to regional-global scales. The derivative properties of VST in conjunction with the NMS's newly designed dynamically constrained core are capable of accurately capturing the deformations of flows by any type of terrain variability. Numerical differencing schemes needed to satisfy critical integral constraints, while also effectively enabling the VST lower boundary, are described. The host model constraints include mass, momentum, energy, vorticity and enstrophy conservation. A quasi-compressible closure cast on multiple-nest rotated spherical grids is the underlying framework used to study the advantages of the VST coordinate system. The principle objective behind the VST formulation is to combine the advantages of both terrain-following and step coordinate systems without suffering either of their disadvantages, while at the same time creating a vertical surface coordinate setting suitable for a scalable, nonhydrostatic model, safeguarded with physically realistic dynamical constraints.

Introducing variable-step topography (VST) coordinates within dynamically constrained nonhydrostatic modeling system (NMS). Part 2: VST performance on orthodox obstacle flows

In this second part of a two-part sequence of papers, the performance metrics and quantitative advantages of a new VST surface coordinate system, implemented within a dynamically constrained, nonhydrostatic, cloud mesoscale atmospheric model, are evaluated in conjunction with seven orthodox obstacle flow problems. [The first part presented a full formulation of the VST model, prefaced by a description of the framework of the newly re-tooled nonhydrostatic modeling system (NMS) operating within integral constraints based on the conservation of the foremost quantities of mass, energy and circulation.] The intent behind VST is to create a vertical surface coordinate system boundary underpinning a nonhydrostatic atmosphere capable of reliable simulations of flows over both smooth and steep terrain without sacrificing dynamical integrity over either type of surface. Model simulation results are analyzed for six classical fluid dynamics problems involving flows relative to obstacles with known analytical or laboratory-simulated solutions, as well as for a seventh noteworthy mountain wave breaking problem that has well-studied numerical solutions. For cases when topography becomes excessively severe or poorly resolved numerically, atmospheric models using transform (terrain-following) coordinates produce noteworthy errors rendering a stable integration only if the topography is smoothed. For cases when topography is slowly varying (smooth or subtle), models using discrete-step coordinates also produce noteworthy errors relative to known solutions. Alternatively, the VST model demonstrates that both limitations of the two conventional approaches, for the entire range of slope severities, can be overcome. This means that VST is ideally suited for a scalable, nonhydrostatic atmospheric model, safeguarded with physically realistic dynamical constraints.

Implementation and Validation of Dynamical Downscaling in a Microscale Simulation of a Lake Michigan Land Breeze

A nested numerical simulation of a land breeze observed over the western Lake Michigan shoreline was completed using the University of Wisconsin Nonhydrostatic Modeling System. This simulation was compared against observations obtained using a variety of instrumentation, most notably a scanning lidar system. The simulated land breeze is shown to closely represent the observations. Most notably, the influence of synoptic scale evolution on the strength and extent of the land breeze is illustrated to be a critical factor in accurately simulating the lifetime of the land breeze. Additionally, model resolution is shown to impact the ability of the simulation to accurately portray the fine scale features of the land breeze circulation, including the depth of the land breeze front and interactions between the land breeze and prevailing synoptic scale flow.

Rain-Rate Characteristics over the Korean Peninsula and Improvement of the Goddard Profiling (GPROF) Database for TMI Rainfall Retrievals

AbstractSummer rainfall characteristics over the Korean Peninsula are examined using six years of Tropical Rainfall Measuring Mission (TRMM) Microwave Imager (TMI) measurements and surface rain measurements from the densely populated rain gauges spread across South Korea. A comparison of the TMI brightness temperature at 85 GHz with the measured surface rain rate reveals that a significant portion of rainfall over the peninsula occurs at warmer brightness temperatures than would be expected from the Goddard profiling (GPROF) database. By incorporating the locally observed rain characteristics into the GPROF algorithm, efforts are made to test whether locally appropriate hydrometeor profiles may be used to improve the retrieved rainfall. Profiles are obtained by simulating rain cases using the cloud-resolving University of Wisconsin Nonhydrostatic Modeling System (UW-NMS) model and matching the calculated radar reflectivities to TRMM precipitation radar (PR) reflectivities. Selected profiles and the corresponding simulated TMI brightness temperatures (limited in this study to values that are larger than 235 K) are added to the GPROF database to form a modified database that is considered to be more suitable for local application over the Korean Peninsula. The rainfall retrieved from the new database demonstrates that heavy-rainfall events—in particular, those associated with warmer clouds—are better captured by the new algorithm as compared with the official TRMM GPROF version-6 retrievals. The results suggest that a more locally suitable rain retrieval algorithm can be developed if locally representative rain characteristics are included in the GPROF algorithm.

Modeling the effects of Southeast Asian monsoon outflow on subtropical anticyclones and midlatitude ozone over the Southern Indian Ocean

This observational and modeling study explores how the Southeast Asian summer monsoon outflow into the Southern Indian Ocean influences the life cycle of local anticyclones and leads to changes in the distribution of column ozone in the Southern Hemisphere (SH). A case study of the evolution of synoptic fields during August 1998 was performed to characterize the circulation leading to the average column ozone distribution known as the “ozone croissant,” with a characteristic ozone maximum south of Australia and a minimum near South America. During this month, five phases of SH circulation were found to lead to distinctive ozone patterns, explaining the monthly location and extent of the SH column ozone maximum. An isentropic trajectory model was used to show the cross‐equatorial flow from the Tibetan High into the SH at the near‐tropopause level of 360 K. Outflow pulses are shown to be responsible for the amplification of observed anticyclones over the SH Indian Ocean and intensification of troughs south of Australia. This couplet establishes an ozone transport pathway from tropical lower stratospheric regions around the edges of anticyclones into the ozone maximum in the amplified troughs. An idealized modeling experiment using the University of Wisconsin‐Madison Nonhydrostatic Modeling System was performed to model strong outflow pulse from tropical convection. Together with model trajectory computations, the modeling study showed a strong anticyclonic response over the Indian Ocean and increased ozone transport into the amplified troughs in a perturbed wave pattern.

Microwave single-scattering properties of randomly oriented soft-ice hydrometeors

Abstract. Large ice hydrometeors are usually present in intense convective clouds and may significantly affect the upwelling radiances that are measured by satellite-borne microwave radiometers – especially, at millimeter-wavelength frequencies. Thus, interpretation of these measurements (e.g., for precipitation retrieval) requires knowledge of the single scattering properties of ice particles. On the other hand, shape and internal structure of these particles (especially, the larger ones) is very complex and variable, and therefore it is necessary to resort to simplifying assumptions in order to compute their single-scattering parameters. In this study, we use the discrete dipole approximation (DDA) to compute the absorption and scattering efficiencies and the asymmetry factor of two kinds of quasi-spherical and non-homogeneous soft-ice particles in the frequency range 50–183 GHz. Particles of the first kind are modeled as quasi-spherical ice particles having randomly distributed spherical air inclusions. Particles of the second kind are modeled as random aggregates of ice spheres having random radii. In both cases, particle densities and dimensions are coherent with the snow hydrometeor category that is utilized by the University of Wisconsin – Non-hydrostatic Modeling System (UW-NMS) cloud-mesoscale model. Then, we compare our single-scattering results for randomly-oriented soft-ice hydrometeors with corresponding ones that make use of: a) effective-medium equivalent spheres, b) solid-ice equivalent spheres, and c) randomly-oriented aggregates of ice cylinders. Finally, we extend to our particles the scattering formulas that have been developed by other authors for randomly-oriented aggregates of ice cylinders.

Resolution enhancement for microwave‐based atmospheric sounding from geostationary orbits

The purpose of this study is to develop and evaluate techniques that improve the spatial resolution of the channels already selected in the preliminary studies for Geostationary Observatory for Microwave Atmospheric Soundings (GOMAS). Reference high resolution multifrequency brightness temperatures scenarios have been derived by applying radiative transfer calculation to the spatially and microphysically detailed output of meteorological events simulated by the University of Wisconsin–Nonhydrostatic Model System. Three approaches, Wiener filter, Super‐resolution and Image‐fusion have been applied to some representative GOMAS frequency channels to enhance the resolution of antenna temperatures. The Wiener filter improved resolution of the largely oversampled images by a factor 1.5–2.0 without introducing any penalty in the radiometric accuracy. Super‐resolution, suitable for not largely oversampled images, improved resolution by a factor ∼1.5 but introducing an increased radiometric noise by a factor 1.4–2.5. The Image‐fusion allows finally to further increase the spatial frequency of the images obtained by the Wiener filter increasing the total resolution up to a factor 5.0 with an increased radiometric noise closely linked to the radiometric frequency and to the examined case study.

Long‐range convective ozone transport during INTEX

Ozone transport across the tropopause near a convective region in the central Pacific is calculated for a 1‐day period using a nested grid, high‐resolution simulation with the University of Wisconsin Nonhydrostatic Modeling System (UWNMS). This ozone is then tracked into the region where DC‐8 lidar ozone observations on 6 July 2004 were taken over the United States during the Intercontinental Chemical Transport Experiment‐North America (INTEX‐NA). On 2 July a convective complex developed north of the Midway Islands and intensified as it propagated eastward. Satellite observations captured this event well, and the UWNMS simulation placed the convection just east of an upper level trough along the subtropical jet. The proximity of the strong convection to the lowered tropopause resulted in significant stratospheric ozone contribution (SOC) to the troposphere. During the 24 h period starting at 1200 UTC 2 July, total net ozone flux into the troposphere was calculated to be about 0.2 Tg across the nested grid domain. Because the convective outflow was feeding the subtropical jet, the SOC was quickly transported across the eastern Pacific and into the INTEX‐NA region. Using an “SOC tracer” we show that a small but significant percentage of observed ozone on the DC8 flight can be attributed to SOC from this one convective event. This supports the idea that the stratosphere over the Pacific is an important source of ozone for air entering North America from the west.

On the potential of sub-mm passive MW observations from geostationary satellites to retrieve heavy precipitation over the Mediterranean Area

Abstract. The general interest in the potential use of the mm and sub-mm frequencies up to 425 GHz resolution from geostationary orbit is increasing due to the fact that the frequent time sampling and the comparable spatial resolution relative to the "classical" (≤89 GHz) microwave frequencies would allow the monitoring of precipitating intense events for the assimilation of rain in now-casting weather prediction models. In this paper, we use the simulation of a heavy precipitating event in front of the coast of Crete island (Greece) performed by the University of Wisconsin - Non-hydrostatic Modeling System (UW-NMS) cloud resolving model in conjunction with a 3D-adjusted plane parallel radiative transfer model to simulate the upwelling brightness temperatures (TB's) at mm and sub-mm frequencies. To study the potential use of high frequencies, we first analyze the relationships of the simulated TB's with the microphysical properties of the UW-NMS simulated precipitating clouds, and then explore the capability of a Bayesian algorithm for the retrieval of surface rain rate, rain and ice water paths at such frequencies.