USGS Science⌕ Search

SEARCH · USGS Science

Results for “Earth”

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,153 records · Page 64Linked to original sources
USGS2009

A global map of rainfed cropland areas (GMRCA) at the end of last millennium using remote sensing

The overarching goal of this study was to produce a global map of rainfed cropland areas (GMRCA) and calculate country-by-country rainfed area statistics using remote sensing data. A suite of spatial datasets, methods and protocols for mapping GMRCA were described. These consist of: (a) data fusion and composition of multi-resolution time-series mega-file data-cube (MFDC), (b) image segmentation based on precipitation, temperature, and elevation zones, (c) spectral correlation similarity (SCS), (d) protocols for class identification and labeling through uses of SCS R 2 -values, bi-spectral plots, space-time spiral curves (ST-SCs), rich source of field-plot data, and zoom-in-views of Google Earth (GE), and (e) techniques for resolving mixed classes by decision tree algorithms, and spatial modeling. The outcome was a 9-class GMRCA from which country-by-country rainfed area statistics were computed for the end of the last millennium. The global rainfed cropland area estimate from the GMRCA 9-class map was 1.13 billion hectares (Bha). The total global cropland areas (rainfed plus irrigated) was 1.53 Bha which was close to national statistics compiled by FAOSTAT (1.51 Bha). The accuracies and errors of GMRCA were assessed using field-plot and Google Earth data points. The accuracy varied between 92 and 98% with kappa value of about 0.76, errors of omission of 2-8%, and the errors of commission of 19-36%. ?? 2008 Elsevier B.V.

International Journal of Applied Earth Observation↗
USGS2013

Next generation of global land cover characterization, mapping, and monitoring

Land cover change is increasingly affecting the biophysics, biogeochemistry, and biogeography of the Earth's surface and the atmosphere, with far-reaching consequences to human well-being. However, our scientific understanding of the distribution and dynamics of land cover and land cover change (LCLCC) is limited. Previous global land cover assessments performed using coarse spatial resolution (300 m–1 km) satellite data did not provide enough thematic detail or change information for global change studies and for resource management. High resolution (∼30 m) land cover characterization and monitoring is needed that permits detection of land change at the scale of most human activity and offers the increased flexibility of environmental model parameterization needed for global change studies. However, there are a number of challenges to overcome before producing such data sets including unavailability of consistent global coverage of satellite data, sheer volume of data, unavailability of timely and accurate training and validation data, difficulties in preparing image mosaics, and high performance computing requirements. Integration of remote sensing and information technology is needed for process automation and high-performance computing needs. Recent developments in these areas have created an opportunity for operational high resolution land cover mapping, and monitoring of the world. Here, we report and discuss these advancements and opportunities in producing the next generations of global land cover characterization, mapping, and monitoring at 30-m spatial resolution primarily in the context of United States, Group on Earth Observations Global 30 m land cover initiative (UGLC).

International Journal of Applied Earth Observation↗
USGS2016

Paleozoic magmatism and porphyry Cu-mineralization in an evolving tectonic setting in the North Qilian Orogenic Belt, NW China

