Benchmark forward gravity schemes: the gravity field of a realistic lithosphere model WINTERC-G

Several alternative gravity forward modelling methodologies and associated numerical codes with their own advantages and limitations are available for the Solid Earth community. With the upcoming state-of-the-art lithosphere density models and accurate global gravity field data sets it is vital to understand the opportunities and limitations of the various approaches. In this paper, we discuss the four widely used techniques: global spherical harmonics (GSH), tesseroid integration (TESS), triangle integration (TRI), and hexahedral integration (HEX). A constant density shell benchmark shows that all four codes can produce similar precise gravitational potential fields. Two additional shell tests were conducted with more complicated density structures: lateral varying density structures and a Moho density interface between crust and mantle. The differences between the four codes were all below 1.5 percent of the modeled gravity signal suitable for reproducing satellite-acquired gravity data. TESS and GSH produced the most similar potential fields (< 0.3 percent). To examine the usability of the forward modelling codes for realistic geological structures, we use the global lithosphere model WINTERC-G, that was constrained, among other data, by satellite gravity field data computed using a spectral forward modeling approach. This spectral code was benchmarked against the GSH and it was confirmed that both approaches produce similar gravity solution with negligible differences between them. In the comparison of the different WINTERC-G-based gravity solutions, again GSH and TESS performed best. Only short-wavelength noise is present between the spectral and tesseroid forward modelling approaches, likely related to the different way in which the spherical harmonic analysis of the varying boundaries of the mass layer is performed. The Spherical harmonic basis functions produces small differences compared to the tesseroid elements especially at sharp interfaces, which introduces mostly short-wavelength differences. Nevertheless, both approaches (GSH and TESS) result in accurate solutions of the potential field with reasonable computational resources. Differences below 0.5 percent are obtained, resulting in residuals of 0.076 mGal standard deviation at 250 km height. The biggest issue for TRI is the characteristic pattern in the residuals that is related to the grid layout. Increasing the resolution and filtering allows for the removal of most of this erroneous pattern, but at the expense of higher computational loads with respect to the other codes. The other spatial forward modelling scheme HEX has more difficulty in reproducing similar gravity field solutions compared to GSH and TESS. These particular approaches need to go to higher resolutions, resulting in enormous computation efforts. The hexahedron-based code performs less than optimal in the forward modelling of the gravity signature, especially of a lateral varying density interface. Care must be taken with any forward modelling software as the approximation of the geometry of the WINTERC-G model may deteriorate the gravity field solution.

Metamorphism of metachert from the Southern Alps, New Zealand. A thermodynamic forward modelling study

