An implicit free surface algorithm for geodynamical simulations

Identifying the dominant controls on Earth’s surface topography is of critical importance to understanding both the short- and long-term evolution of geological processes and past- and present-day dynamics of Earth’s coupled mantle–lithosphere system. The ability to simulate a stress free — or a so-called ‘free surface’ — boundary condition is required to examine such processes via numerical models. However, at present, geodynamical models incorporating a free surface are limited, as most underlying free surface algorithms place severe restrictions on the computational timestep. Consequently, the simulations are often intractable. In this study, we introduce a new approach for incorporating a free surface within geodynamical models: an algorithm, in which free surface elevation is treated as an independent variable and is solved for in conjunction with the momentum and continuity equation, using implicit time integration. We demonstrate that the method is straightforward to implement in existing models and, using a series of analytical and benchmark comparisons, we show that it does not suffer from the timestep constraints of previous schemes. Furthermore, the scheme can be made second order accurate in time, at no additional cost. The method therefore dramatically improves the computational efficiency of geodynamical simulations including a free surface, whilst maintaining solution accuracy

Pressure-Driven Poiseuille Flow: A Major Component of the Torque-Balance Governing Pacific Plate Motion

The Pacific Plate is thought to be driven mainly by slab pull, associated with subduction along the Aleutians-Japan, Marianas-Izu-Bonin, and Tonga-Kermadec trenches. This implies that viscous flow within the sub-Pacific asthenosphere is mainly generated by overlying plate motion (i.e., Couette flow) and that the associated shear stresses at the lithosphere's base are resisting such motion. Recent studies on glacial isostatic adjustment and lithosphere dynamics provide tighter constraints on the viscosity and thickness of Earth's asthenosphere and, therefore, on the amount of shear stress that asthenosphere and lithosphere mutually exchange, by virtue of Newton's third law of motion. In light of these constraints, the notion that subduction is the main driver of present-day Pacific Plate motion becomes somewhat unviable, as the pulling force that would be required by slabs exceeds the maximum available from their negative buoyancy. Here we use coupled global models of mantle and lithosphere dynamics to show that the sub-Pacific asthenosphere features a significant component of pressure-driven (i.e., Poiseuille) flow and that this has driven at least 50% of the Pacific Plate motion since, at least, 15 Ma. A corollary of our models is that a sublithospheric pressure difference as high as ±50 MPa is required across the Pacific domain.I. L. S. acknowledges support from CONICYT Becas-Chile scholarship and from IGN at the University of Copenhagen. D. R. D. acknowledges funding from the Australian Research Council, via grants FT140101262 and DP170100058

Effects of basal drag on subduction dynamics from 2D numerical models

Subducting slabs are an important driver of plate motions, yet the relative importance of different forces in governing subduction motions and styles remains incompletely understood. Basal drag has been proposed to be a minor contributor to subduction forcing because of the lack of correlation between plate size and velocity in observed and reconstructed plate motions. Furthermore, in single subduction system models, low basal drag leads to subduction behaviour most consistent with the observation that trench migration velocities are generally low compared to convergence velocities. By contrast, analytical calculations and global mantle flow models indicate basal drag can be substantial. In this study, we revisit this problem by examining the drag at the base of the lithosphere, for a single subduction system, in 2D models with a free trench and composite non-linear rheology. We compare the behaviour of short and long plates for a range of asthenospheric and lithospheric rheologies. We reproduce results from previous modelling studies, including low ratios of trench over plate motions. However, we also find that any combination of asthenosphere and lithosphere viscosity that produces Earth-like subduction behaviour leads to a correlation of velocities with plate size, due to the role of basal drag. By examining Cenozoic plate motion reconstructions, we find that slab age and plate size are positively correlated: higher slab pull for older plates tends to be offset by higher basal drag below these larger plates. This, in part, explains the lack of plate velocity-size correlation in observations, despite the important role of basal drag in the subduction force balance.Lior Suchoy was supported by the Engineering and Physical Sciences Research Council (EPSRC) (grant no. EP/N509486/1). Ben Maunder and Saskia Goes were supported by the Natural Environment Research Council (NERC) (grant no. NE/K010743/1). Rhodri Davies was supported by the Australian Research Council (grant no. DP170100058)

Global patterns in Earth's dynamic topography since the Jurassic: the role of subducted slabs