The NWW-striking North Qilian Orogenic Belt records the Paleozoic accretion–collision processes in NW China, and hosts Paleozoic Cu–Pb–Zn mineralization that was temporally and spatially related to the closure of the Paleo Qilian-Qinling Ocean. The Wangdian Cu deposit is located in the eastern part of the North Qilian Orogenic Belt, NW China. Copper mineralization is spatially associated with an altered early Paleozoic porphyritic granodiorite, which intruded tonalites and volcaniclastic rocks. Alteration zones surrounding the mineralization progress outward from a potassic to a feldspar-destructive phyllic assemblage. Mineralization consists mainly of quartz-sulfide stockworks and disseminated sulfides, with ore minerals chalcopyrite, pyrite, molybdenite, and minor galena and sphalerite. Gangue minerals include quartz, orthoclase, biotite, sericite, and K-feldspar. Zircon LA-ICPMS U–Pb dating of the ore-bearing porphyritic granodiorite yielded a mean 206 Pb/ 238 U age of 444.6 ± 7.8 Ma, with a group of inherited zircons yielding a mean U–Pb age of 485 ± 12 Ma, consistent with the emplacement age (485.3 ± 6.2 Ma) of the barren precursor tonalite. Rhenium and osmium analyses of molybdenite grains returned model ages of 442.9 ± 6.8 Ma and 443.3 ± 6.2 Ma, indicating mineralization was coeval with the emplacement of the host porphyritic granodiorite. Rhenium concentrations in molybdenite (208.9–213.2 ppm) suggest a mantle Re source. The tonalities are medium-K calc-alkaline. They are characterized by enrichment of light rare-earth elements (LREEs) and large-ion lithophile elements (LILEs), depletion of heavy rare-earth elements (HREEs) and high-field-strength elements (HFSEs), and minor negative Eu anomalies. They have ε Hf ( t ) values in the range of +3.6 to +11.1, with two-stage Hf model ages of 0.67–1.13 Ga, suggesting that the ca. 485 Ma barren tonalites were products of arc magmatism incorporating melts from the mantle wedge and the lithosphere. In contrast, the 40-m.y.-younger ore-bearing porphyritic granodiorite is sub-alkaline and peraluminous. They are enriched in LREEs and LILEs, depleted in HFSEs, and show weak negative Eu anomalies. They display ε Hf ( t ) values of captured or inherited zircons in the range of +8.5 to +10.0, and younger two-stage Hf model ages of 0.78 Ga and 0.86 Ga, similar to those of ca. 485 Ma tonalite. The ca. 445 Ma zircons have ε Hf ( t ) values of −2.1 to +9.9, with two-stage Hf model ages of 0.75–1.27 Ga. Moreover, they have relatively high oxygen fugacity than that of the precursor barren tonalite. The ca. 445 Ma magmas at Wangdian thus formed in a subduction setting, and incorporated melts from the subduction-modified lithosphere that had previously been enriched by additions of chalcophile and siderophile element-rich materials by the earlier magmatism and metasomatism during the Paleo Qilian-Qinling Ocean subduction event.

North Qilian Orogenic Belt↗
USGS2010

Studies from the history of soil science and geology

The United Nations proclaimed the year 2008 as the official International Year of Planet Earth (IYPE), with science and outreach activities spanning 2007–2009. IYPE-sponsored outreach helped focus the attention of the general public on topics such as human health and the environment; ocean and natural resources sustainability; mitigating natural hazards and community resilience; and the effects of climate change. Within the earth science community, the IYPE was a stimulus for retrospection, and for efforts aimed at bridging divides within the community. One such effort was the first joint meeting of the Geological Society of America (GSA) and the Soil Science Society of America (SSSA), held in Houston, Texas, 5–9 October 2008.

Physics and Chemistry of the Earth, Parts A/B/C↗
USGS1998

From polar wander to dynamic planet: A tribute to Keith Runcorn

The evolution of Keith Runcorn's ideas from the static, elastic Earth of Jeffreys to a dynamic, convecting planet are presented based on discussions as his colleague over 25 years from 1963 to his retirement in 1988. Keith reached the concept of a dynamic planet by way of polar wander and continental drift using palaeomagnetism with great zeal. It took some time for Keith to convince himself of the reality of continental drift. Once convinced, he became an evangelist converting others and enthusiastically pursuing possible mechanisms for explaining it, homing in on mantle convection. To establish the nature and history of mantle convection his interests ranged from the world rift system and satellite gravity anomalies to the radiometric age peaks. A great step forward occurred when we realised that a region of uprising convection was not necessary under all the rifts as had been commonly advocated in the 1960s and that the mid-Atlantic, African and Indian Ocean rifts could be equally well explained by one large region of upwelling mantle convection. It was also realised that the plate convergence zones (island arcs, trenches and deep focus earthquakes) were much better correlated with the satellite gravity anomalies and it was much easier to locate the possible regions of downwelling mantle convection. Now, seismic tomography helps to establish the nature of mantle convection and it appears that relations among the Earth's surface features, the geoid anomalies and peturbations of mantle seismic velocities are near to being established. In the next few years a far better and accurate picture of the geometry of mantle convection so enthusiastically advocated by Keith Runcorn is likely to be seen.

