USGS Science⌕ Search

SEARCH · USGS Science

Results for “Earth Observer”

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

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

At least 1,171 records · Page 65Linked to original sources

Paleoseismology along the 1980 surface rupture of the Irpinia Fault: Implications for earthquake recurrence in the southern Apennines, Italy

The Irpinia fault was the source of the M s 6.9 1980 Irpinia earthquake and produced the first unequivocal historical surface faulting in Italy. Trenching of the 1980 fault scarp at Piano di Pecore, a flat intermontane basin about 5 km south of the 1980 instrumental epicenter, provides the first data on earthquake recurrence intervals, slip per event, and slip rate on a major normal fault in the Southern Apennines fault zone. The trenches exposed evidence of four pre-1980 paleoearthquakes that occurred during the past 8600 years. A best estimate average recurrence interval is 2150 years, although the time interval between individual events varies by as much as a factor of 2. Each paleo earthquake is similar to the 1980 surface rupture in amount of slip and style of deformation, which suggests that the 1980 event is characteristic for the Irpinia fault. Slip per event values average 61 cm. The net vertical displacement of 2.12–2.36 m since 8600 cal year B.P. observed in the trenches gives a vertical slip rate of 0.25–0.35 mm/yr, a dip slip rate of 0.29–0.40 mm/yr, and an extension rate of 0.14–0.20 mm/yr. Although fault behavior data are only available for the Irpinia fault they provide a starting point for evaluating earthquake recurrence and rates of deformation in southern Apennines. They suggest that (1) fault specific earthquake recurrence intervals based on the historical seismic record overestimates the occurrence of large magnitude ( M 7) earthquakes and (2) the Holocene rate of extension across the Apennines is ≤1 mm/yr. The 1980 earthquake and the paleoseismologic observations show that repeated and localized surface faulting occurs in southern Apennines and leaves subtle but distinct geomorphic evidence that can be detected with detailed and careful investigation.

Apennine Mountains↗

Observations on basaltic lava streams in tubes from Kilauea Volcano, island of Hawai'i

From 1986 to 1997, the Pu'u 'O'o-Kupaianaha eruption of Kilauea produced a vast pahoehoe flow field fed by lava tubes that extended 10–12 km from vents on the volcano's east rift zone to the ocean. Within a kilometer of the vent, tubes were as much as 20 m high and 10–25 m wide. On steep slopes (4–10°) a little farther away from the vent, some tubes formed by roofing over of lava channels. Lava streams were typically 1–2 m deep flowing within a tube that here was typically 5 m high and 3 m wide. On the coastal plain (<1°), tubes within inflated sheet flows were completely filled, typically 1–2 m high, and several tens of meters wide. Tubes develop as a flow's crust grows on the top, bottom, and sides of the tubes, restricting the size of the fluid core. The tubes start out with nearly elliptical cross-sectional shapes, many times wider than high. Broad, flat sheet flows evolve into elongate tumuli with an axial crack as the flanks of the original flow were progressively buried by breakouts. Temperature measurements and the presence of stalactites in active tubes confirmed that the tube walls were above the solidus and subject to melting. Sometimes, the tubes began downcutting. Progressive downcutting was frequently observed through skylights; a rate of 10 cm/d was measured at one skylight for nearly 2 months.

Journal of Geophysical Research B: Solid Earth↗

The spatial distribution of earthquake stress rotations following large subduction zone earthquakes

Rotations of the principal stress axes due to great subduction zone earthquakes have been used to infer low differential stress and near-complete stress drop. The spatial distribution of coseismic and postseismic stress rotation as a function of depth and along-strike distance is explored for three recent M ≥ 8.8 subduction megathrust earthquakes. In the down-dip direction, the largest coseismic stress rotations are found just above the Moho depth of the overriding plate. This zone has been identified as hosting large patches of large slip in great earthquakes, based on the lack of high-frequency radiated energy. The large continuous slip patches may facilitate near-complete stress drop. There is seismological evidence for high fluid pressures in the subducted slab around the Moho depth of the overriding plate, suggesting low differential stress levels in this zone due to high fluid pressure, also facilitating stress rotations. The coseismic stress rotations have similar along-strike extent as the mainshock rupture. Postseismic stress rotations tend to occur in the same locations as the coseismic stress rotations, probably due to the very low remaining differential stress following the near-complete coseismic stress drop. The spatial complexity of the observed stress changes suggests that an analytical solution for finding the differential stress from the coseismic stress rotation may be overly simplistic, and that modeling of the full spatial distribution of the mainshock static stress changes is necessary.

