SEARCH · USGS Science

Results for “Journal of Geophysical Research - Atmospheres”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5Linked to original sources

USGS2014

Linking channel hydrology with riparian wetland accretion in tidal rivers

The hydrologic processes by which tide affects river channel and riparian morphology within the tidal freshwater zone are poorly understood, yet are fundamental to predicting the fate of coastal rivers and wetlands as sea level rises. We investigated patterns of sediment accretion in riparian wetlands along the non-tidal through oligohaline portion of two coastal plain rivers in Maryland, U.S.A., and how flow velocity, water level, and suspended sediment concentration (SSC) in the channel may have contributed to those patterns. Sediment accretion was measured over a one year period using artificial marker horizons, channel hydrology was measured over a one month period using acoustic Doppler current profilers, and SSC was predicted from acoustic backscatter. Riparian sediment accretion was lowest at the non-tidal sites (mean and standard deviation = 8 ± 8 mm yr -1 ), highest at the upstream tidal freshwater forested wetlands (TFFW) (33 ± 28 mm yr -1 ), low at the midstream TFFW (12 ± 9 mm yr -1 ), and high at the oligohaline (fresh-to-brackish) marshes (19 ± 8 mm yr -1 ). Channel maximum flood and ebb velocity was 2-fold faster at the oligohaline than tidal freshwater zone on both tidal rivers, corresponding with the differences in in-channel SSC: the oligohaline zone's SSC was more than double the tidal freshwater zone's, and was greater than historical SSC at the non-tidal gages. The tidal wave characteristics differed between rivers, leading to significantly greater in-channel SSC during floodplain inundation in the weakly convergent than the strongly convergent tidal river. Overall sediment accretion was higher in the embayed river likely due to a single storm discharge and associated sedimentation.

Maryl↗

USGS2014

Modeling erosion of ice-rich permafrost bluffs along the Alaskan Beaufort Sea coast

The Arctic climate is changing, inducing accelerating retreat of ice-rich permafrost coastal bluﬀs. Along Alaska’s Beaufort Sea coast, erosion rates have increased roughly threefold from 6.8 to 19 m yr −1 since 1955 while the sea ice-free season has increased roughly twofold from 45 to 100 days since 1979. We develop a numerical model of bluﬀ retreat to assess the relative roles of the length of sea ice-free season, sea level, water temperature, nearshore waveﬁeld, and permafrost temperature in controlling erosion rates in this setting. The model captures the processes of erosion observed in short-term monitoring experiments along the Beaufort Sea coast, including evolution of melt notches, topple of ice wedge-bounded blocks, and degradation of these blocks. Model results agree with time-lapse imagery of bluﬀ evolution and time series of ocean-based instrumentation. Erosion is highly episodic with 40% of erosion is accomplished during less than 5% of the sea ice-free season. Among the formulations of the submarine erosion rate we assessed, we advocate those that employ both water temperature and nearshore waveﬁeld. As high water levels are a prerequisite for erosion, any future changes that increase the frequency with which water levels exceed the base of the bluﬀs will increase rates of coastal erosion. The certain increases in sea level and potential changes in storminess will both contribute to this eﬀect. As water temperature also inﬂuences erosion rates, any further expansion of the sea ice-free season into the midsummer period of greatest insolation is likely to result in an additional increase in coastal retreat rates.

Alaska↗

USGS2014

Seismicity within a propagating ice shelf rift: the relationship between icequake locations and ice shelf structure

