Numerical modeling of geologic processes

Numerical modeling of geologic processes

As part of my volcano research I have used numerical tools of solid and fluid mechanics that have helped me to model geodetic observations in active volcanoes. These tools are required to udnerstand to complexity of geological processess. I have worked with a variety of models that include elastic, viscoelastic, magma physics in a variety of software including Pylith, COMSOL Multiphyics, and tools that I have developed in Matlab. I have also developed new analytic models of magma transfer.

Physics-based models of effusive eruptions

Geodetic and extruded volume data sets acquired at both silicic and basaltic volcanoes show quasi-exponential trends that are indicative of coupled conduit flow and magma reservoir pressure drop. I am developing a numerical model that couples the magma extrusion, magma reservoir deflation and including time-dependent gas exsolution and viscosity for the effusive phase of the 2011-2012 rhyodacitic eruption of Cordon Caulle volcano. The major goal is to predict the temporal evolution of one of these rare but powerful and dangerous eruptions, and to provide insights on the magma properties (dissolved and exsolved water, viscosity, crystal content) from geodetic and topographic data only. The latter is a key advance because the vast majority of the models used by volcano geodesists do not provide insights on the magma physical properties, which are key parameters that control the erupted volume, eruption duration, and transitions between effusive and explosive eruptions.

Numerical simulation of a physics-based model with variable amounts of complexity (constant magma properties, constant magma properties with the effect of an increasing load due to lava extrusion over the magma reservoir and time dependent magma properties) compared with the pressure drop of a spheroidal magma reservoir inverted from InSAR data and the total extruded volume calculated from a set of TanDEM-X digital elevation models. The model parameters are the pressure drop, total water content and percetange of crystals in the eruptible magma.

Relevant publications

Trapdoor fault activation: a step towards caldera collapse at Sierra Negra, Galápagos, Ecuador.
Shreve, T., Delgado, F.. Journal of Geophysical Research, 2023, 128, e2023JB026437, doi:10.1029/2023JB026437.

Dynamics of episodic magma injection and migration at Yellowstone caldera: revisiting the 2004-2009 episode of caldera uplift with InSAR and GPS.
Delgado, F., Grandin, R. Journal of Geophysical Research: Solid Earth, 126, e2021JB022341. 2021, doi:10.1029/2021JB022341

Physicochemical models of effusive rhyolitic eruptions constrained with InSAR and DEM data: a case study of the 2011-2012 Cordon Caulle eruption
Delgado, F., Kubanek, J., Anderson, K., Lundgren, P., Pritchard, M. Earth and Planetary Science Letters. 2019, 524, doi:10.1016/j.epsl.2019.115736

Renewed post-eruptive uplift following the 2011-2012 rhyolitic eruption of Cordón Caulle (Southern Andes, Chile): evidence for transient episodes of magma reservoir recharge during 2012-2018
Delgado, F., Pritchard, M., Samsonov, S., Córdova, L.. Journal of Geophysical Research. 2018. doi:10.1029/2018JB016240