Researcher profile

Richard F. Katz

Richard F. Katz contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
13works
0followers
8topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

13 published item(s)

preprint2021arXiv

The role of grain-environment heterogeneity in normal grain growth: a stochastic approach

The size distribution of grains is a fundamental characteristic of polycrystalline solids. In the absence of deformation, the grain-size distribution is controlled by normal grain growth. The canonical model of normal grain growth, developed by Hillert, predicts a grain-size distribution that bears a systematic discrepancy with observed distributions. To address this, we propose a change to the Hillert model that accounts for the influence of heterogeneity in the local environment of grains. In our model, each grain evolves in response to its own local environment of neighbouring grains, rather than to the global population of grains. The local environment of each grain evolves according to an Ornstein-Uhlenbeck stochastic process. Our results are consistent with accepted grain-growth kinetics. Crucially, our model indicates that the size of relatively large grains evolves as a random walk due to the inherent variability in their local environments. This leads to a broader grain-size distribution than the Hillert model and indicates that heterogeneity has a critical influence on the evolution of microstructure.

preprint2016arXiv

The role of volatiles in reactive melt transport in the asthenosphere

Experimental studies of mantle petrology find that small concentrations of water and carbon dioxide have a large effect on the solidus temperature and distribution of melting in the upper mantle. However, it has remained unclear what effect small fractions of deep, volatile-rich melts have on melt transport and reactive melting in the shallow asthenosphere. Here we present theory and computations indicating that low-degree, reactive, volatile-rich melts cause channelisation of magmatic flow at depths approximately corresponding to the anhydrous solidus temperature. These results are obtained with a novel method to simulate the thermochemical evolution of the upper mantle in the presence of volatiles. The method uses a thermodynamically consistent framework for reactive, disequilibrium, multi-component melting. It is coupled with a system of equations representing conservation of mass, momentum, and energy for a partially molten grain aggregate. Application of this method in two-phase, three-component upwelling-column models demonstrates that it reproduces leading-order features of hydrated and carbonated peridotite melting; in particular, it captures the production of low-degree, volatile-rich melt at depths far below the volatile-free solidus. The models predict that segregation of volatile-rich, deep melts promotes a reactive channeling instability that creates fast and chemically isolated pathways of melt extraction. Reactive channeling occurs where volatile-rich melts flux the base of the silicate melting region, enhancing dissolution of fusible components from the ambient mantle. We find this effect to be similarly expressed for models of both hydrated and carbonated mantle melting. These findings indicate that despite their small concentrations, water and carbon dioxide have an important control on the extent and style of magma genesis, as well as on the dynamics of melt transport.

preprint2015arXiv

Mathematical Modelling of Tyndall Star Initiation

The superheating that usually occurs when a solid is melted by volumetric heating can produce irregular solid-liquid interfaces. Such interfaces can be visualised in ice, where they are sometimes known as Tyndall stars. This paper describes some of the experimental observations of Tyndall stars and a mathematical model for the early stages of their evolution. The modelling is complicated by the strong crystalline anisotropy, which results in an anisotropic kinetic undercooling at the interface; it leads to an interesting class of free boundary problems that treat the melt region as infinitesimally thin.

preprint2015arXiv

Melt transport rates in heterogeneous mantle beneath mid-ocean ridges

Recent insights to melt migration beneath ridges suggest that channelized flow is a consequence of melting of a heterogeneous mantle, and that spreading rate modulates the dynamics of the localized flow. A corollary of this finding is that both mantle het- erogeneity and spreading rate have implications for the speed and time scale of melt migration. Here, we investigate these implications using numerical models of magma flow in heterogeneous mantle beneath spreading plates. The models predict that a broad distribution of magma flow speeds is characteristic of the sub-ridge mantle. Within the melting region, magmatic flow is fastest in regions of average fusibility; surprisingly, magmas from sources of above-average fusibility travel to the ridge in a longer time. Spreading rate has comparatively simple consequences, mainly resulting in faster segregation speeds at higher spreading rates. The computed time scales are short enough to preserve deep origin 230Th disequilibria and, under favourable parameter regimes, also 226Ra excesses. An important prediction from the models is that mantle heterogeneity induces significant natural variability into U-series disequilibria, complicating the identification of relationships between disequilibria and ridge properties or chemical signatures of heterogeneity.

