Subduction Zone Earthquakes Research Articles

Migration-Based Local Event-Location Workflow for Ocean-Bottom Seismometer (OBS) Records in Subduction Zones: A Practical Approach for Addressing a Large Number of Events

ABSTRACT An efficient event-location workflow is highly desired to analyze large numbers of local earthquakes recorded by ocean-bottom seismometers (OBSs) in subduction zones. The present study proposes a migration-based event-location approach for evaluating OBS records to examine local subduction-zone earthquakes. This approach can significantly reduce the amount of manual time picks compared with conventional methods. The event-location workflow was designed to detect arrival onsets of both P and S phases. Synthetic tests have shown that the proposed migration-based event-location method is robust against different types of noise, such as environmental noise and local spike noise. This workflow was then applied to real OBS data in the off-Ibaraki region at the southern end of the Japan trench. The results show that this approach is applicable to real data from subduction-zone events: It gives reasonable agreement with manual time picks for both P and S waves and reasonable error bars, and it demonstrates a clear down-dip trend of seismicity. The results also show fair agreement with event distributions from previous studies of the off-Ibaraki region. This proposed workflow can be used to examine the seismicity of local earthquakes around the subduction zone using OBSs. This approach is especially effective when the seismicity is high and/or in cases in which long-term OBS monitoring has recorded a large number of events.

New conditional, scenario-based, and non-conditional cumulative absolute velocity models for subduction tectonic settings

The PEER NGA-Sub ground-motion intensity measure database is used to develop new conditional ground-motion models (CGMMs), a set of scenario-based models, and non-conditional models to estimate the cumulative absolute velocity ([Formula: see text]) of ground motions from subduction zone earthquakes. In the CGMMs, the median estimate of [Formula: see text] is conditioned on the estimated peak ground acceleration ([Formula: see text]), the time-averaged shear-wave velocity in the top 30 m of the soil ([Formula: see text]), the earthquake magnitude ([Formula: see text]), and the spectral acceleration at the period of 1 s ([Formula: see text]). Multiple scenario-based [Formula: see text] models are developed by combining the CGMMs with pseudo-spectral acceleration ([Formula: see text]) ground-motion models (GMMs) for [Formula: see text] and [Formula: see text] to directly estimate [Formula: see text] given an earthquake scenario and site conditions. Scenario-based [Formula: see text] models are capable of capturing the complex ground-motion effects (e.g. soil non-linearity and regionalization effects) included in their underlying [Formula: see text]/[Formula: see text] GMMs. This approach also ensures the consistency of the [Formula: see text] estimates with a [Formula: see text] design spectrum. In addition, two non-conditional [Formula: see text] GMMs are developed using Bayesian hierarchical regressions. Finally, we present comparisons between the developed models. The comparisons show that if non-conditional GMMs are properly constrained, they are consistent with scenario-based GMMs. The [Formula: see text] GMMs developed in this study advance the performance-based earthquake engineering practice in areas affected by subduction zone earthquakes.

Influence of ground motion duration on the seismic vulnerability of aging highway bridges

Despite the potential of long duration earthquakes in subduction zones, little attention has been paid on seismic vulnerability assessment of highway bridges that are simultaneously undergoing prolonged corrosion deterioration. Consequently, this study presents a framework for considering the joint impact of ground motion duration and chloride-induced corrosion deterioration on the lifetime seismic vulnerability of aging highway bridge structures. The proposed framework is then applied to a representative case-study multi-span simply supported highway bridge located in subduction earthquake zone within Seattle, U.S.A. A detailed three-dimensional analytical models are first developed that account for the nonlinear behaviour of different bridge components, including experimentally validated column model capable of capturing cyclic degradation of force–deformation response. Suites of spectrally equivalent short and long duration ground motions are adopted followed by an incremental dynamic analysis-based framework for seismic fragility assessment. Deterministic seismic response assessment reveals a higher seismic demand on multiple bridge components under long duration ground motion record as compared to spectrally equivalent short duration record. Lastly, the bridge components and system fragility curves are developed considering modelling and deterioration uncertainties. Results reveal the significant vulnerability of bridge components and system under the joint impact of long duration earthquake shock and corrosion deterioration.

Keyboard Model of Seismic Cycle of Great Earthquakes in Subduction Zones: Simulation Results and Further Generalization

