scholarly journals Tracers and traceability: implementing the cirrus parameterisation from LACM in the TOMCAT/SLIMCAT chemistry transport model as an example of the application of quality assurance to legacy models

2009 ◽  
Vol 2 (2) ◽  
pp. 1299-1333
Author(s):  
A. M. Horseman ◽  
A. R. MacKenzie ◽  
M. P. Chipperfield

Abstract. A new modelling tool for the investigation of large-scale behaviour of cirrus clouds has been developed. This combines two existing models, the TOMCAT/SLIMCAT chemistry transport model (nupdate library version 0.80, script mpc346_l) and cirrus parameterisation of Ren and MacKenzie (LACM implementation not versioned). The development process employed a subset of best-practice software engineering and quality assurance processes, selected to be viable for small-scale projects whilst maintaining the same traceability objectives. The application of the software engineering and quality control processes during the development has been shown to be not a great overhead, and their use has been of benefit to the developers as well as the end users of the results. We provide a step-by-step guide to the implementation of traceability tailored to the production of geo-scientific research software, as distinct from commercial and operational software. Our recommendations include: maintaining a living "requirements list"; explicit consideration of unit, integration and acceptance testing; and automated revision/configuration control, including control of analysis tool scripts and programs. Initial testing of the resulting model against satellite and in-situ measurements has been promising. The model produces representative results for both spatial distribution of the frequency of occurrence of cirrus ice, and the drying of air as it moves across the tropical tropopause. The model is now ready for more rigorous quantitative testing, but will require the addition of a vertical wind velocity downscaling scheme to better represent extra-tropical continental cirrus.

2010 ◽  
Vol 3 (1) ◽  
pp. 189-203
Author(s):  
A. M. Horseman ◽  
A. R. MacKenzie ◽  
M. P. Chipperfield

Abstract. A new modelling tool for the investigation of large-scale behaviour of cirrus clouds has been developed. This combines two existing models, the TOMCAT/SLIMCAT chemistry transport model (nupdate library version 0.80, script mpc346_l) and cirrus parameterisation of Ren and MacKenzie (LACM implementation not versioned). The development process employed a subset of best-practice software engineering and quality assurance processes, selected to be viable for small-scale projects whilst maintaining the same traceability objectives. The application of the software engineering and quality control processes during the development has been shown to be not a great overhead, and their use has been of benefit to the developers as well as the end users of the results. We provide a step-by-step guide to the implementation of traceability tailored to the production of geo-scientific research software, as distinct from commercial and operational software. Our recommendations include: maintaining a living "requirements list"; explicit consideration of unit, integration and acceptance testing; and automated revision/configuration control, including control of analysis tool scripts and programs. Initial testing of the resulting model against satellite and in-situ measurements has been promising. The model produces representative results for both spatial distribution of the frequency of occurrence of cirrus ice, and the drying of air as it moves across the tropical tropopause. The model is now ready for more rigorous quantitative testing, but will require the addition of a vertical wind velocity downscaling scheme to better represent extra-tropical continental cirrus.


2021 ◽  
Vol 21 (12) ◽  
pp. 9585-9607
Author(s):  
Nuria Pilar Plaza ◽  
Aurélien Podglajen ◽  
Cristina Peña-Ortiz ◽  
Felix Ploeger

Abstract. We investigate the influence of different chemical and physical processes on the water vapour distribution in the lower stratosphere (LS), in particular in the Asian and North American monsoon anticyclones (AMA and NAMA, respectively). Specifically, we use the chemistry transport model CLaMS to analyse the effects of large-scale temperatures, methane oxidation, ice microphysics, and small-scale atmospheric mixing processes in different model experiments. All these processes hydrate the LS and, particularly, the AMA. While ice microphysics has the largest global moistening impact, it is small-scale mixing which dominates the specific signature in the AMA in the model experiments. In particular, the small-scale mixing parameterization strongly contributes to the water vapour transport to this region and improves the simulation of the intra-seasonal variability, resulting in a better agreement with the Aura Microwave Limb Sounder (MLS) observations. Although none of our experiments reproduces the spatial pattern of the NAMA as seen in MLS observations, they all exhibit a realistic annual cycle and intra-seasonal variability, which are mainly controlled by large-scale temperatures. We further analyse the sensitivity of these results to the domain-filling trajectory set-up, here-called Lagrangian trajectory filling (LTF). Compared with MLS observations and with a multiyear reference simulation using the full-blown chemistry transport model version of CLaMS, we find that the LTF schemes result in a drier global LS and in a weaker water vapour signal over the monsoon regions, which is likely related to the specification of the lower boundary condition. Overall, our results emphasize the importance of subgrid-scale mixing and multiple transport pathways from the troposphere in representing water vapour in the AMA.


2007 ◽  
Vol 7 (3) ◽  
pp. 6603-6629 ◽  
Author(s):  
I. Pisso ◽  
B. Legras

