Search NASA⌕ Search

SEARCH · Search NASA

Results for “surface gravity waves”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Footprints of storms on the sea: A view from spaceborne synthetic aperture radar

Synthetic aperture radar (SAR) on board Seasat observed images of stormlike echoes on the sea in 1978. The core of these images is usually an echo-free hole which is attributed to the damping of the short (30-cm) radar detectable gravity waves by the intense rain in the storm core. Although 'the beating down of waves by rain' is consistent with observations by seafarers and with the first scientific explanation of the phenomenon by Reynolds (1875), neither theory nor experiment has provided definitive support. One experiment appears to provide the key; it shows that the kenetic energy of the rain produces sufficient turbulence in a thin fresh water layer to damp 30-cm waves in 10-20 s, thus producing the echo-free hole. A sequence of positive feedbacks then serves to damp the longer waves. The angular dependence of the sea surface echo cross sections seen by Seasat SAR outside the echo-free hole indicates winds diverging from the downdraft induced by the intense rain core. The wind-generated waves and associated echoes extend out to a sharply defined gust front. The sea surface footprint thus mimics the features of a storm microburst. The variations in surface radar cross section due to a combination of rain and wind effects impacts spaceborne measurements of surface winds by scatterometry and rainfall measurements by radar. Portions of this synthesis remain speculative but serve as hypotheses for further research.

Atlas, David↗
DOE OSTI2023

Chemical Mixing Induced by Internal Gravity Waves in Intermediate-mass Stars

Internal gravity waves can cause mixing in the radiative interiors of stars. We study this mixing by introducing tracer particles into 2D hydrodynamic simulations. Following the work of Rogers & McElwaine, we extend our study to different masses (3, 7, and 20 M ⊙ ) and ages (ZAMS, midMS, and TAMS). The diffusion profiles of these models are influenced by various parameters such as the Brunt–Väisälä frequency, density, thermal damping, the geometric effect, and the frequencies of waves contributing to these mixing profiles. We find that the mixing profile changes dramatically across age. In younger stars, we noted that the diffusion coefficient increases toward the surface, whereas in older stars the initial increase in the diffusion profile is followed by a decreasing trend. We also find that mixing is stronger in more massive stars. Hence, future stellar evolution models should include this variation. In order to aid the inclusion of this mixing in 1D stellar evolution models, we determine the dominant waves contributing to these mixing profiles and present a prescription that can be included in 1D models.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1995

The ocean response to operational westerly wind bursts during the 1991 - 1992 El Nino

Numerical simulations of the remotely forced ocean response to westerly wind bursts prior to and during the 1991 - 1992 El Nino are examined; the models are forced by wind stress from the U.S. Navy's atmospheric global operational analysis/forecast system. The study focuses on (1) the relative response of the first and second internal modes to a single episode of westerly bursts; (2) the role of westerly bursts in producing the eastern Pacific sea level variations from October 1990 to February 1992; and (3) the relative importance of the remotely forced sea level responses generated by central and western Pacific wind anomalies. The simulation use the Naval Research Laboratory global multilayer formulation; the suite of experiments includes hydrodynamic simulations that use both one- and three- active-layer reduced gravity models as well as an experiment that also includes thermodynamic effects. The models are executed from January 1, 1990, to March 1, 1992, a period that includes 10 significant westerly wind bursts or burst clusters and the 1991 - 1992 El Nino event. The numerical experiments reveal an ability to accurately simulate the eastern Pacific sea level variations during this period. In response to a single burst, the three-layer hydrodynamic simulation reveals that the second internal mode Kelvin wave yields a sea level change at the eastern boundary that is approximately one third that of the first mode and a surface velocity signature that is equivalent to the first mode. During the onset of the El Nino event the inclusion of higher modes also produces a more realistic representation of the observed eastern boundary sea level signal. Furthermore, by comparing the model response to particular wind bursts with the observed sea level at Baltra, Galapagos, a value of 2.5 - 2.6 m/s is suggested as the most appropriate mean speed for the first internal mode Kelvin wave during the onset phase. A scenario for the onset of the 1991 - 1992 El Nino is given.

Kindle, John C.↗
NASA NTRS1999

Edge-Flames in Von Karman Swirling Flows

