Two-dimensional evaluation of atham-fluidity, a nonhydrostatic atmospheric model using mixed continuous/discontinuous finite elements and anisotropic grid optimization

AbstractThis paper presents the first attempt to apply the compressible nonhydrostatic Active Tracer High-Resolution Atmospheric Model–Fluidity (ATHAM-Fluidity) solver to a series of idealized atmospheric test cases. ATHAM-Fluidity uses a hybrid finite-element discretization where pressure is solved on a continuous second-order grid while momentum and scalars are computed on a first-order discontinuous grid (also known as ). ATHAM-Fluidity operates on two- and three-dimensional unstructured meshes, using triangular or tetrahedral elements, respectively, with the possibility to employ an anisotropic mesh optimization algorithm for automatic grid refinement and coarsening during run time. The solver is evaluated using two-dimensional-only dry idealized test cases covering a wide range of atmospheric applications. The first three cases, representative of atmospheric convection, reveal the ability of ATHAM-Fluidity to accurately simulate the evolution of large-scale flow features in neutral atmospheres at rest. Grid convergence without adaptivity as well as the performances of the Hermite–Weighted Essentially Nonoscillatory (Hermite-WENO) slope limiter are discussed. These cases are also used to test the grid optimization algorithm implemented in ATHAM-Fluidity. Adaptivity can result in up to a sixfold decrease in computational time and a fivefold decrease in total element number for the same finest resolution. However, substantial discrepancies are found between the uniform and adapted grid results, thus suggesting the necessity to improve the reliability of the approach. In the last three cases, corresponding to atmospheric gravity waves with and without orography, the model ability to capture the amplitude and propagation of weak stationary waves is demonstrated. This work constitutes the first step toward the development of a new comprehensive limited area atmospheric model.This research has received funding from the European Union Seventh Framework Program (FP7/2007-2013) under Grant agreement 603663 for the research project PEARL (Preparing for Extreme And Rare events in coastaL regions). The EPSRC multiphase program grant MEMPHIS is also acknowledged.This is the author accepted manuscript. The final version is available from the American Meteorological Society via http://dx.doi.org/10.1175/MWR-D-15-0398.

A general description of entrainment in buoyant cloudy plumes including the effects of mixing-induced evaporation

AbstractIn this work, the steady-state one-dimensional axisymmetric plume theory is revisited and generalized to include the effects of nonhomogeneous updraft velocity and buoyancy profiles across the plume, environmental shear, and, more importantly, evaporative cooling resulting from the mixing between cloudy air and the dry environment. Based on an energy consistency argument, a method is proposed to derive a relationship for the fractional lateral mixing rate (which may here be positive or negative) from the plume’s integral equations, as well as a set of equations for the equivalent plume properties, both of which maintain a high degree of generality by incorporating effects of environmental shear and inhomogeneous radial distributions. In the absence of wind shear, a simpler entrainment-rate closure is proposed, which is then further constrained by systematically varying the plume and environmental conditions and allowing evaporative cooling to occur. The fractional mixing rate is shown to be strongly correlated with the plume buoyancy and, to a lesser extent, to the critical mixing fraction (i.e., the fraction of dry air that needs to be mixed with cloudy air to make the mixture neutrally buoyant). Quantitative estimates of this dependency are given to facilitate implementations of the new model in convection parameterizations. Analyzing the proposed closure suggests that it could capture features observed in recent high-resolution simulations and that it is consistent with the buoyancy-sorting concept. The results therefore support recent findings concerning the parameterization of entrainment for moist atmospheric convection.NERC NE/N013727/

The chance of freezing – a conceptional study to parameterize temperature-dependent freezing by including randomness of ice-nucleating particle concentrations