Abstract. Vertical (cross-isentropic) mixing is produced by small-scale turbulent processes which are still poorly understood and parametrized in numerical models. In this work we provide estimates of local equivalent diffusion in the lower stratosphere by comparing balloon borne high-resolution measurements of chemical tracers with reconstructed mixing ratio from large ensembles of random Lagrangian backward trajectories using European Center for Medium-range Weather Forecasts analysed winds and a chemistry-transport model (REPROBUS). We have investigated cases in subtropical latitudes using data from HIBISCUS campaign. Upper bound on the vertical diffusivity is found to be of the order of 0.5 m2 s−1 in the subtropical region, which is larger than the estimates at higher latitudes. The relation between diffusion and dispersion is studied by estimating Lyapunov exponents and studying their variation according to the presence of active dynamical structures.


2006 ◽  
Vol 41 (1) ◽  
pp. 24-36 ◽  
Author(s):  
Karl-Erich Lindenschmidt ◽  
René Wodrich ◽  
Cornelia Hesse

Abstract A hypothesis stating that more complex descriptions of processes in models simulate reality better (less error) but with more unreliable predictability (more sensitivity) is tested using a river water quality model. This hypothesis was extended stating that applying the model on a domain of smaller scale requires greater complexity to capture the same accuracy as in large-scale model applications which, however, leads to increased model sensitivity. The sediment and pollutant transport model TOXI, a module in the WASP5 package, was applied to two case studies of different scale: a 90-km course of the 5th order (sensu Strahler 1952) lower Saale river, Germany (large scale), and the lock-and-weir system at Calbe (small scale) situated on the same river course. A sensitivity analysis of several parameters relating to the physical and chemical transport processes of suspended solids, chloride, arsenic, iron and zinc shows that the coefficient, which partitions the total heavy metal mass into its dissolved and sorbed fraction, is a very sensitive parameter. Hence, the complexity of the sorptive process was varied to test the hypotheses.


2015 ◽  
Vol 8 (2) ◽  
pp. 205-220
Author(s):  
A. Praga ◽  
D. Cariolle ◽  
L. Giraud

Abstract. To exploit the possibilities of parallel computers, we designed a large-scale bidimensional atmospheric advection model named Pangolin. As the basis for a future chemistry-transport model, a finite-volume approach for advection was chosen to ensure mass preservation and to ease parallelization. To overcome the pole restriction on time steps for a regular latitude–longitude grid, Pangolin uses a quasi-area-preserving reduced latitude–longitude grid. The features of the regular grid are exploited to reduce the memory footprint and enable effective parallel performances. In addition, a custom domain decomposition algorithm is presented. To assess the validity of the advection scheme, its results are compared with state-of-the-art models on algebraic test cases. Finally, parallel performances are shown in terms of strong scaling and confirm the efficient scalability up to a few hundred cores.


2021 ◽  
Author(s):  
Juan Cuesta ◽  
Lorenzo Costantino ◽  
Matthias Beekmann ◽  
Guillaume Siour ◽  
Laurent Menut ◽  
...  

Abstract. We present a comprehensive study integrating satellite observations of ozone pollution, in situ measurements and chemistry transport model simulations for quantifying the role of anthropogenic emission reductions during the COVID-19 lockdown in spring 2020 over Europe. Satellite observations are derived from the IASI+GOME2 multispectral synergism, which provides particularly enhanced sensitivity to near-surface ozone pollution. These observations are first analysed in terms of differences between the average on 1–15 April 2020, when the strictest lockdown restrictions took place, and the same period in 2019. They show clear enhancements of near-surface ozone in Central Europe and Northern Italy, and some other hotspots, which are typically characterized by VOC-limited chemical regimes. An overall reduction of ozone is observed elsewhere, where ozone chemistry is limited by the abundance of NOx. The spatial distribution of positive and negative ozone concentration anomalies observed from space is in relatively good quantitative agreement with surface in situ measurements over the continent (a correlation coefficient of 0.55, a root-mean-squared difference of 11 ppb and the same standard deviation and range of variability). An average bias of ∼8 ppb between the two observational datasets is remarked, which can partly be explained by the fact the satellite approach retrieves partial columns of ozone with a peak sensitivity above the surface (near 2 km of altitude). For assessing the impact of the reduction of anthropogenic emissions during the lockdown, we adjust the satellite and in situ surface observations for withdrawing the influence of meteorological conditions in 2020 and 2019. This adjustment is derived from the chemistry transport model simulations using the meteorological fields of each year and identical emission inventories. This observational estimate of the influence of lockdown emission reduction is consistent for both datasets. They both show lockdown-associated ozone enhancements in hotspots over Central Europe and Northern Italy, with a reduced amplitude with respect to the total changes observed between the two years, and an overall reduction elsewhere over Europe and the ocean. Satellite observations additionally highlight the ozone anomalies in the regions remote from in situ sensors, an enhancement over the Mediterranean likely associated with maritime traffic emissions and a marked large-scale reduction of ozone elsewhere over ocean (particularly over the North Sea), in consistency with previous assessments done with ozonesondes measurements in the free troposphere. These observational assessments are compared with model-only estimations, using the CHIMERE chemistry transport model. For analysing the uncertainty of the model estimates, we perform two sets of simulations with different setups, differing in the emission inventories, their modifications to account for changes in anthropogenic activities during the lockdown and the meteorological fields. Whereas a general qualitative consistency of positive and negative ozone anomalies is remarked between all model and observational estimates, significant changes are seen in their amplitudes. Models underestimate the range of variability of the ozone changes by at least a factor 2 with respect to the two observational data sets, both for enhancements and decreases of ozone, while the large-scale ozone decrease is not simulated. With one of the setups, the model simulates ozone enhancements a factor 3 to 6 smaller than with the other configuration. This is partly linked to the emission inventories of ozone precursors (at least a 30 % difference), but mainly to differences in vertical mixing of atmospheric constituents depending on the choice of the meteorological model.