Iceberg calving is a dominant mass loss mechanism for Antarctic ice shelves, second only to basal melting. An important known process involved in calving is the initiation and propagation of through-penetrating fractures called rifts; however, the mechanisms controlling rift propagation remain poorly understood. To investigate the mechanics of ice-shelf rifting, we analyzed seismicity associated with a propagating rift tip on the Amery Ice Shelf, using data collected during the Austral summers of 2004-2007. We investigated seismicity associated with fracture propagation using a suite of passive seismological techniques including icequake locations, back projection, and moment tensor inversion. We confirm previous results that show that seismicity is characterized by periods of relative quiescence punctuated by swarms of intense seismicity of one to three hours. However, even during periods of quiescence, we find significant seismic deformation around the rift tip. Moment tensors, calculated for a subset of the largest icequakes (M W > -2.0) located near the rift tip, show steeply dipping fault planes, horizontal or shallowly plunging stress orientations, and often have a significant volumetric component. They also reveal that much of the observed seismicity is limited to the upper 50 m of the ice shelf. This suggests a complex system of deformation that involves the propagating rift, the region behind the rift tip, and a system of rift-transverse crevasses. Small-scale variations in the mechanical structure of the ice shelf, especially rift-transverse crevasses and accreted marine ice, play an important role in modulating the rate and location of seismicity associated with propagating ice shelf rifts.

Journal of Geophysical Research F: Earth Surface↗

USGS2014

Predictions of barrier island berm evolution in a time-varying storm climatology

Low-lying barrier islands are ubiquitous features of the world's coastlines, and the processes responsible for their formation, maintenance, and destruction are related to the evolution of smaller, superimposed features including sand dunes, beach berms, and sandbars. The barrier island and its superimposed features interact with oceanographic forces (e.g., overwash) and exchange sediment with each other and other parts of the barrier island system. These interactions are modulated by changes in storminess. An opportunity to study these interactions resulted from the placement and subsequent evolution of a 2 m high sand berm constructed along the northern Chandeleur Islands, LA. We show that observed berm length evolution is well predicted by a model that was fit to the observations by estimating two parameters describing the rate of berm length change. The model evaluates the probability and duration of berm overwash to predict episodic berm erosion. A constant berm length change rate is also predicted that persists even when there is no overwash. The analysis is extended to a 16 year time series that includes both intraannual and interannual variability of overwash events. This analysis predicts that as many as 10 or as few as 1 day of overwash conditions would be expected each year. And an increase in berm elevation from 2 m to 3.5 m above mean sea level would reduce the expected frequency of overwash events from 4 to just 0.5 event-days per year. This approach can be applied to understanding barrier island and berm evolution at other locations using past and future storm climatologies.

Louisiana↗

USGS2014

Transport of fine sediment over a coarse, immobile riverbed

Sediment transport in cobble-boulder rivers consists mostly of fine sediment moving over a coarse, immobile bed. Transport rate depends on several interrelated factors: boundary shear stress, the grain size and volume of fine sediment, and the configuration of fine sediment into interstitial deposits and bed forms. Existing models do not incorporate all of these factors. Approaches that partition stress face a daunting challenge because most of the boundary shear is exerted on immobile grains. We present an alternative approach that divides the bed into sand patches and interstitial deposits and is well constrained by two clear end-member cases: full sand cover and absence of sand. Entrainment from sand patches is a function of their aerial coverage. Entrainment from interstices among immobile grains is a function of sand elevation relative to the size of the immobile grains. The bed-sand coverage function is used to predict the ratio of the rate of entrainment from a partially covered bed to the rate of entrainment from a completely sand-covered bed, which is determined using a standard sand transport model. We implement the bed-sand coverage function in a morphodynamic routing model and test it against observations of sand bed elevation and suspended sand concentration for conditions of nonuniform fine sediment transport in a large flume with steady uniform flow over immobile hemispheres. The results suggest that this approach may provide a simple and robust method for predicting the transport and migration of fine sediment through rivers with coarse, immobile beds.

Journal of Geophysical Research F: Earth Surface↗

USGS2014

Rapid reservoir erosion, hyperconcentrated flow, and downstream deposition triggered by breaching of 38 m tall Condit Dam, White Salmon River, Washington

