scholarly journals Locally-orthogonal unstructured grid-generation for general circulation modelling on the sphere*

2016 ◽  
Author(s):  
Darren Engwirda

Abstract. An algorithm for the generation of non-uniform, locally-orthogonal staggered unstructured spheroidal grids is described. This technique is designed to generate high-quality staggered Voronoi/Delaunay dual meshes appropriate for general circulation modelling on the sphere, including applications to atmospheric simulation, ocean-modelling and numerical weather prediction. Using a recently developed Frontal-Delaunay refinement technique, a method for the construction of unstructured spheroidal Delaunay triangulations is introduced. A locally-orthogonal polygonal grid, derived from the associated Voronoi diagram, is computed as the staggered dual. It is shown that use of the Delaunay-refinement technique allows for the generation of unstructured grids that satisfy a priori constraints on minimum mesh-quality. The initial staggered Voronoi/Delaunay tessellation is iteratively improved through hill-climbing optimisation techniques. Such an approach is shown to produce grids with very high element quality and smooth grading characteristics, while imposing relatively low computational expense. Initial results are presented for a selection of uniform and non-uniform spheroidal grids appropriate for high-resolution, multi-scale general circulation modelling. The use of user-defined mesh-spacing functions to generate smoothly graded, non-uniform grids for multi-resolution type studies is discussed in detail.

2017 ◽  
Vol 10 (6) ◽  
pp. 2117-2140 ◽  
Author(s):  
Darren Engwirda

Abstract. An algorithm for the generation of non-uniform, locally orthogonal staggered unstructured spheroidal grids is described. This technique is designed to generate very high-quality staggered Voronoi–Delaunay meshes appropriate for general circulation modelling on the sphere, including applications to atmospheric simulation, ocean-modelling and numerical weather prediction. Using a recently developed Frontal-Delaunay refinement technique, a method for the construction of high-quality unstructured spheroidal Delaunay triangulations is introduced. A locally orthogonal polygonal grid, derived from the associated Voronoi diagram, is computed as the staggered dual. It is shown that use of the Frontal-Delaunay refinement technique allows for the generation of very high-quality unstructured triangulations, satisfying a priori bounds on element size and shape. Grid quality is further improved through the application of hill-climbing-type optimisation techniques. Overall, the algorithm is shown to produce grids with very high element quality and smooth grading characteristics, while imposing relatively low computational expense. A selection of uniform and non-uniform spheroidal grids appropriate for high-resolution, multi-scale general circulation modelling are presented. These grids are shown to satisfy the geometric constraints associated with contemporary unstructured C-grid-type finite-volume models, including the Model for Prediction Across Scales (MPAS-O). The use of user-defined mesh-spacing functions to generate smoothly graded, non-uniform grids for multi-resolution-type studies is discussed in detail.


2015 ◽  
Vol 72 (1) ◽  
pp. 55-74 ◽  
Author(s):  
Qiang Deng ◽  
Boualem Khouider ◽  
Andrew J. Majda

Abstract The representation of the Madden–Julian oscillation (MJO) is still a challenge for numerical weather prediction and general circulation models (GCMs) because of the inadequate treatment of convection and the associated interactions across scales by the underlying cumulus parameterizations. One new promising direction is the use of the stochastic multicloud model (SMCM) that has been designed specifically to capture the missing variability due to unresolved processes of convection and their impact on the large-scale flow. The SMCM specifically models the area fractions of the three cloud types (congestus, deep, and stratiform) that characterize organized convective systems on all scales. The SMCM captures the stochastic behavior of these three cloud types via a judiciously constructed Markov birth–death process using a particle interacting lattice model. The SMCM has been successfully applied for convectively coupled waves in a simplified primitive equation model and validated against radar data of tropical precipitation. In this work, the authors use for the first time the SMCM in a GCM. The authors build on previous work of coupling the High-Order Methods Modeling Environment (HOMME) NCAR GCM to a simple multicloud model. The authors tested the new SMCM-HOMME model in the parameter regime considered previously and found that the stochastic model drastically improves the results of the deterministic model. Clear MJO-like structures with many realistic features from nature are reproduced by SMCM-HOMME in the physically relevant parameter regime including wave trains of MJOs that organize intermittently in time. Also one of the caveats of the deterministic simulation of requiring a doubling of the moisture background is not required anymore.


2007 ◽  
Vol 7 (13) ◽  
pp. 3519-3536 ◽  
Author(s):  
A. Gobiet ◽  
G. Kirchengast ◽  
G. L. Manney ◽  
M. Borsche ◽  
C. Retscher ◽  
...  

