Disturbances in ionospheric Total Electron Content (dTEC) with frequencies of 1–100 mHz can be driven from above by processes in the magnetosphere and below by processes on the Earth's surface and lower atmosphere. Past studies showed the potential of dTEC as a diagnostic of magnetospheric Ultra Low Frequency (ULF) wave activity and demonstrated that ULF dTEC can impact space weather by, for example, changing ionospheric conductance. However, most past work has focused on single event studies, lacked magnetospheric context, or used sampling rates too low to capture most ULF waves. Here, we perform a statistical study using Time History of Events and Macrsoscale Interactions during Substorms (THEMIS) satellite conjunctions with a ground‐based magnetometer and Global Navigation Satellite System (GNSS) receiver at 65° magnetic latitude. We find that magnetospheric ULF waves generate dTEC variations across the broad range of frequencies examined in this study (2–50 mHz), and that ULF dTEC wave power is correlated with Kp, AE, solar wind speed, and magnetic field wave power observed in the magnetosphere and on the ground. We further find that magnetospheric ULF waves generate dTEC amplitudes up to TECU ( background), with the largest amplitudes occurring during geomagnetically active conditions, at frequencies below 7 mHz, and at local times near midnight. We finally discuss the implications of our results for magnetosphere‐ionosphere coupling and remote sensing techniques related to ULF waves.

The Jet Propulsion Laboratory (JPL) develops JPL-GIM, a software for generating global ionospheric maps (GIMs) of total electron content (TEC) using measurements from multiple Global Navigation Satellite System (GNSS) constellations. Within this overview paper, we delve into the current state and the most recent enhancements of JPL-GIM. An adaptable Kalman filter provides maps with user-defined temporal and spatial resolutions, while concurrently delivering essential covariance data for uncertainty assessment. Furthermore, multiple shell models offer a versatile framework to balance accuracy and computational efficiency. We present the five operational JPL GIM products (JPLG, JPRG, JPLI, JPLD, JPRT), highlighting JPLG and JPRG, our products routinely delivered to the International GNSS Service (IGS), and introduce a new near-real-time product (JPRT). As an added demonstration of JPL-GIM’s capabilities, we present a very-high-resolution (2-minute, multi-GNSS, 1000-station) configuration to showcase JPL-GIM’s ability to resolve long-lasting effects of the 2022 Hunga Tonga-Hunga Ha’apai eruption. Validations using independent datasets confirm the accurate reproduction of ionospheric variations across all latitudinal bands.

Global Navigation Satellite System radio occultation (GNSS-RO) and microwave radiometry (MWR) are two of the most impactful spaceborne remote sensing techniques for numerical weather prediction (NWP). These two techniques provide complementary information about atmospheric temperature and water vapor structure. GNSS-RO provides high vertical resolution measurements with cloud penetration capability, but the temperature and moisture are coupled in the GNSS-RO retrieval process and their separation requires the use of a priori information or auxiliary observations. On the other hand, the MWR measures brightness temperature (Tb) in numerous frequency bands related to the temperature and water vapor structure but is limited by poor vertical resolution (> 2 km) and precipitation. In this study, we combine these two technologies in an optimal estimation approach, 1D variation method (1DVar), to improve the characterization of the complex thermodynamic structures in the lower troposphere. This study employs both simulated and operational observations. GNSS-RO bending angle and MWR Tb observations are used as inputs to the joint retrieval, where bending can be modeled by an Abel integral and Tb can be modeled by a radiative transfer model (RTM) that takes into account atmospheric absorption, as well as surface reflection and emission. By incorporating the forward operators into the 1DVar method, the strength of both techniques can be combined to bridge individual weaknesses. Applying 1DVar to the data simulated from large eddy simulation (LES) is shown to reduce GNSS-RO temperature and water vapor retrieval biases at the lower troposphere while simultaneously capturing the fine-scale variability that MWR cannot resolve. A sensitivity analysis is also conducted to quantify the impact of the a priori information and error covariance used in different retrieval scenarios. The applicability of 1DVar joint retrieval to the actual GNSS-RO and MWR observations is also demonstrated through combining collocated COSMIC-2 and Suomi-NPP (National Polar-orbiting Partnership) measurements.