Condit Dam on the White Salmon River, Washington, a 38 m high dam impounding a large volume (1.8 million m 3 ) of fine-grained sediment (60% sand, 35% silt and clay, and 5% gravel), was rapidly breached in October 2011. This unique dam decommissioning produced dramatic upstream and downstream geomorphic responses in the hours and weeks following breaching. Blasting a 5 m wide hole into the base of the dam resulted in rapid reservoir drawdown, abruptly releasing ~1.6 million m 3 of reservoir water, exposing reservoir sediment to erosion, and triggering mass failures of the thickly accumulated reservoir sediment. Within 90 min of breaching, the reservoir's water and ~10% of its sediment had evacuated. At a gauging station 2.3 km downstream, flow increased briefly by 400 m 3  s −1 during passage of the initial pulse of released reservoir water, followed by a highly concentrated flow phase—up to 32% sediment by volume—as landslide-generated slurries from the reservoir moved downstream. This hyperconcentrated flow, analogous to those following volcanic eruptions or large landslides, draped the downstream river with predominantly fine sand. During the ensuing weeks, suspended-sediment concentration declined and sand and gravel bed load derived from continued reservoir erosion aggraded the channel by >1 m at the gauging station, after which the river incised back to near its initial elevation at this site. Within 15 weeks after breaching, over 1 million m 3 of suspended load is estimated to have passed the gauging station, consistent with estimates that >60% of the reservoir's sediment had eroded. This dam removal highlights the influence of interactions among reservoir erosion processes, sediment composition, and style of decommissioning on rate of reservoir erosion and consequent downstream behavior of released sediment.

Washington↗

USGS2013

Seasonal changes in peatland surface elevation recorded at GPS stations in the Red Lake Peatlands, northern Minnesota, USA

Northern peatlands appear to hold large volumes of free‐phase gas (e.g., CH 4 and CO 2 ), which has been detected by surface deformations, pore pressure profiles, and electromagnetic surveys. Determining the gas content and its impact in peat is challenging because gas storage depends on both the elastic properties of the peat matrix and the buoyant forces exerted by pore fluids. We therefore used a viscoelastic deformation model to estimate these variables by adjusting model runs to reproduce observed changes in peat surface elevation within a 1300 km 2 peatland. A local GPS network documented significant changes in surface elevations throughout the year with the greatest vertical displacements associated with rapid changes in peat water content and unloadings due to melting of the winter snowpack. These changes were coherent with changes in water table elevation and also abnormal pore pressure changes measured by nests of instrumented piezometers. The deformation model reproduced these changes when the gas content was adjusted to 10% of peat volume, and Young's modulus was varied between 5 and 100 kPa as the peat profile shifted from tension to compression. In contrast, the model predicted little peat deformation when the gas content was 3% or lower. These model simulations are consistent with previous estimates of gas volume in northern peatlands and suggest an upper limit of gas storage controlled by the elastic moduli of the peat fabric.

Minnesota↗

USGS2014

Floodplain biogeochemical processing of floodwaters in the Atchafalaya River Basin during the Mississippi River flood of 2011

The 2011 flood in the Lower Mississippi resulted in the second highest recorded river flow diverted into the Atchafalaya River Basin (ARB). The higher water levels during the flood peak resulted in high hydrologic connectivity between the Atchafalaya River and floodplain, with up to 50% of the Atchafalaya River water moving off channel. Water quality samples were collected throughout the ARB over the course of the flood event. Significant nitrate (NO 3 -) reduction (75%) occurred within the floodplain, resulting in a total NO 3 - reduction of 16.6% over the flood. The floodplain was a small but measurable source of dissolved reactive phosphorus (SRP) and ammonium (NH 4 +). Collectively, these results from this large flood event suggest that enhancing river-floodplain connectivity through freshwater diversions will reduce NO 3 - loads to the Gulf of Mexico during large annual floods.

Louisiana↗

USGS2014

The impact of climate and reservoirs on longitudinal riverine carbon ﬂuxes from two major watersheds in the Central and Intermontane West