Physics and Chemistry of the Earth↗
USGS1999

Climatic impact of glacial cycle polar motion: Coupled oscillations of ice sheet mass and rotation pole position

Precessional motion of Earth's rotation axis relative to its orbit is a well-known source of long-period climatic variation. It is less well appreciated that growth and decay of polar ice sheets perturb the symmetry of the global mass distribution enough that the geographic location of the rotation axis will change by at least 15 km and possibly as much as 100 km during a single glacial cycle. This motion of the pole will change the seasonal and latitudinal pattern of temperatures. We present calculations, based on a diurnal average energy balance, which compare the summer and winter temperature anomalies due to a 1° decrease in obliquity with those due to a 1° motion of the rotation pole toward Hudson Bay. Both effects result in peak temperature perturbations of about 1° Celsius. The obliquity change primarily influences the amplitude of the seasonal cycle, while the polar motion primarily changes the annual mean temperatures. The polar motion induced temperature anomaly is such that it will act as a powerful negative feedback on ice sheet growth. We also explore the evolution of the coupled system composed of ice sheet mass and pole position. Oscillatory solutions result from the conflicting constraints of rotational and thermal stability. A positive mass anomaly on an otherwise featureless Earth is in rotational equilibrium only at the poles or the equator. The two polar equilibria are rotationally unstable, and the equatorial equilibrium, though rotationally stable, is thermally unstable. We find that with a plausible choice for the strength of coupling between the thermal and rotational systems, relatively modest external forcing can produce significant response at periods of 10 4 –10 6 years, but it strongly attenuates polar motion at longer periods. We suggest that these coupled oscillations may contribute to the observed dominance of 100 kyr glacial cycles since the mid-Pleistocene and will tend to stabilize geographic patterns that are suitable to glaciations.

Journal of Geophysical Research B: Solid Earth↗
USGS2000

Geochemical evidence for a magmatic CO 2 degassing event at Mammoth Mountain, California, September-December 1997

Recent time series soil CO 2 concentration data from monitoring stations in the vicinity of Mammoth Mountain, California, reveal strong evidence for a magmatic degassing event during the fall of 1997 lasting more than 2 months. Two sensors at Horseshoe Lake first recorded the episode on September 23, 1997, followed 10 days later by a sensor on the north flank of Mammoth Mountain. Direct degassing from shallow intruding magma seems an implausible cause of the degassing event, since the gas released at Horseshoe Lake continued to be cold and barren of other magmatic gases, except for He. We suggest that an increase in compressional strain on the area south of Mammoth Mountain driven by movement of major fault blocks in Long Valley caldera may have triggered an episode of increased degassing by squeezing additional accumulated CO 2 from a shallow gas reservoir to the surface along faults and other structures where it could be detected by the CO 2 monitoring network. Recharge of the gas reservoir by CO 2 emanating from the deep intrusions that probably triggered deep long-period earthquakes may also have contributed to the degassing event. The nature of CO 2 discharge at the soil-air interface is influenced by the porous character of High Sierra soils and by meteorological processes. Solar insolation is the primary source of energy for the Earth atmosphere and plays a significant role in most diurnal processes at the Earth surface. Data from this study suggest that external forcing due largely to local orographic winds influences the fine structure of the recorded CO 2 signals.

California↗
USGS2001

Resolution analysis of finite fault source inversion using one- and three-dimensional Green's functions 2. Combining seismic and geodetic data

