Search

Gravimetry is a discipline of geophysics that deals with observation and interpretation of the earth gravity field. The acquired gravity data serve the study of the earth interior, be it the deep or the near surface one, by means of the inferred subsurface structural density distribution. The subsurface density structure is resolved by solving the gravimetric inverse problem. Diverse methods and approaches exist for solving this non-unique and ill-posed inverse problem. Here, we focused on those methods that do not pre-constrain the number or geometries of the density sources. We reviewed the historical development and the basic principles of the Growth inversion methodology, which belong to the methods based on the growth of the model density structure throughout an iterative exploration process. The process was based on testing and filling the cells of a subsurface domain partition with density contrasts through an iterative mixed weighted adjustment procedure. The procedure iteratively minimized the data misfit residuals jointly with minimizing the total anomalous mass of the model, which facilitated obtaining compact meaningful source bodies of the solution. The applicability of the Growth inversion approach in structural geophysical studies, in geodynamic studies, and in near surface gravimetric studies was reviewed and illustrated. This work also presented the first application of the Growth inversion tool to near surface microgravimetric data with the goal of seeking very shallow cavities in archeological prospection and environmental geophysics.

Gravimetric observations were carried out in 2015, to image the uppermost portion of the volcanic plumbing system of Mt. Etna (Italy). Gravity measurements were performed using two relative gravimeters, along a profile that crosses the summit craters area (elevations between 2,820 and 3,280 m a.s.l.). Accurate positioning of the gravity observation points was determined through GPS measurements. After applying elevation and terrain corrections, the reduced gravity data were used to build a 2D density model of the uppermost part of the volcano edifice. This model was constrained using to-date knowledge of the structural setting of the area and the available volcanological data. We highlighted the presence of low-density material below the summit craters, down to the depth of about 2.1 km, interpreted as highly altered, fumarolized and structurally weakened material. It is also likely that the close presence of the conduits feeding the summit craters of the volcano contributes to the gravity low in the SW half of the measurement profile. Conversely, the gravity low observed at the northern edge of the profile could reflect the high concentration of faults and eruptive fissures in the Pizzi Deneri area, in correspondence of the Ellittico caldera rim.

Introduction: We present the first case of appendiceal intussusception associated with myeloid sarcoma in a young patient. Minimally invasive techniques used along the clinical course are highlighted. Case description: A 2.5-year-old boy was admitted after three weeks of COVID-19 infection with ongoing symptoms of MIS-C. Due to constipation, distended belly and vomiting, US was done which showed ileocolic intussusception. After unsuccessful hydrostatic reduction laparoscopic exploration was performed, where the vermiform appendix was found to be thickened and partially intussuscepted into the coecum. The ileocecal region was exteriorized transumbilically. After manual reduction of the intussusception, a long, thickened, fragile appendix was removed. Histopathology revealed myeloid sarcoma. Bone marrow investigation identified acute myeloid leukemia. During the oncological treatment, laparoscopic cholecystectomy was necessary due to cholecystitis and cholelithiasis. The child recovered uneventfully in terms of surgical complications, with good cosmetic result. Conclusion: No similar case in childhood was found in the English literature. Unusual symptoms and radiological findings of intussusception can conceal unexpected disorders. Minimally invasive technique offered advantages in the treatment of the presented patient and can be recommended to treat intussusception or cholelithiasis, if applicable, during an ongoing oncological treatment as well.

Abstract Globally there is abundant terrestrial surface gravity data used to study the time variation of gravity related to subsurface mass and density changes in different geological, geodynamical and geotechnical environments. We present here a tool for analysing existing and newly acquired, 4D gravity data, which creates new findings from its reuse. Our method calculates in an almost automatic way the possible sources of density change responsible for the observed gravity variations. The specifics of the new methodology are: use of a low number of observation points, relatively small source structures, low signal/noise ratio in the data, and a free 3D source geometry without initial hypothesis. The process is based on the non-linear adjustment of structures defined by aggregation of small cells corresponding to a 3D section of the sub-floor volume. This methodology is implemented in a software tool, named GROWTH-dg, which can be freely downloaded for immediate use, together with a user manual and application examples.