Catastrophic megaearthquakes (M > 8) occurring in the subduction zones are among the most devastating hazards on the planet. In this paper we discuss the seismic cycles of the megathrust earthquakes and propose a blockwise geomechanical model explaining certain features of the stress-deformation cycle revealed in recent decades from seismological and satellite geodesy (GNSS) observations. Starting with an overview of the so-called keyboard model of the seismic cycle by L. Lobkovsky, we outline mathematical formalism describing the motion of seismogenic block system assuming viscous rheology beneath and between the neighboring elastic blocks sitting on top of the subducting slab. By summarizing the GNSS-based evidence from our previous studies concerning the transient motions associated with the 2006–2007 Simushir earthquakes, 2010 Maule earthquake, and 2011 Tohoku earthquake, we demonstrate that those data support the keyboard model and reveal specific effect of the postseismic oceanward motion. However, since the seismogenic blocks in subduction systems are mostly located offshore, the direct analysis of GNSS-measured displacements and velocities is hardly possible in terms of the original keyboard model. Hence, the generalized two-segment keyboard model is introduced, containing both frontal offshore blocks and rear onshore blocks, which allows for direct interpretation of the onshore-collected GNSS data. We present a numerical computation scheme and a series of simulated data, which exhibits the consistency with measured motions and enables estimating the seismic cycle characteristics, important for the long-term earthquake forecasting.

Tsunami Source of the 2021 MW 8.1 Raoul Island Earthquake From DART and Tide‐Gauge Data Inversion

AbstractThe tsunami source of the 2021 MW 8.1 Raoul Island earthquake in the Kermadec subduction zone was estimated by inverting the tsunami signals recorded by Deep‐ocean Assessment and Reporting of Tsunamis (DART) bottom pressure sensors and coastal tide‐gauges. The main asperity of up to 5 m of slip is located northeastward from the hypocenter, with features compatible with the aftershock distribution and rapid back‐projection analysis. Three earthquakes of MW ∼8 or larger which also produced moderate tsunamis happened in the 20th century in the same portion of the subduction zone. This is the first great tsunamigenic event captured by the new New Zealand DART network in the South West Pacific, which proved valuable to estimate a robust image of the tsunami source. We also show a first proof of concept of the capability of this network to reduce the uncertainty associated with tsunami forecasting and to increase the lead time available for evacuation for future alerts.

Effects of the Georgia Sedimentary Basin on the Response of Modern Tall RC Shear-Wall Buildings to M9 Cascadia Subduction Zone Earthquakes

Tall residential RC shear wall buildings (RCSW), which are predominant in Metro Vancouver, have the potential to experience large magnitude earthquakes generated by the Cascadia Subduction Zone (CSZ). Furthermore, the region lies above the Georgia sedimentary basin, which can amplify the intensity of ground motions at medium to long periods and the resulting damage in tall structures. This study provides insights into the effects of the Georgia sedimentary basin amplification on (1) spectral accelerations associated with magnitude 9 (M9) CSZ earthquakes, (2) resulting force- and deformation-controlled actions in modern tall RCSW buildings, and (3) ensuing earthquake-induced repair costs and times. To this end, we leveraged a suite of physics-based ground motion simulations of a range of M9 CSZ earthquake scenarios, which explicitly consider basin effects, and benchmarked these scenarios against a range of seismic hazard intensities, which neglects basin effects. While the M9 simulations have an estimated 500-year return period, at deep basin sites their spectra exceed the 2,475-year hazard in the 1–3 s period range. Nonlinear dynamic analysis results under probabilistic seismic hazard estimates result in negligible collapse risk. In contrast, collapse risk conditioned on the occurrence of the M9 motions results in probabilities as high as 15%. Additionally, seismic demands from the M9 simulations at deep basin sites result in earthquake-induced repair costs and times that exceed those associated with the 2,475-year hazard level.

Spectral damping scaling factors for horizontal components of ground motions from subduction earthquakes using NGA-Subduction data

