scholarly journals A N-dimensional Fortran Interpolation Program (NterGeo.v2020a) for Geophysics Sciences – Application to a back-trajectory program (BACKPLUMES.v2020r1) using CHIMERE or WRF outputs

2020 ◽  
Author(s):  
Bertrand Bessagnet ◽  
Laurent Menut ◽  
Maxime Beauchamp

Abstract. An interpolation program coded in Fortran for irregular N-dimensional cases is presented and freely available. Needs of interpolation procedure over irregular meshes or matrixes with interdependent input data dimensions is frequent in geophysical models. Also, these models often embed look-up tables of physics/chemistry modules. Fortran is a powerful and fast language, highly portable and easy to interface with other existing Fortran models. Our program does not need any libraries and can be compiled with any Fortran compiler. The program is fast and competitive compared to current Python libraries. A novel optional parameter (normalisation option) is provided when considering different types of units on each dimension. For the general program, the inverse distance is used for the weight calculation with a distance defined as a p-distance. Some tests and examples are provided and available in the code package. Moreover, a real case of geophysics application embedding this interpolation program is provided and discussed, it consists in determining back-trajectories using atmospheric dispersion or mesoscale meteorological model outputs, respectively from the widely used models CHIMERE and WRF.

2021 ◽  
Vol 14 (1) ◽  
pp. 91-106
Author(s):  
Bertrand Bessagnet ◽  
Laurent Menut ◽  
Maxime Beauchamp

Abstract. An interpolation programme coded in Fortran for irregular N-dimensional cases is presented and freely available. The need for interpolation procedures over irregular meshes or matrixes with interdependent input data dimensions is frequent in geophysical models. Also, these models often embed look-up tables of physics or chemistry modules. Fortran is a fast and powerful language and is highly portable. It is easy to interface models written in Fortran with each other. Our programme does not need any libraries; it is written in standard Fortran and tested with two usual compilers. The programme is fast and competitive compared to current Python libraries. A normalization option parameter is provided when considering different types of units on each dimension. Some tests and examples are provided and available in the code package. Moreover, a geophysical application embedding this interpolation programme is provided and discussed; it consists in determining back trajectories using chemistry-transport or mesoscale meteorological model outputs, respectively, from the widely used CHIMERE and Weather Research and Forecasting (WRF) models.


2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Zhan-Ning Liu ◽  
Xiao-Yan Yu ◽  
Li-Feng Jia ◽  
Yuan-Sheng Wang ◽  
Yu-Chen Song ◽  
...  

AbstractIn order to study the influence of distance weight on ore-grade estimation, the inverse distance weighted (IDW) is used to estimate the Ni grade and MgO grade of serpentinite ore based on a three-dimensional ore body model and related block models. Manhattan distance, Euclidean distance, Chebyshev distance, and multiple forms of the Minkowski distance are used to calculate distance weight of IDW. Results show that using the Minkowski distance for the distance weight calculation is feasible. The law of the estimated results along with the distance weight is given. The study expands the distance weight calculation method in the IDW method, and a new method for improving estimation accuracy is given. Researchers can choose different weight calculation methods according to their needs. In this study, the estimated effect is best when the power of the Minkowski distance is 3 for a 10 m × 10 m × 10 m block model. For a 20 m × 20 m × 20 m block model, the estimated effect is best when the power of the Minkowski distance is 9.


Atmosphere ◽  
2021 ◽  
Vol 12 (3) ◽  
pp. 388
Author(s):  
Andrea Bazzano ◽  
Stefano Bertinetti ◽  
Francisco Ardini ◽  
David Cappelletti ◽  
Marco Grotti

