Search NASA⌕ Search

SEARCH · Search NASA

Results for “nonhydrostatic”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

89 records · Page 5

GMS-based"Future Time" Rainfall Data Assimilation for Mesoscale Weather Prediction over Korean Peninsula and Future Prospects with GPM Satellite Measurements

This study examines the use of satellite-derived nowcasted (short-term forecasted) rainfall over 3-hour time periods to gain an equivalent time increment in initializing a nonhydrostatic mesoscale model used for predicting convective rainfall events over the Korean peninsula. Infrared (IR) window measurements from the Japanese Geostationary Meteorological Satellite (GMS) are used to specify latent heating for a spinup period of the model - but in future time -- thus initializing in advance of actual time in the framework of a prediction scenario. The main scientific objective of the study is to investigate the strengths and weaknesses of this approach insofar as data assimilation, in which the nowcasted assimilation data are derived independently of the prognostic model itself. Although there have been various recent improvements in formulating the dynamics, thermodynamics, and microphysics of mesoscale models, as well as computer advances which allow the use of high resolution cloud-resolving grids and explicit latent heating over regional domains, spinup remains at the forefront of unresolved mesoscale modeling problems. In general, non-realistic spinup limits the skill in predicting the spatial-temporal distribution of convection and precipitation, primarily in the early hours of a. forecast, stemming from standard prognostic variables not representing the initial diabatic heating field produced by the ambient convection and cloud fields. The long-term goal of this research is to improve short-range (12-hour) quantitative precipitation forecasting (QPF) over the Korean peninsula through the use of innovative data assimilation methods based on geosynchronous satellite measurements. As a step in th~s direction, a non-standard data assimilation experiment in conjunction with GMS-retrieved nowcasted rainfall information introduced to the mesoscale model is conducted. The 3-hourly precipitation forecast information is assimilated through nudging the associated diabatic heating during the early stages of a forecast period. This procedure is expected to enhance details in the moisture field during model integration, and thus improve spinup performance, assuming the errors in the future time latent heating data ate less than intrinsic model background errors.

Smith, Eric A.↗

Comparisons and Contrasts Between Mediterranean Flood-Producing Orographic Storms and Tropical Cyclones Modified by Complex Orography

Often, late-year, flood-producing storms in the Mediterranean Alps region arise from baroclinic disturbances developing on westerly longwaves, with the incipient storm circulations appearing well west of the flood zone. Regardless of the genesis locations and before the storms begin their flow interactions with the orography of the Mediterranean Alps, a process which can produce intense long-lived rainfall over the mountainous terrain of Italy and France, most of their moisture is drawn from warm Mediterranean waters while important elements of their thermodynamic structure are derived from airmass properties originating in north Africa and the Mediterranean basin itself. In fact, because of strong thermodynamic contrasts between the baroclinic zone producing the incipient disturbance and high potential temperature air from the south advecting into and wrapping around the storm center, a Mediterranean storm whose origins are from an evolving shortwave baroclinic perturbation, can become a quasi-warm core hybrid in term of its energetics composition, a process somewhat akin to the development of polar lows. This process begs the question of whether there are similarities between the autumn-winter season flood producing storms of the Mediterranean Alps and warm core tropical cyclones which become altered dynamically and thermodynamically by landfall over steep orographic terrain. This study investigates, from a comparison and contrast perspective, how the similarities and differences between the two types of storms manifest themselves over their respective life cycles. Before the orographically-induced precipitation-release stage in which both types of storm have evolved to the same hydrometeorlogical conclusion, the two storm categories are transforming from entirely distinctive origins, i.e., the classic eastward traveling, baroclinic mid-latitude cyclone forming in western Europe, and the classic eastward-traveling, warm core tropical cyclone developing in the western Pacific. The analyses use high resolution, nonhydrostatic model simulations, satellite data, and various types of ground data to conduct the study. The flood events of Piemonte- 2000 in northern Italy and Typhoon Rusa-2002 in South Korea are selected to focus the analysis on well documented meteorological events.

Smith, Eric A.↗

An Investigation of the Characterization of Cloud Contamination in Hyperspectral Radiances