We analyse spatiotemporal gravity changes observed on the Ischia island (Italy) accompanying the destructive earthquake of 21 August 2017. The 29 May 2016 to 22 September 2017 time-lapse gravity changes observed at 18 benchmarks of the Ischia gravimetric network are first corrected for the gravitational effect of the surface deformation using the deformation-induced topographic effect (DITE) correction. The co-seismic DITE is computed by Newtonian volumetric integration using the Toposk software, a high-resolution LiDAR DEM and the co-seismic vertical displacement field derived from Sentinel-1 InSAR data. We compare numerically the DITE field with its commonly used Bouguer approximation over the island of Ischia with the outcome that the Bouguer approximation of DITE is adequate and accurate in this case. The residual gravity changes are then computed at gravity benchmarks by correcting the observed gravity changes for the planar Bouguer effect of the elevation changes at benchmarks over the same period. The residual gravity changes are then inverted using an inversion approach based on model exploration and growing source bodies, making use of the Growth-dg inversion tool. The found inversion model, given as subsurface time-lapse density changes, is then interpreted as mainly due to a co-seismic or post-seismic disturbance of the hydrothermal system of the island. Pros and weak points of such interpretation are discussed.

Gravitational effect of surface deformation is in 4D microgravimetry treated as the deformation-induced topographic effect (DITE). The DITE field is computed using Newtonian volumetric integration which requires high resolution digital elevation model (DEM) and vertical displacement field in areal form. If only elevation changes on benchmarks of the gravimetric network are available, instead of the vertical displacement field, the DITE on benchmarks can be evaluated only approximately, using a planar Bouguer or a normal free-air-effect (nFAE) approximation. Here we analyze the adequacy and accuracy of these two approximations in a case study for the December 2018 fissure eruption on Etna accompanied by significant surface deformation caused primarily by a relatively shallow dyke. The outcome is that in volcanic areas of similar morphology as that over the Etna summit area, and for surface deformation fields due to relatively shallow dykes, neither the Bouguer nor the nFAE approximation of the DITE is accurate enough. In such situations the residual gravity changes should be computed with both the Bouguer and nFAE corrections and interpreted as two marginal cases. In addition we analyze also a correction for the effect of benchmark elevation change based on the topographically modelled (predicted) vertical gradient of gravity (VGG) meant to approximate the in-situ VGG values at benchmarks. This correction does not appear suitable to approximate the DITE in conditions of our case study or in broader sense.

Results from a detailed gravity survey realized along the planned highway tunnel in the karstic area of Slovak Karst in the eastern Slovakia are presented. Detailed gravity profiles crossed an area of rugged topography, therefore the terrain corrections played a crucial role in the gravity data processing. The airborne laser scanning technique (LiDAR) was used in order to compile a high-resolution digital terrain model (DTM) of the surrounding area and to calculate terrain corrections properly. The difference between the Bouguer anomalies calculated with an available nationwide DTM and those with new LiDAR-based model can be significant in some places as it is presented in the paper. A new method for Bouguer correction density analysis based on surface data is presented. Special underground gravity measurements in the existing nearby railway tunnel were also conducted in order to determine the mean density of the topographic rocks. The Bouguer anomalies were used to interpret lithological contacts and tectonic/karstic discontinuities.

Predicted values of the vertical gradient of gravity (VGG) on benchmarks of Etna’s monitoring system, based on calculation of the topographic contribution to the theoretical free-air gradient, are compared with VGG values observed in situ. The verification campaign indicated that improvements are required when predicting the VGGs at such networks. Our work identified the following factors to be resolved: (a) accuracy of the benchmark position; (b) gravitational effect of buildings and roadside walls adjacent to benchmarks; (c) accuracy of the digital elevation model (DEM) in the proximity of benchmarks. Benchmark positions were refined using precise geodetic methods. The gravitational effects of the benchmark-adjacent walls and buildings were modeled and accounted for in the prediction. New high-resolution DEMs were produced in the innermost zone at some benchmarks based on drone-flown photogrammetry to improve the VGG prediction at those benchmarks. The three described refinements in the VGG prediction improved the match between predicted and in situ observed VGGs at the network considerably. The standard deviation of differences between the measured and predicted VGG values decreased from 36 to 13 μGal/m.