<p>Scattered, scarce occurrences of garnet- and quartz-rich metamorphic rock, probably derived from Mn- and Fe-rich chert, occur within metamorphosed greywacke sequences worldwide. The metamorphism of such garnetiferous metacherts has not previously been investigated using modern thermodynamic forward modelling techniques due to the lack of appropriate, internally-consistent activity-composition (a–x) models for Mn-bearing minerals. The present study applies thermodynamic forward modelling using the recently-proposed a–x models of White et al. (2014) to investigate the metamorphism of garnetiferous metachert samples from the Southern Alps, New Zealand. Pressure-temperature (P–T) pseudosections are used in combination with results from petrography, element composition mapping using micro X-ray fluorescence (µXRF) and scanning electron microscope (SEM) methods, and garnet composition data from analytical transects by electron probe microanalysis (EPMA), to study metachert metamorphism. All the samples are compositionally layered, so the possibility exists that an input bulk rock composition might not match the effective bulk composition at the site of garnet growth. If a mineral assemblage stability field in a calculated P–T pseudosection matched the mineral assemblage in the rock, this was taken as an initial indication of a permissible input bulk rock composition. In that case, refined constraints on the P–T conditions were sought by comparing calculated and measured garnet compositions. The studied rocks include samples that are carbonate-bearing, which require consideration of the effects of fluid composition in mixed H₂O–CO₂ fluids, as well as a sample in which the garnet is strongly zoned, texturally-complex, and inferred to be of polymetamorphic origin. The effects of element fractionation by that garnet were investigated by recalculating the P–T pseudosection using a new bulk rock composition with the garnet core content removed. In none of the samples did the calculated and observed composition isopleths for the garnet cores match, suggesting that initial garnet nucleation in these Mn-rich rocks was locally controlled. For most samples in which the calculated and observed mineral assemblages matched, successful estimates of the peak metamorphic conditions were obtained. A garnet chert (A12E) from the mylonite zone of the Alpine Fault at Vine Creek, near Hokitika, gave a tight intersection of composition isopleths, indicating peak metamorphic conditions of 510 °C/5.5 kbar, after recalculation to correct for element fractionation by the strongly-zoned garnet. This tight, modern constraint is within error of previously-reported results from traditional geothermobarometry (420–600 °C/5.9–13 kbar) and Raman spectroscopy of carbonaceous material (RSCM T = 556 °C) from nearby sites. A peak metamorphic estimate of 520–550 °C/7–10 kbar was obtained from a dolomite-bearing sample from the garnet zone near Fox Glacier (J34), in good comparison with published temperatures from Raman spectroscopy of carbonaceous material in nearby metagreywacke samples (526–546 °C). The prograde metamorphic P–T path was probably steep, based on growth of the garnet core at ~475535 °C/5–9 kbar. The successful results for these garnet chert samples show that the new a-x models for Mn-bearing minerals extend the range of rock types that are amenable to pseudosection modelling. Results obtained in this study also serve to highlight several possible concerns: a) garnet nucleation and initial growth in very Mn-rich rocks may be subject to local compositional or kinetic controls; b) bulk rock compositions may not always mimic the effective bulk composition; c) the existing a–x models for Mn-bearing minerals and white micas may need refining; and d) some rocks may simply be ill-suited to thermodynamic forward modelling. Items a) and b) may be indicated by the common observation of a mismatch between predicted and measured garnet composition isopleths for garnet cores, and by a mismatch between garnet composition isopleths and the appropriate mineral assemblage field for sample AMS01, from the mylonite zone, Hari Hari, Southern Alps. For item c) every P–T pseudosection calculated using the new a–x models for Mn-bearing minerals showed garnet stable to very low temperatures below 300 °C. In addition, the P–T pseudosection for an oligoclase-zone metachporphyroblasts of Fe-Ti oxides (magnetitert (Sample J36) from Hari Mare stream, Franz Josef - Fox Glacier, indicated that the white mica margarite should be present instead of plagioclase (oligoclase), for a rock in which oligoclase is present and margarite is absent, a problem previously noted elsewhere. Item d) is exemplified by a very garnet-rich ferruginous metachert sample (J35, garnet zone, headwater region, Moeraki River, South Westland) which proved impossible to model successfully due to its complex mineral growth and deformation history. This sample contained multiple generations of carbonate with differing compositions, amphibole (not incorporated for modelling with the new a–x models for Mn-bearing minerals), large e associated with smaller, possibly later-formed ilmenite), and the garnet bands were offset by late deformation. The garnetiferous metachert samples studied here preserve in their textures and compositions clues to their growth mechanism and metamorphic history. The textures in at least two of the samples are consistent with the diffusion controlled nucleation and growth model for garnet. This research has successfully used state of the art thermodynamic modelling techniques in combination with the latest internally consistent a-x models on Mn-rich metachert, for the first time, extracting P–T conditions of the metamorphism of garnetiferous metachert from the Southern Alps.</p>

Metamorphism of metachert from the Southern Alps, New Zealand. A thermodynamic forward modelling study