This article develops global models of damping scaling factors (DSFs) for subduction zone earthquakes that are functions of the damping ratio, spectral period, earthquake magnitude, and distance. The Next Generation Attenuation for subduction earthquakes (NGA-Sub) project has developed the largest uniformly processed database of recorded ground motions to date from seven subduction regions: Alaska, Cascadia, Central America and Mexico, South America, Japan, Taiwan, and New Zealand. NGA-Sub used this database to develop new ground motion models (GMMs) at a reference 5% damping ratio. We worked with the NGA-Sub project team to develop an extended database that includes pseudo-spectral accelerations (PSA) for 11 damping ratios between 0.5% and 30%. We use this database to develop parametric models of DSF for both interface and intraslab subduction earthquakes that can be used to adjust any subduction GMM from a reference 5% damping ratio to other damping ratios. The DSF is strongly influenced by the response spectral shape and the duration of motion; therefore, in addition to the damping ratio, the median DSF model uses spectral period, magnitude, and distance as surrogate predictor variables to capture the effects of the spectral shape and the duration of motion. We also develop parametric models for the standard deviation of DSF. The models presented in this article are for the RotD50 horizontal component of PSA and are compared with the models for shallow crustal earthquakes in active tectonic regions. Some noticeable differences arise from the considerably longer duration of interface records for very large magnitude events and the enriched high-frequency content of intraslab records, compared with shallow crustal earthquakes. Regional differences are discussed by comparing the proposed global models with the data from each subduction region along with recommendations on the applicability of the models.

The Earth's Surface Controls the Depth‐Dependent Seismic Radiation of Megathrust Earthquakes

AbstractMegathrust earthquakes exhibit a ubiquitous seismic radiation style: low‐frequency (LF) seismic energy is efficiently emitted from the shallowest portion of the fault, whereas high‐frequency (HF) seismic energy is efficiently emitted from the deepest part of the fault. Although this is observed in many case‐specific studies, we show that it is ubiquitous in global megathrust earthquakes between 1995 and 2021. Previous studies have interpreted this as an effect of systematic depth variation in either the plate interface frictional properties (Lay et al., 2012) or the P wavespeeds (Sallarès & Ranero, 2019). This work suggests an alternative hypothesis: the interaction between waves and ruptures due to the Earth's free surface is the leading mechanism that generates this behavior. Two‐dimensional dynamic rupture simulations of subduction zone earthquakes support this hypothesis. Our simulations show that the interaction between the seismic waves reflected at the Earth's free surface and the updip propagating rupture results in LF radiation at the source. In contrast, the downdip propagation of rupture is less affected by the free surface and is thus dominated by HF radiation typical of buried faults. To a second degree, the presence of a realistic Earth structure derived from P‐wave velocity (VP) tomographic images and realistic VP/VS ratio estimated in boreholes further enhances the contrast in source radiation. We conclude that the Earth's free surface is necessary to explain the observed megathrust earthquake radiation style, and the realistic structure of subduction zone is necessary to better predict earthquake ground motion and tsunami potential.

Collapse risk of pre-Northridge tall steel moment-resisting frames in the Seattle basin during large-magnitude subduction earthquakes

The Pacific Northwest has the potential to experience large-magnitude earthquakes generated by the Cascadia Subduction Zone, which is located approximately 100 km from the city of Seattle. Tall buildings in Seattle are particularly vulnerable to these earthquakes because the city lies above a deep sedimentary basin, which can amplify the intensity of earthquake ground motions at long periods. This paper evaluates the response of an archetype 1970s 50-story steel moment-resisting frame office building in Seattle under 30 simulated scenarios of a magnitude-9 (M9) Cascadia Subduction Zone earthquake with explicit consideration of basin effects. The archetype is assumed to have fracture-prone welded connections consistent with those identified following the 1994 Northridge earthquake. The probability of collapse conditioned on the occurrence of the M9 scenarios, which have a return period of approximately 500 years, is 30%. The annualized collapse risk of the archetype building is also assessed considering all sources (i.e., crustral, intraslab and interface earthquakes) that contribute to the seismic hazard through a multiple stripe analysis procedure. The results indicate a 50-year collapse risk of 6.9% when basin effects are neglected and 10.5% when basin effects are considered. These results exceed the 1% in 50-year target implicit in modern seismic design standards by a factor of 10. These high collapse risks are largely driven by (i) deep sedimentary basin effects and (ii) the expected brittle behavior of fracture-prone welded beam-to-column connections. The simulations of the performance of the building under the M9 scenarios outside of the basin or with ductile beam-to-column connections result in a negligible probability of collapse.