In regions lacking direct observations, the assimilation of radiances from infrared and microwave sounders is the primary method for characterizing the atmosphere in the analysis process. In recent years, technological advances have led to the launching of more advanced sounders, particularly in the thermal infrared spectrum. With the advent of these hyperspectral sounders, the amount of data available for the analysis process has and will continue to be dramatically increased. However, the utilization of infrared radiances in variational assimilation can be problematic in the presence of clouds; specifically the assessment of the presence of clouds in an instantaneous field of view (IFOV) and the contamination in the individual channels within the IFOV. Various techniques have been developed to determine if a channel is contaminated by clouds. The work presented in this paper and subsequent presentation will investigate traditional techniques and compare them to a new technique, the C02 sorting technique, which utilizes the high spectral resolution of the Atmospheric Infrared Sounder (AIRS) within the framework of the Gridpoint Statistical Interpolation (GSI) 3DVAR system. Ultimately, this work is done in preparation for the assessment of short-term forecast impacts with the regional assimilation of AIRS radiances within the analysis fields of the Weather Research and Forecast Nonhydrostatic Mesoscale Model (WRF-NMM) at the NASA Short-term Prediction Research and Transition (SPORT) Center.

McCarty, William↗

Internal Structure of Rhea

We model the interior of Rhea on the basis of observational constraints and the results from geodynamical models available in the literature. Ten main types of models are defined, depending on the presence or absence of a high-pressure ice layer (ice II), and the extent of separation of the rock component from the volatiles. The degree-two gravity coefficients are computed for each of these models in order to assess which properties of the interior are likely to be inferred from Cassini radio science measurements scheduled on 26 November 2005. C22 greater than 2.5 x 10(exp -4) indicates that the satellite is undifferentiated, except for a slight increase in density with depth resulting from material self-compression. C22 between 1.67 x 10(exp -4)(lower bound) and 1.90 x 10(exp -4) indicates the presence of a rocky core, whose radius can be determined from the satellite's mass and ices densities, for a given temperature profile. For other values, most of the ten models cannot be distinguished from each other. However, assumptions on the density of the rock phase, presence or absence of ice II, and the degree of differentiation could allow a unique model to be determined in many cases. While the calculation presented in this work assumes that Rhea is in hydrostatic equilibrium, it is likely that Rhea' gravity field is partly affected by nonhydrostatic anomalies.

gravity science↗

Impact of MODIS High-Resolution Sea-Surface Temperatures on WRF Forecasts at NWS Miami, FL

Over the past few years,studies at the Short-term Prediction Research and Transition (SPoRT) Center have suggested that the use of Moderate Resolution Imaging Spectroradiometer (MODIS) composite sea-surface temperature (SST) products in regional weather forecast models can have a significant positive impact on short-term numerical weather prediction in coastal regions. The recent paper by LaCasse et al. (2007, Monthly Weather Review) highlights lower atmospheric differences in regional numerical simulations over the Florida offshore waters using 2-km SST composites derived from the MODIS instrument aboard the polar-orbiting Aqua and Terra Earth Observing System satellites. To help quantify the value of this impact on NWS Weather Forecast Offices (WFOs), the SPoRT Center and the NWS WFO at Miami, FL (MIA) are collaborating on a project to investigate the impact of using the high-resolution MODIS SST fields within the Weather Research and Forecasting (WRF) prediction system. The scientific hypothesis being tested is: More accurate specification of the lower-boundary forcing within WRF will result in improved land/sea fluxes and hence, more accurate evolution of coastal mesoscale circulations and the associated sensible weather elements. The NWS MIA is currently running the WRF system in real-time to support daily forecast operations, using the National Centers for Environmental Prediction Nonhydrostatic Mesoscale Model dynamical core within the NWS Science and Training Resource Center's Environmental Modeling System (EMS) software; The EMS is a standalone modeling system capable of downloading the necessary daily datasets, and initializing, running and displaying WRF forecasts in the NWS Advanced Weather Interactive Processing System (AWIPS) with little intervention required by forecasters. Twenty-seven hour forecasts are run daily with start times of 0300,0900, 1500, and 2100 UTC on a domain with 4-km grid spacing covering the southern half of Florida and the far western portions of the Bahamas, the Florida Keys, the Straights of Florida, and adjacent waters of the Gulf of Mexico and Atlantic Ocean. Each model run is initialized using the Local Analysis and Prediction System (LAPS) analyses available in AWIPS, invoking the diabatic. "hot-start" capability. In this WRF model "hot-start", the LAPS-analyzed cloud and precipitation features are converted into model microphysics fields with enhanced vertical velocity profiles, effectively reducing the model spin-up time required to predict precipitation systems. The SSTs are initialized with the NCEP Real-Time Global (RTG) analyses at l/12 degree resolution (approx. 9 km); however, the RTG product does not exhibit fine-scale details consistent with its grid resolution. SPoRT is conducting parallel WRF EMS runs identical to the operational runs at NWS MIA in every respect except for the use of MODIS SST composites in place of the RTG product as the initial and boundary conditions over water. The MODIS SST composites for initializing the SPoRT WRF runs are generated on a 2-km grid four times daily at 0400, 0700, 1600, and 1900 UTC, based on the times of the overhead passes of the Aqua and Terra satellites. The incorporation of the MODIS SST composites into the SPoRTWRF runs is staggered such that the 0400UTC composite initializes the 0900 UTC WRF, the 0700 UTC composite initializes the 1500 UTC WRF, the 1600 UTC composite initializes the 2100 UTC WRF, and the 1900 UTC composite initializes the 0300 UTC WRF. A comparison of the SPoRT and Miami forecasts is underway in 2007, and includes quantitative verification of near-surface temperature, dewpoint, and wind forecasts at surface observation locations. In addition, particular days of interest are being analyzed to determine the impact of the MODIS SST data on the development and evolution of predicted sea/land-breeze circulations, clouds, and precipitation. This paper will present verification results comparing the NWS MIA forecasts the SPoRT experimental WRF forecasts, and highlight any substantial differences noted in the predicted mesoscale phenomena.