We report on an 11-year-old girl with cystic fibrosis who presented with thoracic pain and an extensive subcutaneous emphysema and subsequently developed progressive respiratory distress. The chest computed tomography revealed a huge pneumomediastinum. Due to the development of severe respiratory failure, urgent needle thoracocentesis was necessary that resulted in only temporary improvement. Therefore, under general anesthesia two mediastinal drains were introduced. Using active suction, the size of the pneumomediastinum decreased gradually and the drains were removed after 3 weeks. Here, we describe an extremely rare situation, when acute surgical intervention was necessary in a child with spontaneous pneumomediastinum.

Urothelial neoplasms of the bladder (UNB) are considerably rare throughout the pediatric population. UNB develops from the urothelial tissue in the form of a benign disease, generally favoring a successful prognosis in the majority of cases. The authors present the diagnosis and treatment regarding two medical case reports in which urothelial papilloma was diagnosed and effectively treated. Case 1: A 15-year-old male patient was presented to our clinic complaining of a painless yet distinctive, macroscopic form of hematuria. Following a routine examination, which included ultrasound (US) and intravenous pyelography, the urethrocystoscopy revealed an intravesical solitary lesion positioned in the vicinity of the left ureteral orifice. Additionally, histology confirmed urothelial papilloma. During the follow-up, laboratory, urinary control tests, and US results all proved negative. Case 2: A 13-year-old male patient was admitted to our clinic and examined, in regard to complaints associated with recurrent abdominal pain. The pathology was discovered incidentally on abdominal US. Preoperative US and magnetic resonance imaging (MRI) studies ensued, resulting in a scheduled MRI, followed by urethrocystoscopy, which confirmed an intravesical solitary lesion positioned near the right ureteral orifice. Histology revealed urothelial papilloma. During the follow-up control cystoscopy, one resection was repeated due to the presence of a residual tumor. Today, 10 years since the presence of uroepithelial papilloma, both patients are asymptomatic and tumor-free. If there is likely suspicion of recurrence, cystoscopy is recommended.

In pediatric burns the use of systemic antibiotic prophylaxis is a standard procedure in some burn centers, though its beneficial effect on the infectious complications is debated. The present meta-analysis aimed at determining whether systemic antibiotic prophylaxis prevents infectious complications in pediatric patients with burn injuries. We searched the PubMed, EMBASE, and Cochrane Library databases from inception to August 2019. We included 6 studies, in which event rates of infectious complications were reported in children with burn injuries receiving or not receiving systemic antibiotic prophylaxis. We found that the overall odds ratio (OR) of developing an infection (including local and systemic) was not different between the groups (OR = 1.35; 95% CI, 0.44, 4.18). The chances for systemic infectious complications alone were also not different between antibiotic-treated and non-treated patients (OR = 0.74; 95% CI, 0.38, 1.45). Based on the age, affected total body surface area, and country income level, we did not find any subgroup that benefited from the prophylaxis. Our findings provide quantitative evidence for the inefficacy of systemic antibiotic prophylaxis in preventing infections in pediatric burns. To validate our conclusion, multinational, randomized trials in a diverse population of children with burn injuries are warranted.

A full-term male neonate presented with a left sided cervical lump at the level of the thyroid gland. Magnetic resonance imaging (MRI) showed a benign heterogeneous solid mass with lobulated margins. The tumor underwent complete excision. Histology revealed the diagnosis of chondromesenchymal hamartoma in ectopic thyroid tissue.