A nested sampling network on the Colorado (CR) and Missouri Rivers (MR) provided data to assess impacts of large-scale reservoir systems and climate on carbon export. The Load Estimator (LOADEST) model was used to estimate both dissolved inorganic and organic carbon (DIC and DOC) fluxes for a total of 22 sites along the main stems of the CR and MR. Both the upper CR and MR DIC and DOC fluxes increased longitudinally, but the lower CR fluxes decreased while the lower MRs continued to increase. We examined multiple factors through space and time that help explain these flux patterns. Seasonal variability in precipitation and temperature, along with site-level concentration versus discharge relationships proved to be significant factors explaining much of the difference among sites located below reservoirs as compared to sites located in more free-flowing segments of the river. The characterization of variability in carbon exports over space and time provides a basis for understanding carbon cycling and transport within river basins affected by large reservoir systems, particular in arid-to semi-arid ecosystems.

Journal of Geophysical Research: Biogeosciences↗

USGS2014

Spatial variability and landscape controls of near-surface permafrost within the Alaskan Yukon River Basin

The distribution of permafrost is important to understand because of permafrost's influence on high-latitude ecosystem structure and functions. Moreover, near-surface (defined here as within 1 m of the Earth's surface) permafrost is particularly susceptible to a warming climate and is generally poorly mapped at regional scales. Subsequently, our objectives were to (1) develop the first-known binary and probabilistic maps of near-surface permafrost distributions at a 30 m resolution in the Alaskan Yukon River Basin by employing decision tree models, field measurements, and remotely sensed and mapped biophysical data; (2) evaluate the relative contribution of 39 biophysical variables used in the models; and (3) assess the landscape-scale factors controlling spatial variations in permafrost extent. Areas estimated to be present and absent of near-surface permafrost occupy approximately 46% and 45% of the Alaskan Yukon River Basin, respectively; masked areas (e.g., water and developed) account for the remaining 9% of the landscape. Strong predictors of near-surface permafrost include climatic indices, land cover, topography, and Landsat 7 Enhanced Thematic Mapper Plus spectral information. Our quantitative modeling approach enabled us to generate regional near-surface permafrost maps and provide essential information for resource managers and modelers to better understand near-surface permafrost distribution and how it relates to environmental factors and conditions.

Alaska↗

USGS2014

Seismicity, faulting, and structure of the Koyna-Warna seismic region, Western India from local earthquake tomography and hypocenter locations

Although seismicity near Koyna Reservoir (India) has persisted for ~50 years and includes the largest induced earthquake ( M 6.3) reported worldwide, the seismotectonic framework of the area is not well understood. We recorded ~1800 earthquakes from 6 January 2010 to 28 May 2010 and located a subset of 343 of the highest-quality earthquakes using the tomoDD code of Zhang and Thurber (2003) to better understand the framework. We also inverted first arrivals for 3-D Vp , Vs , and Vp / Vs and Poisson's ratio tomography models of the upper 12 km of the crust. Epicenters for the recorded earthquakes are located south of the Koyna River, including a high-density cluster that coincides with a shallow depth (<1.5 km) zone of relatively high Vp and low Vs (also high Vp / Vs and Poisson's ratios) near Warna Reservoir. This anomalous zone, which extends near vertically to at least 8 km depth and laterally northward at least 15 km, is likely a water-saturated zone of faults under high pore pressures. Because many of the earthquakes occur on the periphery of the fault zone, rather than near its center, the observed seismicity-velocity correlations are consistent with the concept that many of the earthquakes nucleate in fractures adjacent to the main fault zone due to high pore pressure. We interpret our velocity images as showing a series of northwest trending faults locally near the central part of Warna Reservoir and a major northward trending fault zone north of Warna Reservoir.

Maharashtra↗

USGS2014

Slip rates and spatially variable creep on faults of the northern San Andreas system inferred through Bayesian inversion of Global Positioning System data