Case, Jonathan L.↗

A Multi-Season Study of the Effects of MODIS Sea-Surface Temperatures on Operational WRF Forecasts at NWS Miami, FL

Studies at the Short-term Prediction Research and Transition (SPORT) Center have suggested that the use of Moderate Resolution Imaging Spectroradiometer (MODIS) sea-surface temperature (SST) composites in regional weather forecast models can have a significant positive impact on short-term numerical weather prediction in coastal regions. Recent work by LaCasse et al (2007, Monthly Weather Review) highlights lower atmospheric differences in regional numerical simulations over the Florida offshore waters using 2-km SST composites derived from the MODIS instrument aboard the polar-orbiting Aqua and Terra Earth Observing System satellites. To help quantify the value of this impact on NWS Weather Forecast Offices (WFOs), the SPORT Center and the NWS WFO at Miami, FL (MIA) are collaborating on a project to investigate the impact of using the high-resolution MODIS SST fields within the Weather Research and Forecasting (WRF) prediction system. The project's goal is to determine whether more accurate specification of the lower-boundary forcing within WRF will result in improved land/sea fluxes and hence, more accurate evolution of coastal mesoscale circulations and the associated sensible weather elements. The NWS MIA is currently running WRF in real-time to support daily forecast operations, using the National Centers for Environmental Prediction Nonhydrostatic Mesoscale Model dynamical core within the NWS Science and Training Resource Center's Environmental Modeling System (EMS) software. Twenty-seven hour forecasts are run dally initialized at 0300, 0900, 1500, and 2100 UTC on a domain with 4-km grid spacing covering the southern half of Florida and adjacent waters of the Gulf of Mexico and Atlantic Ocean. Each model run is initialized using the Local Analysis and Prediction System (LAPS) analyses available in AWIPS. The SSTs are initialized with the NCEP Real-Time Global (RTG) analyses at 1/12deg resolution (approx.9 km); however, the RTG product does not exhibit fine-scale details consistent with its grid resolution. SPORT is conducting parallel WRF EMS runs identical to the operational runs at NWS MIA except for the use of MODIS SST composites in place of the RTG product as the initial and boundary conditions over water, The MODIS SST composites for initializing the SPORT WRF runs are generated on a 2-km grid four times daily at 0400, 0700, 1600, and 1900 UTC, based on the times of the overhead passes of the Aqua and Terra satellites. The incorporation of the MODIS SST data into the SPORT WRF runs is staggered such that SSTs are updated with a new composite every six hours in each of the WRF runs. From mid-February to July 2007, over 500 parallel WRF simulations have been collected for analysis and verification. This paper will present verification results comparing the NWS MIA operational WRF runs to the SPORT experimental runs, and highlight any substantial differences noted in the predicted mesoscale phenomena for specific cases.