Earth, Planets and Space↗

Evaluation of remote mapping techniques for earthquake-triggered landslide inventories in an urban subarctic environment: A case study of the 2018 Anchorage, Alaska Earthquake

Earthquake-induced landslide inventories can be generated using field observations but doing so can be challenging if the affected landscape is large or inaccessible after an earthquake. Remote sensing data can be used to help overcome these limitations. The effectiveness of remotely sensed data to produce landslide inventories, however, is dependent on a variety of factors, such as the extent of coverage, timing, and data quality, as well as environmental factors such as atmospheric interference (e.g., clouds, water vapor) or snow and vegetation cover. With these challenges in mind, we use a combination of field observations and remote sensing data from multispectral, light detection and ranging (lidar), and synthetic aperture radar (SAR) sensors to produce a ground failure inventory for the urban areas affected by the 2018 magnitude (M w ) 7.1 Anchorage, Alaska earthquake. The earthquake occurred during late November at high latitude (∼61°N), and the lack of sunlight, persistent cloud cover, and snow cover that occurred after the earthquake made remote mapping challenging for this event. Despite these challenges, 43 landslides were manually mapped and classified using a combination of the datasets mentioned previously. Using this manually compiled inventory, we investigate the individual performance and reliability of three remote sensing techniques in this environment not typically hospitable to remotely sensed mapping. We found that differencing pre- and post-event normalized difference vegetation index maps and lidar worked best for identifying soil slumps and rapid soil flows, but not as well for small soil slides, soil block slides and rock falls. The SAR-based methods did not work well for identifying any landslide types because of high noise levels likely related to snow. Some landslides, especially those that resulted in minor surface displacement, were identifiable only from the field observations. This work highlights the importance of the rapid collection of field observations and provides guidance for future mappers on which techniques, or combination of techniques, will be most effective at remotely mapping landslides in a subarctic and urban environment.

Alaska↗

Formation of harzburgite by pervasive melt/rock reaction in the upper mantle

Many mantle peridotite samples are too rich in SiO2 (in the form of orthopyroxene) and have ratios of light to heavy rare earth elements that are too high to be consistent with an origin as the residuum of partial melting of the primitive mantle. Trace element studies of melt/rock reaction zones in the Trinity peridotite provide evidence for reaction of the mantle lithosphere with ascending melts, which dissolved calcium-pyroxene and precipitated orthopyroxene as magma mass decreased. This process can account for the observed major and trace element compositions of lithospheric mantle samples, and may accordingly be prevalent in the upper mantle.

Nature↗

Evaluation of ISCCP multisatellite radiance calibration for geostationary imager visible channels using the moon

Since 1983, the International Satellite Cloud Climatology Project (ISCCP) has collected Earth radiance data from the succession of geostationary and polar-orbiting meteorological satellites operated by weather agencies worldwide. Meeting the ISCCP goals of global coverage and decade-length time scales requires consistent and stable calibration of the participating satellites. For the geostationary imager visible channels, ISCCP calibration provides regular periodic updates from regressions of radiances measured from coincident and collocated observations taken by Advanced Very High Resolution Radiometer instruments. As an independent check of the temporal stability and intersatellite consistency of ISCCP calibrations, we have applied lunar calibration techniques to geostationary imager visible channels using images of the Moon found in the ISCCP data archive. Lunar calibration enables using the reflected light from the Moon as a stable and consistent radiometric reference. Although the technique has general applicability, limitations of the archived image data have restricted the current study to Geostationary Operational Environmental Satellite and Geostationary Meteorological Satellite series. The results of this lunar analysis confirm that ISCCP calibration exhibits negligible temporal trends in sensor response but have revealed apparent relative biases between the satellites at various levels. However, these biases amount to differences of only a few percent in measured absolute reflectances. Since the lunar analysis examines only the lower end of the radiance range, the results suggest that the ISCCP calibration regression approach does not precisely determine the intercept or the zero-radiance response level. We discuss the impact of these findings on the development of consistent calibration for multisatellite global data sets.

IEEE Transactions on Geoscience and Remote Sensing↗

Depths of channels in the area of the San Juan Basin Regional Uranium Study, New Mexico, Colorado, Arizona, and Utah