Classical understanding of diffusion flames dictates that they, unlike the premixed flames, do not possess a characteristic propagation velocity and are constrained by stoichiometric requirements at the flame surface. However, it has been commonly observed that when local extinction occurs within a diffusion flame sheet, the edges that are formed propagate with distinct speeds. In general, the propagation speed of these edges depend on their geometrical shape (concave, convex, or straight) among other factors. Recently, Buckmaster investigated the dynamics of straight diffusion flame edges separating burning and quenched regions using simplified one-dimensional models. He showed that these flame edges can have positive, negative, or zero velocity depending on the Damkoehler number of the equilibrium diffusion flame that support them. It was also shown that this unsteady flame-edge behavior is intrinsically linked to S-curve behavior of the diffusion flame with varying Damkoehler number. When the system Damkoehler number lies between the extinction and ignition limits, flame edges can propagate as an "ignition wave" or as a "failure wave," and for a critical Damkoehler number remain as a stationary flame-edge. We have extend Buckmaster's 1-d model to more general edge-flame configurations where the edges appear as "flame holes" or as "flame disks". These two configurations along with the straight-edge case cover the entire range of possible edge-flame geometry observable in planar diffusion-flame sheets. A generalized map of edge-flame propagation velocities as a function of the system Damkoehler number and the edge-flame radius is presented. Experimentally we show that edge flames can be created using diffusion flames embedded in von Karman boundary layers. In a von Karman boundary layer, the flow is generated by spinning a solid (fuel) disk in a quiescent ambient gas. Under normal gravity we were able to produce "flame disks" over a range of fuel-disk rotational velocities varying from 0 to 20 revolutions per second, by orienting the burning surface of the fuel disk facing downward.

Nayagam, Vedha↗
NASA NTRS1992

Sea surface mean square slope from Ku-band backscatter data

A surface mean-square-slope parameter analysis is conducted for 14-GHz airborne radar altimeter near-nadir, quasi-specular backscatter data, which in raw form obtained by least-squares fitting of an optical scattering model to the return waveform show an approximately linear dependence over the 7-15 m/sec wind speed range. Slope data are used to draw inferences on the structure of the high-wavenumber portion of the spectrum. A directionally-integrated model height spectrum that encompasses wind speed-dependent k exp -5/2 and classical Phillips k exp -3 power laws subranges in the range of gravity waves is supported by the data.

Jackson, F. C.↗
NASA NTRS1977

Forced baroclinic ocean motions. II - The linear equatorial bounded case

Results obtained in the first part of this paper (1976), which studied forced baroclinic ocean motions in the linear equatorial unbounded case, are extended to the case of an ocean bounded by two meridians. A complete solution for the linear spin-up in response to a switched-on, x-dependent wind stress is obtained in terms of five components: the unbounded nonoscillatory response to the wind stress forcing; the inertia-gravity waves generated when the forcing switches on, together with their reflections; the quasi-geostrophic Rossby modes forming the eastern boundary response and a western boundary layer; and the equatorial Kelvin wave component of the western boundary response. Kelvin waves are generated if the forcings include a zonal wind component that is symmetric about the equator or an antisymmetric meridional wind stress. For non-Kelvin symmetries, spin-up occurs entirely by the effects of Rossby waves emanating from the eastern boundary. The state approaches the steady solution as increasing numbers of these reach an interior point. In Kelvin symmetries, the Rossby waves act to bring the sea surface tilt to the steady value and the spin-up of the tilt proceeds as in non-Kelvin cases.

Cane, M. A.↗
NASA NTRS1980

Optical image and laser slope meter intercomparisons of high-frequency waves

Spectral analyses of optical images of the ocean surface, obtained by a digital video system, are presented and compared with wave data measured simultaneously by the JPL Waverider-mounted laser slope meter. The image analyses, which incorporate several new ideas, provide two-dimensional wave number spectra of slope, covering wavelengths from 10 cm to 10 m. These slope spectra are converted to wave height spectra by a new technique which includes the effects of sky radiance gradients. Space-time spectra are also presented for waves whose frequencies are less than 2 Hz. The JPL slope frequency spectra are compared with image wave number spectra which have been converted to frequency spectra by use of the gravity wave dispersion relation. Results of comparisons between the frequency spectra obtained from the two different measurements show reasonable agreement for frequencies less than 3 Hz.