We present a detailed calculation of the topographic contribution to the vertical gradient of gravity (VGG) based on high-resolution digital elevation model (DEM) and new developed software (Toposk) for the purpose of predicting the actual VGGs in the field. The calculations presented here were performed for the Central Volcanic Complex (CVC) of Tenerife. We aimed at identifying the most extreme VGGs within the CVC, as well as predicting the VGGs at benchmarks of the former microgravity/deformation network set up to monitor the 2004/5 unrest. We have carried out an observational campaign in June 2016 to verify the predicted VGG values, both the extreme ones and those at the benchmarks. The comparison between the predicted and the in-situ verified VGGs is presented here. We demonstrate the sensitivity of the VGG prediction to the choice of the topo-rock density, which is inherent to the volcanic areas with high variability of rock densities. We illustrate the significance of the use of actual VGG in volcano monitoring microgravimetric surveys on a couple of benchmarks of the CVC network.

We test here the feasibility of ground-penetrating radar (GPR) and microgravity methods in identifying underground voids, such as cellars, tunnels, abandoned mine-workings, etc., in complex urban conditions. For this purpose, we selected a cellar located under a private lot in a residential quarter of the town of Senec in Western Slovakia, which was discovered by chance when a small sinkhole developed on the yard just two meters away from the house. The size of our survey area was limited 1) by the presence of a technical room built at the back of the yard with a staircase leading to the garden, and 2) by the small width of the lot. Therefore the geophysical survey was carried out only in the backyard of the lot as we were not permitted to measure on neighbouring estates. The results from the GPR measurements obtained by the GSSI SIR-3000 system with 400 MHz antenna were visualized in the form of 2D radargrams with the corresponding transformed velocity model of studied cross-sections. Only the profiles running over the pavement next to the house yielded interpretable data because the local geological situation and the regular watering of the lawn covering prevailingly the backyard caused significant attenuation of the emitted GPR signal. The Bouguer gravity map is dominated by a distinctive negative anomaly indicating the presence of a shallow underground void. The quantitative interpretation by means of Euler deconvolution was utilized to validate the depth of the center and location of the cellar. Comparison with the gravitational effect of the cellar model calculated in the in-house program Polygrav shows a quite good correlation between the modelled and observed fields. Only a part of the aerial extent of the anomaly could be traced by the used geophysical methods due to accessibility issues.

We review here the gravitational effects on the temporal (time-lapse) gravity changes induced by the surface deformation (vertical displacements). We focus on two terms, one induced by the displacement of the benchmark (gravity station) in the ambient gravity field, and the other imposed by the attraction of the masses within the topographic deformation rind. The first term, coined often the Free Air Effect (FAE), is the product of the vertical gradient of gravity (VGG) and the vertical displacement of the benchmark. We examine the use of the vertical gradient of normal gravity, typically called the theoretical or normal Free Air Gradient (normal FAG), as a replacement for the true VGG in the FAE, as well as the contribution of the topography to the VGG. We compute a topographic correction to the normal FAG, to offer a better approximation of the VGG, and evaluate its size and shape (spatial behavior) for a volcanic study area selected as the Central Volcanic Complex (CVC) on Tenerife, where this correction reaches 77% of the normal FAG and varies rapidly with terrain. The second term, imposed by the attraction of the vertically displaced topo-masses, referred to here as the Topographic Deformation Effect (TDE) must be computed by numerical evaluation of the Newton volumetric integral. As the effect wanes off quickly with distance, a high resolution DEM is required for its evaluation. In practice this effect is often approximated by the planar or spherical Bouguer deformation effect (BDE). By a synthetic simulation at the CVC of Tenerife we show the difference between the rigorously evaluated TDE and its approximation by the planar BDE. The complete effect, coined here the Deformation Induced Topographic Effect (DITE) is the sum of FAE and TDE.