During December 1977 and January 1978 about 280 measurements were made of the depths of channels (arroyos) more than 6 feet deep in the San Juan Basin area. More than half of the measurements were made at sites where channel depths had been previously measured Between 1964 and 1969. Some channels in the western part of the basin had Been re-measured in 1969 and in 1971. The principal areas Being dissected by arroyos are near highlands along the margins of the Basin and in uplands in the northeastern part of the Basin. The most severe dissection by arroyos and the deepest arroyos--commonly Between 40 and 60 feet deep--are in the southeastern part of the Basin. Dissection By arroyos is least in the central part of the Basin near the Chaco River where most arroyos are less than 10 feet deep. Elsewhere, moderate dissection predominates with most arroyos Between 12 and 40 feet deep. Comparison of measurements made from 1964-71 with those made in 1977-78 shows that more channels in the western San Juan Basin were filling than were downcutting. Downcutting or filling was generally less than 2 feet. About two-fifths of the sites measured showed less than half a foot of downcutting or filling. Maximum downcutting was 4 feet along the Rio San Jose in the southeastern part of the basin. Maximum filling of 7 feet was along the Chaco River at the Chaco Canyon National Monument. Along ii other streams elsewhere in the western part of the basin, channels were filled 3 to 4.5 feet. The few measurements made in the southeastern San Juan Basin indicate that since 1964 downcutting has predominated over filling. Large floods during the summer of 1977 caused some change in channel depths in the southwestern part of the San Juan Basin. Some of the channels appeared to have been filled during the years prior to the cutting that occurred from the 1977 floods. At other places, flood flows aggraded (filled) channels. The rate of erosion and arroyo formation in the entire San Juan Basin is effected by man. The southeastern part of the basin--having been occupied by man for several centuries--shows the greatest effects of man on the rate of arroyo formation. Recent urban developments, particularly near Gallup, also appear to have affected the rate of erosion and arroyo formation. In contrast, arroyos appear to be aggrading below many earth-fill dams. In general, the effects of the petroleum, coal, and uranium exploration and development on arroyo formation have been minimal because the main trenching of the arroyos predates oil and mining operations. Some modification--degradation or aggradation--of the arroyos and local trenching of new arroyos have taken place in the area of some of the mines. Most of the observed effects from the mining operations on erosion and aggradation relate to the discharge of mine and mill water into the nearby streams.

Open-File Report↗

The sedimentology and dynamics of crater-affiliated wind streaks in western Arabia Terra, Mars and Patagonia, Argentina

Wind streaks comprise recent aeolian deposits that have been extensively documented on Venus, Earth and Mars. Martian wind streaks are among the most abundant surface features on the planet and commonly extend from the downwind margins of impact craters. Previous studies of wind streaks emerging from crater interior deposits suggested that the mode of emplacement was primarily related to the deposition of silt-sized particles as these settled from plumes. We have performed geologic investigations of two wind streaks clusters; one situated in western Arabia Terra, a region in the northern hemisphere of Mars, and another in an analogous terrestrial site located in southern Patagonia, Argentina, where occurrences of wind streaks emanate from playas within maar craters. In both these regions we have identified bedforms in sedimentary deposits on crater floors, along wind-facing interior crater margins, and along wind streaks. These observations indicate that these deposits contain sand-sized particles and that sediment migration has occurred via saltation from crater interior deposits to wind streaks. In Arabia Terra and in Patagonia wind streaks initiate from crater floors that contain lithic and evaporitic sedimentary deposits, suggesting that the composition of wind streak source materials has played an important role in development. Spatial and topographic analyses suggest that regional clustering of wind streaks in the studied regions directly correlates to the areal density of craters with interior deposits, the degree of proximity of these deposits, and the craters' rim-to-floor depths. In addition, some (but not all) wind streaks within the studied clusters have propagated at comparable yearly (Earth years) rates. Extensive saltation is inferred to have been involved in its propagation based on the studied terrestrial wind streak that shows ripples and dunes on its surface and the Martian counterpart changes orientation toward the downslope direction where it extends into an impact crater. ?? 2009 Elsevier B.V.

Geomorphology↗

Borehole strainmeter measurements spanning the 2014, M w 6.0 South Napa Earthquake, California: The effect from instrument calibration