Assessing local impacts of the 1700 CE Cascadia earthquake and tsunami using tree-ring growth histories: a case study in South Beach, Oregon, USA

Abstract. We present an investigation of the disturbance history of an old-growth Douglas-fir (Pseudotsuga menziesii) stand in South Beach, Oregon, for possible growth changes due to tsunami inundation caused by the 1700 CE Cascadia Subduction Zone (CSZ) earthquake. A high-resolution model of the 1700 tsunami run-up heights at South Beach, assuming an “L”-sized earthquake, is also presented to better estimate the inundation levels several kilometers inland at the old-growth site. This tsunami model indicates the South Beach fir stand would have been subjected to local inundation depths from 0 to 10 m. Growth chronologies collected from the Douglas-fir stand shows that trees experienced a significant growth reductions in the year 1700 relative to nearby Douglas-fir stands, consistent with the tsunami inundation estimates. The ±1–3-year timing of the South Beach disturbances are also consistent with disturbances previously observed at a Washington state coastal forest ∼220 km to the north. Moreover, the 1700 South Beach growth reductions were not the largest over the >321-year tree chronology at this location, with other disturbances likely caused by climate drivers (e.g., drought or windstorms). Our study represents a first step in using tree growth history to ground truth tsunami inundation models by providing site-specific physical evidence.

An EPIC Tikhonov Regularization: Application to Quasi‐Static Fault Slip Inversion

AbstractImaging subsurface fault slip from surface observations is essential to improving our understanding of the physics of earthquakes and tsunamis. As the estimation of subsurface fault slip is inherently ill‐posed, common inversion methods usually require a regularization term to counteract instabilities. Such regularization introduces biases in inferred slip estimates. Here, we discuss the effects that prior information, implied by a given regularization scheme, can have on fault slip estimates. We propose a novel Equal Posterior Information Condition (EPIC)—based Tikhonov regularization that generalizes the concept of prior information. The EPIC determines variances of prior information based on a chosen form of the structure of the posterior covariance matrix. In the context of subduction zone earthquakes, use of the EPIC counterbalances the spatial heterogeneity of observational constraints on fault slip, improving stability, quality and interpretability of fault slip estimates. We investigate the efficiency of the EPIC in the context of various synthetic fault slip distributions. We also demonstrate the methodology by inferring co‐seismic slip from geodetic data for the 2011 (Mw 9.0) Tohoku‐Oki, Japan, earthquake, obtaining robust slip estimates that are similar to those inferred using a unregularized fully Bayesian approach.

A Significant Increase in Interplate Seismicity near Major Historical Earthquakes Offshore Martinique (FWI)

ABSTRACT Understanding the processes that may be at the origin of major earthquakes in subduction zones is highly challenging, especially in the case of slowly converging areas such as the Lesser Antilles subduction zone. Our study reveals a recorded increase in seismicity rate and cumulative seismic moment over the last two decades offshore Martinique island and, particularly, in the presumed rupture area of the major historical 1839 earthquake. This sustained seismicity is shared between extensive intermediate depth activity since the 2007 Mw 7.4 earthquake that occurred about 20 km north of Martinique and a compressive seismic cluster located in the seismogenic zone of the subduction zone. We also observe a downward migration of the seismicity along the plate interface up to a 60–65 km depth and a recorded increase in the magnitudes of the earthquakes. All these observations may indicate ongoing changes along the plate boundary interface near the area of the historical major rupture of 1839, although more detailed statistical studies and additional data are necessary to confirm this trend.

About Earthquakes in Subduction Zones with the Potential to Cause a Tsunami

The problem of occurrence of starting earthquakes in subduction zones is considered. Subduction is the phenomenon of movement of the oceanic lithospheric plate under the continental one. The oceanic lithospheric plate at a certain depth melts from below and can slide. The paper considers the occurrence of starting earthquakes under the assumption that lithospheric plates have different contact conditions, being on a rigid base in the subduction zone. A molten lithospheric plate has no tangential contact stresses, while the other, oceanic, is rigidly connected to the base. The block element method is used to study the occurrence of the starting earthquake and the peculiarity of its consequences. The conditions to generate of tsunamis as a result of such earthquakes are being studied. Solutions to boundary value problems that are constructed precisely, rather than approximatively, allow us to reveal the mechanisms of destruction of the environment that were not previously known. In particular, the results obtained allowed us to detect a new type of crack that was not previously described. They destroy the environment in a different way than Griffiths cracks, which is demonstrated in this paper and is important in engineering practice.