Fault creep, depending on its rate and spatial extent, is thought to reduce earthquake hazard by releasing tectonic strain aseismically. We use Bayesian inversion and a newly expanded GPS data set to infer the deep slip rates below assigned locking depths on the San Andreas, Maacama, and Bartlett Springs Faults of Northern California and, for the latter two, the spatially variable interseismic creep rate above the locking depth. We estimate deep slip rates of 21.5 ± 0.5, 13.1 ± 0.8, and 7.5 ± 0.7 mm/yr below 16 km, 9 km, and 13 km on the San Andreas, Maacama, and Bartlett Springs Faults, respectively. We infer that on average the Bartlett Springs fault creeps from the Earth's surface to 13 km depth, and below 5 km the creep rate approaches the deep slip rate. This implies that microseismicity may extend below the locking depth; however, we cannot rule out the presence of locked patches in the seismogenic zone that could generate moderate earthquakes. Our estimated Maacama creep rate, while comparable to the inferred deep slip rate at the Earth's surface, decreases with depth, implying a slip deficit exists. The Maacama deep slip rate estimate, 13.1 mm/yr, exceeds long-term geologic slip rate estimates, perhaps due to distributed off-fault strain or the presence of multiple active fault strands. While our creep rate estimates are relatively insensitive to choice of model locking depth, insufficient independent information regarding locking depths is a source of epistemic uncertainty that impacts deep slip rate estimates.

California↗

USGS2014

Latest Quaternary paleoseismology and evidence of distributed dextral shear along the Mohawk Valley fault zone, northern Walker Lane, California

The dextral-slip Mohawk Valley fault zone (MVFZ) strikes northwestward along the eastern margin of the Sierra Nevada in the northern Walker Lane. Geodetic block modeling indicates that the MVFZ may accommodate ~3 mm/yr of regional dextral strain, implying that it is the highest slip-rate strike-slip fault in the region; however, only limited geologic data are available to constrain the system’s slip rate and earthquake history. We mapped the MVFZ using airborne lidar data and field observations and identified a site near Sulphur Creek for paleoseismic investigation. At this site, oblique dextral-normal faulting on the steep valley margin has created a closed depression that floods annually during spring snowmelt to form an ephemeral pond. We excavated three fault-perpendicular trenches at the site and exposed pond sediment that interfingers with multiple colluvial packages eroded from the scarp that bounds the eastern side of the pond. We documented evidence for four surface-rupturing earthquakes on this strand of the MVFZ. OxCal modeling of radiocarbon and luminescence ages indicates that these earthquakes occurred at 14.0 ka, 12.8 ka, 5.7 ka, and 1.9 ka. The mean ~4 kyr recurrence interval is inconsistent with slip rates of ~3 mm/yr; these rates imply surface ruptures of more than 10 m per event, which is geologically implausible for the subdued geomorphic expression and 60 km length of the MVFZ. We propose that unidentified structures not yet incorporated into geodetic models may accommodate significant dextral shear across the northern Walker Lane, highlighting the role of distributed deformation in this region.

California↗

USGS2014

Frequency-dependent effects of rupture for the 2004 Parkfield mainshock, results from UPSAR