Abstract. This study describes and evaluates a Global Navigation Satellite System (GNSS) radio occultation (RO) retrieval scheme particularly aimed at delivering bias-free atmospheric parameters for climate monitoring and research. The focus of the retrieval is on the sensible use of a priori information for careful high-altitude initialisation in order to maximise the usable altitude range. The RO retrieval scheme has been meanwhile applied to more than five years of data (September 2001 to present) from the German CHAllenging Minisatellite Payload for geoscientific research (CHAMP) satellite. In this study it was validated against various correlative datasets including the Michelson Interferometer for Passive Atmospheric Sounding (MIPAS) and the Global Ozone Monitoring for Occultation of Stars (GOMOS) sensors on Envisat, five different atmospheric analyses, and the operational CHAMP retrieval product from GeoForschungsZentrum (GFZ) Potsdam. In the global mean within 10 to 30 km altitude we find that the present validation observationally constrains the potential RO temperature bias to be <0.2 K. Latitudinally resolved analyses show biases to be observationally constrained to <0.2–0.5 K up to 35 km in most cases, and up to 30 km in any case, even if severely biased (about 10 K or more) a priori information is used in the high altitude initialisation of the retrieval. No evidence is found for the 10–35 km altitude range of residual RO bias sources other than those potentially propagated downward from initialisation, indicating that the widely quoted RO promise of "unbiasedness and long-term stability due to intrinsic self-calibration" can indeed be realised given care in the data processing to strictly limit structural uncertainty. The results thus reinforce that adequate high-altitude initialisation is crucial for accurate stratospheric RO retrievals. The common method of initialising, at some altitude in the upper stratosphere, the hydrostatic integral with an upper boundary temperature or pressure value derived from meteorological analyses is prone to introduce biases from the upper boundary down to below 25 km. Also above 30 to 35 km, GNSS RO delivers a considerable amount of observed information up to around 40 km, which is particularly interesting for numerical weather prediction (NWP) systems, where direct assimilation of non-initialised observed RO bending angles (free of a priori) is thus the method of choice. The results underline the value of RO for climate applications.


2007 ◽  
Vol 64 (11) ◽  
pp. 3766-3784 ◽  
Author(s):  
Philippe Lopez

Abstract This paper first reviews the current status, issues, and limitations of the parameterizations of atmospheric large-scale and convective moist processes that are used in numerical weather prediction and climate general circulation models. Both large-scale (resolved) and convective (subgrid scale) moist processes are dealt with. Then, the general question of the inclusion of diabatic processes in variational data assimilation systems is addressed. The focus is put on linearity and resolution issues, the specification of model and observation error statistics, the formulation of the control vector, and the problems specific to the assimilation of observations directly affected by clouds and precipitation.


2021 ◽  
Author(s):  
Florence Matutini ◽  
Jacques Baudry ◽  
Marie-Josée Fortin ◽  
Guillaume Pain ◽  
Joséphine Pithon

Abstract Context – Species distribution modelling is a common tool in conservation biology but two main criticisms remain: (1) the use of simplistic variables that do not account for species movements and/or connectivity and (2) poor consideration of multi-scale processes driving species distributions. Objectives – We aimed to determine if including multi-scale and fine-scale movement processes in SDM predictors would improve accuracy of SDM for low-mobility amphibian species over species-level analysis.Methods – We tested and compared different SDMs for nine amphibian species with four different sets of predictors: (1) simple distance-based predictors; (2) single-scale compositional predictors; (3) multi-scale compositional predictors with a priori selection of scale based on knowledge of species mobility and scale-of-effect (4) multi-scale compositional predictors calculated using a friction-based functional grain to account for resource accessibility with landscape resistance to movement.Results - Using friction-based functional grain predictors produced slight to moderate improvements of SDM performance at large scale. The multi-scale approach, with a priori scale selection led to ambiguous results depending on the species studied, in particular for generalist species.Conclusion - We underline the potential of using a friction-based functional grain to improve SDM predictions for species-level analysis.


Author(s):  
Chaojian Chen ◽  
Mikhail Kruglyakov ◽  
Alexey Kuvshinov

Summary Most of the existing three-dimensional (3-D) electromagnetic (EM) modeling solvers based on the integral equation (IE) method exploit fast Fourier transform (FFT) to accelerate the matrix-vector multiplications. This in turn requires a laterally-uniform discretization of the modeling domain. However, there is often a need for multi-scale modeling and inversion, for instance, to properly account for the effects of non-uniform distant structures, and at the same time, to accurately model the effects from local anomalies. In such scenarios, the usage of laterally-uniform grids leads to excessive computational loads, both in terms of memory and time. To alleviate this problem, we developed an efficient 3-D EM modeling tool based on a multi-nested IE approach. Within this approach, the IE modeling is first performed at a large domain and on a (laterally-uniform) coarse grid, and then the results are refined in the region of interest by performing modeling at a smaller domain and on a (laterally-uniform) denser grid. At the latter stage, the modeling results obtained at the previous stage are exploited. The lateral uniformity of the grids at each stage allows us to keep using the FFT for the matrix-vector multiplications. An important novelty of the paper is a development of a “rim domain” concept which further improves the performance of the multi-nested IE approach. We verify the developed tool on both idealized and realistic 3-D conductivity models, and demonstrate its efficiency and accuracy.