High Fluid‐Pressure Patches Beneath the Décollement: A Potential Source of Slow Earthquakes in the Nankai Trough off Cape Muroto

AbstractPore pressure plays a key role in the generation of earthquakes in subduction zones. However, quantitative constraints for its determination are quite limited. Here, we estimate the subsurface pore pressure by analyzing the transient upwelling flow of drilling mud from borehole C0023A of the International Ocean Discovery Program (IODP) Expedition 370, in the Nankai Trough off Cape Muroto. This upward flow provided the first direct evidence of an overpressured aquifer in the underthrust sediments off Cape Muroto. To estimate the pre‐drilling pore pressure in the overpressured aquifer around a depth of 950–1,050 m below sea floor, we examined the measured porosities of core samples retrieved from nearby IODP wells; we then proceeded to explain the observed time evolution of the flow rate of the upwelling flow by modeling various sized aquifers through solving a radial diffusion equation. It was observed that for a permeability of 10−13 m2, the aquifer possessed an initial excess pore pressure of ∼5–10 MPa above the hydrostatic pressure, with a lateral dimension of several hundred meters and thickness of several tens of meters. The overpressure estimates from the porosity‐depth profile at Site C0023 differ from those at other drill sites in the region, suggesting the possible existence of multiple overpressured aquifers with a patchy distribution in the underthrust sediments of the Nankai Trough. As pore pressure is relevant in maintaining fault stability, the overpressured aquifers may be the source of slow earthquakes that have been observed around the drilling site.

How cementation and fluid flow influence slip behavior at the subduction interface

Abstract Much of the complexity of subduction-zone earthquake size and temporal patterns owes to linkages among fluid flow, stress, and fault healing. To investigate these linkages, we introduce a novel numerical model that tracks cementation and fluid flow within the framework of an earthquake simulator. In the model, there are interseismic increases in cohesion across the plate boundary and decreases in porosity and permeability caused by cementation along the interface. Seismogenic slip is sensitive to the effective stress and therefore fluid pressure; in turn, slip events increase porosity by fracturing. The model therefore accounts for positive and negative feedbacks that modify slip behavior through the seismic cycle. The model produces temporal clustering of earthquakes in the seismic record of the Aleutian margin, which has well-documented along-strike variations in locking characteristics. Model results illustrate how physical, geochemical, and hydraulic linkages can affect natural slip behavior. Specifically, coseismic drops in fluid pressure steal energy from large ruptures, suppress slip, moderate the magnitudes of large earthquakes, and lead to aftershocks.

Seismicity Rate Change as a Tool to Investigate Delayed and Remote Triggering of the 2010–2011 Canterbury Earthquake Sequence, New Zealand

ABSTRACTCrustal earthquakes in low-strain-rate regions are rare in the human life span but can generate disastrous consequences when they occur. Such was the case in the Canterbury earthquake sequence that began in 2010 and eventually led to almost 200 fatalities. Our study explores this earthquake sequence’s origins by producing an enhanced earthquake catalog in the Canterbury Plains and Otago, South Island, New Zealand. We investigate seismicity rate changes from 2005 to before the 2010 Mw 7.2 Darfield earthquake. During this time, major subduction-zone earthquakes, such as the 2009 Mw 7.8 Dusky Sound earthquake, created measurable coseismic and postseismic strain in the region. We use template matching to expand the catalog of earthquakes in the region, and use a support vector machine classifier to remove false positives and poor detections. We then compare the newly obtained seismicity rates with the coseismic and postseismic crustal strain fields, and find that seismicity rate and crustal strain are positively correlated in the low-stress, low-seismicity region of the northern Canterbury Plains. In contrast, near fast-moving plate-boundary faults, the seismicity rate changes rise without much change in the strain rate. Our analysis reveals a substantial seismicity rate decrease in the western rupture area of the Darfield earthquake, which we infer to be an effect of coseismic and postseismic deformation caused by the Dusky Sound earthquake. We show in low-strain-rate regions, stress perturbation of a few kPas creates substantial seismicity rate change. However, the implication that such seismic quiescence is responsible for the nucleation of the Darfield earthquake requires further studies.