The frequency-dependent effects of rupture propagation of the Parkfield, California earthquake (Sept. 28, 2004, M6) to the northwest along the San Andreas fault can be seen in acceleration records at UPSAR (USGS Parkfield Seismic Array) in at least two ways. First, we can see the effects of directivity in the acceleration traces at UPSAR, which is about 11.5 km from the epicenter. Directivity or the seismic equivalent of a Doppler shift has been documented in many cases by comparing short duration, high-amplitude pulses (P or S) in the forward direction with longer duration body waves in the backward direction. In this case we detect a change from a relatively large amplitude, coherent, high-frequency signal at the start of rupture to a low-amplitude, low-coherent, low-frequency signal at about the time the rupture front transfers from the forward azimuth to the back azimuth at about 34-36 s (time is UTC and are the seconds after day 272 and 17 hours and 15 minutes. S arrival is just after 30s) for rays leaving the fault and propagating to UPSAR. The frequency change is obvious in the band about 5 to 30 Hz, which is significantly above the corner frequency of the earthquake (about 0.11Hz). From kinematic source models, the duration of faulting is about 9.2 s and the change in frequency is during faulting as the rupture extends to the northwest. Understanding the systematic change in frequency and amplitude of seismic waves in relation to the propagation of the rupture front is important for predicting strong ground motion. Second, we can filter the acceleration records from the array to determine if the low frequency energy emerges from the same part of the fault as the high frequency signal (e.g. has the same back azimuth and apparent velocity at UPSAR) an important clue to the dynamics of rupture. Analysis of sources of strong motion (characterized by relatively high frequencies) compared to kinematic slip models (relatively low frequency) for the March 11, 2011 Tohoku earthquake as well as Maule (Feb. 27, 2010) and Chi-Chi (Sept. 20, 1999) earthquakes show that high- and low-frequency sources do not have the same locations on the fault. In this paper we filter the accelerograms from UPSAR for the 2004 mainshock in various passbands and then re-compute the cross correlations to determine the vector slowness of the incoming waves. At Parkfield, it appears that for seismic waves with frequencies above 1 Hz there is no discernible frequency-dependent difference in source position (up to 8 Hz) based on estimates of back azimuth and apparent velocity. However at lower frequencies, sources appear to be from shallower depths and trail the high frequencies as the rupture proceeds down the fault. This result is greater than one standard deviation of an estimate of error, based on a new method of estimating error that is a measure of how broad the peak in correlation is and an estimate of the variance of the correlation values. These observations can be understood in terms of a rupture front that is more energetic and coherent near the front of rupture (radiating higher frequencies) and less coherent and less energetic (radiating in a lower frequency band) behind the initial rupture front. This result is a qualitative assessment of changes in azimuth and apparent velocity with frequency and time and does not include corrections to find the source location on the fault.

California↗

USGS2014

High-frequency imaging of elastic contrast and contact area with implications for naturally observed changes in fault properties

During localized slip of a laboratory fault we simultaneously measure the contact area and the dynamic fault normal elastic stiffness. One objective is to determine conditions where stiffness may be used to infer changes in area of contact during sliding on nontransparent fault surfaces. Slip speeds between 0.01 and 10 µm/s and normal stresses between 1 and 2.5 MPa were imposed during velocity step, normal stress step, and slide-hold-slide tests. Stiffness and contact area have a linear interdependence during rate stepping tests and during the hold portion of slide-hold-slide tests. So long as linearity holds, measured fault stiffness can be used on nontransparent materials to infer changes in contact area. However, there are conditions where relations between contact area and stiffness are nonlinear and nonunique. A second objective is to make comparisons between the laboratory- and field-measured changes in fault properties. Time-dependent changes in fault zone normal stiffness made in stress relaxation tests imply postseismic wave speed changes on the order of 0.3% to 0.8% per year in the two or more years following an earthquake; these are smaller than postseismic increases seen within natural damage zones. Based on scaling of the experimental observations, natural postseismic fault normal contraction could be accommodated within a few decimeter wide fault core. Changes in the stiffness of laboratory shear zones exceed 10% per decade and might be detectable in the field postseismically.

Journal of Geophysical Research↗

USGS2014

Continuous uplift near the seaward edge of the Prince William Sound megathrust: Middleton Island, Alaska