Author(s):  
Benjamin Metha ◽  
Michele Trenti ◽  
Tingjin Chu

Abstract Thanks to recent advances in integral field spectroscopy (IFS), modern surveys of nearby galaxies are capable of resolving metallicity maps of H ii regions down to scales of ∼50pc. However, statistical analysis of these metallicity maps has seldom gone beyond fitting basic linear regressions and comparing parameters to global galaxy properties. In this paper (the first of a series), we introduce techniques from spatial statistics that are well suited for detailed analysis of both small- and large-scale metallicity variations within the interstellar media (ISMs) of local galaxies. As a first application, we compare the observed structure of small-scale metallicity fluctuations within 7 local galaxies observed by the PHANGS collaboration to predictions from a stochastic, physically motivated, analytical model developed by Krumholz & Ting. We show that while the theoretical model underestimates the amount of correlated scatter in the galactic metallicity distributions by 3 − 4 orders of magnitude, it provides good estimates of the physical scale of metallicity correlations. We conclude that the ISM of local spiral galaxies is far from homogeneous, with regions of size ∼1 kpc showing significant departures from the mean metallicity at each galactocentric radius.


2006 ◽  
Vol 6 (7) ◽  
pp. 1853-1864 ◽  
Author(s):  
A. Hodzic ◽  
R. Vautard ◽  
H. Chepfer ◽  
P. Goloub ◽  
L. Menut ◽  
...  

Abstract. This study describes the atmospheric aerosol load encountered during the large-scale pollution episode that occurred in August 2003, by means of the aerosol optical thicknesses (AOTs) measured at 865 nm by the Polarization and Directionality of the Earth's Reflectances (POLDER) sensor and the simulation by the CHIMERE chemistry-transport model. During this period many processes (stagnation, photochemistry, forest fires) led to unusually high particle concentrations and optical thicknesses. The observed/simulated AOT comparison helps understanding the ability of the model to reproduce most of the gross AOT features observed in satellite data, with a general agreement within a factor 2 and correlations in the 0.4–0.6 range. However some important aerosol features are missed when using regular anthropogenic sources. Additional simulations including emissions and high-altitude transport of smoke from wildfires that occurred in Portugal indicate that these processes could dominate the AOT signal in some areas. Our results also highlight the difficulties of comparing simulated and POLDER-derived AOTs due to large uncertainties in both cases. Observed AOT values are significantly lower than the simulated ones (30–50%). Their comparison with the ground-based Sun photometer Aerosol Robotic Network (AERONET) measurements suggests, for the European sites considered here, an underestimation of POLDER-derived aerosol levels with a factor between 1 and 2. AERONET AOTs compare better with simulations (no particular bias) than POLDER AOTs.


2005 ◽  
Vol 5 (7) ◽  
pp. 1815-1834 ◽  
Author(s):  
L. J. Labrador ◽  
R. von Kuhlmann ◽  
M. G. Lawrence

Abstract. The impact of different assumptions concerning the source magnitude as well as the vertical placement of lightning-produced nitrogen oxides is studied using the global chemistry transport model MATCH-MPIC. The responses of NOx, O3, OH, HNO3 and peroxyacetyl-nitrate (PAN) are investigated. A marked sensitivity to both parameters is found. NOx burdens globally can be enhanced by up to 100% depending on the vertical placement and source magnitude strength. In all cases, the largest enhancements occur in the tropical upper troposphere, where lifetimes of most trace gases are longer and where they thus become more susceptible to long-range transport by large-scale circulation patterns. Comparison with observations indicate that 0 and 20 Tg(N)/yr production rates of NOx from lightning are too low and too high, respectively. However, no single intermediate production rate or vertical distribution can be singled out as best fitting the observations, due to the large scatter in the datasets. This underscores the need for further measurement campaigns in key regions, such as the tropical continents.


Sign in / Sign up

Export Citation Format

Share Document