Non-hydrostatic Model Research Articles

Прогноз толщины ледяного покрова в прибрежных областях Карского и Охотского морей

Based on the CICE viscoplastic sea ice model, the ice cover of the Gulf of Ob and the northwestern Okhotsk Sea was simulated for the winter seasons of 2021/2022, 2022/2023, and 2023/2024. The description of the seasonal variability of sea ice characteristics in these areas and examples of their numerical modeling are given. Using atmospheric forcing based on the WRF-ARW nonhydrostatic atmospheric model, the sea ice thickness in the northwestern Okhotsk Sea was predicted for April 2024. A comparison of the forecasts with the ice maps of the Hydrometcentre of Russia and sea ice thickness measurements at the Ayan and Bolshoy Shantar stations gave satisfactory results.

Efficient spectral element method for the Euler equations on unbounded domains

Mitigating the impact of waves leaving a numerical domain has been a persistent challenge in numerical modeling. Reducing wave reflection at the domain boundary is crucial for accurate simulations. Absorbing layers, while common, often incur significant computational costs. This paper introduces an efficient application of a Legendre-Laguerre basis for absorbing layers for two-dimensional non-linear compressible Euler equations. The method couples a spectral-element bounded domain with a semi-infinite region, employing a tensor product of Lagrange and scaled Laguerre basis functions. Semi-infinite elements are used in the absorbing layer with Rayleigh damping. In comparison to existing methods with similar absorbing layer extensions, this approach, a pioneering application to the Euler equations of compressible and stratified flows, demonstrates substantial computational savings. The study marks the first application of semi-infinite elements to mitigate wave reflection in the solution of the Euler equations, particularly in nonhydrostatic atmospheric modeling. A comprehensive set of tests demonstrates the method's versatility for general systems of conservation laws, with a focus on its effectiveness in damping vertically propagating mountain gravity waves, a benchmark for atmospheric models. Across all tests, the model presented in this paper consistently exhibits notable performance improvements compared to a traditional Rayleigh damping approach.

Robustness of the relationship between tropical high-cloud cover and large-scale circulations

High clouds play an important role in climate variations by shading the sun and contributing to the greenhouse effect. It has been suggested that these variations are strongly controlled by atmospheric circulations; however, the relationship between high clouds and circulations in global simulations has been confirmed only for a limited number of general circulation models (GCMs). Increasing our understanding of this relationship requires a more complete dataset of state-of-the-art GCMs. This study quantifies the relationship between high clouds and circulations by spatial and temporal correlation using data from 28 atmospheric GCMs along with a global nonhydrostatic model, a new-type of GCM, and reveals a robust and strong relationship in most models. This study also finds that the sensitivity, which is defined as the slope of the relationship between high clouds and mass fluxes of circulations, is highly model-dependent. This suggests that a similar change of circulations does not guarantee that of high clouds, which could be one reason for the complexity in projecting high-cloud changes due to warming. Moreover, this study confirms in both observation and models that the relationship with circulations is stronger for medium-thickness and thick high clouds than for thin clouds. Furthermore, the relationship between high clouds and circulations is more convincing in near-equatorial regions between 15°S–15°N, the ascending branch of the Hadley circulation, where deep convection is especially active.

Computational study on complex wave hydrodynamics of multidirectional extreme waves at fringing reef

Due to the special topography of coral reefs, most of the energy of incoming waves can be reduced through the wave breakings at the reef edge, thus playing a vital role in protecting the coastal regions behind coral reefs. In severe weather conditions, wind-generated waves tend to propagate in multiple directions in real ocean environment. The abrupt multidirectional wave-focusing mechanisms may result in the formation of huge waves in both deep and shallow waters, threatening the safety of coastal infrastructures and well-being of coastal communities at the reef islands. Therefore, it is crucial to thoroughly examine the complex hydrodynamics of multidirectional extreme waves at the fringing reefs. Previous studies on reef wave hydrodynamics predominantly focus on unidirectional waves. Only a small number of studies have considered the wave hydrodynamics of fringing reefs under multidirectional waves, let alone multidirectional extreme waves. To address the knowledge gap in previous studies, this study analyzes the complex hydrodynamic processes of multidirectional extreme waves at the fringing reef by applying a nonhydrostatic numerical wave model (NHWAVE). Influences of several key factors are analyzed, such as significant wave height, peak wave period, water depth, incident wave angle, and reef topography. It is expected that the research findings of this study will contribute to a deeper understanding of the hydrodynamics of multidirectional extreme waves at the fringing reef.