Case, Jonathan L.↗

Cloud-System Resolving Models: Status and Prospects

Cloud-system resolving models (CRM), which are based on the nonhydrostatic equations of motion and typically have a grid-spacing of about a kilometer, originated as cloud-process models in the 1970s. This paper reviews the status and prospects of CRMs across a wide range of issues, such as microphysics and precipitation; interaction between clouds and radiation; and the effects of boundary-layer and surface-processes on cloud systems. Since CRMs resolve organized convection, tropical waves and the large-scale circulation, there is the prospect for several advances in both basic knowledge of scale-interaction requisite to parameterizing mesoscale processes in climate models. In superparameterization, CRMs represent convection, explicitly replacing many of the assumptions necessary in contemporary parameterization. Global CRMs have been run on an experimental basis, giving prospect to a new generation of climate weather prediction in a decade, and climate models due course. CRMs play a major role in the retrieval of surface-rain and latent heating from satellite measurements. Finally, enormous wide dynamic ranges of CRM simulations present new challenges for model validation against observations.

Tao, Wei-Kuo↗

Cloud-Resolved Simulations of Lightning NOx in an Observed Hector Thunderstorm

Cloud chemistry simulations are being performed for a "Hector" storm observed on 16 November 2005 during the SCOUT-03lACTIVE campaigns based in Darwin, Australia. The primary objective of these simulations is to estimate the average production of NO per lightning flash during the storm. The 3-D WRF-AqChem model is being used for these calculations. This modeling package contains the WRF nonhydrostatic cloud-resolving model, online gas- and aqueous-phase chemistry, and a lightning algorithm (Barth et al., 2007). Early morning soundings of temperature, water vapor and winds are used to initialize the model. Surface heating of the Tiwi Islands is simulated in the model to induce convection. Observations from the Egrett, Falcon, Geophysica, and Dornier aircraft in air undisturbed by the storm are used to construct composite initial condition chemical profiles. Convective transport in the model is tested using tracer species such as CO and O3. Lightning flashes observed by the LINET network are input to the model and a lightning placement scheme is used to inject the resulting NO into the simulated cloud. Various scenarios of NO production per flash are used for cloud-to-ground and intracloud flashes in a series of simulations for the storm. Resulting NO, mixing ratios from each simulation are compared with upper tropospheric anvil observations (from the Geophysica and Egrett aircraft) to determine the best fit with the mean NOx at anvil altitudes, the profile shape, and the frequency distribution of NOx values. We will compare the results for lightning NO production from this tropical thunderstorm with similar analyses conducted for several midlatitude and subtropical convective events.

Pickering, Kenneth↗

Impact of Lake Okeechobee Sea Surface Temperatures on Numerical Predictions of Summertime Convective Systems over South Florida

The NASA Short-term Prediction Research and Transition (SPoRT) Center, the Florida Institute of Technology, and the NOAA/NWS Weather Forecast Office at Miami, FL (MFL) are collaborating on a project to investigate the impact of using high-resolution, 2-km Moderate Resolution Imaging Spectroradiometer (MODIS) sea surface temperature (SST) composites within the Weather Research and Forecasting (WRF) prediction system. The NWS MFL is currently running WRF in real-time to support daily forecast operations, using the National Centers for Environmental Prediction Nonhydrostatic Mesoscale Model dynamical core within the NWS Science and Training Resource Center's Environmental Modeling System (EMS) software. Twenty-seven hour forecasts are run daily initialized at 0300, 0900, 1500, and 2100 UTC on a domain with 4-km grid spacing covering the southern half of Florida and adjacent waters of the Gulf of Mexico and Atlantic Ocean. The SSTs are initialized with the NCEP Real-Time Global (RTG) analyses at 1/12deg resolution. The project objective is to determine whether more accurate specification of the lower-boundary forcing over water using the MODIS SST composites within the 4-km WRF runs will result in improved sea fluxes and hence, more accurate e\~olutiono f coastal mesoscale circulations and the associated sensible weather elements. SPoRT conducted parallel WRF EMS runs from February to August 2007 identical to the operational runs at NWS MFL except for the use of MODIS SST composites in place of the RTG product as the initial and boundary conditions over water. During the course of this evaluation, an intriguing case was examined from 6 May 2007, in which lake breezes and convection around Lake Okeechobee evolved quite differently when using the high-resolution SPoRT MODIS SST composites versus the lower-resolution RTG SSTs. This paper will analyze the differences in the 6 May simulations, as well as examine other cases from the summer 2007 in which the WRF-simulated Lake Okeechobee breezes evolved differently due to the SST initialization. The effects on wind fields and precipitation systems will be emphasized, including validation against surface mesonet observations and Stage IV precipitation grids.