2021 ◽  
Author(s):  
Natacha Galmiche ◽  
Nello Blaser ◽  
Morten Brun ◽  
Helwig Hauser ◽  
Thomas Spengler ◽  
...  

&lt;p&gt;Probability distributions based on ensemble forecasts are commonly used to assess uncertainty in weather prediction. However, interpreting these distributions is not trivial, especially in the case of multimodality with distinct likely outcomes. The conventional summary employs mean and standard deviation across ensemble members, which works well for unimodal, Gaussian-like distributions. In the case of multimodality this misleads, discarding crucial information.&amp;#160;&lt;/p&gt;&lt;p&gt;We aim at combining previously developed clustering algorithms in machine learning and topological data analysis to extract useful information such as the number of clusters in an ensemble. Given the chaotic behaviour of the atmosphere, machine learning techniques can provide relevant results even if no, or very little, a priori information about the data is available. In addition, topological methods that analyse the shape of the data can make results explainable.&lt;/p&gt;&lt;p&gt;Given an ensemble of univariate time series, a graph is generated whose edges and vertices represent clusters of members, including additional information for each cluster such as the members belonging to them, their uncertainty, and their relevance according to the graph. In the case of multimodality, this approach provides relevant and quantitative information beyond the commonly used mean and standard deviation approach that helps to further characterise the predictability.&lt;/p&gt;


2012 ◽  
Vol 5 (1) ◽  
pp. 87-110 ◽  
Author(s):  
A. Kerkweg ◽  
P. Jöckel

Abstract. The numerical weather prediction model of the Consortium for Small Scale Modelling (COSMO), maintained by the German weather service (DWD), is connected with the Modular Earth Submodel System (MESSy). This effort is undertaken in preparation of a new, limited-area atmospheric chemistry model. Limited-area models require lateral boundary conditions for all prognostic variables. Therefore the quality of a regional chemistry model is expected to improve, if boundary conditions for the chemical constituents are provided by the driving model in consistence with the meteorological boundary conditions. The new developed model is as consistent as possible, with respect to atmospheric chemistry and related processes, with a previously developed global atmospheric chemistry general circulation model: the ECHAM/MESSy Atmospheric Chemistry (EMAC) model. The combined system constitutes a new research tool, bridging the global to the meso-γ scale for atmospheric chemistry research. MESSy provides the infrastructure and includes, among others, the process and diagnostic submodels for atmospheric chemistry simulations. Furthermore, MESSy is highly flexible allowing model setups with tailor made complexity, depending on the scientific question. Here, the connection of the MESSy infrastructure to the COSMO model is documented and also the code changes required for the generalisation of regular MESSy submodels. Moreover, previously published prototype submodels for simplified tracer studies are generalised to be plugged-in and used in the global and the limited-area model. They are used to evaluate the TRACER interface implementation in the new COSMO/MESSy model system and the tracer transport characteristics, an important prerequisite for future atmospheric chemistry applications. A supplementary document with further details on the technical implementation of the MESSy interface into COSMO with a complete list of modifications to the COSMO code is provided.


2020 ◽  
Vol 375 (1807) ◽  
pp. 20190383 ◽  
Author(s):  
Sara Bernardi ◽  
Marco Scianna

Collective dynamics in animal groups is a challenging theme for the modelling community, being treated with a wide range of approaches. This topic is here tackled by a discrete model. Entering in more details, each agent, represented by a material point, is assumed to move following a first-order Newtonian law, which distinguishes speed and orientation. In particular, the latter results from the balance of a given set of behavioural stimuli, each of them defined by a direction and a weight, that quantifies its relative importance. A constraint on the sum of the weights then avoids implausible simultaneous maximization/minimization of all movement traits. Our framework is based on a minimal set of rules and parameters and is able to capture and classify a number of collective group dynamics emerging from different individual preferred behaviour, which possibly includes attractive, repulsive and alignment stimuli. In the case of a system of animals subjected only to the first two behavioural inputs, we also show how analytical arguments allow us to a priori relate the equilibrium interparticle spacing to critical model coefficients. Our approach is then extended to account for the presence of predators with different hunting strategies, which impact on the behaviour of a prey population. Hints for model refinement and applications are finally given in the conclusive part of the article. This article is part of the theme issue ‘Multi-scale analysis and modelling of collective migration in biological systems’.


Sign in / Sign up

Export Citation Format

Share Document