We analyzed Synthetic Aperture Radar (SAR) images from Copernicus Sentinel-1A and 1B satellites operated by the European Space Agency and the Advanced Land Observation Satellite-2 (ALOS-2) satellite operated by the Japan Aerospace Exploration Agency and Global Navigation Satellite System (GNSS) data from the Network of the Americas for the 4 July 2019 Mw 6.4 and 5 July (local; 6 July UTC) Mw 7.1 Ridgecrest earthquakes. We integrated geodetic measurements for the 3D vector field of coseismic surface deformation for the two events, using SAR data from Sentinel-1 and ALOS-2 satellites. We combined less precise large-scale displacements from SAR images by pixel offset tracking or matching, including the along-track component, with the more precise SAR interferometry (Interferometric Synthetic Aperture Radar [InSAR]) measurements in the radar line of sight (LoS) direction and intermediate-precision along-track InSAR to estimate all three components of the surface displacement for the two events together. We also estimated the coseismic deformation for the two earthquakes from time-series processing of continuous Global Navigation Satellite System data stations in the area. InSAR coherence and coherence change maps the surface disruptions due to fault ruptures reaching the surface. Large slip in the Mw 6.4 earthquake was on a NE-striking fault that intersects with the NW-striking fault that was the main rupture in the Mw 7.1 earthquake. The main fault bifurcates towards the southeast ending 3 km from the Garlock Fault. The Garlock fault had triggered slip of about 20 mm in the radar LoS along a short section directly south of the main rupture. About 3 km northwest of the Mw 7.1 epicenter, the surface fault separates into two strands that form a pull-apart with about 1 m of down-drop. Further northwest is a wide zone of complex deformation.

Abstract This paper describes the methods used to estimate positions, velocities, breaks, and seasonal terms from daily Global Navigation Satellite System (GNSS) measurements. Break detection and outlier removal have been automated so that decades of daily measurements from thousands of stations can be processed in a few hours. New measurements are added, and parameters are updated every week. Model parameters allow separation of interseismic, annual, coseismic, and postseismic signals. Tools available through GeoGateway (http://geo-gateway.org) allow rapid visualization and analysis of these terms for results that can be subsetted in time or space. Results show highly variable and nonlinear motion for GPS stations in southern California. The variable motion is related to seasonal motions, distributed tectonic motion, earthquakes, and postseismic motions that can continue for years. In some areas results suggest that additional processes are responsible for the observed motions. In general, following earthquakes, stations return to their long‐term motions after 2–3 years, though some exceptions occur. The use of the tools shows nonlinear motion in the Salton Trough of southern California related to the 2010 M7.2 El Mayor‐Cucapah earthquake, 2012 Brawley earthquake swarm, and a creep event on the Superstition Hills fault in 2017.

The Global Positioning System (GPS) has revolutionized the ability to monitor Earth-system processes, including Earth’s water cycle. Several analysis centers process GPS data to estimate ground-antenna positions at daily temporal resolution. Differences in processing strategies can lead to inconsistencies in coordinate-position estimates and therefore influence the analysis of crustal displacement associated with variations in atmospheric and hydrologic mass loading. Here, we compare five GPS data products produced by three processing centers: the Nevada Geodetic Laboratory, Jet Propulsion Laboratory, and UNAVCO Consortium. We find that 5 to 30% of the scatter in residual GPS time series (commonly considered noise) can be explained by atmospheric loading in the contiguous USA and Alaska, but that the percentages vary widely by data product. Positions derived using high-resolution troposphere models (e.g., ECMWF) exhibit significantly lower scatter after correcting for atmospheric loading than positions estimated using constant or slowly varying troposphere models (e.g., GPT2w). The data products also exhibit differences in seasonal deformation (commonly attributed, in large part, to fluctuations in hydrologic mass loading): median vector differences in estimated seasonal amplitude range from 0.4–1.0 mm in the vertical component and 0.1–0.3 mm in the horizontal components, or about 10–40% of the mean amplitudes of seasonal oscillation. Newer products exhibit lower total scatter and stronger correlations than older products. Network-coherent differences in estimates of seasonal deformation reveal reference-frame inconsistencies between data products. We also cross-check two independent models of atmospheric pressure loading: ESMGFZ and LoadDef.

GipsyX/RTGx is the Jet Propulsion Laboratory?s (JPL) next generation software package for positioning, navigation, timing, and Earth science using measurements from three geodetic techniques: Global Navigation Satellite Systems (GNSS), Satellite Laser Ranging (SLR), and Doppler Orbitography and Radiopositioning Integrated by Satellite (DORIS); with Very Long Baseline Interferometry (VLBI) under development. The software facilitates combined estimation of geodetic and geophysical parameters using a Kalman filter approach on real or simulated data in both post -processing and in real-time. The estimated parameters include station coordinates and velocities, satellite orbits and clocks, Earth orientation, ionospheric and tropospheric delays. The software is also capable of full realiza- tion of a dynamic terrestrial reference through analysis and combination of time series of ground station coordinates. Applying lessons learned from its predecessors, GIPSY -OASIS and Real Time GIPSY (RTG), GipsyX/RTGx was re -designed from the ground up to offer improved precision, accuracy, usability, and operational flexibility. We present some key aspects of its new archi- tecture, and describe some of its major applications, including Real-time orbit determination and ephemeris predictions in the U.S. Air Force Next Generation GPS Operational Control Segment (OCX), as well as in JPL?s Global Differential GPS (GDGPS) System, sup- porting User Range Error (URE) of <5 cm RMS; precision post -processing GNSS orbit determination, including JPL?s contributions to the International GNSS Service (IGS) with URE in the 2 cm RMS range; Precise point positioning (PPP) with ambiguity resolution, both statically and kinematically, for geodetic applications with 2 mm horizontal, and 6.5 mm vertical repeatability for static positioning; Operational orbit and clock determination for Low Earth Orbiting (LEO) satellites, such as NASA?s Gravity Recovery and Climate Experiment (GRACE) mission with GRACE relative clock alignment at the 20 ps level; calibration of radio occultation data from LEO satellites for weather forecasting and climate studies; Satellite Laser Ranging (SLR) to GNSS and LEO satellites, DORIS-based and multi -technique orbit determination for LEO; production of terrestrial reference frames and Earth rotation parameters in support of JPL?s contribution to the International Terrestrial Reference Frame (ITRF). (c) 2020 COSPAR. Published by Elsevier Ltd. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/).