Case, Jonathan L.↗

The High-Resolution Wave-Propagation Method Applied to Meso- and Micro-Scale Flows

The high-resolution wave-propagation method for computing the nonhydrostatic atmospheric flows on meso- and micro-scales is described. The design and implementation of the Riemann solver used for computing the Godunov fluxes is discussed in detail. The method uses a flux-based wave decomposition in which the flux differences are written directly as the linear combination of the right eigenvectors of the hyperbolic system. The two advantages of the technique are: 1) the need for an explicit definition of the Roe matrix is eliminated and, 2) the inclusion of source term due to gravity does not result in discretization errors. The resulting flow solver is conservative and able to resolve regions of large gradients without introducing dispersion errors. The methodology is validated against exact analytical solutions and benchmark cases for non-hydrostatic atmospheric flows.

Ahmad, Nashat N.↗

Vertical Structure of Ice Cloud Layers From CloudSat and CALIPSO Measurements and Comparison to NICAM Simulations

The shape of the vertical profile of ice cloud layers is examined using 4 months of CloudSat and Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO) global measurements taken on January, April, July, and October 2007. Ice clouds are selected using temperature profiles when the cloud base is located above the 253K temperature level. The obtained ice water content (IWC), effective radius, or extinction coefficient profiles are normalized by their layer mean values and are expressed in the normalized vertical coordinate, which is defined as 0 and 1 at the cloud base and top heights, respectively. Both CloudSat and CALIPSO observations show that the maximum in the IWC and extinction profiles shifts toward the cloud bottom, as the cloud depth increases. In addition, clouds with a base reaching the surface in a high-latitude region show that the maximum peak of the IWC and extinction profiles occurs near the surface, which is presumably due to snow precipitation. CloudSat measurements show that the seasonal difference in normalized cloud vertical profiles is not significant, whereas the normalized cloud vertical profile significantly varies depending on the cloud type and the presence of precipitation. It is further examined if the 7 day Nonhydrostatic Icosahedral Atmospheric Model (NICAM) simulation results from 25 December 2006 to 1 January 2007 generate similar cloud profile shapes. NICAM IWC profiles also show maximum peaks near the cloud bottom for thick cloud layers and maximum peaks at the cloud bottom for low-level clouds near the surface. It is inferred that oversized snow particles in the NICAM cloud scheme produce a more vertically inhomogeneous IWC profile than observations due to quick sedimentation.

Ham, Seung-Hee↗

Toward Improved Land Surface Initialization in Support of Regional WRF Forecasts at the Kenya Meteorological Department

Flooding and drought are two key forecasting challenges for the Kenya Meteorological Department (KMD). Atmospheric processes leading to excessive precipitation and/or prolonged drought can be quite sensitive to the state of the land surface, which interacts with the boundary layer of the atmosphere providing a source of heat and moisture. The development and evolution of precipitation systems are affected by heat and moisture fluxes from the land surface within weakly-sheared environments, such as in the tropics and sub-tropics. These heat and moisture fluxes during the day can be strongly influenced by land cover, vegetation, and soil moisture content. Therefore, it is important to represent the land surface state as accurately as possible in numerical weather prediction models. Enhanced regional modeling capabilities have the potential to improve forecast guidance in support of daily operations and high-end events over east Africa. KMD currently runs a configuration of the Weather Research and Forecasting (WRF) model in real time to support its daily forecasting operations, invoking the Nonhydrostatic Mesoscale Model (NMM) dynamical core. They make use of the National Oceanic and Atmospheric Administration / National Weather Service Science and Training Resource Center's Environmental Modeling System (EMS) to manage and produce the WRF-NMM model runs on a 7-km regional grid over eastern Africa. Two organizations at the National Aeronautics and Space Administration Marshall Space Flight Center in Huntsville, AL, SERVIR and the Short-term Prediction Research and Transition (SPoRT) Center, have established a working partnership with KMD for enhancing its regional modeling capabilities. To accomplish this goal, SPoRT and SERVIR will provide experimental land surface initialization datasets and model verification capabilities to KMD. To produce a land-surface initialization more consistent with the resolution of the KMD-WRF runs, the NASA Land Information System (LIS) will be run at a comparable resolution to provide real-time, daily soil initialization data in place of interpolated Global Forecast System soil moisture and temperature data. Additionally, real-time green vegetation fraction data from the Visible Infrared Imaging Radiometer Suite will be incorporated into the KMD-WRF runs, once it becomes publicly available from the National Environmental Satellite Data and Information Service. Finally, model verification capabilities will be transitioned to KMD using the Model Evaluation Tools (MET) package, in order to quantify possible improvements in simulated temperature, moisture and precipitation resulting from the experimental land surface initialization. The transition of these MET tools will enable KMD to monitor model forecast accuracy in near real time. This presentation will highlight preliminary verification results of WRF runs over east Africa using the LIS land surface initialization.