Using numerical tests for a prescribed heterogeneous earthquake slip distribution, we examine the importance of accurate Green's functions (GF) for finite fault source inversions which rely on coseismic GPS displacements and leveling line uplift alone and in combination with near-source strong ground motions. The static displacements, while sensitive to the three-dimensional (3-D) structure, are less so than seismic waveforms and thus are an important contribution, particularly when used in conjunction with waveform inversions. For numerical tests of an earthquake source and data distribution modeled after the 1994 Northridge earthquake, a joint geodetic and seismic inversion allows for reasonable recovery of the heterogeneous slip distribution on the fault. In contrast, inaccurate 3-D GFs or multiple 1-D GFs allow only partial recovery of the slip distribution given strong motion data alone. Likewise, using just the GPS and leveling line data requires significant smoothing for inversion stability, and hence, only a blurred vision of the prescribed slip is recovered. Although the half-space approximation for computing the surface static deformation field is no longer justifiable based on the high level of accuracy for current GPS data acquisition and the computed differences between 3-D and half-space surface displacements, a layered 1-D approximation to 3-D Earth structure provides adequate representation of the surface displacement field. However, even with the half-space approximation, geodetic data can provide additional slip resolution in the joint seismic and geodetic inversion provided a priori fault location and geometry are correct. Nevertheless, the sensitivity of the static displacements to the Earth structure begs caution for interpretation of surface displacements, particularly those recorded at monuments located in or near basin environments.

Journal of Geophysical Research B: Solid Earth↗
USGS2001

Resolution analysis of finite fault source inversion using one- and three-dimensional Green's functions 1. Strong motions

We develop a methodology to perform finite fault source inversions from strong motion data using Green's functions (GFs) calculated for a three-dimensional (3-D) velocity structure. The 3-D GFs are calculated numerically by inserting body forces at each of the strong motion sites and then recording the resulting strains along the target fault surface. Using reciprocity, these GFs can be recombined to represent the ground motion at each site for any (heterogeneous) slip distribution on the fault. The reciprocal formulation significantly reduces the required number of 3-D finite difference computations to at most 3NS, where NS is the number of strong motion sites used in the inversion. Using controlled numerical resolution tests, we have examined the relative importance of accurate GFs for finite fault source inversions which rely on near-source ground motions. These experiments use both 1-D and 3-D GFs in inversions for hypothetical rupture models in order (1) to analyze the ability of the 3-D methodology to resolve trade-offs between complex source phenomena and 3-D path effects, (2) to address the sensitivity of the inversion results to uncertainties in the 3-D velocity structure, and (3) to test the adequacy of the 1-D GF method when propagation effects are known to be three-dimensional. We find that given "data" from a prescribed 3-D Earth structure, the use of well-calibrated 3-D GFs in the inversion provides very good resolution of the assumed slip distribution, thus adequately separating source and 3-D propagation effects. In contrast, using a set of inexact 3-D GFs or a set of hybrid 1-D GFs allows only partial recovery of the slip distribution. These findings suggest that in regions of complex geology the use of well-calibrated 3-D GFs has the potential for increased resolution of the rupture process relative to 1-D GFs. However, realizing this full potential requires that the 3-D velocity model and associated GFs should be carefully validated against the true 3-D Earth structure before performing the inverse problem with actual data.

Journal of Geophysical Research B: Solid Earth↗
USGS2011

Seismically observed seiching in the Panama Canal

A large portion of the seismic noise spectrum is dominated by water wave energy coupled into the solid Earth. Distinct mechanisms of water wave induced ground motions are distinguished by their spectral content. For example, cultural noise is generally <1 s period, microseisms dominate the seismic spectrum from periods of 2 to 20 s, and the Earth's "hum" is in the range of 50 to 600 s. We show that in a large lake in the Panama Canal there is an additional source of long-period noise generated by standing water waves, seiches, induced by disturbances such as passing ships and wind pressure. We compare seismic waveforms to water level records and relate these observations to changes in local tilt and gravity due to an oscillating seiche. The methods and observations discussed in this paper provide a first step toward quantifying the impact of water inundation as recorded by seismometers. This type of quantified understanding of water inundation will help in future estimates of similar phenomena such as the seismic observations of tsunami impact. Copyright 2011 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2016

