Continuous GPS Measurements Research Articles

GPS estimates of ocean tide loading in NW-France: determination of ocean tide loading constituents and comparison with a recent ocean tide model

SUMMARY Ocean loading effects cause 3-D displacements large enough to affect space geodesy measurements either at the subdiurnal periods or at longer time scales by the means of spurious signals. GPS measurements, in turn, could provide local improvements of the models in coastal areas if their ability to precisely monitor such rapid displacements is assessed. In this paper, we use 105 days of continuous GPS measurements collected in 2004 in the French Brittany and Cotentin region to investigate: (1) the precision achieved by the GPS analysis on measuring 3-D subdiurnal displacements and amplitude and phase of the tidal constituents, (2) the quality of the most recent ocean tide model FES2004 in such a complex coastal context. Indeed, in this area, tide amplitudes are among the largest in the world (up to 16 cm of loading displacements on the vertical component) and are believed to show strong shallow-water tides. From a state of the art GPS analysis using the scientific GAMIT software over 2-h sessions, we test two independent strategies for the realization of the reference frame. The position time series are then compared with the displacements predicted by the FES2004 model applied on an elastic Earth model. The two sets of results are consistent with each other at the same level of agreement than with the predicted displacements, namely 3‐5 mm on the horizontal components, 10 mm on the vertical. This assesses the capability of this technique for measuring 3-D ocean tide loading deformation. We validate the FES2004 model in the Brittany area, even though it slightly (2‐7 mm) underestimates the three components amplitudes. With a harmonic analysis of the observed position time series, we obtain nevertheless an agreement at a submillimetre level for the M2, N2, O1, Q1 tidal constituents and at a millimetre level for the K1 and S2 tidal constituents. Moreover, we can extract a significant M4 load signal at the 95 per cent confidence level from the GPS time series at the stations located in the Mont St Michel area. The detection of other shallow-water constituents would be helpful to understand the amplitude deficit between the FES2004 predicted and GPS observed displacements.

Tide-induced perturbations of glacier velocities

Recent observations showing substantial diurnal changes in velocities of glaciers flowing into the ocean, measured at locations far inland of glacier grounding lines, add fuel to the ongoing debate concerning the ability of glaciers to transmit longitudinal-stress perturbations over large distances. Resolution of this debate has major implications for the prediction of glacier mass balance, because it determines how rapidly a glacier can respond dynamically to changes such as weakening or removal of an ice shelf. Current IPCC assessment of sea-level rise takes little account of such changes, on the assumption that dynamic responses would be too slow to have any appreciable effect on ice discharge fluxes. However, this assumption must be questioned in view of observations showing massive increases in glacier velocities following removal of parts of the Larsen Ice Shelf, Antarctica, and of others showing diurnal velocity changes apparently linked to the tides. Here, I use a simple force-perturbation model to calculate the response of glacier strain rates to tidal rise and fall, assuming associated longitudinal-force perturbations are transmitted swiftly far inland of the glacier grounding line. Results show reasonable agreement with observations from an Alaskan glacier, where the velocity changes extended only a short distance up-glacier. However, for larger Antarctic glaciers, big velocity changes extending far upstream cannot be explained by this mechanism, unless ice-shelf “back forces” change substantially with the tides. Additional insight will require continuous measurement of velocity and strain-rate profiles along flow lines of glaciers and ice shelves. An example is suggested, involving continuous GPS measurements at a series of locations along the centre line of Glaciar San Rafael, Chile, extending from near the calving front to perhaps 20 km inland. Tidal range here is about ± 0.8 m, which should be sufficient to cause a variation in ice-front velocity of ± 2 cm h − 1 about its average value of 75 cm h − 1 , assuming local seawater depth of 150 m and glacier thickness of 200–400 m.

Coseismic Slip and Afterslip of the GreatMw 9.15 Sumatra–Andaman Earthquake of 2004