Case. Jonathan↗

Global Cloud-Resolving Models

Global cloud-resolving models (GCRMs) are a new category of atmospheric global models designed to solve different flavors of the nonhydrostatic equations through the use of kilometer-scale global meshes. GCRMs make it possible to explicitly simulate deep convection, thereby avoiding the need for cumulus parameterization and allowing for clouds to be “resolved” by microphysical models responding to grid-scale forcing. GCRMs require high-resolution discretization over the globe, for which a variety of mesh structures have been proposed and employed. The first GCRM was constructed 15 years ago, and in recent years, other groups have also begun adopting this approach, enabling the first intercomparison studies of such models. Because conventional general circulation models (GCMs) suffer from large biases associated with cumulus parameterization, GCRMs are attractive tools for researchers studying global weather and climate. In this review, GCRMs are described, with some emphasis on their historical development and the associated literature documenting their use. The advantages of GCRMs are presented, and currently existing GCRMs are listed and described. Future prospects for GCRMs are also presented in the final section.

Satoh, Masaki↗

High-resolution Wave Propagation Method for Stratified Flows

The implementation of the multidimensional f-waves Riemann solver for the time-dependent, three-dimensional, nonhydrostatic, meso- and microscale atmospheric flows is described in detail. The Riemann solver employs flux-based wave decomposition (f-waves) for the calculation of Godunov fluxes in which the flux differences are written directly as the linear combination of the right eigenvectors of the hyperbolic system. The scheme incorporates the source term due to gravity without introducing discretization errors which is an important property in the context of atmospheric flows. The resulting flow solver is conservative, accurate, stable, and well-balanced. The implementation of the solver is evaluated using benchmark test cases for atmospheric dynamics.

Riemann problem↗

The Diurnal Cycle of Precipitation and Organized Convection in a Set of Global Mesoscale Simulations with the NASA GEOS AGCM

The NASA Global Modeling and Assimilation Office (GMAO) has conducted a series of 40-day nonhydrostatic global simulations with horizontal grid spacing ranging from 200 km to 3 km, as part of the DYAMOND model intercomparison project. The Goddard Earth Observing System (GEOS) model was run with the Grell-Freitas scale-aware convection scheme, which smoothly reduces parameterized deep convection with increasing resolution. Here we evaluate the diurnal cycle and other statistics of precipitation and organized convection as a function of resolution. For validation we use the 0.1 degree IMERG precipitation and 4 km Merged IR brightness temperature datasets, focusing on four regions: the continental United States, Amazonia, the equatorial Indian ocean, and the Maritime Continent. Early results indicate good phase agreement but excessive magnitude of the continental diurnal cycle of precipitation at coarser resolutions, with improved magnitude as resolution increases and the role of parameterization is reduced. Convective cloud clusters are identified with a brightness temperature threshold, and we find realistic numbers of the largest clusters (>10^4 km sq) at all resolutions, while the number of smaller clusters increases with resolution, approaching observations when dx=3 km. The observed diurnal cycle in the cluster size distribution is also reproduced, with realistic magnitude in the highest resolution runs. Precipitation characteristics across cluster sizes are also examined. The results show the potential for global mesoscale simulations as a community science resource, and we invite collaboration to explore these runs in greater detail.