Lead content, enrichment factors, and isotopic composition (208Pb/206Pb and 207Pb/206Pb) measured in atmospheric particulate matter (PM10) samples collected for nine years at Ny-Ålesund (Svalbard islands, Norwegian Arctic) during spring and summer are presented and discussed. The possible source areas (PSA) for particulate inferred from Pb isotope ratio values were compared to cluster analysis of back-trajectories. Results show that anthropogenic Pb dominates over natural crustal Pb, with a recurring higher influence in spring, compared to summer. Crustal Pb accounted for 5–16% of the measured Pb concentration. Anthropogenic Pb was affected by (i) a Central Asian PSA with Pb isotope signature compatible with ores smelted in the Rudny Altai region, at the Russian and Kazakhstan border, which accounted for 85% of the anthropogenic Pb concentration, and (ii) a weaker North American PSA, contributing for the remaining 15%. Central Asian PSA exerted an influence on 71–86% of spring samples, without any significant interannual variation. On the contrary, 59–87% of summer samples were influenced by the North American PSA, with higher contributions during 2015 and 2018. Back-trajectory analysis agreed on the seasonal difference in PSA and highlighted a possible increased influence for North American air masses during summer 2010 and 2018, but not for summer 2015.


2019 ◽  
Vol 19 (6) ◽  
pp. 3557-3578 ◽  
Author(s):  
Martin Osborne ◽  
Florent F. Malavelle ◽  
Mariana Adam ◽  
Joelle Buxmann ◽  
Jaqueline Sugier ◽  
...  

Abstract. On 15–16 October 2017, ex-hurricane Ophelia passed to the west of the British Isles, bringing dust from the Sahara and smoke from Portuguese forest fires that was observable to the naked eye and reported in the UK's national press. We report here detailed observations of this event using the UK operational lidar and sun-photometer network, established for the early detection of aviation hazards, including volcanic ash. We also use ECMWF ERA5 wind field data and MODIS imagery to examine the aerosol transport. The observations, taken continuously over a period of 30 h, show a complex picture, dominated by several different aerosol layers at different times and clearly correlated with the passage of different air masses associated with the intense cyclonic system. A similar evolution was observed at several sites, with a time delay between them explained by their different location with respect to the storm and associated meteorological features. The event commenced with a shallow dust layer at 1–2 km in altitude and culminated in a deep and complex structure that lasted ∼12 h at each site over the UK, correlated with the storm's warm sector. For most of the time, the aerosol detected was dominated by mineral dust mixtures, as highlighted by depolarisation measurements, but an intense biomass burning aerosol (BBA) layer was observed towards the end of the event, lasting around 3 h at each site. The aerosol optical depth at 355 nm (AOD355) during the whole event ranged from 0.2 to 2.9, with the larger AOD correlated to the intense BBA layer. Such a large AOD is unprecedented in the UK according to AERONET records for the last 20 years. The Raman lidars permitted the measurement of the aerosol extinction coefficient at 355 nm, the particle linear depolarisation ratio (PLDR), and the lidar ratio (LR) and made the separation of the dust (depolarising) aerosol from other aerosol types possible. A specific extinction has also been computed to provide an estimate of the atmospheric concentration of both aerosol types separately, which peaked at 420±200 µg m−3 for the dust and 558±232 µg m−3 for the biomass burning aerosols. Back trajectories computed using the Numerical Atmospheric-dispersion Modelling Environment (NAME) were used to identify the sources and strengthen the conclusions drawn from the observations. The UK network represents a significant expansion of the observing capability in northern Europe, with instruments evenly distributed across Great Britain, from Camborne in Cornwall to Lerwick in the Shetland Islands, and this study represents the first attempt to demonstrate its capability and validate the methods in use. Its ultimate purpose will be the detection and quantification of volcanic plumes, but the present study clearly demonstrates the advanced capabilities of the network.


2014 ◽  
Vol 2014 ◽  
pp. 1-19 ◽  
Author(s):  
Kristi A. Gebhart ◽  
William C. Malm ◽  
Marco A. Rodriguez ◽  
Michael G. Barna ◽  
Bret A. Schichtel ◽  
...  

The Rocky Mountain Atmospheric Nitrogen and Sulfur (RoMANS II) study with field operations during November 2008 through November 2009 was designed to evaluate the composition and sources of reactive nitrogen in Rocky Mountain National Park, Colorado, USA. As part of RoMANS II, a mesoscale meteorological model was utilized to provide input for back trajectory and chemical transport models. Evaluation of the model's ability to capture important transport patterns in this region of complex terrain is discussed. Previous source-receptor studies of nitrogen in this region are also reviewed. Finally, results of several back trajectory analyses for RoMANS II are presented. The trajectory mass balance (TrMB) model, a receptor-based linear regression technique, was used to estimate mean source attributions of airborne ammonia concentrations during RoMANS II. Though ammonia concentrations are usually higher when there is transport from the east, the TrMB model estimates that, on average, areas to the west contribute a larger mean fraction of the ammonia. Possible reasons for this are discussed and include the greater frequency of westerly versus easterly winds, the possibility that ammonia is transported long distances as ammonium nitrate, and the difficulty of correctly modeling the transport winds in this area.