Investigating the sign of stratocumulus adjustments to aerosols in the ICON global storm-resolving model

Abstract. Aerosols can cause brightening of stratocumulus clouds, thereby cooling the climate. Observations and models disagree on the magnitude of this cooling, partly because of the aerosol-induced liquid water path (LWP) adjustment, with climate models predicting an increase in the LWP and satellites observing a weak decrease in response to increasing aerosols. With higher-resolution global climate models, which allow the simulation of mesoscale circulations in which stratocumulus clouds are embedded, there is hope to start bridging this gap. In this study, we present boreal summertime simulations conducted with the ICOsahedral Non-hydrostatic (ICON) global storm-resolving model (GSRM). Compared to geostationary satellite data, ICON produces realistic cloud coverage in the stratocumulus regions; however, these clouds look cumuliform, and the sign of LWP adjustments disagrees with observations. We investigate this disagreement with a causal approach, which combines time series with knowledge of cloud processes, allowing us to diagnose the sources of observation–model discrepancies. The positive ICON LWP adjustment results from a superposition of processes, with an overestimated positive response due to (1) precipitation suppression, (2) a lack of wet scavenging, and (3) cloud deepening under a weak inversion, despite (4) small negative influences from cloud-top entrainment enhancement. We also find that precipitation suppression and entrainment enhancement occur at different intensities during the day and the night, implying that daytime satellite studies suffer from selection bias. This causal methodology can guide modelers on how to modify model parameterizations and setups to reconcile conflicting studies concerning the sign and magnitude of LWP adjustments across different data sources.

Updating the radiation infrastructure in MESSy (based on MESSy version 2.55)

Abstract. The calculation of the radiative transfer is a key component of global circulation models. In this article, we describe the most recent updates of the radiation infrastructure in the Modular Earth Submodel System (MESSy). These updates include the implementation of the PSrad radiation scheme within the RAD submodel. Furthermore, the radiation-related submodels CLOUDOPT (for the calculation of cloud optical properties) and AEROPT (for the calculation of aerosol optical properties) have been updated and are now more flexible in order to deal with different sets of shortwave and longwave bands of radiation schemes. In the wake of these updates, a new submodel (ALBEDO), which features solar-zenith-angle-dependent albedos and a new satellite-based background (white sky) albedo, was created. All of these developments are backward compatible, and previous features of the MESSy radiation infrastructure remain available. Moreover, these developments mark an important step in the use of the ECHAM/MESSy Atmospheric Chemistry (EMAC) model, as the update of the radiation scheme was a key aspect in the development of the sixth generation of the European Centre for Medium-Range Weather Forecasts – HAMburg (ECHAM6) model from ECHAM5. The developments presented here are also aimed towards using the MESSy infrastructure with the ICOsahedral Non-hydrostatic (ICON) model as a base model. The improved infrastructure will also aid in the implementation of additional radiation schemes once this should be needed. We have optimized the set of free parameters for two general circulation model-type (GCM-type) setups for pre-industrial and present-day conditions: one with the radiation scheme that was used to date (i.e. the radiation scheme of ECHAM5) and one with the newly implemented PSrad radiation scheme. After this parameter optimization, we performed four model simulations and evaluated the corresponding model results using reanalysis and observational data. The most apparent improvements related to the updated radiation scheme are the reduced cold biases in the tropical upper troposphere and lower stratosphere and the extratropical lower stratosphere and a strengthened polar vortex. The former is also related to improved stratospheric humidity and its variability if the new radiation scheme is employed. Using the multiple radiation call capability of MESSy, we have applied the two model configurations to calculate instantaneous and stratospheric-adjusted radiative forcings related to changes in greenhouse gases. Overall, we find that for many forcing experiments the simulations with the new radiation scheme show improved radiative forcing values. This is in particular the case for methane radiative forcings, which are considerably higher when assessed with the new radiation scheme and thus in better agreement with reference values.

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.