AbstractWe determine coseismic and the first-month postseismic deformation associated with the Sumatra–Andaman earthquake of 26 December 2004 from near- field Global Positioning System (gps) surveys in northwestern Sumatra and along the Nicobar-Andaman islands, continuous and campaign gps measurements from Thailand and Malaysia, and in situ and remotely sensed observations of the vertical motion of coral reefs. The coseismic model shows that the Sunda subduction megathrust ruptured over a distance of about 1500 km and a width of less than 150 km, releasing a total moment of 6.7–7.0 × 1022 N m, equivalent to a magnitude Mw ∼9.15. The latitudinal distribution of released moment in our model has three distinct peaks at about 4° N, 7° N, and 9° N, which compares well to the latitudinal variations seen in the seismic inversion and of the analysis of radiated T waves. Our coseismic model is also consistent with interpretation of normal modes and with the amplitude of very-long-period surface waves. The tsunami predicted from this model fits relatively well the altimetric measurements made by the jason and topex satellites. Neither slow nor delayed slip is needed to explain the normal modes and the tsunami wave. The near-field geodetic data that encompass both coseismic deformation and up to 40 days of postseismic deformation require that slip must have continued on the plate interface after the 500-sec-long seismic rupture. The postseismic geodetic moment of about 2.4 × 1022 N m (Mw ∼8.8) is equal to about 30 ± 5% of the coseismic moment release. Evolution of postseismic deformation is consistent with rate-strengthening frictional afterslip.Online material: Summary of geodetic data used in this study.

Space geodetic determination of spatial variability in relative sea level change, Los Angeles basin

We combine Synthetic Aperture Radar Interferometry (InSAR), tide gauge, and continuous GPS measurements to determine the spatial variation in vertical land motion (VLM) along the coast of the Los Angeles basin over the past decade, and to examine the impact of spatially variable VLM on relative sea level trends. By identifying radar scattering targets with long‐term coherence we make height corrections which allow interferogram creation for nearly the entire ERS‐1 catalog and permit estimation of average deformation rates with minimal temporal aliasing. Between Los Angeles Harbor and Newport Beach, mean VLM trends range from ∼3.4 to −4.3 mm/yr, reflecting the high level of ground water and oil extraction activity in the region. West of Los Angeles Harbor, VLM rates and spatial variability are roughly half as large. The 8‐year VLM trends exceed the long‐term sea level trend (0.8 mm/yr) determined from the 80 year Los Angeles Harbor tide gauge. The high degree of observed VLM variability emphasizes the need for the spatially continuous measurements provided by InSAR; a single tide gauge assessment of regional RSL would otherwise have limited applicability.

Current plate movements across the Mid‐Atlantic Ridge determined from 5 years of continuous GPS measurements in Iceland

We analyze data spanning up to 5 years from 18 continuous GPS stations in Iceland, computing daily positions of the stations with three different high‐level geodetic processing software packages. We observe large‐scale crustal deformation due to plate spreading across Iceland. The observed plate divergence between the North American and the Eurasian plates is in general agreement with existing models of plate motion. Spreading is taken up within a ∼100–150 km wide plate boundary zone that runs through the island. Of the two parallel branches of the plate boundary in south Iceland, the eastern volcanic zone is currently taking up the majority of the spreading and little is left for the western volcanic zone. The plate boundary deformation field has been locally and temporarily affected in south Iceland by two Mw = 6.5 earthquakes in June 2000, inflation at Katla volcano during 2000 to 2004, and an eruption of Hekla volcano in February 2000. All stations with significant vertical velocities are moving up relative to the reference station REYK, with the highest velocity exceeding 20 mm/yr in the center of the island.

Seismicity, deformation and seismic hazard in the western rift of Corinth: New insights from the Corinth Rift Laboratory (CRL)