2022 ◽  
Vol 3 (1) ◽  
pp. 21-44
Author(s):  
Sonja Murto ◽  
Rodrigo Caballero ◽  
Gunilla Svensson ◽  
Lukas Papritz

Abstract. Atmospheric blocking can influence Arctic weather by diverting the mean westerly flow and steering cyclones polewards, bringing warm, moist air to high latitudes. Recent studies have shown that diabatic heating processes in the ascending warm conveyor belt branch of extratropical cyclones are relevant to blocking dynamics. This leads to the question of the extent to which diabatic heating associated with mid-latitude cyclones may influence high-latitude blocking and drive Arctic warm events. In this study we investigate the dynamics behind 50 extreme warm events of wintertime high-Arctic temperature anomalies during 1979–2016. Classifying the warm events based on blocking occurrence within three selected sectors, we find that 30 of these events are associated with a block over the Urals, featuring negative upper-level potential vorticity (PV) anomalies over central Siberia north of the Ural Mountains. Lagrangian back-trajectory calculations show that almost 60 % of the air parcels making up these negative PV anomalies experience lifting and diabatic heating (median 11 K) in the 6 d prior to the block. Further, almost 70 % of the heated trajectories undergo maximum heating in a compact region of the mid-latitude North Atlantic, temporally taking place between 6 and 1 d before arriving in the blocking region. We also find anomalously high cyclone activity (on average five cyclones within this 5 d heating window) within a sector northwest of the main heating domain. In addition, 10 of the 50 warm events are associated with blocking over Scandinavia. Around 60 % of the 6 d back trajectories started from these blocks experience diabatic heating, of which 60 % undergo maximum heating over the North Atlantic but generally closer to the time of arrival in the block and further upstream relative to heated trajectories associated with Ural blocking. This study suggests that, in addition to the ability of blocks to guide cyclones northwards, Atlantic cyclones play a significant role in the dynamics of high-latitude blocking by providing low-PV air via moist-diabatic processes. This emphasizes the importance of the mutual interactions between mid-latitude cyclones and Eurasian blocking for wintertime Arctic warm extremes.


2020 ◽  
Author(s):  
Catherine Hayer ◽  
Mike Burton

<p>The use of polar-orbiting satellite instruments to monitor volcanoes has been an established technique for decades. However, a major limitation is the temporal resolution provided by these satellite platforms. For UV instruments, one or occasionally two observations per day are possible for tropical latitudes, though an improved temporal resolution is seen at high latitudes. The SO<sub>2</sub> altitude within the atmospheric column is usually highly unconstrained and is one of the largest sources of uncertainty within the SO<sub>2</sub> retrieval. This method assigns a best-fit altitude to each pixel, instead of using a single value for the whole plume.</p><p>TROPOMI is an UV spectrometer, launched on the Sentinel-5P platform in October 2017. The instrument has a swath of 2600 km and a spatial resolution of 5.5x7.5 km (improving to 3.5x7.5 km from August 2019). Sentinel-5P flies with the A-Train constellation, with an equatorial overpass time of 13:30 local time.</p><p>Applying the NOAA HYbrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) back trajectory model, the injection time, injection and measurement altitudes of the SO<sub>2</sub> in each pixel within the satellite image is derived. Back trajectories are run for each pixel at a range of altitudes. The natural variability in the wind field at different altitudes (wind shear) means that only some of those trajectories will return to the volcano, constraining the measurement altitude to those trajectories. The SO<sub>2</sub> concentration is interpolated to this altitude. Finding the point in the trajectory when it most closely approaches the volcano provides the time and altitude of injection.</p><p>Combining the corrected SO<sub>2</sub> concentrations with the injection time produces the SO<sub>2</sub> flux that generated the observed SO<sub>2</sub> cloud, and with the injection altitude to calculate the mass eruption rate. These parameters can also be used to improve eruption plume modelling by improving the constraints on the eruption column characteristics.</p><p>The method is applied to the December 2019 eruption of White Island, New Zealand.</p>