Skill assessment of a total water level and coastal change forecast during the landfall of a hurricane

The Total Water Level and Coastal Change Forecast (TWL&CC Forecast) provides coastal communities with 6-day notice of potential elevated water levels and coastal change (i.e., dune erosion, overwash, or inundation) on sandy beaches that threatens safety, infrastructure, or resources. This continuously operating model provides hourly information for select regions along U.S. Gulf of Mexico and Atlantic Ocean coastlines. The objective of this work is to assess the skill of forecasts during a period of elevated water levels along the coasts of North Carolina (NC) and South Carolina, USA caused by Hurricane Isaias in August 2020, using a combination of observations and model hindcasts. Water levels and waves were observed throughout the storm at three locations near Wrightsville Beach, NC, which provided information to assess forecast skill; a wave buoy offshore, a tide gage at a local pier, and a pressure sensor deployed at the pier. In addition to observations, the non-hydrostatic phase-resolving model SWASH (Simulating WAves till SHore) was forced with hourly wave energy spectra derived from a coupled Delft3D-SWAN simulation during the peak of Isaias, to complement observations by computing nearshore wave height and wave-induced setup and runup at the shoreline. During the storm peak, SWASH-simulated water levels at the sensor position were comparable to those at the maximum landward extent (bias = −0.05 m; gain = 0.26; r2 = 0.99), suggesting that observations at the USGS sensor location were a useful proxy for total water level (TWL; sum of tide, surge and wave runup) at the shoreline that are predicted by the TWL&CC Forecast. The TWL forecast at Wrightsville Beach was consistent with observations from the USGS sensor (bias = −0.38 m and −0.74 m, scatter index = 0.22 and 0.28 for the two forecast model grids considered, respectively; weighted regression considering model uncertainty explained 95 percent of variability in observed TWL). Observed TWL was within the confidence interval of the TWL&CC Forecast for the 5 h at the storm peak. Forecast mean water levels (MWL; sum of tide, surge and wave setup) and tide gage observations were also consistent (bias = 0.07 m and 0.02 m for the forecast model grids; scatter index = 0.46; r2 = 0.80). Forecast MWL at the storm peak was within 0.06 m of the observed MWL from the tide gage for both sites. In the region where Isaias made landfall, eight additional pressure sensors were compared to the peak TWL forecast (bias = 0.14 m; scatter index = 0.18). Forecast TWL explained 90 percent of observed variability in TWL when considering uncertainty of the forecast with a weighted regression. The results demonstrate that wave-driven water levels contributed a significant portion of the forecast TWL during Isaias (52 percent during the three peak hours of the storm), and that TWL were represented using the forecast model. Mean absolute error of the coastal change forecast and observed overwash is 0.4 and 0.14 for the two forecast model grids considered. The skill demonstrated by this computationally efficient method indicates that the forecasting system can provide fast and reliable predictions of TWL across hundreds of km of coastline at sub-km resolution, days to hours in advance of when storms threaten coastal regions.

Multi-layer non-hydrostatic free-surface flow model with kinematic seafloor for seismic tsunami generation

Earthquake induced tsunamis pose devastating threat to coastal communities worldwide. Accurate description of tsunami generation by kinematic fault rupture is of importance to investigate earthquake events and evaluate tsunami impacts. The paper develops a multi-layer non-hydrostatic model with a kinematic bottom boundary condition and conducts comprehensive validation to examine the model’s capability in resolving seismic tsunami generation. The two-dimensional governing equations of the multi-layer non-hydrostatic free-surface flow system equipped with kinematic seafloor displacement is first introduced in Cartesian coordinates. A combined finite difference and finite volume scheme is utilized to discretize the governing equations with flow variables arranged on a staggered Arakawa C-grid. Application of pressure correction technique to the discretized formulations yields Poisson-type equations from which the non-hydrostatic pressure is solved for the next time step to complete the temporal integration. Based on existing analytical solutions, four groups of tsunami generation cases considering broad ranges of source parameters, including horizontal scale, rise time, and rupture velocity, are designed to demonstrate performance of the proposed non-hydrostatic model with one-, two-, and three-layers as well as the hydrostatic one. Comparison of 59 generation cases including extreme scenarios indicates the non-hydrostatic model performs better than the hydrostatic model in reproducing the entire waveform and predicting the maximum wave amplitude. High modeling accuracy can be achieved through incorporation of more layers. The proposed multi-layer non-hydrostatic model is a powerful tool for investigating earthquake source mechanisms and evaluating coastal tsunami hazards.