The 24 August 2014 M w 6.0 South Napa, California earthquake produced significant offsets on 12 borehole strainmeters in the San Francisco Bay area. These strainmeters are located between 24 and 80 km from the source and the observed offsets ranged up to 400 parts-per-billion (ppb), which exceeds their nominal precision by a factor of 100. However, the observed offsets of tidally calibrated strains differ by up to 130 ppb from predictions based on a moment tensor derived from seismic data. The large misfit can be attributed to a combination of poor instrument calibration and better modeling of the strain fit from the earthquake. Borehole strainmeters require in-situ calibration, which historically has been accomplished by comparing their measurements of Earth tides with the strain-tides predicted by a model. Although the borehole strainmeter accurately measure the deformation within the borehole, the long-wavelength strain signals from tides or other tectonic processes recorded in the borehole are modified by the presence of the borehole and the elastic properties of the grout and the instrument. Previous analyses of surface-mounted, strainmeter data and their relationship with the predicted tides suggest that tidal models could be in error by 30%. The poor fit of the borehole strainmeter data from this earthquake can be improved by simultaneously varying the components of the model tides up to 30% and making small adjustments to the point-source model of the earthquake, which reduces the RMS misfit from 130 ppb to 18 ppb. This suggests that relying on tidal models to calibrate borehole strainmeters significantly reduces their accuracy.

California↗

Source spectral properties of small-to-moderate earthquakes in southern Kansas

The source spectral properties of injection-induced earthquakes give insight into their nucleation, rupture processes, and influence on ground motion. Here we apply a spectral decomposition approach to analyze P -wave spectra and estimate Brune-type stress drop for more than 2000 M L 1.5–5.2 earthquakes occurring in southern Kansas from 2014 to 2016. We find that these earthquakes are characterized by low stress drop values (median ∼0.4MPa) compared to natural seismicity in California. We observe a significant increase in stress drop as a function of depth, but the shallow depth distribution of these events is not by itself sufficient to explain their lower stress drop. Stress drop increases with magnitude from M1.5–M3.5, but this scaling trend may weaken above M4 and also depends on the assumed source model. Although we observe a nonstationary, sequence-specific temporal evolution in stress drop, we find no clear systematic relation with the activity of nearby injection wells.

Kansas↗

The Anemomilos prediction methodology for Dst

This paper describes new capabilities for operational geomagnetic D isturbance s torm t ime (Dst) index forecasts. We present a data‐driven, deterministic algorithm called Anemomilos for forecasting Dst out to a maximum of 6 days for large, medium, and small storms, depending upon transit time to the Earth. This capability is used for operational satellite management and debris avoidance in Low Earth Orbit (LEO). Anemomilos has a 15 min cadence, 1 h time granularity, 144 h prediction window (+6 days), and up to 1 h latency. A new finding is that nearly all flare events above a certain irradiance threshold, occurring within a defined solar longitude/latitude region and having sufficient estimated liftoff velocity of ejected material, will produce a geoeffective Dst perturbation. Three solar observables are used for operational Dst forecasting: flare magnitude, integrated flare irradiance through time, and event location. Magnitude is a proxy for ejecta quantity or mass and, combined with speed derived from the integrated flare irradiance, represents the kinetic energy. Speed is estimated as the line‐of‐sight velocity for events within 45° radial of solar disk center. Storms resulting from high‐speed streams emanating from coronal holes are not modeled or predicted. A new result is that solar disk, not limb, observable features are used for predictive techniques. Comparisons between Anemomilos predicted and measured Dst for every hour over 25 months in three continuous time frames between 2001 (high solar activity), 2005 (low solar activity), and 2012 (rising solar activity) are shown. The Anemomilos operational algorithm was developed for a specific customer use related to thermospheric mass density forecasting. It is an operational space weather technology breakthrough using solar disk observables to predict geomagnetically effective Dst up to several days at 1 h time granularity. Real‐time forecasts are presented at http://sol.spacenvironment.net/~sam_ops/index.html?

Space Weather↗

Hydrometeor-enhanced tephra sedimentation: Constraints from the 18 May 1980 eruption of Mount St. Helens