This paper presents the main recent results obtained by the seismological and geophysical monitoring arrays in operation in the rift of Corinth, Greece. The Corinth Rift Laboratory (CRL) is set up near the western end of the rift, where instrumental seismicity and strain rate is highest. The seismicity is clustered between 5 and 10 km, defining an active layer, gently dipping north, on which the main normal faults, mostly dipping north, are rooting. It may be interpreted as a detachment zone, possibly related to the Phyllade thrust nappe. Young, active normal faults connecting the Aigion to the Psathopyrgos faults seem to control the spatial distribution of the microseismicity. This seismic activity is interpreted as a seismic creep from GPS measurements, which shows evidence for fast continuous slip on the deepest part on the detachment zone. Offshore, either the shallowest part of the faults is creeping, or the strain is relaxed in the shallow sediments, as inferred from the large NS strain gradient reported by GPS. The predicted subsidence of the central part of the rift is well fitted by the new continuous GPS measurements. The location of shallow earthquakes (between 5 and 3.5 km in depth) recorded on the on-shore Helike and Aigion faults are compatible with 50° and 60° mean dip angles, respectively. The offshore faults also show indirect evidence for high dip angles. This strongly differs from the low dip values reported for active faults more to the east of the rift, suggesting a significant structural or rheological change, possibly related to the hypothetical presence of the Phyllade nappe. Large seismic swarms, lasting weeks to months, seem to activate recent synrift as well as pre-rift faults. Most of the faults of the investigated area are in their latest part of cycle, so that the probability of at least one moderate to large earthquake ( M = 6 to 6.7) is very high within a few decades. Furthermore, the region west to Aigion is likely to be in an accelerated state of extension, possibly 2 to 3 times its mean interseismic value. High resolution strain measurement, with a borehole dilatometer and long base hydrostatic tiltmeters, started end of 2002. A transient strain has been recorded by the dilatometer, lasting one hour, coincident with a local magnitude 3.7 earthquake. It is most probably associated with a slow slip event of magnitude around 5 ± 0.5. The pore pressure data from the 1 km deep AIG10 borehole, crossing the Aigion fault at depth, shows a 1 MPa overpressure and a large sensitivity to crustal strain changes.

Deformation of the North American plate interior from a decade of continuous GPS measurements

A combination of two independent geodetic solutions using data from close to 300 continuous GPS stations covering the central and eastern United States shows that surface deformation in the North American plate interior is best fit by a model that includes rigid rotation of North America with respect to ITRF2000 and a component of strain qualitatively consistent with that expected from glacial isostatic adjustment (GIA). After correcting for the North American plate motion, residual horizontal velocities show a north‐to‐south deformation gradient of ∼1 mm yr−1, mostly localized between 1000 and 2200 km from the GIA center, corresponding to strain rates of about 10−9 yr−1. At distances farther than 2100 km from the GIA center, horizontal residual velocities are random with no evidence for regions of elevated strain rates. In particular, we find no detectable residual motion at the 95% confidence level in the New Madrid Seismic Zone, where the average weighted misfit of 0.7 mm yr−1 is the same as the weighted misfit of our rigid plate model. Vertical velocities show (1) a maximum uplift rate of 10 mm yr−1 at the assumed GIA center, (2) a hinge line located 1500 km from that center, and (3) a subsidence rate up to 1.4 mm yr−1 in the forebulge, with a maximum located about 2000 km from the GIA center. Our results have the potential to better constrain glacial isostatic adjustment models and contribute to a better definition of stable North America for tectonic and geodetic applications.

Contributions to the Deformation Analysis in Germany Based on Precise and Continuous GPS Measurements

Modern space geodetic techniques enable deformation monitoring of continental plate interiors with high spatial and temporal coverage. Resolving data and results are currently evaluated for their application for the integrated assessment of seismic hazard and risk in Germany. This goes especially for regions where earthquakes are generally rare but high magnitudes are still not unrealistic while vulnerability of today's society is steadily growing. The present contribution deals with the continuous monitoring of tectonic fracture systems in Germany using the GPS. The estimation of the station velocities with GPS and the resulting geodetic strain is supposed to provide additional input to the earthquake hazard assessment. Unfortunately, the low expected and currently seen velocities (<1-2 mm/year) make it extremely difficult to distinguish between noise and a tectonic signal. Because of the short observation interval the velocity uncertainties are about 2 mm/year in the horizontal components. The essential goal of this program is to provide and model highly precise deformation data and to discuss its needs for a better assessment of geological hazard, especially for the most active tectonic regions in Germany, the Rhine-Graben, the Swabian Alb, the Alpine foreland, and the Vogtland. Here we present preliminary results from 2 years of measurements at currently 150 GPS stations throughout Germany. The time span of this program has proven to be too short and the density of the station network to be not dense enough yet for reliable significant horizontal station velocities and supporting the earthquake hazard assessment.

Earthquake triggering by migrating slow slip initiated by M8 earthquake along Kuril Trench, Japan