Dynamical downscaling CMIP6 models over New Zealand: added value of climatology and extremes

Dynamical downscaling provides physics-based high-resolution climate change projections across regional and local scales. This is particularly important for island nations characterized by complex terrain, where the coarse resolution of global climate model (GCM) output often prohibits direct use. One of the main motivations for dynamical downscaling is to reduce biases relative to the host GCM at the local scale, which can be quantified through assessing ‘added value’. However, added value from downscaling is not guaranteed; quantifying this can help users make informed decisions about how best to use available climate projection data. Here we describe the experiment design of the updated national climate projections for New Zealand based on dynamical downscaling. The global non-hydrostatic Conformal Cubic Atmospheric Model (CCAM) is primarily used for downscaling, with a global stretched grid targeting high resolution over New Zealand (12-km) and the wider South Pacific region (12–35-km). Focusing on the historical simulations, we assess added value for a range of metrics, climatological fields, extreme indices, and tropical cyclones. The main strengths of the downscaling include generally large improvements relative to the host GCM for temperature and orographic precipitation. Inter-annual variability in temperature is well captured across New Zealand, and several temperature and precipitation-based extreme indices show large improvements. The representation of tropical cyclones reaching at least category 2 intensity is generally improved relative to the large consistent under-representation in the host GCMs. The remaining biases are explored and discussed forming the basis for ongoing bias-correction work.

A performance study of horizontally explicit vertically implicit (HEVI) time-integrators for non-hydrostatic atmospheric models

We conduct a thorough study of different forms of horizontally explicit and vertically implicit (HEVI) time-integration strategies for the compressible Euler equations on spherical domains typical of nonhydrostatic global atmospheric applications. We compare the computational time and complexity of two nonlinear variants (NHEVI-GMRES and NHEVI-LU) and a linear variant (LHEVI). We report on the performance of these three variants for a number of additive Runge-Kutta methods ranging in order of accuracy from second through fifth, and confirm the expected order of accuracy of the HEVI methods for each time-integrator. To gauge the maximum usable time-step of each HEVI method, we run simulations of a nonhydrostatic baroclinic instability for 100 days and then use this time-step to compare the time-to-solution of each method. The results show that NHEVI-LU is 3x faster than NHEVI-GMRES, and LHEVI is 8x faster than NHEVI-GMRES, for the idealized cases tested. The baroclinic instability and inertia-gravity wave simulations indicate that the optimal choice of HEVI method is LHEVI. For the fastest time-to-solution, the second and third order methods are best although the better accuracy of the high-order methods (particularly the fourth order method) should be considered.In the future, we will report on whether these results hold for more complex problems using, e.g., real atmospheric data and/or a higher model top typical of space weather applications.

The impact of mesh size, turbulence parameterization, and land‐surface‐exchange scheme on simulations of the mountain boundary layer in the hectometric range

AbstractThe horizontal grid spacing of numerical weather prediction models keeps decreasing towards the hectometric range. We perform limited‐area simulations with the Icosahedral Nonhydrostatic (ICON) model across horizontal grid spacings (1 km, 500 m, 250 m, 125 m) in the Inn Valley, Austria, and evaluate the model with observations from the Cross‐Valley Flow in the Inn Valley Investigated by Dual‐Doppler LIDAR Measurements (CROSSINN) measurement campaign. This allows us to investigate whether increasing the horizontal resolution automatically improves the representation of the flow structure, surface exchange, and common meteorological variables. Increasing the horizontal resolution results in an improved simulation of the thermally induced circulation. However, the model still faces challenges with scale interactions and the evening transition of the up‐valley flow. Differences between two turbulence schemes (1D turbulence kinetic energy (TKE) and 3D Smagorinsky) emerge due to their different surface transfer formulations, yielding a delayed evening transition in the 3D Smagorinsky scheme. Generally speaking, the correct simulation of the mountain boundary layer depends mostly on the representation of model topography and surface exchange, and the choice of turbulence parameterization is secondary.

