1207 research outputs found
Sort by
Method to observe Jupiter’s radio emissions at high resolution using multiple LOFAR stations: a first case study of the Io-decametric emission using the Irish IE613, French FR606 and German DE604 stations
Get PDFThe Low Frequency Array (LOFAR) is an international radio telescope array, consisting of 38 stations in the Netherlands and 14 international stations spread over Europe. Here we present an observation method to study the jovian decametric radio emissions from several LOFAR stations (here Birr Castle in Ireland, Nançay in France and Postdam in Germany), at high temporal and spectral resolution. This method is based on prediction tools, such as radio emission simulations and probability maps, and data processing. We report an observation of Io-induced decametric emission from June 2021, and a first case study of the substructures that compose the macroscopic emissions (called millisecond bursts). The study of these bursts make it possible to determine the electron populations at the origin of these emissions. We then present several possible future avenues for study based on these observations. The methodology and study perspectives described in this paper can be applied to new observations of jovian radio emissions induced by Io, but also by Ganymede or Europa, or jovian auroral radio emissions
Caldera resurgence during the 2018 eruption of Sierra Negra volcano, Galápagos Islands
Get PDFRecent large basaltic eruptions began after only minor surface uplift and seismicity, and resulted in caldera subsidence. In contrast, some eruptions at Galápagos Island volcanoes are preceded by prolonged, large amplitude uplift and elevated seismicity. These systems also display long-term intra-caldera uplift, or resurgence. However, a scarcity of observations has obscured the mechanisms underpinning such behaviour. Here we combine a unique multiparametric dataset to show how the 2018 eruption of Sierra Negra contributed to caldera resurgence. Magma supply to a shallow reservoir drove 6.5m of pre-eruptive uplift and
seismicity over thirteen years, including an Mw5.4 earthquake that triggered the eruption.
Although co-eruptive magma withdrawal resulted in 8.5m of subsidence, net uplift of the inner-caldera on a trapdoor fault resulted in 1.5m of permanent resurgence. These
observations reveal the importance of intra-caldera faulting in affecting resurgence, and the mechanisms of eruption in the absence of well-developed rift systems
Retrieving reflection arrivals from passive seismic data using Radon correlation
Get PDFSince explosive and impulsive seismic sources such as dynamite, air guns, gas guns or even vibroseis can have a big impact on the environment, some companies have decided to record ambient seismic noise and use it to estimate the physical properties of the subsurface. Big challenges arise when the aim is extracting body waves from recorded passive signals, especially in the presence of strong surface waves. In passive seismic signals, such body waves are usually
weak in comparison to surface waves that are much more prominent. To understand the characteristics of passive signals and the effect of natural source locations, three simple synthetic models were created. To extract body waves from simulated passive signals we propose and test a
Radon-correlation method. This is a time-spatial correlation of amplitudes with a train of time-shifted Dirac delta functions through different hyperbolic paths. It is tested on a two-layer horizontal model, a three-layer model that includes a dipping layer (with and without lateral
heterogeneity) and also on synthetic Marmousi model data sets. Synthetic tests show that the introduced method is able to reconstruct reflection events at the correct time-offset positions that are hidden in results obtained by the general cross-correlation method. Also, a depth migrated section shows a good match between imaged horizons and the true model. It is possible to generate off-end virtual gathers by applying the method to a linear array of receivers and to construct a velocity model by semblance velocity analysis of individually extracted gathers
Wind-envelope interaction as the origin of the slow cyclic brightness variations of luminous blue variables
Get PDFLuminous blue variables (LBVs) are hot, very luminous massive stars displaying large quasi-periodic variations in brightness, radius, and photospheric temperature on timescales of years to decades. The physical origin of this variability, called S Doradus cycle after its prototype, has remained elusive. We study the feedback of stellar wind mass-loss on the envelope structure in stars near the Eddington limit. We calculated a time-dependent hydrodynamic stellar evolution, applying a stellar wind mass-loss prescription with a temperature dependence inspired by the predicted systematic increase in mass-loss rates below 25 kK. We find that when the wind mass-loss rate crosses a well-defined threshold, a discontinuous change in the wind base conditions leads to a restructuring of the stellar envelope. The induced drastic radius and temperature changes, which occur on the thermal timescale of the inflated envelope, in turn impose mass-loss variations that reverse the initial changes, leading to a cycle that lacks a stationary equilibrium configuration. Our proof-of-concept model broadly reproduces the typical observational phenomenology of the S Doradus variability. We identify three key physical ingredients that are required to trigger the instability: inflated envelopes in close proximity to the Eddington limit, a temperature range where decreasing opacities do not lead to an accelerating outflow, and a mass-loss rate that increases with decreasing temperature, crossing a critical threshold value within this temperature range. Our scenario and model provide testable predictions, and open the door for a consistent theoretical treatment of the LBV phase in stellar evolution, with consequences for their further evolution as single stars or in binary systems
PION: simulating bow shocks and circumstellar nebulae
Get PDFExpanding nebulae are produced by mass-loss from stars, especially during late stages of evolution. Multidimensional simulation of these nebulae requires high resolution near the star and permits resolution that decreases with distance from the star, ideally with adaptive time-steps. We report the implementation and testing of static mesh-refinement in the radiation-magnetohydrodynamics (R-MHD) code PION, and document its performance for 2D and 3D calculations. The bow shock produced by a hot, magnetized, slowly rotating star as it moves through the magnetized ISM is simulated in 3D, highlighting differences compared with 2D calculations. Latitude-dependent, time-varying magnetized winds are modelled and compared with simulations of ring nebulae around blue supergiants from the literature. A 3D simulation of the expansion of a fast wind from a Wolf-Rayet star into the slow wind from a previous red supergiant phase of evolution is presented, with results compared with results in the literature and analytic theory. Finally, the wind-wind collision from a binary star system is modelled with 3D MHD, and the results compared with previous 2D hydrodynamic calculations. A PYTHON library is provided for reading and plotting simulation snapshots, and the generation of synthetic infrared emission maps using TORUS is also demonstrated. It is shown that state-of-the-art 3D MHD simulations of wind-driven nebulae can be performed using PION with reasonable computational resources. The source code and user documentation is made available for the community under a BSD3 licence
North Atlantic Oscillation (NAO) climate index hidden in ocean generated secondary microseisms
Get PDFWith the role of surface winds in the generation of ocean waves, secondary microseisms induced by ocean wave-wave interactions represent a unique interconnection between the solid Earth, the ocean and global atmospheric circulation patterns. In this study, temporal changes in the ocean generated seismic wavefield in Northeast Atlantic are monitored offshore West of Ireland using ocean bottom seismometers located on top of a thick sedimentary basin. Comparisons with numerical seismic simulations and ocean wave model hindcast data suggest those variations are correlated with changing patterns in ocean wave interactions closely linked to secondary microseism generation areas. Here we show how those changes, accentuated by the specific structure of the Irish offshore margin, reveal the signature of the North Atlantic Oscillation (NAO) in the Earth’s background seismic wavefield. In the North Atlantic, secondary microseism sources likely fluctuate with the changing storm track in response to variations in the NAO
Pre-migration diffraction separation using generative adversarial networks
Get PDFDiffraction imaging is the process of separating diffraction events from the seismic wavefield and imaging them independently, highlighting subsurface discontinuities.
While there are many analytic-based methods for diffraction imaging which use kinematic, dynamic or both, properties of the diffracted wavefield, they can be slow and require parameterization. Here, we propose an image-to-image generative adversarial network to automatically separate diffraction events on pre-migrated seismic data in a fraction of the time of conventional methods. To train the generative adversarial network, plane-wave destruction was applied to a range of synthetic and real images from field data to create training data. These training data were screened and any areas where the plane-wave destruction did not perform well, such as synclines and areas of complex dip, were removed to prevent bias in the neural network. A total of 14,132 screened images were used to train the final generative adversarial network.
The trained network has been applied across several geologically distinct field datasets, including a 3D example. Here, generative adversarial network separation
is shown to be comparable to a benchmark separation created with plane-wave destruction, and up to 12 times faster. This demonstrates the clear potential in generative adversarial networks for fast and accurate diffraction separation
Optimized workflows for high-frequency seismic interferometry using dense arrays
Get PDFHigh-frequency seismic surface waves sample the top few tens of meters to the top few kilometres of the subsurface. They can be used to determine 3-D distributions of shear-wave
velocities and to map the depths of discontinuities (interfaces) within the crust. Passive seismic imaging, using ambient noise as the source of signal, can thus be an effective tool of exploration for mineral, geothermal and other resources, provided that sufficient high-frequency signal is available in the ambient noise wavefield and that accurate, high-frequency measurements can be performed on this signal. Ambient noise imaging using the ocean-generated noise at 5–30 s periods is now a standard method, but less signal is available at frequencies high enough for deposit-scale imaging (0.2–30 Hz), and few studies have reported successful measurements in broad frequency bands. Here, we develop a workflow for the measurement of high-frequency, surface wave phase velocities in very broad frequency ranges. Our workflow comprises (1) a new noise cross-correlation procedure that accounts for the non-stationary
properties of the high-frequency noise sources, removes bandpass filtering, replaces temporal normalization with short time window stacking, and drops the explicit spectral normalization by adopting cross-coherence; (2) a new phase-velocity measurement method that extends the bandwidth of reliable measurements by exploiting the (resolved) 2π ambiguity of phasevelocity measurements and (3) interstation-distance-dependent quality control that uses the similarity of subgroups of dispersion curves to reject outliers and identify the frequency ranges with accurate measurements. The workflow is highly automated and applicable to large arrays. Applying our method to data from a large-N array that operated for one month near Marathon, Ontario, Canada, we use rectangular subarrays with 150-m station spacing and, typically, 1 hr of data and obtain Rayleigh-wave phase-velocity measurements in a 0.5–30 Hz
frequency range, spanning over 5.9 octaves, twice the typical frequency range of 1.5–3 octaves in previous studies. Phase-velocity maps and the subregion-average 1-D velocity models they constrain show a high-velocity anomaly consistent with the known, west-dipping gabbro intrusions beneath the area. The new structural information can improve our understanding of the geometry of the gabbro intrusions, hosting the Cu-PGE Marathon deposi
Retrieving reflection arrivals from passive seismic data using Radon correlation
Get PDFSince explosive and impulsive seismic sources such as dynamite, air guns, gas guns or even vibroseis can have a big impact on the environment, some companies have decided to record ambient seismic noise and use it to estimate the physical properties of the subsurface. Big challenges arise when the aim is extracting body waves from recorded passive signals, especially in the presence of strong surface waves. In passive seismic signals, such body waves are usually weak in comparison to surface waves that are much more prominent. To understand the characteristics of passive signals and the effect of natural source locations, three simple synthetic models were created. To extract body waves from simulated passive signals we propose and test a Radon-correlation method. This is a time-spatial correlation of amplitudes with a train of time-shifted Dirac delta functions through different hyperbolic paths. It is tested on a two-layer horizontal model, a three-layer model that includes a dipping layer (with and without lateral heterogeneity) and also on synthetic Marmousi model data sets. Synthetic tests show that the introduced method is able to reconstruct reflection events at the correct time-offset positions that are hidden in results obtained by the general cross-correlation method. Also, a depth migrated section shows a good match between imaged horizons and the true model. It is possible to generate off-end virtual gathers by applying the method to a linear array of receivers and to construct a velocity model by semblance velocity analysis of individually extracted gathers
Observing Jupiter's radio emissions using multiple LOFAR stations: a first case study of the Io-decametric emission using the Irish IE613, French FR606 and German DE604 stations
Get PDFThe Low Frequency Array (LOFAR) is an international radio telescope array, consisting of 38 stations in the Netherlands and 14 international stations spread over Europe. Here we present an observation method to study the jovian decametric radio emissions from several LOFAR stations (here DE604, FR606 and IE613), at high temporal and spectral resolution. This method is based on prediction tools, such as radio emission simulations and probability maps, and data processing. We report an observation of Io-induced decametric emission from June 2021, and a first case study of the substructures that compose the macroscopic emissions (called millisecond bursts). The study of these bursts make it possible to determine the electron populations at the origin of these emissions. We then present several possible future avenues for study based on these observations. The methodology and study perspectives described in this paper can be applied to new observations of jovian radio emissions induced by Io, but also by Ganymede or Europa, or jovian auroral radio emissions