The geomagnetic blitz of September 1941

Seventy-five years ago, on 18–19 September 1941, the Earth experienced a great magnetic storm, one of the most intense ever recorded. It arrived at a poignant moment in history, when radio and electrical technology was emerging as a central part of daily life and when much of the world was embroiled in World War II, which the United States had not yet officially entered. Auroras danced across the night sky as voltage surged in power grid lines. A radio blackout interrupted fan enjoyment of a baseball game, while another radio program was interrupted by private phone conversations. Citizens, already on edge, wondered if neon lights were some sort of antiaircraft signal. And far away in the North Atlantic, the illuminated night sky exposed an Allied convoy to German attack. These effects raised awareness within the scientific community and among the public of the societal significance of the effects that the Sun and outer space can have on the Earth—what we now call space weather .

Eos, Earth and Space Science News↗
USGS2020

Timescales and processes of methane hydrate formation and breakdown, with application to geologic systems

Gas hydrate is an ice-like form of water and low molecular weight gas stable at temperatures of roughly -10ºC to 25ºC and pressures of ~3 to 30 MPa in geologic systems. Natural gas hydrates sequester an estimated one-sixth of Earth’s methane and are found primarily in deepwater marine sediments on continental margins, but also in permafrost areas and under continental ice sheets. When gas hydrate is removed from its stability field, its breakdown has implications for the global carbon cycle, ocean chemistry, marine geohazards, and interactions between the geosphere and the ocean-atmosphere system. Gas hydrate breakdown can also be artificially driven as a component of studies assessing the resource potential of these deposits. Furthermore, geologic processes and perturbations to the ocean-atmosphere system (e.g., warming temperatures) can cause not only dissociation, but also more widespread dissolution of hydrate or even formation of new hydrate in reservoirs. Linkages between gas hydrate and disparate aspects of Earth’s near-surface physical, chemical, and biological systems render an assessment of the rates and processes affecting the persistence of gas hydrate an appropriate Centennial Grand Challenge. This paper reviews the thermodynamic controls on methane hydrate stability and then describes the relative importance of kinetic, mass transfer, and heat transfer processes in the formation and breakdown (dissociation and dissolution) of gas hydrate. Results from numerical modeling, laboratory, and some fields studies are used to summarize the rates of hydrate formation and breakdown, followed by an extensive treatment of hydrate dynamics in marine and cryospheric gas hydrate systems.

Journal of Geophysical Research-Solid Earth↗
USGS2024

GRACE and GRACE Follow-On gravity observations of intermediate-depth earthquakes contrasted with those of shallow events

Earthquakes involve mass redistribution within the solid Earth and the ocean, and as a result, perturb the Earth's gravitational field. For most of the shallow (<60 km) earthquakes with M w > 8.0, the GRACE satellite gravity measurements suggest considerable volumetric disturbance of rocks. At a spatial scale of hundreds of km, the effect of volumetric change exceeds gravity change by vertical deformation; for example, negative gravity anomalies associated with volumetric expansion are characteristic patterns after shallow thrust events. In this study, however, we report contrasting observations of gravity change from two intermediate-depth (100–150 km) earthquakes of 2016 & 2017 Mw 8.0 (two combined) Papua New Guinea thrust faulting events and 2019 Mw 8.0 Peru normal faulting and highlight the importance of compressibility in earthquake deformation. The combined 2016/17 thrust events resulted in a positive gravity anomaly of 5–6 microGal around the epicenter, while the 2019 normal faulting produced a negative gravity anomaly of 3–4 microGal. Our modeling found that these gravity changes are manifestation of vertical deformation with limited volumetric change, distinct from gravity changes after the shallow earthquakes. The stronger resistance of rocks to volume change at intermediate-depth results in largely incompressible deformation and thus in a gravity change dominated by vertical deformation. In addition, malleable rocks under high pressure and temperature at depth facilitated substantial afterslip and/or fast viscoelastic relaxation causing additional vertical deformation and gravity change equivalent to the coseismic change. For the Papua New Guinea events, this means that postseismic relaxation enhanced coseismic uplift and relative sea level decrease.