Implementation and tests of the NLT3D scheme in the ICON model

AbstractThe main goal of this article is to test the long‐term performance of the three‐dimensional non‐local turbulence (NLT) parameterization scheme at different grid sizes in the so‐called gray zone between classical mesoscale modeling ( several km) and large eddy simulations (LES: several 100 m). For this, NLT has been implemented in the numerical weather prediction Icosahedral Nonhydrostatic model (ICON) of Deutscher Wetterdienst (DWD). Results are compared with a one‐dimensional version of NLT (NLT) and with two operational turbulence schemes available in ICON. Comparisons with observations from radiosondes, the operational surface synoptic (SYNOP) station network, and RAdar‐OnLine‐ANeichung (RADOLAN) radar data of DWD indicate that all turbulence schemes investigated perform reasonably well. Nonetheless, a more detailed study of the model results reveals several interesting differences between the turbulence parameterizations to be discussed in detail. Median absolute errors (MAE) from point‐to‐point comparisons between numerical results and SYNOP observations tend to be smaller than those from comparisons with averaging simulated fields over an environment around each station location. This behavior indicates an information loss caused by the averaging process. For the 2‐m temperature () and the hourly precipitation sums (), MAEs decrease with decreasing grid sizes, thus suggesting an information gain for finer grids. The nighttime MAEs of and obtained with NLT and NLT are similar to or lower than those of the operational turbulence schemes of ICON. Moreover, during a shallow warm‐air intrusion, NLT and especially NLT yield a more realistic representation of the horizontal structures of and, during nighttime stable boundary‐layer situations, also . Radiosonde profiles of the potential temperature confirm a reasonable vertical mixing as obtained with NLT and NLT.

Lagrangian transport based on the winds of the icosahedral nonhydrostatic model (ICON)

Lagrangian transport based on the winds of the icosahedral nonhydrostatic model (ICON)

The application of a non-hydrostatic RANS model for simulating irregular wave breaking on a barred and sloping beach

In the current study, a non-hydrostatic two-dimensional vertical (2DV) numerical model with a shock-capturing technique is developed to simulate irregular wave breaking on a barred beach. The complete form of the 2DV Reynolds-averaged Navier-Stokes (RANS) equations is discretized using the finite volume approach. In the spatial discretization phase, a flux limiter function is employed for the velocity advection terms to prevent non-physical oscillations in regions with steep gradients. A novel method has been introduced during the discretization stage, leading to a significant reduction in computational expenses. For temporal discretization, the Leapfrog method, known for its second-order accuracy in time, is employed to address wave-damping issues. Model validation involves a comparison between simulation results and experimental data pertaining to irregular wave breaking, affirming the satisfactory performance of the model. To evaluate the accuracy of the model, the Root Mean Square Error (RMSE) was calculated. The assessment showed that the model is capable of estimating the significant wave height reported in the experiments used with an RMSE between 0.002 and 0.089, and the mean wave periods with an RMSE between 0.061 and 0.252.

High-resolution tsunami hazard assessment for the Guangdong-Hong Kong-Macao Greater Bay Area based on a non-hydrostatic tsunami model

High-resolution tsunami hazard assessment for the Guangdong-Hong Kong-Macao Greater Bay Area based on a non-hydrostatic tsunami model

An evaluation of microphysics in a numerical model using Doppler velocity measured by ground-based radar for application to the EarthCARE satellite

Abstract. The Cloud Profiling Radar (CPR) of the Earth Cloud, Aerosol and Radiation Explorer (EarthCARE) has a new capability to observe the Doppler velocity related to the vertical air motion of the terminal velocity of hydrometeors. The new observation from space will be used to evaluate and improve the model. Before the launch of EarthCARE, we need to develop a methodology for using the CPR data for model evaluations. In this study, we evaluated simulated data by a stretched version of the global non-hydrostatic model over Japan with a ground-based CPR using an instrument design similar to the EarthCARE CPR. We chose two cases with different precipitation events in September 2019 using two cloud microphysics schemes. We introduced the categorization method for evaluating microphysics using Doppler velocity. The results show that the liquid and solid phases of hydrometeors are divided in Doppler velocity, and the model's terminal velocities of rain, snow, and graupel categories can be evaluated with the observation. The results also show that the choice of microphysics scheme has a more significant impact than the dependence on precipitation cases. We discussed the application of the EarthCARE-like simulation results using a satellite simulator.