Large earthquakes in subduction zones result from tectonic strain that has accumulated over decades to centuries along coupled regions of the plate interface. Additional elastic energy can also be released through aseismic slip during the post-seismic and inter-seismic phase. Assuming a fully elastic crust (ie that cannot accumulate permanent deformation), measuring and comparing strain increase and release often leads to a balanced budget. However, we know that part of the tectonic contraction we evidence with geodesy translates into permanent deformation resulting in mountain building or vertical motion of the coastal domain in subduction zones over millions of years. Although increments of permanent strain should be small compared to interseismic loading strain, we do not know yet how much of this contraction is stored permanently and how this deformation is partitioned between the different phases of …

A nearly 20-year hiatus in major seismic activity in southern California ended on 4 July 2019 with a sequence of intersecting earthquakes near the city of Ridgecrest, California. This sequence included a foreshock with a moment magnitude (M w) of 6.4 followed by a M w 7.1 mainshock nearly 34 hours later. Geodetic, seismic, and seismicity data provided an integrative view of this sequence, which ruptured an unmapped multiscale network of interlaced orthogonal faults. This complex fault geometry persists over the entire seismogenic depth range. The rupture of the mainshock terminated only a few kilometers from the major regional Garlock fault, triggering shallow creep and a substantial earthquake swarm. The repeated occurrence of multifault ruptures, as revealed by modern instrumentation and analysis techniques, poses a formidable challenge in quantifying regional seismic hazards.