<p>Scattered, scarce occurrences of garnet- and quartz-rich metamorphic rock, probably derived from Mn- and Fe-rich chert, occur within metamorphosed greywacke sequences worldwide. The metamorphism of such garnetiferous metacherts has not previously been investigated using modern thermodynamic forward modelling techniques due to the lack of appropriate, internally-consistent activity-composition (a–x) models for Mn-bearing minerals. The present study applies thermodynamic forward modelling using the recently-proposed a–x models of White et al. (2014) to investigate the metamorphism of garnetiferous metachert samples from the Southern Alps, New Zealand. Pressure-temperature (P–T) pseudosections are used in combination with results from petrography, element composition mapping using micro X-ray fluorescence (µXRF) and scanning electron microscope (SEM) methods, and garnet composition data from analytical transects by electron probe microanalysis (EPMA), to study metachert metamorphism. All the samples are compositionally layered, so the possibility exists that an input bulk rock composition might not match the effective bulk composition at the site of garnet growth. If a mineral assemblage stability field in a calculated P–T pseudosection matched the mineral assemblage in the rock, this was taken as an initial indication of a permissible input bulk rock composition. In that case, refined constraints on the P–T conditions were sought by comparing calculated and measured garnet compositions. The studied rocks include samples that are carbonate-bearing, which require consideration of the effects of fluid composition in mixed H₂O–CO₂ fluids, as well as a sample in which the garnet is strongly zoned, texturally-complex, and inferred to be of polymetamorphic origin. The effects of element fractionation by that garnet were investigated by recalculating the P–T pseudosection using a new bulk rock composition with the garnet core content removed. In none of the samples did the calculated and observed composition isopleths for the garnet cores match, suggesting that initial garnet nucleation in these Mn-rich rocks was locally controlled. For most samples in which the calculated and observed mineral assemblages matched, successful estimates of the peak metamorphic conditions were obtained. A garnet chert (A12E) from the mylonite zone of the Alpine Fault at Vine Creek, near Hokitika, gave a tight intersection of composition isopleths, indicating peak metamorphic conditions of 510 °C/5.5 kbar, after recalculation to correct for element fractionation by the strongly-zoned garnet. This tight, modern constraint is within error of previously-reported results from traditional geothermobarometry (420–600 °C/5.9–13 kbar) and Raman spectroscopy of carbonaceous material (RSCM T = 556 °C) from nearby sites. A peak metamorphic estimate of 520–550 °C/7–10 kbar was obtained from a dolomite-bearing sample from the garnet zone near Fox Glacier (J34), in good comparison with published temperatures from Raman spectroscopy of carbonaceous material in nearby metagreywacke samples (526–546 °C). The prograde metamorphic P–T path was probably steep, based on growth of the garnet core at ~475535 °C/5–9 kbar. The successful results for these garnet chert samples show that the new a-x models for Mn-bearing minerals extend the range of rock types that are amenable to pseudosection modelling. Results obtained in this study also serve to highlight several possible concerns: a) garnet nucleation and initial growth in very Mn-rich rocks may be subject to local compositional or kinetic controls; b) bulk rock compositions may not always mimic the effective bulk composition; c) the existing a–x models for Mn-bearing minerals and white micas may need refining; and d) some rocks may simply be ill-suited to thermodynamic forward modelling. Items a) and b) may be indicated by the common observation of a mismatch between predicted and measured garnet composition isopleths for garnet cores, and by a mismatch between garnet composition isopleths and the appropriate mineral assemblage field for sample AMS01, from the mylonite zone, Hari Hari, Southern Alps. For item c) every P–T pseudosection calculated using the new a–x models for Mn-bearing minerals showed garnet stable to very low temperatures below 300 °C. In addition, the P–T pseudosection for an oligoclase-zone metachporphyroblasts of Fe-Ti oxides (magnetitert (Sample J36) from Hari Mare stream, Franz Josef - Fox Glacier, indicated that the white mica margarite should be present instead of plagioclase (oligoclase), for a rock in which oligoclase is present and margarite is absent, a problem previously noted elsewhere. Item d) is exemplified by a very garnet-rich ferruginous metachert sample (J35, garnet zone, headwater region, Moeraki River, South Westland) which proved impossible to model successfully due to its complex mineral growth and deformation history. This sample contained multiple generations of carbonate with differing compositions, amphibole (not incorporated for modelling with the new a–x models for Mn-bearing minerals), large e associated with smaller, possibly later-formed ilmenite), and the garnet bands were offset by late deformation. The garnetiferous metachert samples studied here preserve in their textures and compositions clues to their growth mechanism and metamorphic history. The textures in at least two of the samples are consistent with the diffusion controlled nucleation and growth model for garnet. This research has successfully used state of the art thermodynamic modelling techniques in combination with the latest internally consistent a-x models on Mn-rich metachert, for the first time, extracting P–T conditions of the metamorphism of garnetiferous metachert from the Southern Alps.</p>