2011 ◽  
Vol 11 (12) ◽  
pp. 5655-5672 ◽  
Author(s):  
I. Xueref-Remy ◽  
C. Messager ◽  
D. Filippi ◽  
M. Pastel ◽  
P. Nedelec ◽  
...  

Abstract. Atmospheric airborne measurements of CO2 are very well suited for estimating the time-varying distribution of carbon sources and sinks at the regional scale due to the large geographical area covered over a short time. We present here an analysis of two cross-European airborne campaigns carried out on 23–26 May 2001 (CAATER-1) and 2–3 October 2002 (CAATER-2) over Western Europe. The area covered during CAATER-1 and CAATER-2 was 4° W to 14° E long; 44° N to 52° N lat and 1° E to 17° E long; 46° N to 52° N lat respectively. High precision in situ CO2, CO and Radon 222 measurements were recorded. Flask samples were collected during both campaigns to cross-validate the in situ data. During CAATER-1 and CAATER-2, the mean CO2 concentration was 370.1 ± 4.0 (1-σ standard deviation) ppm and 371.7 ± 5.0 (1-σ) ppm respectively. A HYSPLIT back-trajectories analysis shows that during CAATER 1, northwesterly winds prevailed. In the planetary boundary layer (PBL) air masses became contaminated over Benelux and Western Germany by emissions from these highly urbanized areas, reaching about 380 ppm. Air masses passing over rural areas were depleted in CO2 because of the photosynthesis activity of the vegetation, with observations as low as 355 ppm. During CAATER-2, the back-trajectory analysis showed that air masses were distributed among the 4 sectors. Air masses were enriched in CO2 and CO over anthropogenic emission spots in Germany but also in Poland, as these countries have part of the most CO2-emitting coal-based plants in Europe. Simultaneous measurements of in situ CO2 and CO combined with back-trajectories helped us to distinguish between fossil fuel emissions and other CO2 sources. The ΔCO/ΔCO2 ratios (R2 = 0.33 to 0.88, slopes = 2.42 to 10.37), calculated for anthropogenic-influenced air masses over different countries/regions matched national inventories quite well, showing that airborne measurements can help to identify the origin of fossil fuel emissions in the PBL even when distanced by several days/hundreds of kms from their sources. We have compared airborne CO2 observations to nearby ground station measurements and thereby, confirmed that measurements taken in the lower few meters of the PBL (low-level ground stations) are representative of the local scale, while those located in the free troposphere (FT) (moutain stations) are representative of atmospheric CO2 regionally on a scale of a few hundred kilometers. Stations located several 100 km away from each other differ from a few ppm in their measurements indicating the existence of a gradient within the free troposphere. Observations at stations located on top of small mountains may match the airborne data if the sampled air comes from the FT rather than coming up from the valley. Finally, the analysis of the CO2 vertical variability conducted on the 14 profiles recorded in each campaign shows a variability at least 5 to 8 times higher in the PBL (the 1-σ standard deviation associated to the CO2 mean of all profiles within the PBL is 4.0 ppm and 5.7 ppm for CAATER-1 and CAATER-2, respectively) than in the FT (within the FT, 1-σ is 0.5 ppm and 1.1 ppm for CAATER-1 and CAATER-2, respectively). The CO2 jump between the PBL and the FT equals 3.7 ppm for the first campaign and −0.3 ppm for the second campaign. A very striking zonal CO2 gradient of about 11 ppm was observed in the mid-PBL during CAATER-2, with higher concentrations in the west than in the east. This gradient may originate from differences in atmospheric mixing, ground emission rates or Autumn's earlier start in the west. More airborne campaigns are currently under analysis in the framework of the CARBOEUROPE-IP project to better assess the likelihood of these different hypotheses. In a companion paper (Xueref-Remy et al., 2011, Part 2), a comparison of vertical profiles from observations and several modeling frameworks was conducted for both campaigns.