Here we investigate the applicability of the harmonic inversion method to time-lapse gravity changes observed in volcanic areas. We carry out our study on gravity changes occured over the period of 2004–2005 during the unrest of the Central Volcanic Complex on Tenerife, Canary Islands. The harmonic inversion method is unique in that it calculates the solution of the form of compact homogeneous source bodies via the mediating 3-harmonic function called quasigravitation. The latter is defined in the whole subsurface domain and it is a linear integral transformation of the surface gravity field. At the beginning the seeds of the future source bodies are introduced: these are quasispherical bodies located at the extrema of the quasigravitation (calculated from the input gravity data) and their differential densities are free parameters preselected by the interpreter. In the following automatic iterative process the source bodies change their size and shape according to the local values of quasigravitation (calculated in each iterative step from the residual surface gravity field); the process stops when the residual surface gravity field is sufficiently small. In the case of inverting temporal gravity changes, the source bodies represent the volumetric domains of temporal mass-density changes. The focus of the presented work is to investigate the dependence of the size and shape of the found source bodies on their differential densities. We do not aim here (yet) at interpreting the found solutions in terms of volcanic processes associated with intruding or rejuvenating magma and/or migrating volatiles.

The empirical model of the depth-dependent density change within the upper continental mantle is derived in this study. The density of the upper(most) mantle underlying the continental crust is obtained from the estimated values of the crust–mantle (Moho) density contrast. Since the continental crustal thickness varies significantly, these upper mantle density values to a large extent reflect the density changes with depth. The estimation of the Moho density contrast is done through solving Moritz’s generalization of the Vening-Meinesz inverse problem of isostasy. The solution combines gravity and seismic data in the least-squares estimation model. The estimated upper mantle density (beneath the continental crust) varies between 2770 and 3649 kg/m3. The upper mantle density increases almost proportionally with depth at a rate of 13±2 kg/m3 per 1 km at the investigated depth interval from 6 to 58 km.

The robustness of the gravimetric forward modelling is investigated by applying the harmonic inversion procedure at the input which is the difference of the calculated and measured surface gravity. The gravity data are taken from two profiles in the Carpathian-Pannonian Basin region. The result of the inversion are density models obtained from the original two-layer models with various horizontal boundary depth and density contrast. The deformation of the originally planar boundary is the measure of the mismatch between calculated and measured data. The calculated deformation has reached up to tens of kilometers and thus the uncertainties in determining the geometry of disturbing bodies by the forward modelling are substantial.

We introduce the generic expressions for computing the gravitational field (potential and its radial derivative) generated by an arbitrary density (contrast) layer with a variable depth and thickness having a laterally-distributed radial density variation. The information on the geometry and density distribution of a volumetric mass layer is described by means of spherical harmonics. These generic expressions can uniformly be applied to model all major known density structures within the Earth’s interior using methods for a spherical harmonic analysis and synthesis of the gravitational field. This is demonstrated on specific examples given for various density models commonly adopted for the approximation of crust density structures.

We derive the expressions for computing the ice density contrast stripping corrections to the topography corrected gravity field quantities by means of the spherical harmonics. The expressions in the spectral representation utilize two types of the spherical functions, namely the spherical height functions and the newly introduced lower-bound ice functions. The spherical height functions describe the global geometry of the upper topographic bound. The spherical lower-bound ice functions combined with the spherical height functions describe the global thickness of the continental ice sheet. The newly derived formulas are utilized in the forward modelling of the gravitational field quantities generated by the ice density contrast. The 30×30 arc-sec global elevation data from GTOPO30 are used to generate the global elevation model (GEM) coefficients. The spatially averaged global elevation data from GTOPO30 and the 2×2 arc-deg ice-thickness data from the CRUST 2.0 global crustal model are used to generate the global lower-bound ice model (GIM) coefficients. The mean value of the ice density contrast 1753 kg/m3 (i.e., difference of the reference constant density of the continental upper crust 2670 kg/m3and the density of glacial ice 917 kg/m3) is adopted. The numerical examples are given for the gravitational potential and attraction generated by the ice density contrast computed globally with a low-degree spectral resolution complete to degree and order 90 of the GEM and GIM coefficients.