Arnold, Nathan↗

Dispersion and Aging of Volcanic Aerosols after the La Soufriere Eruption in April 2021

Volcanic aerosols change the atmospheric composition and thereby affect weather and climate. Aerosol dynamic processes such as nucleation, condensation, and coagulation modify the shape, size, and mass of aerosol particles, which influence their atmospheric lifetime and radiative properties. Nevertheless, most models omit these processes for ash particles. In this work, we explore the ash aerosol aging and sulfate production during the first 4 days following the 2021 La Soufrière (St. Vincent) eruption with the ICON-ART model (ICOsahedral Nonhydrostatic model with Aerosol and Reactive Trace gases). Online coupling of ICON-ART with a one-dimensional volcanic plume model calculates volcanic emission, which makes it possible to resolve the different eruption phases of the noncontinuous La Soufrière eruption. We compared our simulated aerosol distribution and composition with observations from the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) instrument, the Multiangle Imaging SpectroRadiometer (MISR) Research Aerosol (RA) Algorithm, and the Barbados Cloud Observatory (BCO). We show that online coupling is essential to adequately model the emissions and plume development close to the volcano. The modeled aerosol aging is in very good agreement with observations from MISR near the emission source and with CALIOP at larger distances. Furthermore, particle aging occurs faster in the troposphere than in the stratosphere due to the availability of water vapor and OH, but a layer of coated ash appears at the plume top due to faster oxidation of SO2 and lofting by aerosol-radiation interaction. This paper gives the first direct comparison of aerosol aging in volcanic eruption plumes between simulations and observations

volcanic aerosols↗

Modeling Dust Mineralogical Composition: Sensitivity to Soil Mineralogy Atlases and Their Expected Climate Impacts

Soil dust aerosols are a key component of the climate system, as they interact with short- and long-wave radiation, alter cloud formation processes, affect atmospheric chemistry and play a role in biogeochemical cycles by providing nutrient inputs such as iron and phosphorus. The influence of dust on these processes depends on its physicochemical properties, which, far from being homogeneous, are shaped by its regionally varying mineral composition. The relative amount of minerals in dust depends on the source region and shows a large geographical variability. However, many state-of-the-art Earth system models (ESMs), upon which climate analyses and projections rely, still consider dust mineralogy to be invariant. The explicit representation of minerals in ESMs is more hindered by our limited knowledge of the global soil composition along with the resulting size-resolved airborne mineralogy than by computational constraints. In this work we introduce an explicit mineralogy representation within the state-of-the-art Multiscale Online Nonhydrostatic AtmospheRe CHemistry (MONARCH) model. We review and compare two existing soil mineralogy datasets, which remain a source of uncertainty for dust mineralogy modeling and provide an evaluation of multiannual simulations against available mineralogy observations. Soil mineralogy datasets are based on measurements performed after wet sieving, which breaks the aggregates found in the parent soil. Our model predicts the emitted particle size distribution (PSD) in terms of its constituent minerals based on brittle fragmentation theory (BFT), which reconstructs the emitted mineral aggregates destroyed by wet sieving. Our simulations broadly reproduce the most abundant mineral fractions independently of the soil composition data used. Feldspars and calcite are highly sensitive to the soil mineralogy map, mainly due to the different assumptions made in each soil dataset to extrapolate a handful of soil measurements to arid and semi-arid regions worldwide. For the least abundant or more difficult-to-determine minerals, such as iron oxides, uncertainties in soil mineralogy yield differences in annual mean aerosol mass fractions of up to ∼ 100 %. Although BFT restores coarse aggregates including phyllosilicates that usually break during soil analysis, we still identify an overestimation of coarse quartz mass fractions (above 2 µm in diameter). In a dedicated experiment, we estimate the fraction of dust with undetermined composition as given by a soil map, which makes up ∼ 10 % of the emitted dust mass at the global scale and can be regionally larger. Changes in the underlying soil mineralogy impact our estimates of climate-relevant variables, particularly affecting the regional variability of the single-scattering albedo at solar wavelengths or the total iron deposited over oceans. All in all, this assessment represents a baseline for future model experiments including new mineralogical maps constrained by high-quality spaceborne hyperspectral measurements, such as those arising from the NASA Earth Surface Mineral Dust Source Investigation (EMIT) mission.

Soil dust↗