Evidence of a slow slip temporally ahead of the M7 earthquake doublet that occurred in November/December 2004 offshore of eastern Hokkaido, Japan, was registered in continuous GPS measurements. The GPS data suggest that a slow slip around a forthcoming rupture zone on the plate boundary triggered these earthquakes. This slow slip episode was itself triggered by the M8 Tokachi earthquake of September 2003. This is a demonstration of the geodetic visualization of the mechanism of a triggering of a subsequent earthquake caused by a preceding large earthquake. The geodetic data will improve the reliability of the warning of an occurrence of a subsequent earthquake after a large one that is currently dependent on only co‐seismic static stress increase.

A quantitative assessment of the SBAS algorithm performance for surface deformation retrieval from DInSAR data

We investigate in this work the performance of the Small BAseline Subset (SBAS) approach that is a Differential Synthetic Aperture Radar Interferometry (DInSAR) algorithm allowing the generation of mean deformation velocity maps and displacement time series from a data set of subsequently acquired SAR images. In particular, we have carried out a quantitative assessment of the SBAS procedure performance by processing SAR data acquired by the European Remote Sensing Satellite (ERS) sensors and comparing the achieved results with geodetic measurements that are assumed as reference. The analysis has been focused on the Napoli bay (Italy) and Los Angeles (California) test areas where different deformation phenomena are present and, at the same time, a large amount of ERS SAR data is available as well as geometric leveling (in the Napoli zone) and continuous GPS (in the Los Angeles zone) measurements, to be used for our performance analysis. Moreover, due to the presence of large urbanized zones, the selected test sites are also characterized by extended, highly coherent areas in the DInSAR maps. The presented study shows that the SBAS technique provides an estimate of the mean deformation velocity with a standard deviation of about 1 mm/year for a typical ERS data set including between 40 and 60 images. Moreover, the single displacement measurements, computed with respect to a reference point of known motion, show a sub-centimetric accuracy with a standard deviation of about 5 mm, consistently in both the SAR/leveling and SAR/GPS comparisons; we also show that there is an increase of this standard deviation value as we move away from the reference SAR pixel, with an estimated spatial variation value of about 0.05 mm/km.

Gravity tide and seasonal gravity variation at Ny-Ålesund, Svalbard in Arctic

Abstract We have analyzed 4 years of gravity data (September 1999–August 2003) obtained from the superconducting gravimeter at Ny-Alesund, Svalbard in Arctic. The tidal analysis results indicate that the time variation in amplitude and phase of the short-period gravity tides may be due to the seasonal variation of the ocean tides, especially in the semidiurnal band. The gravity residuals, which were reduced from the following three effects from the original data: (1) the observed short- and long-period tides (up to the annual Sa wave) including the loading and attraction of ocean tide, (2) the air pressure and (3) the polar motion, show clear correlation with the computed hydrological effects (i.e. the effects of soil moisture and snow) in both the amplitude and phase. However, we found that the gravity data reveal a problem in the snow model at Ny-Alesund that was used for the computation. The vertical motion obtained from the continuous GPS measurements carried out near the gravity station shows seasonal variations that are consistent with the gravity data at least in their sign.

Geodetic constraints for the mechanism of Anatahan eruption of May 2003

Anatahan Island is located at the southern end of the Mariana volcanic chain. On May 10, 2003, the eastern crater of the island erupted for the first time in recorded history. The Plinian eruption column reached an altitude as high as 13 km on May 11, and a thick layer of ash covered the island. The eruptive activity continued to June 2003, but most of the erupted material was expelled during the first week. The volcanic activity declined in the second half of 2003, but resumed in April 2004. In order to determine crustal deformation associated with the eruption, we conducted GPS measurements in July 2003 at a benchmark (ANAT) located approximately 7 km west–northwest of the active crater, where GPS campaign measurements had been repeated four times since 1992. In the period from January to July 2003 during the eruption, significant subsidence – as much as 21 cm – was detected, but horizontal movement was negligible. We began taking continuous GPS measurements at the same site in July 2003 to monitor the transient deformation that was probably associated with magma migration. To assess the spatial extent of the deformation more accurately, we established another permanent GPS site (ANA2) at a site approximately 3 km from the active crater in the northeastern part of the island in January 2004. The coordinates of this time series at ANAT probably show a change in trends at the beginning of 2004. Another subsidence of 2.8 cm and a westward motion of 2.1 cm were estimated to have occurred in the period from July to December 2003. This was followed by an uplift of 5.2 cm and movement in an eastward direction of 1.0 cm in the period from January to June 2004. We developed three preliminarily models of inflation/deflation sources for three different time periods. During the period from January to July 2003, a deformation source was located beneath the ANAT site and acted as a deflation source. Considering the gap in the GPS time series and errors in data (especially after July 2003), we expected that the deformation sources were located beneath the western part of Anatahan Island and not below the active crater.