Uncertainty remains on the origin of distal mass deposition maxima observed in many recent tephra fall deposits. In this study the link between ash aggregation and the formation of distal mass deposition maxima is investigated through reanalysis of tephra fallout from the Mount St. Helens 18 May 1980 (MSH80) eruption. In addition, we collate all the data needed to model distal ash sedimentation from the MSH80 eruption cloud. Four particle size subpopulations were present in distal fallout with modes at 2.2 Φ, 4.2 Φ, 5.9 Φ, and 8.3 Φ. Settling rates of the coarsest subpopulation closely matched predicted single‐particle terminal fall velocities. Sedimentation of particles <100 μ m was greatly enhanced, predominantly through aggregation of a particle subpopulation with modal diameter 5.9 ± 0.2 Φ (19 ± 3 μ m). Mammatus on the MSH80 cloud provided a mechanism to transport very fine ash particles, with predicted atmospheric lifetimes of days to weeks, from the upper troposphere to the surface in a matter of hours. In this mechanism, ash particles initiate ice hydrometeor formation high in the troposphere. Subsequently, the volcanic cloud rapidly subsides as mammatus develop from increased particle loading and cloud base sublimation. Rapid fallout occurs as the cloud passes through the melting level in a process analogous to snowflake aggregation. Aggregates sediment en masse and form the distal mass deposition maxima observed in many recent volcanic ash fall deposits. This work provides a data resource that will facilitate tephra sedimentation modeling and allow model intercomparisons.

Washington↗

Interseismic ground deformation and fault slip rates in the greater San Francisco Bay Area from two decades of space geodetic data

The detailed spatial variations of strain accumulation and creep on major faults in the northern San Francisco Bay Area (North Bay), which are important for seismic potential and evaluation of natural hazards, remain poorly understood. Here we combine interferometric synthetic aperture radar data from the ERS‐1/2 and Envisat satellites between 1992 and 2010 with continuous and campaign GPS data to obtain a high spatial and temporal coverage of ground deformation of the North Bay. The SAR data from both ascending and descending orbits are combined to separate horizontal and vertical components of the deformation. We jointly invert the horizontal component of the mean velocities derived from these data to infer the deep strike‐slip rates on major locked faults. We use the estimated deep rates to simulate the long‐wavelength deformation due to interseismic elastic strain accumulation along these locked faults. After removing the long‐wavelength signal from the InSAR horizontal mean velocity field, we estimate fault‐parallel surface creep rates of up to 2 mm/year along the central section of the Rodgers Creek fault and surface creep rates ranging between 2 and 4 mm/year along the Concord fault. No surface creep is geodetically resolved along the West Napa and Green Valley fault zones. We identified characteristically repeating earthquakes on the Rodgers Creek fault, the West Napa fault, the Green Valley fault, and the Concord fault. Nontectonic deformation in the Geysers geothermal field and in Late Cenozoic basins (Rohnert Park and Sonoma basins) are also observed, likely due to hydrological and sediment‐compaction processes, respectively.

Journal of Geophysical Research B: Solid Earth↗

Hydrous carbonates on Mars?: Evidence from Mariner 6/7 infrared spectrometer and ground‐based telescopic spectra

Absorption features at 2.28 and 5.4 μm identified in Mariner 6/7 infrared spectrometer and terrestrial telescopic spectra are consistent with the spectra of hydrous magnesium carbonates such as hydromagnesite and artinite. Spectral characteristics of these hydrous carbonates are different from those of the anhydrous carbonates, as the former do not have the strong spectral features typically associated with anhydrous carbonates such as calcite and siderite. Theoretical mixing indicates that, depending on the type of hydrous carbonate, 10–20 wt % can be incorporated into the regolith without contradicting the spectral observations or the Viking x ray fluorescence chemical analysis. Hydrous carbonates form as weathering products of mafic minerals in the presence of H 2 O and CO 2 , even in the Antarctic. Their formation as evaporite minerals from either original magmas or hydrothermally altered rocks is consistent with the Martian environment, provided liquid water is or has been at least transiently present. On Earth, formation of hydrous Mg carbonates is associated with the production of amorphous iron oxides, which is consistent with both the environment and the inferred surface mineralogy of Mars. These minerals are about 60 wt % H 2 O, CO 3 , and OH; if they are abundant everywhere at the 10% level, then about 6% of the surface weight could be volatiles bound in this type of mineral. Although the stability of hydrous carbonates in a Martian environment is uncertain, there may be kinetic factors inhibiting the dehydration of these minerals, which may persist metastably in the current environment. Although the spectroscopic evidence for anhydrous carbonates is scant, the possible presence of hydrous carbonates provides an appealing mechanism for the existence of carbonates on Mars.

Journal of Geophysical Research↗