Journal of Geophysical Research: Solid Earth↗
USGS2006

Speculations on continental crustal evolution

The evolution of the continental crust is a topic that has challenged Earth scientists since the earliest hypotheses of crustal evolution were put forth by such luminaries as Hutton, the 18th century Scottish scientist, and later by Stille (Germany), Argand (France), and Dana (United States). Recent geophysical observations provide important constraints on hypotheses of crustal evolution, and the most important of these observations are reviewed in a companion paper [ Mooney and Meissner , 1991], henceforth referred to as Paper 1. In this article we briefly speculate on crustal evolution using both geological and geophysical data as guidelines. For the past 25 years, the basic framework for models of the evolution of the Earth's crust and lithosphere has been plate tectonics. This framework has been particularly successful in explaining the processes that form and modify the oceanic crust but has had somewhat more limited success in its application to the continental crust. Some of the basic tenets of plate tectonics, as applied to the continental crust, are listed below. Continental masses coalesce and disperse as ocean basins open and close in the Wilson cycle.

Eos, Earth and Space Science News↗
USGS1997

Correlated errors in geodetic time series: Implications for time-dependent deformation

Analysis of frequent trilateration observations from the two-color electronic distance measuring networks in California demonstrate that the noise power spectra are dominated by white noise at higher frequencies and power law behavior at lower frequencies. In contrast, Earth scientists typically have assumed that only white noise is present in a geodetic time series, since a combination of infrequent measurements and low precision usually preclude identifying the time-correlated signature in such data. After removing a linear trend from the two-color data, it becomes evident that there are primarily two recognizable types of time-correlated noise present in the residuals. The first type is a seasonal variation in displacement which is probably a result of measuring to shallow surface monuments installed in clayey soil which responds to seasonally occurring rainfall; this noise is significant only for a small fraction of the sites analyzed. The second type of correlated noise becomes evident only after spectral analysis of line length changes and shows a functional relation at long periods between power and frequency of 1/ƒ α , where ƒ is frequency and α≈2. With α=2, this type of correlated noise is termed random-walk noise, and its source is mainly thought to be small random motions of geodetic monuments with respect to the Earth's crust, though other sources are possible. Because the line length changes in the two-color networks are measured at irregular intervals, power spectral techniques cannot reliably estimate the level of 1/ƒ α noise. Rather, we also use here a maximum likelihood estimation technique which assumes that there are only two sources of noise in the residual time series (white noise and random-walk noise) and estimates the amount of each. From this analysis we find that the random-walk noise level averages about 1.3 mm/√yr and that our estimates of the white noise component confirm theoretical limitations of the measurement technique. In addition, the seasonal noise can be as large as 3 mm in amplitude but typically is less than 0.5 mm. Because of the presence of random-walk noise in these time series, modeling and interpretation of the geodetic data must account for this source of error. By way of example we show that estimating the time-varying strain tensor (a form of spatial averaging) from geodetic data having both random-walk and white noise error components results in seemingly significant variations in the rate of strain accumulation; spatial averaging does reduce the size of both noise components but not their relative influence on the resulting strain accumulation model.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Photogrammetric portrayal of Mars topography