Extent and duration of the 2003 Cascadia slow earthquake

Inversion of continuous GPS measurements from the Pacific Northwest show the 2003 Cascadia slow earthquake to be among the largest of ten transients recognized here. Twelve stations bracketing slow slip indicate transient slip propagated bi‐directionally from initiation in the southern Puget basin, reaching 300 km along‐strike over a period of seven weeks. This event produced, for the first time, resolvable vertical subsidence, and horizontal displacement reaching six mm in southern Washington State. Inverted for non‐negative thrust slip, a maximum of 3.8 cm of slip is inferred, centered at 28 km depth near the sharp arch in the subducting Juan de Fuca plate. Nearly all slip lies shallower than 38 km. Inverted slip shows a total moment release equal to Mw = 6.6 and a high degree of spatial localization rather than near‐uniform slip. This suggests rupture concentrated along asperities holds for slow earthquakes as well as conventional events.

New slip rate estimates for the creeping segment of the San Andreas fault, California

New continuous and differential global positioning system (GPS) measurements of recent slip rates and 30 yr alignment-array offsets from the central creeping segment of the San Andreas fault yield a maximum right-lateral slip rate of 25 ′ 1 mm/yr. This slip rate is 20% slower than the 30 mm/yr slip rate accepted for this segment of the fault and 35% slower than the 39 mm/yr slip rate predicted between the Sierra Nevada-Great Valley block and the Pacific plate. New continuous GPS measurements between pairs of sites that flank the creeping segment at intersite distances of 1.0 km and 70 km give relative fault-parallel slip rates of 28 ′ 2 and 30 ′ 2 mm/yr, respectively. These observations indicate that right-lateral deformation rates increase with distance from the fault. Possible explanations for the gradient observed in the geodetic data are elastic strain accumulation along the creeping segment or significant distributed deformation on off-fault structures.

Spatial and temporal evolution of rapid basal sliding on Bench Glacier, Alaska, USA

AbstractWe measured the surface velocity field during the summers of 1999 and 2000 on the 7 km long, 185 m thick Bench Glacier, Alaska, USA. In the spring of both years, a short-lived pulse of surface velocity, 2-4 times the annual mean velocity, propagated up-glacier from the terminus at a rate of ~200-250md-1. Displacement attributable to rapid sliding is ~5-10% of the annual surface motion, while the high-velocity event comprised 60-95% of annual basal motion. Sliding during the propagating speed-up event peaked at 6-14 cm d-1, with the highest rates in mid-glacier. Continuous horizontal and vertical GPS measurements at one stake showed divergence and then convergence of the ice surface with the bed as the velocity wave passed, with maximum surface uplift of 8-16 cm. High divergence rates coincided with high horizontal velocities, suggesting rapid sliding on the up-glacier side of bedrock steps. Initiation of the annual speed-up event occurred during the peak in englacial water storage, while the glacier was entirely snow-covered. Basal motion during the propagating speed-up event enlarges cavities and connections among them, driving a transition from a poorly connected hydrologic system to a well-connected linked-cavity system. Sliding is probably halted by the development of a conduit system.

P08 GPSによる富士山周辺の歪速度分布 : 富士山の側噴火火口配列の謎

P08 GPSによる富士山周辺の歪速度分布 : 富士山の側噴火火口配列の謎

Crustal deformation measurements in Guerrero, Mexico

GPS measurements of crustal deformation in Guerrero, southern Mexico, include surveys collected between 1992 and 2001 as well as continuous GPS measurements at a few sites. These geodetic observations are used to calculate interseismic deformation rates and assess the presence and possible location of transient deformation during the period encompassing 1992.25 to 2001.75. The data are used to examine transient deformation in 1998 previously described from data at a single site by Lowry et al. [2001]. Survey measurements and continuous data from a site near Popocatépetl volcano confirm the 1998 transient, and survey data also suggest another transient occurred following the 14 September 1995 (Mw = 7.3) Copala earthquake. All of the available GPS position estimates have been inverted for a combined model of slip during each event plus the steady state slip on the plate interface. Modeling of the steady state deformation rates confirms that the Guerrero seismic gap is partially frictionally locked at depths shallower than about 25 km and accumulating strain that may eventually be released in a great earthquake. The data also suggest that there is frictional coupling to much greater (&gt;40 km) depths, which releases more frequently in aseismic slip events. The locations and sizes of the transient events are only partially constrained by the available data. However, the transient models which best fit the GPS coordinate time series suggest that aseismic slip was centered downdip of the seismogenic portion of the plate‐bounding thrust in both events, and the moment release had equivalent magnitudes Mw = 7.1 + 1.3/−1.0 in 1995–1996 and 7.1 + 0.4/−0.1 in 1998.