Forests ◽  
2021 ◽  
Vol 12 (1) ◽  
pp. 64
Author(s):  
Ana Galveias ◽  
Ana R. Costa ◽  
Daniele Bortoli ◽  
Russell Alpizar-Jara ◽  
Rui Salgado ◽  
...  

Research Highlights: Daily airborne Cupressaceae pollen disruption ranged from 20 to 90%; relative humidity (RH), rainfall and atmospheric pressure (AtP) were the major meteorological determinants of this phenomenon. Background and Objectives: Cupressaceae family includes several species that are widely used as ornamental plants pollinating in late winter-early spring and might be responsible for allergic outbreaks. Cupressaceae pollen disruption may favour allergen dissemination, potentiating its allergenicity. The aim of this work was to characterize the Cupressaceae pollen aerobiology in Évora, South of Portugal, in 2017 and 2018, particularly the pollen disruption, and to identify the meteorological parameters contributing to this phenomenon. Materials and Methods: Pollen was collected using a Hirst type 7-day pollen trap and was identified following the standard methodology. Temperature, RH, rainfall, global solar radiation (Global Srad), AtP, wind speed and direction were obtained from a weather station installed side-by-side to the Hirst platform. Back trajectories (12-h) of air masses arriving at Évora were calculated using the HYSPLIT model. Results: Cupressaceae pollen index was higher in 2017 compared to 2018 (>5994 and 3175 pollen/m3, respectively) and 36 ± 19% (2017) and 59 ± 23% (2018) of the pollen was disrupted. Higher levels of disrupted pollen coincided with RH > 60% and rainfall. Temperature, Global Srad and AtP correlated negatively with pollen disruption. Wind speed and wind direction did not significantly correlate with pollen disruption. Intra-diurnal pollen pattern peaked between 9:00 am–2:00 pm, suggesting local origin, confirmed by the back trajectory analysis. Intra-diurnal pollen disruption profile followed hourly pollen pattern and it negatively correlated with AtP, temperature and Global Srad but was uncorrelated with RH. Conclusions: The results suggest that RH, rainfall and AtP are the main factors affecting airborne Cupressaceae pollen integrity and in conjunction with daily pollen concentration may be used to predict the risk of allergy outbreaks to this pollen type.


2014 ◽  
Vol 7 (1) ◽  
pp. 107-128 ◽  
Author(s):  
S. Freitag ◽  
A. D. Clarke ◽  
S. G. Howell ◽  
V. N. Kapustin ◽  
T. Campos ◽  
...  

Abstract. The history of air masses is often investigated using backward trajectories to gain knowledge about processes along the air parcel path as well as possible source regions. Here, we describe a refined approach that incorporates airborne gas, aerosol, and environmental data into back trajectories and show how this technique allows for simultaneous evaluation of air mass history and back trajectory reliability without the need to calculate trajectory errors. We use the HYbrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) model and add a simple semi-automated computing routine to facilitate high-frequency coverage of back trajectories initiated along free tropospheric (FT) flight tracks and profiles every 10 s. We integrate our in situ physiochemical data by color-coding each of these trajectories with its corresponding in situ tracer values measured at the back trajectory start points along the flight path. The unique color for each trajectory aids assessment of trajectory reliability through the visual clustering of air mass pathways of similar coloration. Moreover, marked changes in trajectories associated with marked changes evident in measured physiochemical or thermodynamic properties of an air mass add credence to trajectories. This is particularly true when these air mass properties are linked to trajectory features characteristic of recognized sources or processes. This visual clustering of air mass pathways is of particular value for large-scale 3-D flight tracks common to aircraft experiments where air mass features of interest are often spatially distributed and temporally separated. The cluster-visualization tool used here reveals that most FT back trajectories with pollution signatures measured in the central equatorial Pacific reach back to sources on the South American continent over 10 000 km away and 12 days back in time, e.g., the Amazonian basin. We also demonstrate the distinctions in air mass properties between these and trajectories that penetrate deep convection in the Inter-Tropical Convergence Zone. Additionally, for the first time we show consistency of modeled precipitation along back trajectories with scavenging signatures in the aerosol measured for these trajectories.


Sign in / Sign up

Export Citation Format

Share Document