Deficit of Large Aftershocks as an Indicator of Afterslip at the Sources of Earthquakes in Subduction Zones

Deficit of Large Aftershocks as an Indicator of Afterslip at the Sources of Earthquakes in Subduction Zones

Ground-Motion Intensity Measure Correlations on Interface and Intraslab Subduction Zone Earthquakes Using the NGA-Sub Database

ABSTRACTThe NGA-Sub (subduction zone earthquake) database developed by the Pacific Earthquake Engineering Research center is used to derive new correlation coefficients for a number of ground-motion intensity measure (IM) parameters from ground motions in subduction zone earthquakes, considering both interface and intraslab tectonic settings. The IMs include peak ground acceleration, pseudospectral accelerations with periods from 0.01 to 10 s, Arias intensity, cumulative absolute velocity, peak ground velocity, and significant duration. Comparisons of the estimated correlation coefficients for ground motions from the interface and intraslab tectonic settings generally show a good agreement. Our estimations are also in good agreement with correlation coefficients estimated in previous studies that used ground motions from shallow crustal earthquakes, supporting the concept that any variation in correlation coefficients comes from spectral shape (i.e., the distribution of peaks and troughs) rather than tectonic region. This study also explores the influence of parameters such as magnitude, distance, and site conditions on the estimated correlation coefficients. We did not find apparent trends of the correlation coefficients with respect to these parameters. Finally, we propose analytical models to estimate correlation coefficients between the IMs explored in this study, considering both subduction interface and intraslab tectonic settings.

Subslab heterogeneity and giant megathrust earthquakes

The nucleation and rupture processes of giant megathrust earthquakes (M ≥ 9.0) in subduction zones are still controversial. Most previous studies have focused on the subducting plate interface, and the structure beneath the subducting slab and its influence on earthquake generation remain unclear. Here, we present high-resolution seismic velocity tomography beneath six regions where giant earthquakes have occurred. Subslab low-velocity (slow) anomalies are revealed, which may reflect hot mantle upwelling. The giant earthquake hypocentres are generally located above the edges of the slow anomalies or above the gaps between them. Large coseismic slips of the giant earthquakes mainly occurred above gaps between the slow anomalies. We suggest that differential buoyancy force between the slow anomalies and their gaps may be an important factor for earthquake nucleation, and the rupture extent of a giant earthquake may be constrained by the slow anomalies. Hence, it is necessary to conduct seismic tomography to investigate the detailed subslab structure, which may help to pinpoint the potential location and damage zone of a future giant earthquake. Mantle heterogeneity beneath subducting plates may influence giant megathrust earthquakes, according to seismic tomography of the subslab structure beneath six megathrusts that have ruptured in M ≥ 9.0 earthquakes.

The Preservation of Climate‐Driven Landslide Dams in Western Oregon

AbstractBedrock landsliding, including the formation of landslide dams, is a predominant geomorphic process in steep landscapes. Clarifying the importance of hydrologic and seismic mechanisms for triggering deep‐seated landslides remains an ongoing effort, and formulation of geomorphic metrics that predict dam preservation is crucial for quantifying secondary landslide hazards. Here, we identify >200 landslide‐dammed lakes in western Oregon and utilize dendrochronology and enhanced 14C dating (“wiggle matching”) of “ghost forests” to establish slope failure timing at 20 sites. Our dated landslide dataset reveals bedrock landsliding has been common since the last Cascadia Subduction Zone earthquake in January 1700 AD. Our study does not reveal landslides that date to 1700 AD. Rather, we observe temporal clustering of at least four landslides in the winter of 1889/1890 AD, coincident with a series of atmospheric rivers that generated one of the largest regionally recorded floods. We use topographic and field analyses to assess the relation between dam preservation and topographic characteristics of the impounded valleys. In contrast to previous studies, we do not observe systematic scaling between dam size and upstream drainage area, though dam stability indices for our sites correspond with “stable” dams elsewhere. Notably, we observe that dams are preferentially preserved at drainage areas of ∼1.5 to 13 km2 and valley widths of ∼25 to 80 m, which may reflect the reduced downstream influence of debris flows and the accumulation of mature conifer trees upstream from landslide‐dammed lake outlets. We suggest that wood accumulation upstream of landslide dams tempers large stream discharges, thus inhibiting dam incision.