Ice-nucleating particle concentrations (INPCs) can spread over several orders of magnitude at any given temperature. However, this variability is rarely accounted for in heterogeneous ice-nucleation parameterizations. In this paper, we present an approach to incorporate the random variation in the INPC into the parameterization of immersion freezing and analyze this novel concept with various sensitivity tests. In the new scheme, the INPC is drawn from a relative frequency distribution of cumulative INPCs. At each temperature, this distribution describing the INPCs is expressed as a lognormal frequency distribution. The new parameterization scheme does not require aerosol information from the driving model to represent the heterogeneity of INPCs. The scheme's performance is tested in a large-eddy simulation of a relatively warm Arctic mixed-phase stratocumulus. We find that it leads to reasonable ice masses in the cloud, especially when compared to immersion freezing schemes that yield one fixed INPC per temperature and lead to almost no ice production in the simulated cloud. The scheme is sensitive to the median of the frequency distribution and highly sensitive to the standard deviation of the distribution, as well as to the frequency of drawing a new INPC and the resolution of the model. Generally, a higher probability of drawing large INPCs leads to substantially more ice in the simulated cloud. We expose inherent challenges to introducing such a parameterization and explore possible solutions and potential developments.</p

Intercomparison of Large-Eddy Simulations of Arctic Mixed-Phase Clouds: Importance of Ice Size Distribution Assumptions

Large-eddy simulations of mixed-phase Arctic clouds by 11 different models are analyzed with the goal of improving understanding and model representation of processes controlling the evolution of these clouds. In a case based on observations from the Indirect and Semi-Direct Aerosol Campaign (ISDAC), it is found that ice number concentration, Ni, exerts significant influence on the cloud structure. Increasing Ni leads to a substantial reduction in liquid water path (LWP), in agreement with earlier studies. In contrast to previous intercomparison studies, all models here use the same ice particle properties (i.e., mass-size, mass-fall speed, and mass-capacitance relationships) and a common radiation parameterization. The constrained setup exposes the importance of ice particle size distributions (PSDs) in influencing cloud evolution. A clear separation in LWP and IWP predicted by models with bin and bulk microphysical treatments is documented and attributed primarily to the assumed shape of ice PSD used in bulk schemes. Compared to the bin schemes that explicitly predict the PSD, schemes assuming exponential ice PSD underestimate ice growth by vapor deposition and overestimate mass-weighted fall speed leading to an underprediction of IWP by a factor of two in the considered case. Sensitivity tests indicate LWP and IWP are much closer to the bin model simulations when a modified shape factor which is similar to that predicted by bin model simulation is used in bulk scheme. These results demonstrate the importance of representation of ice PSD in determining the partitioning of liquid and ice and the longevity of mixed-phase clouds

The free troposphere as a potential source of arctic boundary layer aerosol particles

This study investigates aerosol particle transport from the free troposphere to the boundary layer in the summertime high Arctic. Observations from the Arctic Summer Cloud Ocean Study field campaign show several occurrences of high aerosol particle concentrations above the boundary layer top. Large-eddy simulations suggest that when these enhanced aerosol concentrations are present, they can be an important source of aerosol particles for the boundary layer. Most particles are transported to the boundary layer by entrainment. However, it is found that mixed-phase stratocumulus clouds, which often extend into the inversion layer, also can mediate the transport of particles into the boundary layer by activation at cloud top and evaporation below cloud base. Finally, the simulations also suggest that aerosol properties at the surface sometimes may not be good indicators of aerosol properties in the cloud layer

The free troposphere as a potential source of arctic boundary layer aerosol particles

This study investigates aerosol particle transport from the free troposphere to the boundary layer in the summertime high Arctic. Observations from the Arctic Summer Cloud Ocean Study field campaign show several occurrences of high aerosol particle concentrations above the boundary layer top. Large-eddy simulations suggest that when these enhanced aerosol concentrations are present, they can be an important source of aerosol particles for the boundary layer. Most particles are transported to the boundary layer by entrainment. However, it is found that mixed-phase stratocumulus clouds, which often extend into the inversion layer, also can mediate the transport of particles into the boundary layer by activation at cloud top and evaporation below cloud base. Finally, the simulations also suggest that aerosol properties at the surface sometimes may not be good indicators of aerosol properties in the cloud layer