Middleton Island, located at the seaward edge of the continental shelf 50 km from the base of the inner wall of the Aleutian Trench, affords an opportunity to make land-based measurements of uplift near the toe of the Prince William Sound megathrust, site of the 1964, M  = 9.2, Alaska earthquake. Leveling surveys (1973–1993) on Middleton Island indicate roughly uniform tilting (~1 µrad/a down to the northwest) of the island, and GPS surveys (1993–2012) show an uplift rate of 14 mm/a of the island relative to fixed North America. The data are consistent with a combined (coseismic and postseismic) uplift (in meters) due to the 1964 earthquake as a function of time τ (years after the earthquake) u ( τ ) = (3.5 + 1.21 log 10  [1 + 1.67  τ ]) H ( τ ) where 3.5 is the coseismic uplift and H ( τ ) is 0 for τ  < 0 and 1 otherwise. The current uplift on Middleton Island is attributed to continuous slip on a fault splaying off from the megathrust, and the long-term uplift is the superposition of the effects of past earthquakes, each earthquake being similar to the 1964 event. Then, the predicted uplift at time t due to a sequence of earthquakes at times t i would be . From studies of strandlines associated with the uplifted terraces on Middleton Island, Plafker et al. (1992) estimated the occurrence times of the last six earthquakes and measured the present-day elevations of those strandlines. The predicted uplift is in rough agreement with those measurements. About half of the predicted uplift is due to postseismic relaxation from previous earthquakes.

Alaska↗

USGS2014

Time-averaged discharge rate of subaerial lava at Kīlauea Volcano, Hawai‘i, measured from TanDEM-X interferometry: Implications for magma supply and storage during 2011-2013

Differencing digital elevation models (DEMs) derived from TerraSAR add-on for Digital Elevation Measurements (TanDEM-X) synthetic aperture radar imagery provides a measurement of elevation change over time. On the East Rift Zone (EZR) of Kīlauea Volcano, Hawai‘i, the effusion of lava causes changes in topography. When these elevation changes are summed over the area of an active lava flow, it is possible to quantify the volume of lava emplaced at the surface during the time spanned by the TanDEM-X data—a parameter that can be difficult to measure across the entirety of an ~100 km 2 lava flow field using ground-based techniques or optical remote sensing data. Based on the differences between multiple TanDEM-X-derived DEMs collected days to weeks apart, the mean dense-rock equivalent time-averaged discharge rate of lava at Kīlauea between mid-2011 and mid-2013 was approximately 2 m 3 /s, which is about half the long-term average rate over the course of Kīlauea's 1983–present ERZ eruption. This result implies that there was an increase in the proportion of lava stored versus erupted, a decrease in the rate of magma supply to the volcano, or some combination of both during this time period. In addition to constraining the time-averaged discharge rate of lava and the rates of magma supply and storage, topographic change maps derived from space-based TanDEM-X data provide insights into the four-dimensional evolution of Kīlauea's ERZ lava flow field. TanDEM-X data are a valuable complement to other space-, air-, and ground-based observations of eruptive activity at Kīlauea and offer great promise at locations around the world for aiding with monitoring not just volcanic eruptions but any hazardous activity that results in surface change, including landslides, floods, earthquakes, and other natural and anthropogenic processes.

Hawai'i↗

USGS2014

The 2010 slow slip event and secular motion at Kilauea, Hawai`i inferred from TerraSAR-X InSAR data

We present here an Small BAseline Subset (SBAS) algorithm to extract both transient and secular ground deformations on the order of millimeters in the presence of tropospheric noise on the order of centimeters, when the transient is of short duration and known time, and the background deformation is smooth in time. We applied this algorithm to study the 2010 slow slip event as well as the secular motion of Kīlauea's south flank using 49 TerraSAR-X images. We also estimate the tropospheric delay variation relative to a given reference pixel using an InSAR SBAS approach. We compare the InSAR SBAS solution for both ground deformation and tropospheric delays with existing GPS measurements and confirm that the ground deformation signal andtropospheric noise in InSAR data are successfully separated. We observe that the coastal region on the south side of the Hilina Pali moves at a higher background rate than the region north side of the Pali. We also conclude that the 2010 SSE displacement is mainly horizontal and the maximum magnitude of the 2010 SSE vertical component is less than 5 mm.

Hawaii↗