Uninhabited aerial vehicle synthetic aperture radar (UAVSAR) observations 2009–2017 of the Yuha Desert area and Global Positioning System (GPS) time series encompassing the region reveal a northward migrating pattern of deformation following the 4 April 2010 Mw7.2 El Mayor-Cucapah (EMC) earthquake. The north end of the EMC rupture exhibits an asymmetric pattern of deformation that is substantial and smooth northeast of the rupture and limited but with surface fracturing slip northwest. The earthquake triggered ~1 cm of surface coseismic slip at the Yuha fault, which continued to slip postseismically. 2.5 cm of Yuha fault slip occurred by the time of the 15 June 2010 Mw5.7 Ocotillo aftershock and 5 cm of slip occurred by 2017 following a logarithmic afterslip decay 16-day timescale. The Ocotillo aftershock triggered 1.4 cm of slip on a northwest trend extending to the Elsinore fault and by 7 years after the EMC earthquake 2.4 cm of slip had accumulated with a distribution following an afterslip function with a 16-day timescale consistent with other earthquakes and a rate strengthening upper crustal sedimentary layer. GPS data show broad coseismic uplift of the Salton Trough and delayed postseismic motion that may be indicative of fluid migration there and subsidence west of the rupture extension, which continues following the earthquake. The data indicate that the Elsinore, Laguna Salada, and EMC ruptures are part of the same fault system. The results also suggest that north-south shortening and east-west extension across the region drove fracture advancing step tectonics north of the EMC earthquake rupture. Plain Language Summary Airborne radar and Global Positioning System observations of the Yuha Desert and Salton Trough suggest that the 2010 M7.2 El Mayor Cucapah earthquake rupture, the Laguna Salada fault in Baja California, Mexico, and the Elsinore fault in California are part of the same fault system. The M5.7 Ocotillo aftershock links these faults and is a result of stress propagation from afterslip on the El Mayor Cucapah rupture. Further east in the Salton Trough fluid migration also contributes to the observed results.

On 26 August 2017, Hurricane Harvey struck the Gulf Coast as a category four cyclone depositing ~95 km3 of water, making it the wettest cyclone in U.S.HISTORY:Water left in Harvey's wake should cause elastic loading and subsidence of Earth's crust, and uplift as it drains into the ocean and evaporates. To track daily changes of transient water storage, we use Global Positioning System (GPS) measurements, finding a clear migration of subsidence (up to 21 mm) and horizontal motion (up to 4 mm) across the Gulf Coast, followed by gradual uplift over a 5-week period. Inversion of these data shows that a third of Harvey's total stormwater was captured on land (25.7 ± 3.0 km3), indicating that the rest drained rapidly into the ocean at a rate of 8.2 km3/day, with the remaining stored water gradually lost over the following 5 weeks at ~1 km3/day, primarily by evapotranspiration. These results indicate that GPS networks can remotely track the spatial extent and daily evolution of terrestrial water storage following transient, extreme precipitation events, with implications for improving operational flood forecasts and understanding the response of drainage systems to large influxes of water.

Uninhabited aerial vehicle synthetic aperture radar (UAVSAR) observations 2009–2017 of the Yuha Desert area and Global Positioning System (GPS) time series encompassing the region reveal a northward migrating pattern of deformation following the 4 April 2010 M w 7.2 El Mayor‐Cucapah (EMC) earthquake. The north end of the EMC rupture exhibits an asymmetric pattern of deformation that is substantial and smooth northeast of the rupture and limited but with surface fracturing slip northwest. The earthquake triggered ~1 cm of surface coseismic slip at the Yuha fault, which continued to slip postseismically. 2.5 cm of Yuha fault slip occurred by the time of the 15 June 2010 M w 5.7 Ocotillo aftershock and 5 cm of slip occurred by 2017 following a logarithmic afterslip decay 16‐day timescale. The Ocotillo aftershock triggered 1.4 cm of slip on a northwest trend extending to the Elsinore fault and by 7 years after the EMC earthquake 2.4 cm of slip had accumulated with a distribution following an afterslip function with a 16‐day timescale consistent with other earthquakes and a rate strengthening upper crustal sedimentary layer. GPS data show broad coseismic uplift of the Salton Trough and delayed postseismic motion that may be indicative of fluid migration there and subsidence west of the rupture extension, which continues following the earthquake. The data indicate that the Elsinore, Laguna Salada, and EMC ruptures are part of the same fault system. The results also suggest that north‐south shortening and east‐west extension across the region drove fracture advancing step tectonics north of the EMC earthquake rupture.