Continuous GPS measurements of postglacial adjustment in Fennoscandia: 2. Modeling results

Data collected under the auspices of the BIFROST GPS project yield a geographically dense suite of estimates of present‐day, three‐dimensional (3‐D) crustal deformation rates in Fennoscandia [Johansson et al., 2002]. A preliminary forward analysis of these estimates [Milne et al., 2001] has indicated that models of ongoing glacial isostatic adjustment (GIA) in response to the final deglaciation event of the current ice age are able to provide an excellent fit to the observed 3‐D velocity field. In this study we revisit our previous GIA analysis by considering a more extensive suite of forward calculations and by performing the first formal joint inversion of the BIFROST rate estimates. To establish insight into the physics of the GIA response in the region, we begin by decomposing a forward prediction into the three contributions associated with the ice, ocean, and rotational forcings. From this analysis we demonstrate that recent advances in postglacial sea level theory, in particular the inclusion of rotational effects and improvements in the treatment of the ocean load in the vicinity of an evolving continental margin, involve peak signals that are larger than the observational uncertainties in the BIFROST network. The forward analysis is completed by presenting predictions for a pair of Fennoscandian ice histories and an extensive suite of viscoelastic Earth models. The former indicates that the BIFROST data set provides a powerful discriminant of such histories. The latter yields bounds on the (assumed constant) upper and lower mantle viscosity (νUM, νLM); specifically, we derive a 95% confidence interval of 5 × 1020 ≤ νUM ≤ 1021 Pa s and 5 × 1021 ≤ νLM ≤ 5 × 1022 Pa s, with some preference for (elastic) lithospheric thickness in excess of 100 km. The main goal of the (Bayesian) inverse analysis is to estimate the radial resolving power of the BIFROST GPS data as a function of depth in the mantle. Assuming a reasonably accurate ice history, we demonstrate that this resolving power varies from ∼200 km near the base of the upper mantle to ∼700 km in the top portion of the lower mantle. We conclude that the BIFROST data are able to resolve structure on radial length scale significantly smaller than a single upper mantle layer. However, these data provide little constraint on viscosity in the bottom half of the mantle. Finally, elements of both the forward and inverse analyses indicate that radial and horizontal velocity estimates provide distinct constraints on mantle viscosity.

GPS連続観測による日本列島上下地殻変動とその意義

GPS連続観測による日本列島上下地殻変動とその意義

Dynamics of active magmatic and hydrothermal systems at Taal Volcano, Philippines, from continuous GPS measurements

A dense network of continuous single‐ and dual‐frequency GPS receivers at Taal Volcano, Philippines, reveals four major stages of volcanogenic deformation: deflation, from installation in June 1998 to December 1998; inflation, from January to March 1999; deflation, from April 1999 to February 2000; and inflation, from February to November 2000. The largest displacements occurred during the February–November 2000 period of inflation, which was characterized by ∼120 mm of uplift of the center of Volcano Island relative to the northern caldera rim at average rates up to 216 mm/yr. Deformation sources were modeled using a constrained least squares inversion algorithm. The source of 1999 deflation and inflation in 2000 were modeled as contractional and dilatational Mogi point sources centered at 4.2 and 5.2 km depth, respectively, beneath Volcano Island. The locations of the inflationary and deflationary sources are indistinguishable within the 95% confidence estimates. Modeling using a running 4‐month time window from June 1999 to March 2001 reveals little evidence for source migration. We suggest that the two periods of inflation observed at Taal result from episodic intrusions of magma into a shallow reservoir centered beneath Volcano Island. Subsequent deflation may result from exsolution of magmatic fluids and/or gases into an overlying, unconfined hydrothermal system.