Lubard, S. C.↗
NASA NTRS1987

Atmosphere dynamics in the equatorial meteor zone

The study of the atmospheric circulation of the Earth from its surface to the altitudes of 100 to 110 km is essential for establishing atmospheric motion regularities with a view toward perfecting weather forecasting. The main results of the Soviet equatorial meteor expedition (SEME) are presented. A continuous cycle of measurements was carried out. Considerable interdiurnal variation of the zonal component was observed. Importantly, in the meridional component, the prevalence of a two day component was established in the equatorial meteor zone for the first time. The pronounced westward motion of the atmosphere over the equator is noted. The SEME data analysis has shown that the meteor zone is characterized by flashes of intensity of the internal gravity waves and turbulence at highest instability moments of atmosphere due to tidal motion.

Kascheev, B. L.↗
NASA NTRS1981

Numerical prediction of the monsoon depression of 5-7 July 1979

A well defined monsoon depression was used for two assimilation and forecast experiments: (1) using conventional surface and upper air data, (2) using these data plus Monex data. The data sets were assimilated and used with a general circulation model to make numerical predictions. The model, the analysis and assimilation procedure, the differences in the analyses due to different data inputs, and the differences in the numerical predictions are described. The MONEX data have a positive impact, although the differences after 24 hr are not significant. The MONEX assimilation does not agree with manual analysis location of depression center. The 2.5 x 3 deg horizontal resolution of the prediction model is too coarse. The assimilation of geopotential height data derived from satellite soundings generated gravity waves with amplitudes similar to the meteorologically significant features investigated.

Shukla, J.↗
NASA NTRS1995

Assimilation of TOPEX Sea Level Measurements with a Reduced-Gravity, Shallow Water Model of the Tropical Pacific Ocean

Sea surface height variability measured by TOPEX is analyzed in the tropical Pacific Ocean by way of assimilation into a wind-driven, reduced-gravity, shallow water model using an approximate Kalman filter and smoother. The analysis results in an optimal fit of the dynamic model to the observations, providing it dynamically consistent interpolation of sea level and estimation of the circulation. Nearly 80% of the expected signal variance is accounted for by the model within 20 deg of the equator, and estimation uncertainty is substantially reduced by the voluminous observation. Notable features resolved by the analysis include seasonal changes associated with the North Equatorial Countercurrent and equatorial Kelvin and Rossby waves. Significant discrepancies are also found between the estimate and TOPEX measurements, especially near the eastern boundary. Improvements in the estimate made by the assimilation are validated by comparisons with independent tide gauge and current meter observations. The employed filter and smoother are based on approximately computed estimation error covariance matrices, utilizing a spatial transformation and an symptotic approximation. The analysis demonstrates the practical utility of a quasi-optimal filter and smoother.

Fukumori, Ichiro↗
DOE OSTI2018

Membrane-induced hydroelastic migration of a particle surfing its own wave

While coupling between fluid flow and soft elastic surfaces is common in biology and engineering, an analytical description is challenging as it often involves non-linear dynamics. Here we show using theory and experiments that a small particle moving along an elastic membrane through a viscous fluid is repelled from the membrane due to hydroelastic forces. The flow field produces an elastic disturbance in the membrane leading to particle–wave coupling. We derive an analytic expression for the particle trajectory and find that the normal migration velocity of the particle is quadratic in its speed and depends on a combination of the tension and bending resistance of the membrane. Experimentally, we measure the normal displacement of spheres sedimenting under gravity along a suspended elastic membrane and find quantitative agreement with the theoretical predictions with no fitting parameters. We experimentally demonstrate that the effect is strong enough for separation and sorting of particles on the basis of both their size and density. We discuss the significance of our results for particles interacting with biological membranes, and propose the use of our model for membrane elasticity measurements.

catalysis (homogeneous), solar (photovoltaic), bio↗
NASA NTRS1978

Modulation of centimetric waves by long gravity waves - Progress report on field and laboratory results