preprint2015arXiv

Melt-preferred orientation, anisotropic permeability, and melt-band formation in a deforming, partially molten aggregate

Shear deformation of partially molten rock in laboratory experiments causes the emergence of melt-enriched sheets (bands in cross-section) that are aligned at about 15-20 degrees to the shear plane. Deformation and deviatoric stress also cause the coherent alignment of pores at the grain scale. This leads to a melt-preferred orientation that may, in turn, give rise to an anisotropic permeability. Here we develop a simple, general model of anisotropic permeability in partially molten rocks. We use linearised analysis and nonlinear numerical solutions to investigate its behaviour under simple-shear deformation. In particular, we consider implications of the model for the emergence and angle of melt-rich bands. Anisotropic permeability affects the angle of bands and, in a certain parameter regime, it can give rise to low angles consistent with experiments. However, the conditions required for this regime have a narrow range and seem unlikely to be entirely met by experiments. Anisotropic permeability may nonetheless affect melt transport and the behaviour of partially molten rocks in Earth's mantle.

preprint2015arXiv

Three-field block-preconditioners for models of coupled magma/mantle dynamics

For a prescribed porosity, the coupled magma/mantle flow equations can be formulated as a two-field system of equations with velocity and pressure as unknowns. Previous work has shown that while optimal preconditioners for the two-field formulation can be obtained, the construction of preconditioners that are uniform with respect to model parameters is difficult. This limits the applicability of two-field preconditioners in certain regimes of practical interest. We address this issue by reformulating the governing equations as a three-field problem, which removes a term that was problematic in the two-field formulation in favour of an additional equation for a pressure-like field. For the three-field problem, we develop and analyse new preconditioners and we show numerically that they are optimal in terms of problem size and less sensitive to model parameters, compared to the two-field preconditioner. This extends the applicability of optimal preconditioners for coupled mantle/magma dynamics into parameter regimes of physical interest.

preprint2015arXiv

Torsion of a cylinder of partially molten rock with a spherical inclusion: theory and simulation

The processes that are involved in migration and extraction of melt from the mantle are not yet fully understood. Gaining a better understanding of material properties of partially molten rock could help shed light on the behavior of melt on larger scales in the mantle. In this study, we simulate three-dimensional torsional deformation of a partially molten rock that contains a rigid, spherical inclusion. We compare the computed porosity patterns to those found in recent laboratory experiments. The laboratory experiments show emergence of melt-rich bands throughout the rock sample, and pressure shadows around the inclusion. The numerical model displays similar melt-rich bands only for a small bulk-to-shear-viscosity ratio (five or less). The results are consistent with earlier two-dimensional numerical simulations; however, we show that it is easier to form melt-rich bands in three dimensions compared to two. The addition of strain-rate dependence of the viscosity causes a distinct change in the shape of pressure shadows around the inclusion. This change in shape presents an opportunity for experimentalists to identify the strain-rate dependence and therefore the dominant deformation mechanism in torsion experiments with inclusions.

preprint2015arXiv

Variations in mid-ocean ridge CO2 emissions driven by glacial cycles

The geological record shows links between glacial cycles and volcanic productivity, both subaerially and at mid-ocean ridges. Sea-level-driven pressure changes could also affect chemical properties of mid-ocean ridge volcanism. We consider how changing sea-level could alter the CO2 emissions rate from mid-ocean ridges, on both the segment and global scale. We develop a simplified transport model for a highly incompatible element through a homogenous mantle; variations in the melt concentration the emission rate of the element are created by changes in the depth of first silicate melting. The model predicts an average global mid-ocean ridge CO2 emissions-rate of 53 Mt/yr, in line with other estimates. We show that falling sea level would cause an increase in ridge CO2 emissions with a lag of about 100 kyrs after the causative sea level change. The lag and amplitude of the response are sensitive to mantle permeability and plate spreading rate. For a reconstructed sea-level time series of the past million years, we predict variations of up to 12% (7 Mt/yr) in global mid-ocean ridge CO2 emissions. The magnitude and timing of the predicted variations in CO2 emissions suggests a potential role for ridge carbon emissions in glacial cycles.

preprint2014arXiv

An experimental test of the viscous anisotropy hypothesis for partially molten rocks