On April 25th 2015, the Mw 7.8 Gorkha (Nepal) earthquake ruptured a portion of the Main Himalayan Thrust underlying Kathmandu and surrounding regions. We develop kinematic slip models of the Gorkha earthquake using both a regularized multi-time-window (MTW) approach and an unsmoothed Bayesian formulation, constrained by static and high rate GPS observations, synthetic aperture radar (SAR) offset images, interferometric SAR (InSAR), and teleseismic body wave records. These models indicate that Kathmandu is located near the updip limit of fault slip and approximately 20km south of the centroid of fault slip. Fault slip propagated unilaterally along-strike in an ESE direction for approximately 140km with a 60km cross-strike extent. The deeper portions of the fault are characterized by a larger ratio of high frequency (0.03–0.2Hz) to low frequency slip than the shallower portions. From both the MTW and Bayesian results, we can resolve depth variations in slip characteristics, with higher slip roughness, higher rupture velocity, longer rise time and higher complexity of subfault source time functions in the deeper extents of the rupture. The depth varying nature of rupture characteristics suggests that the up-dip portions are characterized by relatively continuous rupture, while the down-dip portions may be better characterized by a cascaded rupture. The rupture behavior and the tectonic setting indicate that the earthquake may have ruptured both fully seismically locked and a deeper transitional portions of the collision interface, analogous to what has been seen in major subduction zone earthquakes.

Drought struck California during 7 of the 9years from 2007 to 2015, reducing the state's available water resources. Pumping of Central Valley groundwater has produced spectacular land subsidence. Uplift of the adjacent Sierra Nevada mountains has been proposed to be either tectonic uplift or solid Earth's elastic response to unloading of Central Valley groundwater. We find that of the 24mm of uplift of the Sierra Nevada from October 2011 to October 2015, just 5mm is produced by Central Valley groundwater loss, less than 2mm is tectonic uplift, and 17mm is solid Earth's elastic response to water loss in the Sierra Nevada. We invert GPS vertical displacements recording solid Earth's elastic response to infer changes in water storage across the western U.S. from January 2006 to October 2017. We find water changes to be sustained over periods of drought or heavy precipitation: the Sierra Nevada lost 1519km(3) of water during drought from October 2006 to October 2009, gained 1814km(3) of water during heavy precipitation from October 2009 to October 2011, and lost 4521km(3) of water during severe drought from October 2011 to October 2015 (95% confidence limits). Such large changes are not in hydrology models: snow accumulation in October is negligible and long-term soil moisture change is small. We infer that there must be large loss of either deep soil moisture or groundwater in river alluvium and in crystalline basement in the Sierra Nevada. The results suggest there to be parching of water in the ground during the summer of years of drought and seeping of melting snow into the Sierra Nevada in the spring of years of heavy precipitation.

The 2016 Kumamoto earthquake sequence occurred on the Futagawa–Hinagu fault zone near the Aso volcano on Kyushu island. The sequence was initiated with two major (Mw ≥ 6.0) foreshocks, and the mainshock (Mw = 7.0) occurred 25 h after the second major foreshock. We combine GPS, strong motion, synthetic aperture radar images, and surface offset data in a joint inversion to resolve the kinematic rupture process of the mainshock and coseismic displacement of the foreshocks. The joint inversion results reveal a unilateral rupture process for the mainshock involving sequential rupture of four major asperities. The slip area of the foreshocks and mainshock and the aftershock loci form a detailed complementary pattern. The mainshock rupture terminates near the rim of the caldera, leaving a ~10 km long gap of aftershocks. This area is characterized by high temperature and low shear wave velocity, density, and resistivity, which may be related to the partially melted geothermal condition. Ductile material property near the volcano may act as a “material barrier” to the dynamic rupture. Topographic weight of the caldera increases compressional normal stress on the fault plane, which may behave as a “stress barrier.” Long‐term seismic hazard and deformation behaviors related to these two types of barriers are discussed in terms of the associated frictional mechanism. Significant postseismic creeps observed near the volcano area indicates a velocity strengthening frictional behavior near the rupture termination, which confirms that the “material barrier” mechanism is likely the dominant rupture termination mechanism.