Results are presented from laboratory and field experiments on the modulation of short waves by long waves. The field study employed a wave follower capable of tracking ocean waves with frequencies less than 1.0 Hz and heights less than 2.0 m. A high-response laser-optical system was used to detect upwind-downwind and cross-wind slopes of short waves. The laboratory study was conducted with wind over periodic long waves. The laboratory findings are discussed and compared with laboratory radar measurements and also short wave measurements obtained in the field. It is found that long waves significantly modulate the short wave dispersion by their orbital velocity, and that demodulation is necessary if the data collected by remote sensors are to be compared to surface penetrating devices. The modulation level is weak for wavelengths in the range 2.76-3.30 cm. Other relevant results are also presented.

Shemdin, O. H.↗
NASA NTRS2004

High-Order Shock-Capturing Methods for Modeling Dynamics of the Solar Atmosphere

We use one-dimensional high-order central shock capturing numerical methods to study the response of various model solar atmospheres to forcing at the solar surface. The dynamics of the atmosphere is modeled with the Euler equations in a variable-sized flux tube in the presence of gravity. We study dynamics of the atmosphere suggestive of spicule formation and coronal oscillations. These studies are performed on observationally-derived model atmospheres above the quiet sun and above sunspots. To perform these simulations, we provide a new extension of existing second- and third- order shock-capturing methods to irregular grids. We also solve the problem of numerically maintaining initial hydrostatic balance via the introduction of new variables in the model equations and a careful initialization mechanism. We find several striking results: all model atmospheres respond to a single impulsive perturbation with several strong shock waves consistent with the rebound-shock model. These shock waves lift material and the transition region well into the initial corona, and the sensitivity of this lift to the initial impulse depends non-linearly on the details of the atmosphere model. We also reproduce an observed 3-minute coronal oscillation above sunspots compared to 5-minute oscillations above the quiet sun.

Bryson, Steve↗
NASA NTRS1993

The NASA/GISS Mars general circulation model: Preliminary experiments

The NASA/GISS Mars General Circulation Model (GCM) is an adapted version of the GISS Global Climate/Middle Atmosphere Model, specifically developed for the diagnostic validation and objective analysis of measured atmospheric temperatures from the Mars Observer Pressure Modulator Infrared Radiometer (PMIRR) experiment. The GISS Mars GCM has 23 vertical layers extending from the surface to approximately 80 km altitude, representing a vertical resolution of about 0.3 scale heights. The primitive (vertically hydrostatic) equations are solved in finite difference form on the Krakawa B grid, with a horizontal resolution of 8 deg x 10 deg (latitude-longitude). The model includes a diurnal solar cycle, heat transport within a two-layer ground, and a high-order 'slopes-scheme' for the advection of heat in the upper atmosphere. The radiative transfer scheme is based on the correlated k distribution method for the treatment of nongray gaseous absorption thermal emission, and multiple scattering, including options for suspended dust. A special feature of the model of particular importance for Mars is a parameterization of gravity-wave-induced drag incorporating orographic forcing, wind shear, convection, and radiative damping. The implementation of the GISS Mars model includes global maps of topography, roughness, and albedo.

Allison, Michael↗
NASA NTRSDate not supplied

Response of the QBO to Abrupt Increases in CO2 Using Three Atmospheric Chemistry Configurations

Long-term projections of the Quasi-Biennial Oscillation (QBO) remain highly uncertain. This is partly due to the paucity of models which are able to properly simulate that phenomenon. Only 5 of the 47 CMIP5 models are capable of spontaneously generating a realistic QBO (Butchart et al., 2018), and even those models exhibit large biases in key QBO characteristics (e.g. amplitude, period, vertical extent) when compared with observations. Furthermore, only 1 of these 5 employed interactive atmospheric chemistry, which is known to modulate QBO dynamics. We here investigate the QBO response to increased greenhouse gases using the NASA Goddard Institute for Space Studies Middle Atmosphere Model E2.2. Compared to lower vertical resolution versions of Model E, version 2.2 has a higher model top (0.002 hPa), and additional interactive non-orographic gravity wave drag sources from convection and shear, which produce a sufficiently realistic QBO, thus rendering it suitable for use in climate change studies. Steady-state responses to doubled and quadrupled CO2 from a pre-industrial control are analyzed, as well as the transient response to a 1% per year CO2 increase. In addition, we systematically explore the impact of interactive chemistry in modulating the QBO response to increased CO2 by contrasting interactive, prescribed, and linearized ozone chemistry configurations of the model. Overall, in response to increase CO2 concentrations the QBO is seen to increase in frequency and weaken in amplitude, consistent with previous results, but the memory of the tropical stratosphere may complicate assessments of trends in chemistry and surface impacts. We also discuss implications for the trade-off between ensemble size and the complexity of the chemistry scheme in the model.

