scholarly journals Novel statistical emulator construction for volcanic ash transport model Ash3d with physically motivated measures

Author(s):  
Qingyuan Yang ◽  
E. Bruce Pitman ◽  
Elaine Spiller ◽  
Marcus Bursik ◽  
Andrea Bevilacqua

Statistical emulators are a key tool for rapidly producing probabilistic hazard analysis of geophysical processes. Given output data computed for a relatively small number of parameter inputs, an emulator interpolates the data, providing the expected value of the output at untried inputs and an estimate of error at that point. In this work, we propose to fit Gaussian Process emulators to the output from a volcanic ash transport model, Ash3d. Our goal is to predict the simulated volcanic ash thickness from Ash3d at a location of interest using the emulator. Our approach is motivated by two challenges to fitting emulators—characterizing the input wind field and interactions between that wind field and variable grain sizes. We resolve these challenges by using physical knowledge on tephra dispersal. We propose new physically motivated variables as inputs and use normalized output as the response for fitting the emulator. Subsetting based on the initial conditions is also critical in our emulator construction. Simulation studies characterize the accuracy and efficiency of our emulator construction and also reveal its current limitations. Our work represents the first emulator construction for volcanic ash transport models with considerations of the simulated physical process.

2016 ◽  
Vol 16 (14) ◽  
pp. 9189-9200 ◽  
Author(s):  
Guangliang Fu ◽  
Arnold Heemink ◽  
Sha Lu ◽  
Arjo Segers ◽  
Konradin Weber ◽  
...  

Abstract. The forecast accuracy of distal volcanic ash clouds is important for providing valid aviation advice during volcanic ash eruption. However, because the distal part of volcanic ash plume is far from the volcano, the influence of eruption information on this part becomes rather indirect and uncertain, resulting in inaccurate volcanic ash forecasts in these distal areas. In our approach, we use real-life aircraft in situ observations, measured in the northwestern part of Germany during the 2010 Eyjafjallajökull eruption, in an ensemble-based data assimilation system combined with a volcanic ash transport model to investigate the potential improvement on the forecast accuracy with regard to the distal volcanic ash plume. We show that the error of the analyzed volcanic ash state can be significantly reduced through assimilating real-life in situ measurements. After a continuous assimilation, it is shown that the aviation advice for Germany, the Netherlands and Luxembourg can be significantly improved. We suggest that with suitable aircrafts measuring once per day across the distal volcanic ash plume, the description and prediction of volcanic ash clouds in these areas can be greatly improved.


1985 ◽  
Vol 90 (D6) ◽  
pp. 10620-10630 ◽  
Author(s):  
Arthur T. Hopkins ◽  
Charles J. Bridgman

2016 ◽  
Vol 144 (2) ◽  
pp. 575-589 ◽  
Author(s):  
S. Lu ◽  
H. X. Lin ◽  
A. W. Heemink ◽  
G. Fu ◽  
A. J. Segers

Abstract Volcanic ash forecasting is a crucial tool in hazard assessment and operational volcano monitoring. Emission parameters such as plume height, total emission mass, and vertical distribution of the emission plume rate are essential and important in the implementation of volcanic ash models. Therefore, estimation of emission parameters using available observations through data assimilation could help to increase the accuracy of forecasts and provide reliable advisory information. This paper focuses on the use of satellite total-ash-column data in 4D-Var based assimilations. Experiments show that it is very difficult to estimate the vertical distribution of effective volcanic ash injection rates from satellite-observed ash columns using a standard 4D-Var assimilation approach. This paper addresses the ill-posed nature of the assimilation problem from the perspective of a spurious relationship. To reduce the influence of a spurious relationship created by a radiate observation operator, an adjoint-free trajectory-based 4D-Var assimilation method is proposed, which is more accurate to estimate the vertical profile of volcanic ash from volcanic eruptions. The method seeks the optimal vertical distribution of emission rates of a reformulated cost function that computes the total difference between simulated and observed ash columns. A 3D simplified aerosol transport model and synthetic satellite observations are used to compare the results of both the standard method and the new method.


2011 ◽  
Vol 11 (5) ◽  
pp. 1961-1977 ◽  
Author(s):  
J. Flemming ◽  
A. Inness ◽  
L. Jones ◽  
H. J. Eskes ◽  
V. Huijnen ◽  
...  