Role of rotational and divergence effects induced by wind turbines wakes on near-ground air warming

Various studies, involving numerical simulations, satellite measurements, and field measurement campaigns, have demonstrated that wind turbine wakes can lead to changes in near-ground air temperature. To investigate these phenomena an Actuator Disk with Rotation (ADR) is implemented within the LES non-hydrostatic atmospheric model Meso-NH and validated using the NewMEXICO wind tunnel measurements and the VERTEX field measurement campaign. The NewMEXICO case allows the aerodynamic validation of the approach, especially the effect of the wake on wind speed components. On the VERTEX case, meteorological effects of wind turbines are studied under various atmospheric stability conditions through the application of homogeneous sensible heat fluxes at the ground level. Results show wake-related warming and cooling near the ground for stable and unstable conditions respectively, and no significant change for neutral conditions. These results show a good tendency agreement with the VERTEX campaign’s measurements. The comparison between an Actuator Disk with No Rotation (ADNR) and an ADR simulation show the dual contribution of the rotation of the wake and the growth of the stream tube on the temperature change in the wake but also on the near-ground air temperature.

A multi-faceted methodology for calibration of coastal vegetation drag coefficient

The accurate prediction of wave height attenuation due to vegetation is crucial for designing effective and efficient natural and nature-based solutions for flood mitigation, shoreline protection, and coastal ecosystem preservation. Central to these predictions is the estimation of the vegetation drag coefficient (Cd). The present study undertakes a comprehensive evaluation of three distinct methodologies for estimating the drag coefficient: traditional manual calibration, calibration using a novel application of state-of-the-art metaheuristic optimization algorithms, and the integration of an empirical bulk drag coefficient formula (Tanino and Nepf, 2008) into the XBeach non-hydrostatic wave model. These methodologies were tested using a series of existing laboratory experiments involving nearshore vegetation on a sloping beach. A key innovation of the study is the first application of metaheuristic optimization algorithms for calibrating the drag coefficient, which enables efficient automated searches to identify optimal values aligning with measurements. We found that the optimization algorithms rapidly converge to precise drag coefficients, enhancing accuracy and overcoming limitations in manual calibration which can be laborious and inconsistent. While the integrated empirical formula also demonstrates reasonable performance, the optimization approach exemplifies the potential of computational techniques to transform traditional practices of model calibration. Comparing these strategies provides a framework to determine effective methodology based on constraints in determining the vegetation drag coefficient.

MieAI: a neural network for calculating optical properties of internally mixed aerosol in atmospheric models

Aerosols influence weather and climate by interacting with radiation through absorption and scattering. These effects heavily rely on the optical properties of aerosols, which are mainly governed by attributes such as morphology, size distribution, and chemical composition. These attributes undergo continuous changes due to chemical reactions and aerosol micro-physics, resulting in significant spatio-temporal variations. Most atmospheric models struggle to incorporate this variability because they use pre-calculated tables to handle aerosol optics. This offline approach often leads to substantial errors in estimating the radiative impacts of aerosols along with posing significant computational burdens. To address this challenge, we introduce a computationally efficient and robust machine learning approach called MieAI. It allows for relatively inexpensive calculation of the optical properties of internally mixed aerosols with a log-normal size distribution. Importantly, MieAI fully incorporates the variability in aerosol chemistry and microphysics. Our evaluation of MieAI against traditional Mie calculations, using number concentrations from the ICOsahedral Nonhydrostatic model with Aerosol and Reactive Trace gases (ICON-ART) simulations, demonstrates that MieAI exhibits excellent predictive accuracy for aerosol optical properties. MieAI achieves this with errors well within 10%, and it operates more than 1000 times faster than the benchmark approach of Mie calculations. Due to its generalized nature, the MieAI approach can be implemented in any chemistry transport model which represents aerosol size distribution in the form of log-normally distributed internally mixed modes. This advancement has the potential to replace frequently employed look-up tables and plays a substantial role in the ongoing attempts to reduce uncertainties in estimating aerosol radiative forcing.