Late Neogene geohistory analysis of the Humboldt basin and its relationship to convergence of the Juan de Fuca plate

Geohistory analysis of Neogene Humboldt basin strata provides important constraints for hypotheses of the tectonic evolution of the southern Cascadia subduction margin, leading up to the arrival of the Mendocino triple junction. This analysis suggests that the tectonic evolution of the Humboldt basin area was dominated by coupling between the downgoing Juan de Fuca plate and the continental margin. This coupling is reflected in the timing of major hiatuses within the basin sedimentary sequence and margin uplift and subsidence which occur during periods of tectonic plate adjustment. Stratigraphic evidence indicates that Humboldt basin originated at the base of the continental slope in early Miocene time. Syndepositional uplift of basin strata began in the late Pliocene and was both thermal isostatic and tectonic in origin. Isostatic uplift was a function of an increasingly more buoyant slab being subducted, whereas tectonic uplift was due to imbricated thrusting of the accretionary complex and underplating of offscraped sediment during subduction. A component of margin uplift is postulated to have been caused by a change in the rate of convergence between the Juan de Fuca and North American plates. Coeval with late Pliocene uplift documented onshore was a sharp decrease in covergence rate ∼3 Ma. A reduction in rate of tectonic uplift, observed in the Eel River section, in early Pleistocene time was coeval with a marked increase in relative motion parallel to the continental margin. This localized subsidence may have been caused by syndepositional folding.

Journal of Geophysical Research Solid Earth↗

Historical record of Landsat global coverage

The long-term, 34+ year record of global Landsat remote sensing data is a critical resource to study the Earth system and human impacts on this system. The National Satellite Land Remote Sensing Data Archive (NSLRSDA) is charged by public law to: “maintain a permanent, comprehensive Government archive of global Landsat and other land remote sensing data for long-term monitoring and study of the changing global environment” (U.S. Congress, 1992). The advisory committee for NSLRSDA requested a detailed analysis of observation coverage within the U.S. Landsat holdings, as well as that acquired and held by International Cooperator (IC) stations. Our analyses, to date, have found gaps of varying magnitude in U.S. holdings of Landsat global coverage data, which appear to reflect technical or administrative variations in mission operations. In many cases it may be possible to partially fill these gaps in U.S. holdings through observations that were acquired and are now being held at International Cooperator stations.

Photogrammetric Engineering and Remote Sensing↗

Assessment of forest degradation in Vietnam using Landsat time series data

Landsat time series data were used to characterize forest degradation in Lam Dong Province, Vietnam. We conducted three types of image change analyses using Landsat time series data to characterize the land cover changes. Our analyses concentrated on the timeframe of 1973–2014, with much emphasis on the latter part of that range. We conducted a field trip through Lam Dong Province to develop a better understanding of the ground conditions of the region, during which we obtained many photographs of representative forest sites with Global Positioning System locations to assist us in our image interpretations. High-resolution Google Earth imagery and Landsat data of the region were used to validate results. In general, our analyses indicated that many land-use changes have occurred throughout Lam Dong Province, including gradual forest to non-forest transitions. Recent changes are most marked along the relatively narrow interfaces between agricultural and forest areas that occur towards the boundaries of the province. One important observation is that the most highly protected national reserves in the region have not changed much over the entire Landsat timeframe (1972–present). Spectral changes within these regions have not occurred at the same levels as those areas adjacent to the reserves.

Forests↗

A new database of giant impacts over a wide range of masses and with material strength: A first analysis of outcomes

In the late stage of terrestrial planet formation, planets are predicted to undergo pairwise collisions known as giant impacts. Here, we present a high-resolution database of giant impacts for differentiated colliding bodies of iron–silicate composition, with target masses ranging from 1 × 10 −4 M ⊕ up to super-Earths (5 M ⊕ ). We vary the impactor-to-target mass ratio, core–mantle (iron–silicate) fraction, impact velocity, and impact angle. Strength in the form of friction is included in all simulations. We find that, due to strength, the collisions with bodies smaller than about 2 × 10 −3 M ⊕ can result in irregular shapes, compound-core structures, and captured binaries. We observe that the characteristic escaping velocity of smaller remnants (debris) is approximately half of the impact velocity, significantly faster than currently assumed in N -body simulations of planet formation. Incorporating these results in N -body planet formation studies would provide more realistic debris–debris and debris–planet interactions.

The Planetary Science Journal↗