Comparison Study of Crustal Structure in Aceh Region based on Volcanic Arc System using Receiver Function Method

Aceh region has a very complex crustal structure from the forearc ridge to the backarc basin. This study aims to determine the velocity model of P and S waves and the depth of Moho discontinuity. This research was conducted using teleseismic earthquake data (30°-90° from the station) with M>6 from four seismic stations belonging to the BMKG in Aceh region. The stations are qualified based on the volcanic arc system zone. Furthermore, the velocity model determined by result of forward modelling, while the depth of the Moho layer estimated by migrated receiver function from time domain to the depth domain. At station SNSI that represented the forearc ridge zone, the depth of Moho is ±28 km, at station TPTI represent the forearc basin is ±16 km, while at zone with higher topography, namely volcanic arc zone represented by station KCSI, the Moho depth was identified at ±38 km, and the backarc basin represented by station LASI with ±40 km depth of Moho. This variation occurs because the composition of the earth’s layers below the station is diverse also different topography for each station.

Image processing to delineate the boundaries of peripheral arterial walls

The analysis of the arterial wall properties is vital in the prediction of stroke events and arterial hypertension in humans. Numerous researchers have experimented with several approaches to model arterial vessels and to analyse their biomechanical behaviour for many years now. Our study is focussed on image processing of peripheral arterial cross sections to detect and isolate the distinct layers. These boundaries will enable the creation of FEM models for further analysis of arterial wall properties. In a clinical setting, it facilitates doctors to identify the optimum pressure that can be applied to the artery for the treatment of stenosis without damaging the morphology of the blood vessels. This paper aims at distinguishing the various layers of arterial walls from images by minimizing human intervention. Cross section images of arteries from various sources were collected[10][11]. The boundaries from the image were obtained using image processing techniques of MATLAB(R2021a). The approach identified was to convert the input RGB images to grayscale, thresholding and applying morphological operators to delineate the Intima, Media, and Adventitia. These regions of interests (ROI) were then superimposed to generate an image with differentiated boundaries and void of unnecessary noise and inhomogeneity. This approach gave us an insight of the differences in various methods of boundary detection and to infer the optimum approach for accurate demarcation of boundaries of the three layers of arterial walls. It paves a pathway for forward modelling and to perform detailed FEM analysis in in-vitro diagnostics. In a nutshell, it was observed that the edge detection procedure implemented could be used for healthy and stenotic arteries. Further studies must be conducted to test the efficiency across a wide range of images and hence generalise its usage. Upon satisfactory boundary detection, forward modelling could be performed using the identified geometric forms.

1D Stochastic Inversion of Airborne Time-Domain Electromagnetic Data with Realistic Prior and Accounting for the Forward Modeling Error

Airborne electromagnetic surveys may consist of hundreds of thousands of soundings. In most cases, this makes 3D inversions unfeasible even when the subsurface is characterized by a high level of heterogeneity. Instead, approaches based on 1D forwards are routinely used because of their computational efficiency. However, it is relatively easy to fit 3D responses with 1D forward modelling and retrieve apparently well-resolved conductivity models. However, those detailed features may simply be caused by fitting the modelling error connected to the approximate forward. In addition, it is, in practice, difficult to identify this kind of artifacts as the modeling error is correlated. The present study demonstrates how to assess the modelling error introduced by the 1D approximation and how to include this additional piece of information into a probabilistic inversion. Not surprisingly, it turns out that this simple modification provides not only much better reconstructions of the targets but, maybe, more importantly, guarantees a correct estimation of the corresponding reliability.