Quasi-Biennial Oscillation (QBO)↗
NASA NTRS2011

NASA Tech Briefs, January 2011

The topics include: 1) Distributed Aerodynamic Sensing and Processing Toolbox; 2) Collaborative Supervised Learning for Sensor Networks; 3) Hazard Detection Software for Lunar Landing; 4) Onboard Nonlinear Engine Sensor and Component Fault Diagnosis and Isolation Scheme; 5) Network-Capable Application Process and Wireless Intelligent Sensors for ISHM; 6) Interface Supports Multiple Broadcast Transceivers for Flight Applications; 7) FPGA Sequencer for Radar Altimeter Applications; 8) Miniature Sapphire Acoustic Resonator - MSAR; 9) Process-Hardened, Multi-Analyte Sensor for Characterizing Rocket Plume Constituents; 10) SAD5 Stereo Correlation Line-Striping in an FPGA; 11) Hybrid Composite Cryogenic Tank Structure; 12) Nanoscale Deformable Optics; 13) Reliability-Based Design Optimization of a Composite Airframe Component; 14) Zinc Oxide Nanowire Interphase for Enhanced Lightweight Polymer Fiber Composites; 15) Plasma Igniter for Reliable Ignition of Combustion in Rocket Engines; 16) Wire Test Grip Fixture; 17) A Sub-Hertz, Low-Frequency Vibration Isolation Platform; 18) Carbon Nanofibers Synthesized on Selective Substrates for Nonvolatile Memory and 3D Electronics; 19) Nanoparticle/Polymer Nanocomposite Bond Coat or Coating; 20) High-Resolution Wind Measurements for Offshore Wind Energy Development; 21) Spring Tire; 22) Marsviewer 2008; 23) Mission Services Evolution Center Message Bus; 24) Major Constituents Analysis for the Vehicle Cabin Atmosphere Monitor; 25) Astronaut Health Participant Summary Application; 26) Adaption of the AMDIS Method to Flight Status on the VCAM Instrument; 27) Natural Language Interface for Safety Certification of Safety-Critical Software; 28) Cryogenic Caging for Science Instrumentation; 29) Wide-Range Neutron Detector for Space Nuclear Applications; 30) In Situ Guided Wave Structural Health Monitoring System; 31) Multiplexed Energy Coupler for Rotating Equipment; 32) Attitude Estimation in Fractionated Spacecraft Cluster Systems; 33) Full Piezoelectric Multilayer-Stacked Hybrid Actuation/Transduction Systems; 34) Active Flow Effectors for Noise and Separation Control; 35) Method and System for Temporal Filtering in Video Compression Systems; 36) Apparatus for Measuring Total Emissivity of Small, Low-Emissivity Samples; 37) Multiple-Zone Diffractive Optic Element for Laser Ranging Applications; 38) Simplified Architecture for Precise Aiming of a Deep-Space Communication Laser Transceiver; 39) Two-Photon-Absorption Scheme for Optical Beam Tracking; 40) High-Sensitivity, Broad-Range Vacuum Gauge Using Nanotubes for Micromachined Cavities; 41) Wide-Field Optic for Autonomous Acquisition of Laser Link; 42) Extracting Zero-Gravity Surface Figure of a Mirror; 43) Modeling Electromagnetic Scattering From Complex Inhomogeneous Objects; 44) Visual Object Recognition and Tracking of Tools; 45) Method for Implementing Optical Phase Adjustment; 46) Visual SLAM Using Variance Grid Maps; 47) Rapid Calculation of Spacecraft Trajectories Using Efficient Taylor Series Integration; 48) Efficient Kriging Algorithms; 49) Predicting Spacecraft Trajectories by the WeavEncke Method; 50) An Augmentation of G-Guidance Algorithms; 51) Comparison of Aircraft Icing Growth Assessment Software; 52) Silicon-Germanium Voltage-Controlled Oscillator at 105 GHz; 53) Estimation of Coriolis Force and Torque Acting on Ares-1; 54) Null Lens Assembly for X-Ray Mirror Segments; and 55) High-Precision Pulse Generator.