Chemical differentiation of rocky planets occurs by melt segregation away from the region of melting. The mechanics of this process, however, are complex and incompletely understood. In partially molten rocks undergoing shear deformation, melt pockets between grains align coherently in the stress field; it has been hypothesized that this anisotropy in microstructure creates an anisotropy in the viscosity of the aggregate. With the inclusion of anisotropic viscosity, continuum, two-phase-flow models reproduce the emergence and angle of melt-enriched bands that form in laboratory experiments. In the same theoretical context, these models also predict sample-scale melt migration due to a gradient in shear stress. Under torsional deformation, melt is expected to segregate radially inward. Here we present new torsional deformation experiments on partially molten rocks that test this prediction. Microstructural analyses of the distribution of melt and solid reveal a radial gradient in melt fraction, with more melt toward the centre of the cylinder. The extent of this radial melt segregation grows with progressive strain, consistent with theory. The agreement between theoretical prediction and experimental observation provides a validation of this theory, which is critical to understanding the large-scale geodynamic and geochemical evolution of Earth.

preprint2014arXiv

Analysis of block-preconditioners for models of coupled magma/mantle dynamics

This article considers the iterative solution of a finite element discretisation of the magma dynamics equations. In simplified form, the magma dynamics equations share some features of the Stokes equations. We therefore formulate, analyse and numerically test a Elman, Silvester and Wathen-type block preconditioner for magma dynamics. We prove analytically and demonstrate numerically the optimality of the preconditioner. The presented analysis highlights the dependence of the preconditioner on parameters in the magma dynamics equations that can affect convergence of iterative linear solvers. The analysis is verified through a range of two- and three-dimensional numerical examples on unstructured grids, from simple illustrative problems through to large problems on subduction zone-like geometries. The computer code to reproduce all numerical examples is freely available as supporting material.

preprint2014arXiv

Glacial cycles drive variations in the production of oceanic crust

Glacial cycles redistribute water between oceans and continents causing pressure changes in the upper mantle, with consequences for melting of Earth's interior. Using Plio-Pleistocene sea-level variations as a forcing function, theoretical models of mid-ocean ridge dynamics that include melt transport predict temporal variations in crustal thickness of hundreds of meters. New bathymetry from the Australian-Antarctic ridge shows significant spectral energy near the Milankovitch periods of 23, 41, and 100 ky, consistent with model predictions. These results suggest that abyssal hills, one of the most common bathymetric features on Earth, record the magmatic response to changes in sea level. The models and data support a link between glacial cycles at the surface and mantle melting at depth, recorded in the bathymetric fabric of the sea floor.

preprint2014arXiv

Grain-size dynamics beneath mid-ocean ridges: Implications for permeability and melt extraction

Grain size is an important control on mantle viscosity and permeability, but is difficult or impossible to measure in situ. We construct a two-dimensional, single phase model for the steady-state mean grain size beneath a mid-ocean ridge. The mantle rheology is modelled as a composite of diffusion creep, dislocation creep, dislocation accommodated grain boundary sliding, and a plastic stress limiter. The mean grain size is calculated by the piezometric relationship of Austin and Evans [2007]. We investigate the sensitivity of our model to global variations in grain growth exponent, potential temperature, spreading-rate, and mantle hydration. We interpret the mean mean grain-size field in the context of permeability. The permeability structure due to mean grain size may be approximated as a high permeability region beneath a low permeability region. The transition between high and low permeability regions forms a boundary that is steeply sloped toward the ridge axis. We hypothesise that such a permeability structure generated from the variability of the mean grain size may be able to focus melt towards the ridge axis, analogous to a Sparks and Parmentier [1991]-type focusing. This focusing may, in turn, constrain the region where significant melt fractions are observed by seismic or magnetotelluric surveys. This interpretation of melt focusing via the grain-size permeability structure is consistent with MT observation of the asthenosphere beneath the East Pacific Rise [Baba et al., 2006, Key et al., 2013].

preprint2011arXiv

Formation of Chimneys in Mushy Layers: Experiment and Simulation

In this fluid dyanmics video, we show experimental images and simulations of chimney formation in mushy layers. A directional solidification apparatus was used to freeze 25 wt % aqueous ammonium chloride solutions at controlled rates in a narrow Hele-Shaw cell (1mm gap). The convective motion is imaged with schlieren. We demonstrate the ability to numerically simulate mushy layer growth for direct comparison with experiments.