Special photogrammetric techniques have been developed to portray Mars topography, using Mariner and Viking imaging and nonimaging topographic information and earth-based radar data. Topography is represented by the compilation of maps at three scales: global, intermediate, and very large scale. The global map is a synthesis of topographic information obtained from Mariner 9 and earth-based radar, compiled at a scale of 1:25,000,000 with a contour interval of 1 km; it gives a broad quantitative view of the planet. At intermediate scales, Viking Orbiter photographs of various resolutions are used to compile detailed contour maps of a broad spectrum of prominent geologic features; a contour interval as small as 20 m has been obtained from very high resolution orbital photography. Imagery from the Viking lander facsimile cameras permits construction of detailed, very large scale (1:10) topographic maps of the terrain surrounding the two landers; these maps have a contour interval of 1 cm. This paper presents several new detailed topographic maps of Mars.

Journal of Geophysical Research Solid Earth↗
USGS2012

Eolian features in the Western Desert of Egypt and some applications to Mars

Relations of landform types to wind regimes, bedrock composition, sediment supply, and topography are shown by field studies and satellite photographs of the Western Desert of Egypt. This desert, which lies at the core of the largest hyperarid region on earth, provides analogs of Martian wind-formed features. These include sand dunes, alternating light and dark streaks, knob ‘shadows,’ and yardangs. Surface particles have been segregated by wind into deposits (dunes, sand sheets, and light streaks) that can be differentiated by their grain size distributions, surface shapes, and colors. Throughgoing sand of mostly fine to medium grain size is migrating southward in longitudinal dune belts and barchan chains whose long axes lie parallel to the prevailing northerly winds, but topographic variations such as scarps and depressions strongly influence the zones of deposition and dune morphology. Sand from the longitudinal dunes on the plains is commonly redistributed into barchans in the depressions. These barchans are generally simple crescents that are morphologically similar to many of the dunes seen on Viking orbiter pictures of the north polar sand sea on mars. Light streaks are depositional features consisting of dune belts and elongate sheets of coarse to medium sand and granules. Intervening dark streaks are erosional features consisting of strips of desert-varnished bedrock and lag gravel surfaces exposed between the sand deposits. The shape of both light and dark streaks is controlled by wind flow around topographic highs. Dark zones (shadows) in the lee of mountains, hills, and knobs are erosional products from the topographic highs; they change shape only in response to movement of the adjacent lighter-colored sand deposits. Streamlined yardangs carved in crystalline limestone constitute one of the largest yardang fields on earth. Yardangs occur also in sandstone of the Nubian Series and in lacustrine sediments. The variables that affect the patterns of wind erosion and deposition in the Western Desert are topographic effects on wind velocities and directions, resistance of the bedrock, sand supply, and climatic change with time; vegetation is essentially absent and is not a controlling factor.

Journal of Geophysical Research Solid Earth↗
USGS2012

Ice sculpture in the Martian outflow channels

Many landforms in Martian outflow channels have characteristics that suggest sculpture by glaciers, ice streams, or ice sheets. Viking Orbiter and terrestrial satellite images were examined at similar resolution to compare features of the Martian outflow channels to features produced by the movement of ice on earth. Many resemblances were found. They include the anastomoses, sinuosities, and U‐shaped cross profiles of valleys; hanging valleys; and linear scour marks on valley walls, grooves and ridges on valley floors, and streamlining of bedrock highs. The question of whether ice could have moved in the Martian environment is investigated. Since gradients on Martian channel floors are extremely low or even reversed, a certain thickness of ice is required to initiate glacial motion. If ice filled the Martian channels to the level of interior plateaus, flow could have occurred, even under present climatic conditions, although warmer climates in the past would have been more favorable. Additionally, if the channel fluids were brines, ice might have been wet based and would have moved readily by slipping over the ground. It is envisioned that springs or small catastrophic outbursts discharged fluids from structural outlets or chaotic terrains. The fluids built icings that may have grown into substantial masses and eventually flowed like glaciers down preexisting valleys. Alternatively, the fluids may have formed rivers or floods that formed ice jams and consolidated into icy masses in places where obstacles blocked their flow. These masses of ice or slush were probably episodic and temporary, but where they occurred within the channels, they may have produced landforms that are characteristic of glacial sculpture on earth.

Journal of Geophysical Research B: Solid Earth↗