Abstract. The 2008 Antarctic ozone hole was one of the largest and most long-lived in recent years. Predictions of the ozone hole were made in near-real time (NRT) and hindcast mode with the Integrated Forecast System (IFS) of the European Centre for Medium-Range Weather Forecasts (ECMWF). The forecasts were carried out both with and without assimilation of satellite observations from multiple instruments to provide more realistic initial conditions. Three different chemistry schemes were applied for the description of stratospheric ozone chemistry: (i) a linearization of the ozone chemistry, (ii) the stratospheric chemical mechanism of the Model of Ozone and Related Chemical Tracers, version 3, (MOZART-3) and (iii) the relaxation to climatology as implemented in the Transport Model, version 5, (TM5). The IFS uses the latter two schemes by means of a two-way coupled system. Without assimilation, the forecasts showed model-specific shortcomings in predicting start time, extent and duration of the ozone hole. The assimilation of satellite observations from the Microwave Limb Sounder (MLS), the Ozone Monitoring Instrument (OMI), the Solar Backscattering Ultraviolet radiometer (SBUV-2) and the SCanning Imaging Absorption spectroMeter for Atmospheric CartograpHY (SCIAMACHY) led to a significant improvement of the forecasts when compared with total columns and vertical profiles from ozone sondes. The combined assimilation of observations from multiple instruments helped to overcome limitations of the ultraviolet (UV) sensors at low solar elevation over Antarctica. The assimilation of data from MLS was crucial to obtain a good agreement with the observed ozone profiles both in the polar stratosphere and troposphere. The ozone analyses by the three model configurations were very similar despite the different underlying chemistry schemes. Using ozone analyses as initial conditions had a very beneficial but variable effect on the predictability of the ozone hole over 15 days. The initialized forecasts with the MOZART-3 chemistry produced the best predictions of the increasing ozone hole whereas the linear scheme showed the best results during the ozonehole closure.


2013 ◽  
Vol 48 (3) ◽  
pp. 232-242
Author(s):  
Ian H. Halket ◽  
Peter F. Rasmussen ◽  
John C. Doering

One-dimensional substance transport models assume that the river reach modelled has a uniform cross-sectional shape which manifests as a constant average velocity in the model equations. Rarely do rivers meet this criterion. Their channels are seldom uniform in shape but rather alternate in a quasi-periodic manner between pool and riffle sections. This bedform sequencing imparts a corresponding variation in the average cross-sectional velocity which is not accounted for in constant velocity transport models. The literature points out that the pool and riffle planform may be the reason for the sometimes poor predictions obtained from these models. This paper presents a new variable velocity transport model and confirms that the fluctuation in average cross-sectional velocity caused by the pool and riffle planform does have a marked effect on transport in rivers. The pool and riffle planform promotes an enhanced decay of a substance when a first-order biochemical reaction is simulated with the new transport equation. Investigation of the analytical solution shows that the enhanced decay is the result of the overall lower velocity experienced in a pool and riffle channel as opposed to a uniform channel. This difference in transport velocity between a pool and riffle channel and a uniform channel becomes more pronounced as flow declines a critical finding for total maximum daily load calculations because these regulatory limits are usually determined for low flow levels by models that do not account for this phenomenon.


Author(s):  
Scott D. Chambers ◽  
Elise-Andree Guérette ◽  
Khalia Monk ◽  
Alan D. Griffiths ◽  
Yang Zhang ◽  
...  

We propose a new technique to prepare statistically-robust benchmarking data for evaluating chemical transport model meteorology and air quality parameters within the urban boundary layer. The approach employs atmospheric class-typing, using nocturnal radon measurements to assign atmospheric mixing classes, and can be applied temporally (across the diurnal cycle), or spatially (to create angular distributions of pollutants as a top-down constraint on emissions inventories). In this study only a short (<1-month) campaign is used, but grouping of the relative mixing classes based on nocturnal mean radon concentrations can be adjusted according to dataset length (i.e., number of days per category), or desired range of within-class variability. Calculating hourly distributions of observed and simulated values across diurnal composites of each class-type helps to: (i) bridge the gap between scales of simulation and observation, (ii) represent the variability associated with spatial and temporal heterogeneity of sources and meteorology without being confused by it, and (iii) provide an objective way to group results over whole diurnal cycles that separates ‘natural complicating factors’ (synoptic non-stationarity, rainfall, mesoscale motions, extreme stability, etc.) from problems related to parameterizations, or between-model differences. We demonstrate the utility of this technique using output from a suite of seven contemporary regional forecast and chemical transport models. Meteorological model skill varied across the diurnal cycle for all models, with an additional dependence on the atmospheric mixing class that varied between models. From an air quality perspective, model skill regarding the duration and magnitude of morning and evening “rush hour” pollution events varied strongly as a function of mixing class. Model skill was typically the lowest when public exposure would have been the highest, which has important implications for assessing potential health risks in new and rapidly evolving urban regions, and also for prioritizing the areas of model improvement for future applications.


