Description: As Earth's atmospheric temperatures and human populations increase, more people are becoming vulnerable to natural and human-induced disasters. This is particularly true in Central America, where the growing human population is experiencing climate extremes (droughts and floods), and the region is susceptible to geological hazards, such as earthquakes and volcanic eruptions, and environmental deterioration in many forms (soil erosion, lake eutrophication, heavy metal contamination, etc.). Instrumental and historical data from the region are insufficient to understand and document past hazards, a necessary first step for mitigating future risks. Long, continuous, well-resolved geological records can, however, provide a window into past climate and environmental changes that can be used to better predict future conditions in the region. The Lake Izabal Basin (LIB), in eastern Guatemala, contains the longest known continental records of tectonics, climate, and environmental change in the northern Neotropics. The basin is a pull-apart depression that developed along the North American and Caribbean plate boundary ∼ 12 Myr ago and contains 〉 4 km of sediment. The sedimentological archive in the LIB records the interplay among several Earth System processes. Consequently, exploration of sediments in the basin can provide key information concerning: (1) tectonic deformation and earthquake history along the plate boundary; (2) the timing and causes of volcanism from the Central American Volcanic Arc; and (3) hydroclimatic, ecologic, and geomicrobiological responses to different climate and environmental states. To evaluate the LIB as a potential site for scientific drilling, 65 scientists from 13 countries and 33 institutions met in Antigua, Guatemala, in August 2022 under the auspices of the International Continental Scientific Drilling Program (ICDP) and the US National Science Foundation (NSF). Several working groups developed scientific questions and overarching hypotheses that could be addressed by drilling the LIB and identified optimal coring sites and instrumentation needed to achieve the project goals. The group also discussed logistical challenges and outreach opportunities. The project is not only an outstanding opportunity to improve our scientific understanding of seismotectonic, volcanic, paleoclimatic, paleoecologic, and paleobiologic processes that operate in the tropics of Central America, but it is also an opportunity to improve understanding of multiple geological hazards and communicate that knowledge to help increase the resilience of at-risk Central American communities.

Description: The potential for future earthquakes on faults is often inferred from inversions of geodetically derived surface velocities for locking on faults using kinematic models such as block models. This can be challenging in complex deforming zones with many closely spaced faults or where deformation is not readily described with block motions. Furthermore, surface strain rates are more directly related to coupling on faults than surface velocities. We present a methodology for estimating slip deficit rate directly from strain rate and apply it to New Zealand for the purpose of incorporating geodetic data in the 2022 revision of the New Zealand National Seismic Hazard Model. The strain rate inversions imply slightly higher slip deficit rates than the preferred geologic slip rates on sections of the major strike‐slip systems including the Alpine Fault, the Marlborough Fault System and the northern part of the North Island Fault System. Slip deficit rates are significantly lower than even the lowest geologic estimates on some strike‐slip faults in the southern North Island Fault System near Wellington. Over the entire plate boundary, geodetic slip deficit rates are systematically higher than geologic slip rates for faults slipping less than one mm/yr but lower on average for faults with slip rates between about 5 and 25 mm/yr. We show that 70%–80% of the total strain rate field can be attributed to elastic strain due to fault coupling. The remaining 20%–30% shows systematic spatial patterns of strain rate style that is often consistent with local geologic style of faulting. Plain Language Summary The potential for future earthquakes on faults is often inferred from velocities of the ground surface derived from satellite geodesy, but this approach can be challenging in complex deforming zones with many closely spaced faults. We present a new methodology for estimating the rate at which energy is accumulating on faults using measurements of surface strain rates. The method is applied to New Zealand for the purpose of incorporating geodetic data in the 2022 revision of the New Zealand National Seismic Hazard Model. We show that 70%–80% of the total deformation field can be attributed to energy accumulation on known active faults while the source of the remaining 20%–30% remains unknown. Along some of the major faults in New Zealand we find some important differences in rates of energy accumulation from what is expected from geologic data. Estimated rates are significantly lower than even the lowest geologic estimates on some faults in the fault system near highly‐populated Wellington. Key Points We develop a method to invert geodetically derived strain rates for slip deficit rates on faults We find small but systematic differences between slip deficit rates and geologic slip rates About 70%–80% of the surface strain can be attributed to elastic strain due to coupling on faults

Description: Geodetic data in plate boundary zones reflect the accrual of tectonic strain and stress, which will ultimately be released in earthquakes, and so they can provide valuable insights into future seismic hazards. To incorporate geodetic measurements of contemporary deformation into the 2022 revision of the New Zealand National Seismic Hazard Model 2022 (NZ NSHM 2022), we derive a range of strain-rate models from published interseismic Global Navigation Satellite Systems velocities for New Zealand. We calculate the uncertainty in strain rate excluding strain from the Taupō rift–Havre trough and Hikurangi subduction zone, which are handled separately, and the corresponding moment rates. A high shear strain rate occurs along the Alpine fault and the North Island dextral fault belt, as well as the eastern coast of the North Island. Dilatation rates are primarily contractional in the South Island and less well constrained in the North Island. Total moment accumulation derived using Kostrov-type summation varies from 0.64 to 2.93×1019 N·m/yr depending on method and parameter choices. To account for both aleatory and epistemic uncertainty in the strain-rate results, we use four different methods for estimating strain rate and calculate various average models and uncertainty metrics. The maximum shear strain rate is similar across all methods, whereas the dilatation rate and overall strain rate style differ more significantly. Each method provides an estimate of its own uncertainty propagated from the data uncertainties, and variability between methods provides an additional estimate of epistemic uncertainty. Epistemic uncertainty in New Zealand tends to be higher than the aleatory uncertainty estimates provided by any single method, and epistemic uncertainty on dilatation rate exceeds the aleatory uncertainty nearly everywhere. These strain-rate models were provided to the NZ NSHM 2022 team and used to develop fault-slip deficit rate models and scaled seismicity rate models.