Source record↗
NASA NTRS1993

Lee waves: Benign and malignant

The flow of an incompressible fluid over an obstacle will produce an oscillation in which buoyancy is the restoring force, called a gravity wave. For disturbances of this scale, the atmosphere may be treated as dynamically incompressible, even though there exists a mean static upward density gradient. Even in the linear approximation - i.e., for small disturbances - this model explains a great many of the flow phenomena observed in the lee of mountains. However, nonlinearities do arise importantly, in three ways: (1) through amplification due to the decrease of mean density with height; (2) through the large (scaled) size of the obstacle, such as a mountain range; and (3) from dynamically singular levels in the fluid field. These effects produce a complicated array of phenomena - large departure of the streamlines from their equilibrium levels, high winds, generation of small scales, turbulence, etc. - that present hazards to aircraft and to lee surface areas. The nonlinear disturbances also interact with the larger-scale flow in such a manner as to impact global weather forecasts and the climatological momentum balance. If there is no dynamic barrier, these waves can penetrate vertically into the middle atmosphere (30-100 km), where recent observations show them to be of a length scale that must involve the coriolis force in any modeling. At these altitudes, the amplitude of the waves is very large, and the phenomena associated with these wave dynamics are being studied with a view to their potential impact on high performance aircraft, including the projected National Aerospace Plane (NASP). The presentation shows the results of analysis and of state-of-the-art numerical simulations, validated where possible by observational data, and illustrated with photographs from nature.

Wurtele, M. G.↗
NASA NTRS1993

Collisional and dynamical history of Gaspra

Interpretation of the impact record on Gaspra requires understanding of the effects of collisions on a target body of Gaspra's size and shape, recognition of impact features that may have different morphologies from craters on larger planets, and models of the geological processes that erase and modify impact features. Crater counts on the 140 sq km of Gaspra imaged at highest resolution by the Galileo spacecraft show a steep size-frequency distribution (cumulative power-law index near -3.5) from the smallest resolvable size (150 m diameter) up through the large feature (1.5 km diameter crater) of familiar crater-like morphology. In addition, there appear to be as many as eight roughly circular concavities with diameters greater than 3 km visible on the asteroid. If we restrict our crater counts to features with traditionally recognized crater morphologies, these concavities would not be included. However, if we define craters to include any concave structures that may represent local or regional damage at an impact size, then the larger features on Gaspra are candidates for consideration. Acceptance of the multi-km features as craters has been cautious for several reasons. First, scaling laws (the physically plausible algorithms for extrapolating from experimental data) indicate that Gaspra could not have sustained such large-crater-forming impacts without being disrupted; second, aside from concavity, the larger structures have no other features (e.g. rims) that can be identified with known impact craters; and third, extrapolation of the power-law size distribution for smaller craters predicts no craters larger than 3 km over the entire surface. On the other hand, recent hydrocode modeling of impacts shows that for given impact (albeit into a sphere), the crater size is much larger than given by scaling laws. Gaspra-size bodies can sustain formation of up to 8-km craters without disruption. Besides allowing larger impact craters, this result doubles the lifetime since the last catastrophic fragmentation event up to one billion years. Events that create multi-km craters also globally damage the material structure, such that regolith is produced, whether or not Gaspra 'initially' had a regolith, contrary to other models in which initial regolith is required in order to allow current regolith. Because the globally destructive shock wave precedes basin formation, crater size is closer to the large size extrapolated from gravity-scaling rather than the strength-scaling that had earlier been assumed for such small bodies. This mechanism may also help explain the existence of Stickney on Phobos. Moreover, rejection of the large concavities as craters based on unfamiliar morphology would be premature, because (aside from Stickney) we have no other data on such large impact structures on such a small, irregular body. The eight candidate concavities cover an area greater than that counted for smaller craters, because they are most apparent where small craters cannot be seen: on low resolution images and at the limb on high resolution images. We estimate that there are at least two with diameter greater than 4 km per 140 sq km, which would have to be accounted for in any model that claims these are impact craters.

Greenberg, R.↗