We evaluate the spatial and temporal evolution of Earth's long-wavelength surface dynamic topography since the Jurassic using a series of high-resolution global mantle convection models. These models are Earth-like in terms of convective vigour, thermal structure, surface heat-flux and the geographic distribution of heterogeneity. The models generate a degree-2-dominated spectrum of dynamic topography with negative amplitudes above subducted slabs (i.e. circum-Pacific regions and southern Eurasia) and positive amplitudes elsewhere (i.e. Africa, north-western Eurasia and the central Pacific). Model predictions are compared with published observations and subsidence patterns from well data, both globally and for the Australian and southern African regions. We find that our models reproduce the long-wavelength component of these observations, although observed smaller-scale variations are not reproduced. We subsequently define geodynamic rules for how different surface tectonic settings are affected by mantle processes: (i) locations in the vicinity of a subduction zone show large negative dynamic topography amplitudes; (ii) regions far away from convergent margins feature long-term positive dynamic topography; and (iii) rapid variations in dynamic support occur along the margins of overriding plates (e.g. the western US) and at points located on a plate that rapidly approaches a subduction zone (e.g. India and the Arabia Peninsula). Our models provide a predictive quantitative framework linking mantle convection with plate tectonics and sedimentary basin evolution, thus improving our understanding of how subduction and mantle convection affect the spatio-temporal evolution of basin architecture.This research was supported by resources provided by the Pawsey Supercomputing Centre with funding from the Australian Government and the Government of Western Australia and with the assistance of resources from the National Computational Infrastructure (NCI), which is supported by the Australian Government. Sascha Brune was funded by the Marie Curie International Outgoing Fellowship 326115 and the Helmholtz Young Investigators Group CRYSTALS. Christian Heine was supported by ARC Linkage Project LP0989312 with Shell E & P and TOTAL. D. Rhodri Davies is funded by an ARC Future Fellowship (FT140101262) and Simon Williams and R. Dietmar Müller are supported by ARC grants DP130101946 and IH13020001

Automated Inline Analysis of Myocardial Perfusion MRI with Deep Learning

Recent development of quantitative myocardial blood flow (MBF) mapping allows direct evaluation of absolute myocardial perfusion, by computing pixel-wise flow maps. Clinical studies suggest quantitative evaluation would be more desirable for objectivity and efficiency. Objective assessment can be further facilitated by segmenting the myocardium and automatically generating reports following the AHA model. This will free user interaction for analysis and lead to a 'one-click' solution to improve workflow. This paper proposes a deep neural network based computational workflow for inline myocardial perfusion analysis. Adenosine stress and rest perfusion scans were acquired from three hospitals. Training set included N=1,825 perfusion series from 1,034 patients. Independent test set included 200 scans from 105 patients. Data were consecutively acquired at each site. A convolution neural net (CNN) model was trained to provide segmentation for LV cavity, myocardium and right ventricular by processing incoming 2D+T perfusion Gd series. Model outputs were compared to manual ground-truth for accuracy of segmentation and flow measures derived on global and per-sector basis. The trained models were integrated onto MR scanners for effective inference. Segmentation accuracy and myocardial flow measures were compared between CNN models and manual ground-truth. The mean Dice ratio of CNN derived myocardium was 0.93 +/- 0.04. Both global flow and per-sector values showed no significant difference, compared to manual results. The AHA 16 segment model was automatically generated and reported on the MR scanner. As a result, the fully automated analysis of perfusion flow mapping was achieved. This solution was integrated on the MR scanner, enabling 'one-click' analysis and reporting of myocardial blood flow.Comment: This work has been submitted to Radiology: Artificial Intelligence for possible publicatio

Compressibility of ferropericlase at high-temperature: evidence for the iron spin crossover in seismic tomography

The iron spin crossover in ferropericlase, the second most abundant mineral in Earth's lower mantle, causes changes in a range of physical properties, including seismic wave velocities. Understanding the effect of temperature on the spin crossover is essential to detect its signature in seismic observations and constrain its occurrence in the mantle. Here, we report the first experimental results on the spin crossover-induced bulk modulus softening at high temperatures, derived directly from time-resolved x-ray diffraction measurements during continuous compression of (Mg0.8Fe0.2)O in a resistive-heated dynamic diamond-anvil cell. We present new theoretical calculations of the spin crossover at mantle temperatures benchmarked by the experiments. Based on our results, we create synthetic seismic tomography models to investigate the signature of the spin crossover in global seismic tomography. A tomographic filter is applied to allow for meaningful comparisons between the synthetic models and data-based seismic tomography models, like SP12RTS. A negative anomaly in the correlation between Vs variations and Vc variations (S-C correlation) is found to be the most suitable measure to detect the presence of the spin crossover in tomographic models. When including the effects of the spin crossover, the misfit between the synthetic model and SP12RTS is reduced by 63%, providing strong evidence for the presence of the spin crossover, and hence ferropericlase, in the lower mantle. Future improvement of seismic resolution may facilitate a detailed mapping of spin state using the S-C correlation, providing constraints on mantle temperatures by taking advantage of the temperature sensitivity of the spin crossover