The April 25, 2015 M7.8 Gorkha earthquake caused more than 8,000 fatalities and widespread building damage in central Nepal. Four days after the earthquake, the Italian Space Agency's (ASI's) COSMO-SkyMed Synthetic Aperture Radar (SAR) satellite acquired data over Kathmandu area. Nine days after the earthquake, the Japan Aerospace Exploration Agency's (JAXA's) ALOS-2 SAR satellite covered larger area. Using these radar observations, we rapidly produced damage proxy maps derived from temporal changes in Interferometric SAR (InSAR) coherence. These maps were qualitatively validated through comparison with independent damage analyses by National Geospatial-Intelligence Agency (NGA) and the UNITAR's (United Nations Institute for Training and Research's) Operational Satellite Applications Programme (UNOSAT), and based on our own visual inspection of DigitalGlobe's WorldView optical pre- vs. post-event imagery. Our maps were quickly released to responding agencies and the public, and used for damage assessment, determining inspection/imaging priorities, and reconnaissance fieldwork.

Continuous GPS monitoring on global and regional scales has become an essential component of geophysical and meteorological infrastructure for studying fundamental Earth processes that drive natural hazards, weather, and climate. The NASA-funded “Solid Earth Science ESDR System (SESES)” project provides long-term Earth Science Data Records (ESDRs), the result of a combined solution of independent GPS analyses by the Jet Propulsion Laboratory and Scripps Institution of Oceanography using a common source of metadata archived at the Scripps Orbit and Permanent Array Center. The project has now produced up to twenty-two years of consistent, calibrated and validated ESDR products for over 3,200 GPS stations in western North America, other plate boundaries, and global networks. We describe the methodology to estimate a single set of time series with 24-h resolution of station displacements in north, east and vertical components. This is followed by a time series analysis for velocities, coseismic offsets, postseismic deformation, seasonal signals and nuisance offsets, primarily due to GPS antenna changes. Realistic one-sigma velocity are on the order of 0.03–0.05 mm/year in horizontal components and 0.1–0.3 mm/year in the vertical based on time series of 10–20 year duration. We present examples of time series that are well modeled by this parameterization and of time series that exhibit residual transient motions exhibiting episodic tremor and slip (ETS) processes. The project also catalogs seismic displacement and velocity waveforms estimated for a set of historical earthquakes in Japan and the U.S. through a seismogeodetic combination of GPS and collocated strong-motion accelerometers.

The subduction zone in northern Chile is a well‐identified seismic gap that last ruptured in 1877. On 1 April 2014, this region was struck by a large earthquake following a two week long series of foreshocks. This study combines a wide range of observations, including geodetic, tsunami, and seismic data, to produce a reliable kinematic slip model of the M w =8.1 main shock and a static slip model of the M w =7.7 aftershock. We use a novel Bayesian modeling approach that accounts for uncertainty in the Green's functions, both static and dynamic, while avoiding nonphysical regularization. The results reveal a sharp slip zone, more compact than previously thought, located downdip of the foreshock sequence and updip of high‐frequency sources inferred by back‐projection analysis. Both the main shock and the M w =7.7 aftershock did not rupture to the trench and left most of the seismic gap unbroken, leaving the possibility of a future large earthquake in the region.

Detailed geodetic imaging of earthquake ruptures enhances our understanding of earthquake physics and associated ground shaking. The 25 April 2015 moment magnitude 7.8 earthquake in Gorkha, Nepal was the first large continental megathrust rupture to have occurred beneath a high-rate (5-hertz) Global Positioning System (GPS) network. We used GPS and interferometric synthetic aperture radar data to model the earthquake rupture as a slip pulse ~20 kilometers in width, ~6 seconds in duration, and with a peak sliding velocity of 1.1 meters per second, which propagated toward the Kathmandu basin at ~3.3 kilometers per second over ~140 kilometers. The smooth slip onset, indicating a large (~5-meter) slip-weakening distance, caused moderate ground shaking at high frequencies (>1 hertz; peak ground acceleration, ~16% of Earth's gravity) and minimized damage to vernacular dwellings. Whole-basin resonance at a period of 4 to 5 seconds caused the collapse of tall structures, including cultural artifacts.