2016 ◽  
Author(s):  
Andreas Ostler ◽  
Ralf Sussmann ◽  
Prabir K. Patra ◽  
Sander Houweling ◽  
Marko De Bruine ◽  
...  

Abstract. The distribution of methane (CH4) in the stratosphere can be a major driver of spatial variability in the dry-air column-averaged CH4 mixing ratio (XCH4), which is being measured increasingly for the assessment of CH4 surface emissions. Chemistry-transport models (CTMs) therefore need to simulate the tropospheric and stratospheric fractional columns of XCH4 accurately for estimating surface emissions from XCH4. Simulations from three CTMs are tested against XCH4 observations from the Total Carbon Column Network (TCCON). We analyze how the model-TCCON agreement in XCH4 depends on the model representation of stratospheric CH4 distributions. Model equivalents of TCCON XCH4 are computed with stratospheric CH4 fields from both the model simulations and from satellite-based CH4 distributions from MIPAS (Michelson Interferometer for Passive Atmospheric Sounding) and MIPAS CH4 fields adjusted to ACE-FTS (Atmospheric Chemistry Experiment Fourier Transform Spectrometer) observations. In comparison to simulated model fields we find an improved model-TCCON XCH4 agreement for all models with MIPAS-based stratospheric CH4 fields. For the Atmospheric Chemistry Transport Model (ACTM) the average XCH4 bias is significantly reduced from 38.1 ppb to 13.7 ppb, whereas small improvements are found for the models TM5 (Transport Model, version 5; from 8.7 ppb to 4.3 ppb), and LMDz (Laboratoire de Météorologie Dynamique model with Zooming capability; from 6.8 ppb to 4.3 ppb), respectively. MIPAS stratospheric CH4 fields adjusted to ACE-FTS reduce the average XCH4 bias for ACTM (3.3 ppb), but increase the average XCH4 bias for TM5 (10.8 ppb) and LMDz (20.0 ppb). These findings imply that the range of satellite-based stratospheric CH4 is insufficient to resolve a possible stratospheric contribution to differences in total column CH4 between TCCON and TM5 or LMDz. Applying transport diagnostics to the models indicates that model-to-model differences in the simulation of stratospheric transport, notably the age of stratospheric air, can largely explain the inter-model spread in stratospheric CH4 and, hence, its contribution to XCH4. This implies that there is a need to better understand the impact of individual model transport components (e.g., physical parameterization, meteorological data sets, model horizontal/vertical resolution) on modeled stratospheric CH4.


2021 ◽  
Vol 9 ◽  
Author(s):  
Allan Fries ◽  
Jonathan Lemus ◽  
Paul A. Jarvis ◽  
Amanda B. Clarke ◽  
Jeremy C. Phillips ◽  
...  

Settling-driven gravitational instabilities observed at the base of volcanic ash clouds have the potential to play a substantial role in volcanic ash sedimentation. They originate from a narrow, gravitationally unstable region called a Particle Boundary Layer (PBL) that forms at the lower cloud-atmosphere interface and generates downward-moving ash fingers that enhance the ash sedimentation rate. We use scaled laboratory experiments in combination with particle imaging and Planar Laser Induced Fluorescence (PLIF) techniques to investigate the effect of particle concentration on PBL and finger formation. Results show that, as particles settle across an initial density interface and are incorporated within the dense underlying fluid, the PBL grows below the interface as a narrow region of small excess density. This detaches upon reaching a critical thickness, that scales with (ν2/g′)1/3, where ν is the kinematic viscosity and g′ is the reduced gravity of the PBL, leading to the formation of fingers. During this process, the fluid above and below the interface remains poorly mixed, with only small quantities of the upper fluid phase being injected through fingers. In addition, our measurements confirm previous findings over a wider set of initial conditions that show that both the number of fingers and their velocity increase with particle concentration. We also quantify how the vertical particle mass flux below the particle suspension evolves with time and with the particle concentration. Finally, we identify a dimensionless number that depends on the measurable cloud mass-loading and thickness, which can be used to assess the potential for settling-driven gravitational instabilities to form. Our results suggest that fingers from volcanic clouds characterised by high ash concentrations not only are more likely to develop, but they are also expected to form more quickly and propagate at higher velocities than fingers associated with ash-poor clouds.


Sign in / Sign up

Export Citation Format

Share Document