The effect of sepsis and its inflammatory response on mechanical clot characteristics: a prospective observational study

Purpose: Sepsis and its progression are known to have a major influence on the coagulation system. Current coagulation tests are of limited use when assessing coagulation in sepsis patients. This study aims to assess the potential for a new functional biomarker of clot microstructure, fractal dimension, to identify changes in the mechanical properties of clot microstructure across the sepsis spectrum (sepsis, severe sepsis and septic shock). Methods: A total of 100 patients that presented acutely to a large teaching hospital were included in this prospective observational study (50 sepsis, 20 severe sepsis and 30 septic shock) against a matched control of 44 healthy volunteers. Fractal analysis was performed, as well as standard markers of coagulation, and six plasma markers of inflammation. Results: Fractal dimension was significantly higher in the sepsis and severe sepsis groups than the healthy control (1.78 ± 0.07 and 1.80 ± 0.05 respectively vs 1.74 ± 0.03) (p < 0.001), indicating a significant increase in mechanical clot strength and elasticity consistent with a hypercoagulable state. Conversely, fractal dimension was significantly lower in septic shock (1.66 ± 0.10, p < 0.001), indicating a significant reduction in mechanical clot strength and functionality consistent with a hypocoagulable state. This corresponded with a significant increase in the inflammatory response. Conclusions: This study confirms that clot microstructure is significantly altered through the various stages of sepsis. Of particular importance was the marked change in clot development between severe sepsis and septic shock, which has not been previously reported

Technical development and feasibility of a reusable vest to integrate cardiovascular magnetic resonance with electrocardiographic imaging

BACKGROUND: Electrocardiographic imaging (ECGI) generates electrophysiological (EP) biomarkers while cardiovascular magnetic resonance (CMR) imaging provides data about myocardial structure, function and tissue substrate. Combining this information in one examination is desirable but requires an affordable, reusable, and high-throughput solution. We therefore developed the CMR-ECGI vest and carried out this technical development study to assess its feasibility and repeatability in vivo. METHODS: CMR was prospectively performed at 3T on participants after collecting surface potentials using the locally designed and fabricated 256-lead ECGI vest. Epicardial maps were reconstructed to generate local EP parameters such as activation time (AT), repolarization time (RT) and activation recovery intervals (ARI). 20 intra- and inter-observer and 8 scan re-scan repeatability tests. RESULTS: 77 participants were recruited: 27 young healthy volunteers (HV, 38.9 ± 8.5 years, 35% male) and 50 older persons (77.0 ± 0.1 years, 52% male). CMR-ECGI was achieved in all participants using the same reusable, washable vest without complications. Intra- and inter-observer variability was low (correlation coefficients [rs] across unipolar electrograms = 0.99 and 0.98 respectively) and scan re-scan repeatability was high (rs between 0.81 and 0.93). Compared to young HV, older persons had significantly longer RT (296.8 vs 289.3 ms, p = 0.002), ARI (249.8 vs 235.1 ms, p = 0.002) and local gradients of AT, RT and ARI (0.40 vs 0.34 ms/mm, p = 0,01; 0.92 vs 0.77 ms/mm, p = 0.03; and 1.12 vs 0.92 ms/mm, p = 0.01 respectively). CONCLUSION: Our high-throughput CMR-ECGI solution is feasible and shows good reproducibility in younger and older participants. This new technology is now scalable for high throughput research to provide novel insights into arrhythmogenesis and potentially pave the way for more personalised risk stratification. CLINICAL TRIAL REGISTRATION: Title: Multimorbidity Life-Course Approach to Myocardial Health-A Cardiac Sub-Study of the MRC National Survey of Health and Development (NSHD) (MyoFit46). National Clinical Trials (NCT) number: NCT05455125. URL: https://clinicaltrials.gov/ct2/show/NCT05455125?term=MyoFit&draw=2&rank=1