The currently available global geopotential models and the global elevation and bathymetry data allow modelling the topography-corrected and bathymetry stripped reference gravity field to a very high spectral resolution (up to degree 2160 of spherical harmonics) using methods for a spherical harmonic analysis and synthesis of the gravity field. When modelling the topography-corrected and crust-density-contrast stripped reference gravity field, additional stripping corrections are applied due to the ice, sediment and other major known density contrasts within the Earth’s crust. The currently available data of global crustal density structures have, however, a very low resolution and accuracy. The compilation of the global crust density contrast stripped gravity field is thus limited to a low spectral resolution, typically up to degree 180 of spherical harmonics. In this study we derive the expressions used in forward modelling of the bathymetry-generated gravitational field quantities and the corresponding bathymetric stripping corrections to gravity field quantities by means of the spherical bathymetric (ocean bottom depth) functions. The expressions for the potential and its radial derivative are formulated for the adopted constant (average) ocean saltwater density contrast and for the spherical approximation of the geoid surface. These newly derived expressions are utilized in numerical examples to compute the gravitational potential and attraction generated by the ocean density contrast. The computation is realized globally on a 1 × 1 arc-deg geographical grid at the Earth’s surface.

We compile the global maps of atmospheric effects on the gravity field quantities using the spherical harmonic representation of the gravitational field. A simple atmospheric density distribution is assumed within a lower atmosphere (< 6 km). Dis- regarding temporal and lateral atmospheric density variations, the radial atmospheric density model is defined as a function of the nominal atmospheric density at the sea level and the height. For elevations above 6 km, the atmospheric density distribution from the United States Standard Atmosphere 1976 is adopted. The 5 × 5 arc-min global elevation data from the ETOPO5 are used to generate the global elevation model coefficients. These coefficients (which represent the geometry of the lower bound of atmospheric masses) are utilized to compute the atmospheric effects with a spectral resolution complete to degree and order 180. The atmospheric effects on gravity disturbances, gravity anomalies and geoid undulations are evaluated globally on a 1 × 1 arc-deg grid.

We investigate globally the correlation of the step-wise consolidated crust- stripped gravity field quantities with the topography, bathymetry, and the Moho bound- ary. Global correlations are quantified in terms of Pearson’s correlation coefficient. The elevation and bathymetry data from the ETOPO5 are used to estimate the correlation of the gravity field quantities with the topography and bathymetry. The 2×2 arc-deg discrete data of the Moho depth from the global crustal model CRUST 2.0 are used to estimate the correlation of the gravity field quantities with the Moho boundary. The results reveal that the topographically corrected gravity field quantities have the highest absolute correlation with the topography. The negative correlation of the topographically corrected gravity disturbances with the topography over the continents reaches –0.97. The ocean, ice and sediment density contrasts stripped and topographically corrected gravity field quantities have the highest correlation with the bathymetry (ocean bottom relief ). The correlation of the ocean, ice and sediment density contrasts stripped and topographically corrected gravity disturbances over the oceans reaches 0.93. The consolidated crust-stripped grav- ity field quantities have the highest absolute correlation with the Moho boundary. In particular, the global correlation of the consolidated crust-stripped gravity disturbances with the Moho boundary is found to be –0.92. Among all the investigated gravity field quantities, the consolidated crust-stripped gravity disturbances are thus the best suited for a refinement of the Moho density interface by means of the gravimetric modeling or inversion.

With the rapid growth of digital photography, sharing of photos with friends and family has become very popular. When people share their photos, they usually organize them into albums according to events or places. To tell the story of some important events in one’s life, it is desirable to have an efficient summarization tool which can help people to receive a quick overview of an album containing large number of photos. In this paper, we present and analyze an approach for photo album summarization through a novel social game “Epitome” as a Facebook application. This social game can collect research data, and, at the same time, it provides a collage or a cover photo of the user’s photo album, while the user enjoys playing the game. The proof of concept of the proposed method is demonstrated through a set of experiments on several photo albums. As a benchmark comparison to this game, we perform automatic visual analysis considering several state-of-the-art features. We also evaluate the usability of the game by making use of a questionnaire on several subjects who played the “Epitome” game. Furthermore, we address privacy issues concerning shared photos in Facebook applications.