Search USGSSearch

SEARCH · Search USGS

Results for “Journal of Maps”

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,045 records · Page 58Linked to original sources

Eruptive history and tectonic setting of Medicine Lake Volcano, a large rear-arc volcano in the southern Cascades

Medicine Lake Volcano (MLV), located in the southern Cascades ∼ 55 km east-northeast of contemporaneous Mount Shasta, has been found by exploratory geothermal drilling to have a surprisingly silicic core mantled by mafic lavas. This unexpected result is very different from the long-held view derived from previous mapping of exposed geology that MLV is a dominantly basaltic shield volcano. Detailed mapping shows that < 6% of the ∼ 2000 km 2 of mapped MLV lavas on this southern Cascade Range shield-shaped edifice are rhyolitic and dacitic, but drill holes on the edifice penetrated more than 30% silicic lava. Argon dating yields ages in the range ∼ 475 to 300 ka for early rhyolites. Dates on the stratigraphically lowest mafic lavas at MLV fall into this time frame as well, indicating that volcanism at MLV began about half a million years ago. Mafic compositions apparently did not dominate until ∼ 300 ka. Rhyolite eruptions were scarce post-300 ka until late Holocene time. However, a dacite episode at ∼ 200 to ∼ 180 ka included the volcano's only ash-flow tuff, which was erupted from within the summit caldera. At ∼ 100 ka, compositionally distinctive high-Na andesite and minor dacite built most of the present caldera rim. Eruption of these lavas was followed soon after by several large basalt flows, such that the combined area covered by eruptions between 100 ka and postglacial time amounts to nearly two-thirds of the volcano's area. Postglacial eruptive activity was strongly episodic and also covered a disproportionate amount of area. The volcano has erupted 9 times in the past 5200 years, one of the highest rates of late Holocene eruptive activity in the Cascades. Estimated volume of MLV is ∼ 600 km 3 , giving an overall effusion rate of ∼ 1.2 km 3 per thousand years, although the rate for the past 100 kyr may be only half that. During much of the volcano's history, both dry HAOT (high-alumina olivine tholeiite) and hydrous calcalkaline basalts erupted together in close temporal and spatial proximity. Petrologic studies indicate that the HAOT magmas were derived by dry melting of spinel peridotite mantle near the crust mantle boundary. Subduction-derived H 2 O-rich fluids played an important role in the generation of calcalkaline magmas. Petrology, geochemistry and proximity indicate that MLV is part of the Cascades magmatic arc and not a Basin and Range volcano, although Basin and Range extension impinges on the volcano and strongly influences its eruptive style. MLV may be analogous to Mount Adams in southern Washington, but not, as sometimes proposed, to the older distributed back-arc Simcoe Mountains volcanic field.

Journal of Volcanology and Geothermal Research

Petrology and geochronology of Mesoproterozoic basement of the Mount Rogers area of southwestern Virginia and northwestern North Carolina: Implications for the Precambrian tectonic evolution of the southern Blue Ridge province

Results from new geologic mapping, SHRIMP U-Pb geochronology, and petrologic studies indicate that Mesoproterozoic basement in the northern French Broad massif near Mount Rogers consists of multiple, mostly granitic plutons, map- and outcrop-scale xenoliths of pre-existing crustal rocks, and remnants of formerly overlying meta-sedimentary lithologies. Zircon and titanite ages demonstrate that these rocks collectively record nearly 350 m.y. of tectonic evolution including periods of igneous intrusion at ca . 1190 to 1130 Ma (Early Magmatic Suite) and ca . 1075 to 1030 Ma (Late Magmatic Suite) and three episodes of regional metamorphism at ca . 1170 to 1140, 1070 to 1020, and 1000 to 970 Ma. The existence of ca . 1.3 Ga age crust is indicated by (1) orthogranofels of ca . 1.32 Ga age in a map-scale xenolith, (2) inherited zircons of ca . 1.33 to 1.29 Ga age in Early Magmatic Suite plutons, and (3) ca . 1.36 to 1.30 Ga age detrital zircons in meta-sedimentary lithologies. Mineral assemblages developed in amphibolites and granofelses indicate that metamorphism during both Mesoproterozoic episodes occurred at upper amphibolite- to lower granulite-facies conditions. Syn-orogenic Early Magmatic Suite plutons emplaced at ca . 1190 to 1145 Ma are characterized by high-K, variably magnesian, dominantly calc-alkalic compositions, and have trace-element characteristics indicative of continental-arc magmatic origin involving melting of thick continental crust. In contrast, ca . 1140 Ma age quartz syenite displays A-type features indicating derivation from depleted crustal sources with increased mantle input during waning stages of regional contraction. Plutons of the compositionally bimodal Late Magmatic Suite include (1) ca . 1060 Ma meta-granite with geochemical characteristics transitional between silicic rocks of arc systems and post-collisional granites of A-type lineage, and (2) ca . 1055 Ma monzodioritic rocks with A-type compositional characteristics that likely reflect derivation from fertile, mafic sources in the lower crust. Collectively, these data suggest that Mesoproterozoic rocks of the study area preserve evidence of multiple orogenic episodes that likely involved continental-arc development and deformation at ca . 1150 Ma followed by crustal thickening at ca . 1060 Ma. Field relations and geochronologic data indicate that regional uplift and sedimentation occurred at ca . 1100 Ma between the two episodes of overlapping magmatism and orogenesis. The nature and timing of Mesoproterozoic events recorded in basement rocks of the study area illustrate significant differences in the lithologic assemblages and geologic history preserved by Mesoproterozoic basement of the adjacent Shenandoah and French Broad massifs, suggesting that the Blue Ridge massifs occupied different locations within the regional Grenville-age orogen until about 1070 Ma when the effects of Ottawan-age tectonics began to affect both areas. The near ubiquity of Ottawan-age orogenic activity recorded in Mesoproterozoic rocks of the Blue Ridge, other Appalachian inliers such as the New Jersey Highlands, and the Grenville province of Canada, including the Adirondacks, suggests that these formerly disparate terranes were amalgamated to form a common, regional orogen by this time.

North Carolina, Virginia

Sea-floor environments within Long Island Sound: A regional overview

Modern sea-floor sedimentary environments within the glaciated, topographically complex Long Island Sound estuary have been interpreted and mapped from an extensive collection of sidescan sonographs, bottom samples, and video-camera observations together with supplemental bathymetric, marine-geologic, and bottom-current data. Four categories of environments are present that reflect the dominant long-term processes of erosion or nondeposition; coarsegrained bedload transport; sediment sorting and reworking; and fine-grained deposition. (1) Environments of erosion or nondeposition contain exposures of glacial drift, coarse lag deposits, and possibly bedrock and include sediments which range from boulder fields to gravelly coarse-to-medium sands. (2) Environments of coarse-grained bedload transport are mantled by sand ribbons and sand waves and contain mostly coarse-to-fine sands with only small amounts of mud. (3) Environments of sediment sorting and reworking comprise both uniform and heterogeneous sediment types and contain variable amounts of fine sand and mud. (4) Environments of fine-grained deposition are blanketed by muds and sandy muds. The patchy distribution of sedimentary environments within Long Island Sound reflects both regional and local changes in bottom processes. Regional changes are primarily the result of a strong, east-to-west decreasing gradient of bottom tidal-current speeds, coupled with the net (westward) estuarine bottom drift. The regional current regime has produced a westward succession of environments along the basin floor beginning with erosion or nondeposition at the narrow eastern entrance to the Sound, changing to an extensive area of coarse-grained bedload transport, passing into a contiguous band of sediment sorting, and ending with broad areas of fine-grained deposition in the central and western Sound. However, local changes in processes are superimposed on the regional conditions within the central and western parts of the basin and along the nearshore margins. Within the central and western basin, localized sedimentary environments are produced where the bottom flow is enhanced by, and interacts with, the bottom topography, whereas along the nearshore margins, they variously reflect wave-produced currents, the irregular bathymetry, the indented shoreline, and the proximal supply of sediments. Results from this study (1) confirm the high trapping efficiency of fine-grained sediments in the Sound, (2) suggest that fine-grained sediments accumulate at an average (regional) rate of 0.08 g/cm 2 /y , and (3) indicate that the postglacial delta in the eastern Sound was a significant source of fine-grained sediments now buried beneath depositional areas.

Connecticut, New York

Upper crustal structure in Puget Lowland, Washington: Results from the 1998 Seismic Hazards Investigation in Puget Sound

A new three-dimensional (3-D) model shows seismic velocities beneath the Puget Lowland to a depth of 11 km. The model is based on a tomographic inversion of nearly one million first-arrival travel times recorded during the 1998 Seismic Hazards Investigation in Puget Sound (SHIPS), allowing higher-resolution mapping of subsurface structures than previously possible. The model allows us to refine the subsurface geometry of previously proposed faults (e.g., Seattle, Hood Canal, southern Whidbey Island, and Devils Mountain fault zones) as well as to identify structures (Tacoma, Lofall, and Sequim fault zones) that warrant additional study. The largest and most important of these newly identified structures lies along the northern boundary of the Tacoma basin; we informally refer to this structure here as the Tacoma fault zone. Although tomography cannot provide information on the recency of motion on any structure, Holocene earthquake activity on the Tacoma fault zone is suggested by seismicity along it and paleoseismic evidence for abrupt uplift of tidal marsh deposits to its north. The tomography reveals four large, west to northwest trending low-velocity basins (Tacoma, Seattle, Everett, and Port Townsend) separated by regions of higher velocity ridges that are coincident with fault-bounded uplifts of Eocene Crescent Formation basalt and pre-Tertiary basement. The shapes of the basins and uplifts are similar to those observed in gravity data; gravity anomalies calculated from the 3-D tomography model are in close agreement with the observed anomalies. In velocity cross sections the Tacoma and Seattle basins are asymmetric: the basin floor dips gently toward a steep boundary with the adjacent high-velocity uplift, locally with a velocity "overhang" that suggests a basin vergent thrust fault boundary. Crustal fault zones grow from minor folds into much larger structures along strike. Inferred structural relief across the Tacoma fault zone increases by several kilometers westward along the fault zone to Lynch Cove, where we interpret it as a zone of south vergent faulting overthrusting Tacoma basin. In contrast, structural relief along the Seattle fault zone decreases west of Seattle, which we interpret as evidence that the N-S directed compression is being accommodated by slip transfer between the Seattle and Tacoma fault zones. Together, the Tacoma and Seattle fault zones raise the Seattle uplift, one of a series of east-west trending, pop-up structures underlying Puget Lowland from the Black Hills to the San Juan Islands.

Washington

Shear wave velocity structure in North America from large-scale waveform inversions of surface waves

A two-step nonlinear and linear inversion is carried out to map the lateral heterogeneity beneath North America using surface wave data. The lateral resolution for most areas of the model is of the order of several hundred kilometers. The most obvious feature in the tomographic images is the rapid transition between low velocities in the tectonically active region west of the Rocky Mountains and high velocities in the stable central and eastern shield of North America. The model also reveals smaller-scale heterogeneous velocity structures. A high-velocity anomaly is imaged beneath the state of Washington that could be explained as the subducting Juan de Fuca plate beneath the Cascades. A large low-velocity structure extends along the coast from the Mendocino to the Rivera triple junction and to the continental interior across the southwestern United States and northwestern Mexico. Its shape changes notably with depth. This anomaly largely coincides with the part of the margin where no lithosphere is consumed since the subduction has been replaced by a transform fault. Evidence for a discontinuous subduction of the Cocos plate along the Middle American Trench is found. In central Mexico a transition is visible from low velocities across the Trans-Mexican Volcanic Belt (TMVB) to high velocities beneath the Yucatan Peninsula. Two elongated low-velocity anomalies beneath the Yellowstone Plateau and the eastern Snake River Plain volcanic system and beneath central Mexico and the TMVB seem to be associated with magmatism and partial melting. Another low-velocity feature is seen at depths of approximately 200 km beneath Florida and the Atlantic Coastal Plain. The inversion technique used is based on a linear surface wave scattering theory, which gives tomographic images of the relative phase velocity perturbations in four period bands ranging from 40 to 150 s. In order to find a smooth reference model a nonlinear inversion based on ray theory is first performed. After correcting for the crustal thickness the phase velocity perturbations obtained from the subsequent linear waveform inversion for the different period bands are converted to a three-layer model of S velocity perturbations (layer 1, 25–100 km; layer 2, 100–200 km; layer 3, 200–300 km). We have applied this method on 275 high-quality Rayleigh waves recorded by a variety of instruments in North America (IRIS/USGS, IRIS/IDA, TERRAscope, RSTN). Sensitivity tests indicate that the lateral resolution is especially good in the densely sampled western continental United States, Mexico, and the Gulf of Mexico.

Journal of Geophysical Research B: Solid Earth

Distribution of late Cenozoic volcanic vents in the Cascade Range: Volcanic arc segmentation and regional tectonic considerations

Spatial, temporal, and compositional distributions of approximately 4000 volcanic vents formed since 16 Ma in Washington, Oregon, northern California, and northwestern Nevada illustrate the evolution of volcanism related to subduction of the Juan de Fuca plate system and extension of the Basin and Range province. Vent data were obtained from published map compilations and include monogenetic and small polygenetic volcanoes in addition to major composite centers. On the basis of the distribution of 2821 vents formed since 5 Ma, the Cascade Range is divided into five segments, with vents of the High Lava Plains along the northern margin of the Basin and Range province in Oregon forming a sixth segment. Some aspects of the Cascade Range segmentation can be related to gross structural features of the subducting Juan de Fuca plate. The orientation of the volcanic front of segments one and two changes from NW in northern Washington to NE in southern Washington, paralleling the strike of the subducting Juan de Fuca plate. Segments one and two are separated by a 90-km volcanic gap between Mount Rainier and Glacier Peak that is landward of the portion of the subducting plate having the least average dip to a depth of 60 km. A narrow, N-S trending belt of predominantly andesitic vents in Oregon constitutes a third segment, which is landward of the seismically quiet portion of the subduction zone. The narrowness of this segment may indicate steep dip of the subducting plate beneath the Cascade arc in Oregon. Vents are sparse between segment four (containing the Mount Shasta and Medicine Lake centers) and segment five (containing Lassen Peak), where the Juan de Fuca and Gorda North plates are characterized by differing age, amounts of subcrustal seismicity, and probably geometry. From the relation between seismicity at depth of 60 km and the position of the volcanic front of vents formed since 5 Ma, transitions between subducting-plate segments of varying geometry likely occur near boundaries between independently defined volcanic segments in northern Oregon and northern California. In the Basin and Range province east of the Cascade arc, volcanism migrated into the region adjacent to the Cascade Range during the interval 5–10 Ma. Since 5 Ma, the impingement of the two provinces is characterized by cessation of basin-range volcanism in southern Oregon, continuation of basaltic volcanism in northeastern California where the impingement process may not yet be complete, and contraction of the area of mafic volcanism around Mount Shasta, Medicine Lake, and Lassen Peak. In central Oregon where the northern margin of basin-range volcanism (the High Lava Plains) intersects the Cascade arc, impingement of basin-range extensional volcanism approximately coincides in time and space with the development of the High Cascade graben between Three Sisters and Mount Jefferson.

Journal of Geophysical Research Solid Earth

Precaldera mafic magmatism at Long Valley, California: Magma-tectonic siting and incubation of the Great Rhyolite System

The iconic volcanic center at Long Valley has released ∼820 km 3 of rhyolite in at least 110 eruptions. From 2.2 Ma until 0.23 Ma, products were exclusively rhyolitic, and ∼ 700 km 3 were high-silica rhyolite severely depleted in Sr, Ba, and Eu. The rhyolitic interval was preceded by an interval from 3.9 to 2.6 Ma with numerous basalt-andesite-dacite eruptions accompanied by no rhyolite at all. We have now mapped the circumcaldera products of this interval, defined 107 eruptive units, characterized them all chemically and petrographically, and dated many by 40 Ar/ 39 Ar. Here we display and describe them by sector around the caldera, interpret the nature of the transcrustal magma system that eventuated in the 35-km-wide Long Valley granite-rhyolite pluton, and analyze regional tectonic factors that did or did not contribute to siting the system. Nine Miocene (12–6 Ma) eruptive units close to Long Valley were followed by a Pliocene flare-up that released >300 mafic eruptions in a SW–NE swath 170 km long, centered across the later site of Long Valley. The basalts and their fractionates are intraplate alkalic products dominated by a continental lithosphere that had long been fluxed by Mesozoic subduction. Tertiary arc volcanism had not impinged on the area of Long Valley. Volumes estimated for the 107 Neogene precaldera eruptive units (only 40 of which exceeded 0.1 km 3 ) total ∼ 27 km 3 ± 50%—only ∼3% of the subsequent volume of rhyolite erupted. Such a volume of high-silica rhyolite with ultra-low Sr and Eu is not a product of partial melting but requires as proximate parent a leucogranitic crystal mush that is itself the upper level of a long-lived plutonic reservoir that extends to the lower crust. The 27 km 3 of Neogene magma that erupted was a small contingent of the mantle-derived basaltic flux needed to energize (and contribute its fractionated melt to) a 30-km-deep compositionally graded crustal column, which culminated in ∼10,000 km 3 of granitoid mush from which 820 km 3 of Quaternary high-silica rhyolitic melt escaped and erupted. Pliocene basaltic eruptions ceased at ∼2.6 Ma, probably because the basaltic flux intensified sufficiently to render the mushy upper crust impenetrable. The 2.6–2.2 Ma quiescent interval represented culmination of thermal activation of the plutonic column and refinement of its leucogranitic mushy upper layer, from which extreme melts escaped for the next 2 Myr. The Pliocene mafic swath crosses the Sierran rangefront fault zone coincident with a left-stepping extensional reentrant that also began developing at ∼3 Ma. Moreover, Long Valley overlies a dextral offset in the initial Sr-isotope 0.706 line, which may reflect the rifted or attenuated edge of Proterozoic crust and mantle lithosphere. Concatenation of these three influences may account for siting of intensified edge-focused magmatism that produced the great Quaternary pluton.

California

Topography of Valles Marineris: Implications for erosional and structural history

Compilation of a simplified geologic/geomorphic map onto digital terrain models of the Valles Marineris permitted an evaluation of elevations in the vicinity of the troughs and the calculation of depth of troughs below surrounding plateaus, thickness of deposits inside the troughs, volumes of void spaces above geologic/geomorphic units, and volumes of deposits. The central troughs north Ophir, north and central Candor, and north Melas Chasmata lie as much as 11 km below the adjacent plateaus. In Ophir and Candor Chasmata, interior layered deposits reach 8 km in elevation. If the deposits are lacustrine and if all troughs were interconnected, lake waters standing 8 km high would have spilled out of Coprates Chasma onto the surrounding plateaus having surface elevations of only 4–5 km. In this case, interior deposits above about 4 km in the central troughs would not be lacustrine. They could be volcanic. On the other hand, the troughs may not have been interconnected at the time of interior‐deposit emplacement; they may have formed isolated ancestral basins. The existence of such basins is supported by independent structural and stratigraphic evidence. The ancestral basins may have eventually merged, perhaps through renewed faulting, to form northern subsidiary troughs in Ophir and Candor Chasmata and the Coprates/north Melas/Ius graben system. The peripheral troughs are only 2–5 km deep, shallower than the central troughs. They may have formed from a combination of erosional collapse and structural activity. Chaotic terrain is seen in the peripheral troughs near a common contour level of about 4 km on the adjacent plateaus, which supports the idea of release of water under artesian pressure from confined aquifers. The layered deposits in the peripheral troughs may have formed in isolated depressions that harbored lakes and predated the formation of the deep outflow channels. If these layered deposits are of volcanic origin, they may have been emplaced beneath ice in the manner of table mountains. Areal and volumetric computations show that erosion widened the troughs by about one‐third and that deposits occupy one‐sixth of the interior space. Even though the volume eroded is larger than the volume deposited, topographic and geologic considerations imply that material eroded from trough walls was probably part of the interior layered deposits but not their sole source. Additional material may have come from subterranean piping, from reworking of local disintegration products on the floors, such as chaotic materials, or from eolian influx. But overall it is likely that the additional material is volcanic and that it forms mostly the upper, more diversely bedded layers of the interior deposits.

Journal of Geophysical Research E: Planets

Origin of the Lyme Dome and implications for the timing of multiple Alleghanian deformational and intrusive events in southern Connecticut

Geologic mapping, structural analysis, and geochronology in the area of the Lyme dome, southern Connecticut provides constraints on the origin of the rocks in the core of the dome, the absolute timing of the principal deformational and thermal events attributed to Alleghanian orogenesis, and the processes that generated the dome. Detrital zircon geochronology in combination with ages on intrusive rocks brackets the deposition of quartzite in the core of the dome sometime between ca. 925 and 620 Ma. Granite and granodiorite intruded the Neoproteorozic metasedimentary rocks in the core of the dome at ca. 620 to 610 Ma. Four major early Permian events associated with the Alleghanian orogeny affected the rocks in the Lyme dome area. Syn-tectonic migmatization and widespread penetrative deformation (D1, ca. 300 - 290 Ma) included emplacement of alaskite at 290 ?? 4 Ma during regional foliation development and aluminosilicate-orthoclase metamorphic conditions. Rocks of the Avalon terrane may have wedged between Gander cover rocks and Gander basement in the core of the Lyme during D1. Limited structural evidence for diapiric uplift of the Lyme dome indicates that diapirism started late in D1 and was completed by D2 (ca. 290 - 280 Ma) when horizontal WNW contractional stresses dominated over vertical stresses. Second sillimanite metamorphism continued and syn-tectonic D2 granite pegmatite (288 ?? 4 Ma) and the Joshua Rock Granite Gniess (284 ?? 3 Ma) intruded at this time. North-northwest extension during D3 (ca. 280 - 275 Ma) led to granitic pegmatite intrusion along S3 cleavage planes and in extensional zones in boudin necks during hydraulic failure and decompression melting. Intrusion of a Westerly Granite dike at 275 ?? 4 Ma suggests that D3 extension was active, and perhaps concluding, by ca. 275 Ma. Late randomly oriented but gently dipping pegmatite dikes record a final stage of intrusion during D4 (ca. 275 - 260 Ma), and a switch from NNW extension to vertical unloading and exhumation. Monazite and metamorphic zircon rim ages record this event at ca. 259 Ma. The evolution of the Lyme dome involved D1 mylonitization, intrusion, and migmatization during north-directed contraction, limited late D1 diapirism, D2 migmatization during WNW contraction with associated flexural flow and fold interference, D3 NNW horizontal extension and decompression melting, and final D4 vertical extension and rapid exhumation. Late regional uplift, extension, and normal faulting at higher crustal levels may have been caused by diapiric rise of the lower crust, below the structural level of the Lyme dome. The rocks record no evidence of Acadian metamorphism or deformation, suggesting that the Gander zone here was not tectonically juxtaposed with Avalon until the Alleghanian orogeny.

American Journal of Science

Automated cropland mapping of continental Africa using Google Earth Engine cloud computing

The automation of agricultural mapping using satellite-derived remotely sensed data remains a challenge in Africa because of the heterogeneous and fragmental landscape, complex crop cycles, and limited access to local knowledge. Currently, consistent, continent-wide routine cropland mapping of Africa does not exist, with most studies focused either on certain portions of the continent or at most a one-time effort at mapping the continent at coarse resolution remote sensing. In this research, we addressed these limitations by applying an automated cropland mapping algorithm (ACMA) that captures extensive knowledge on the croplands of Africa available through: (a) ground-based training samples, (b) very high (sub-meter to five-meter) resolution imagery (VHRI), and (c) local knowledge captured during field visits and/or sourced from country reports and literature. The study used 16-day time-series of Moderate Resolution Imaging Spectroradiometer (MODIS) normalized difference vegetation index (NDVI) composited data at 250-m resolution for the entire African continent. Based on these data, the study first produced accurate reference cropland layers or RCLs (cropland extent/areas, irrigation versus rainfed, cropping intensities, crop dominance, and croplands versus cropland fallows) for the year 2014 that provided an overall accuracy of around 90% for crop extent in different agro-ecological zones (AEZs). The RCLs for the year 2014 (RCL2014) were then used in the development of the ACMA algorithm to create ACMA-derived cropland layers for 2014 (ACL2014). ACL2014 when compared pixel-by-pixel with the RCL2014 had an overall similarity greater than 95%. Based on the ACL2014, the African continent had 296 Mha of net cropland areas (260 Mha cultivated plus 36 Mha fallows) and 330 Mha of gross cropland areas. Of the 260 Mha of net cropland areas cultivated during 2014, 90.6% (236 Mha) was rainfed and just 9.4% (24 Mha) was irrigated. Africa has about 15% of the world’s population, but only about 6% of world’s irrigation. Net cropland area distribution was 95 Mha during season 1, 117 Mha during season 2, and 84 Mha continuous. About 58% of the rainfed and 39% of the irrigated were single crops (net cropland area without cropland fallows) cropped during either season 1 (January-May) or season 2 (June-September). The ACMA algorithm was deployed on Google Earth Engine (GEE) cloud computing platform and applied on MODIS time-series data from 2003 through 2014 to obtain ACMA-derived cropland layers for these years (ACL2003 to ACL2014). The results indicated that over these twelve years, on average: (a) croplands increased by 1 Mha/yr, and (b) cropland fallows decreased by 1 Mha/year. Cropland areas computed from ACL2014 for the 55 African countries were largely underestimated when compared with an independent source of census-based cropland data, with a root-mean-square error (RMSE) of 3.5 Mha. ACMA demonstrated the ability to hind-cast (past years), now-cast (present year), and forecast (future years) cropland products using MODIS 250-m time-series data rapidly, but currently, insufficient reference data exist to rigorously report trends from these results.

ISPRS Journal of Photogrammetry and Remote Sensing

Mars Global Surveyor Thermal Emission Spectrometer experiment: Investigation description and surface science results

The Thermal Emission Spectrometer (TES) investigation on Mars Global Surveyor (MGS) is aimed at determining (1) the composition of surface minerals, rocks, and ices; (2) the temperature and dynamics of the atmosphere; (3) the properties of the atmospheric aerosols and clouds; (4) the nature of the polar regions; and (5) the thermophysical properties of the surface materials. These objectives are met using an infrared (5.8- to 50-μm) interferometric spectrometer, along with broadband thermal (5.1- to 150-μm) and visible/near-IR (0.3- to 2.9-μm) radiometers. The MGS TES instrument weighs 14.47 kg, consumes 10.6 W when operating, and is 23.6×35.5×40.0 cm in size. The TES data are calibrated to a 1-σ precision of 2.5 −6 ×10 −8 W cm −2 sr −1 /cm −1 , 1.6×10 −6 W cm −2 sr −1 , and ∼0.5 K in the spectrometer, visible/near-IR bolometer, and IR bolometer, respectively. These instrument subsections are calibrated to an absolute accuracy of ∼4×10 −8 W cm −2 sr −1 /cm −1 (0.5 K at 280 K), 1–2%, and ∼1–2 K, respectively. Global mapping of surface mineralogy at a spatial resolution of 3 km has shown the following: (1) The mineralogic composition of dark regions varies from basaltic, primarily plagioclase feldspar and clinopyroxene, in the ancient, southern highlands to andesitic, dominated by plagioclase feldspar and volcanic glass, in the younger northern plains. (2) Aqueous mineralization has produced gray, crystalline hematite in limited regions under ambient or hydrothermal conditions; these deposits are interpreted to be in-place sedimentary rock formations and indicate that liquid water was stable near the surface for a long period of time. (3) There is no evidence for large-scale (tens of kilometers) occurrences of moderate-grained (>50-μm) carbonates exposed at the surface at a detection limit of ∼10%. (4) Unweathered volcanic minerals dominate the spectral properties of dark regions, and weathering products, such as clays, have not been observed anywhere above a detection limit of ∼10%; this lack of evidence for chemical weathering indicates a geologic history dominated by a cold, dry climate in which mechanical, rather than chemical, weathering was the significant form of erosion and sediment production. (5) There is no conclusive evidence for sulfate minerals at a detection limit of ∼15%. The polar region has been studied with the following major conclusions: (1) Condensed CO 2 has three distinct end-members, from fine-grained crystals to slab ice. (2) The growth and retreat of the polar caps observed by MGS is virtually the same as observed by Viking 12 Martian years ago. (3) Unique regions have been identified that appear to differ primarily in the grain size of CO 2 ; one south polar region appears to remain as black slab CO 2 ice throughout its sublimation. (4) Regional atmospheric dust is common in localized and regional dust storms around the margin and interior of the southern cap. Analysis of the thermophysical properties of the surface shows that (1) the spatial pattern of albedo has changed since Viking observations, (2) a unique cluster of surface materials with intermediate inertia and albedo occurs that is distinct from the previously identified low-inertia/bright and high-inertia/dark surfaces, and (3) localized patches of high-inertia material have been found in topographic lows and may have been formed by a unique set of aeolian, fluvial, or erosional processes or may be exposed bedrock.

Journal of Geophysical Research E: Planets

Poroelastic rebound along the Landers 1992 earthquake surface rupture

Maps of surface displacement following the 1992 Landers, California, earthquake, generated by interferometric processing of ERS-1 synthetic aperture radar (SAR) images, reveal effects of various postseismic deformation processes along the 1992 surface rupture. The large-scale pattern of the postseismic displacement field includes large lobes, mostly visible on the west side of the fault, comparable in shape with the lobes observed in the coseismic displacement field. This pattern and the steep displacement gradient observed near the Emerson-Camp Rock fault cannot be simply explained by afterslip on deep sections of the 1992 rupture. Models show that horizontal slip occurring on a buried dislocation in a Poisson's material produces a characteristic quadripole pattern in the surface displacement field with several centimeters of vertical motion at distances of 10–20 km from the fault, yet this pattern is not observed in the postseismic interferograms. As previously proposed to explain local strain in the fault step overs [ Peltzer et al ., 1996b], we argue that poroelastic rebound caused by pore fluid flow may also occur over greater distances from the fault, compensating the vertical ground shift produced by fault afterslip. Such a rebound is explained by the gradual change of the crustal rocks' Poisson's ratio value from undrained (coseismic) to drained (postseismic) conditions as pore pressure gradients produced by the earthquake dissipate. Using the Poisson's ratio values of 0.27 and 0.31 for the drained and undrained crustal rocks, respectively, elastic dislocation models show that the combined contributions of afterslip on deep sections of the fault and poroelastic rebound can account for the range change observed in the SAR data and the horizontal displacement measured at Global Positioning System (GPS) sites along a 60-km-long transect across the Emerson fault [ Savage and Svarc , 1997]. Using a detailed surface slip distribution on the Homestead Valley, Kickapoo, and Johnson Valley faults, we modeled the poroelastic rebound in the Homestead Valley pull apart. A Poisson's ratio value of 0.35 for the undrained gouge rocks in the fault zone is required to account for the observed surface uplift in the 3.5 years following the earthquake. This large value implies a seismic velocity ratio V p /V s of 2.1, consistent with the observed low V s values of fault zone guided waves at shallow depth [ Li et al ., 1997]. The SAR data also reveal postseismic creep along shallow patches of the Eureka Peak and Burnt Mountain faults with a characteristic decay time of 0.8 years. Coseismic, dilatant hardening (locking process) followed by post-seismic, pore pressure controlled fault creep provide a plausible mechanism to account for the decay time of the observed slip rate along this section of the fault.

Journal of Geophysical Research B: Solid Earth

Geophysical modeling of the northern Appalachian Brompton-Cameron, Central Maine, and Avalon terranes under the New Jersey Coastal Plain

A regional terrane map of the New Jersey Coastal Plain basement was constructed using seismic, drilling, gravity and magnetic data. The Brompton-Cameron and Central Maine terranes were coalesced as one volcanic island arc terrane before obducting onto Laurentian, Grenville age, continental crust in the Taconian orogeny [Rankin, D.W., 1994. Continental margin of the eastern United States: past and present. In: Speed, R.C., (Ed.), Phanerozoic Evolution of North American Continent-Ocean Transitions. DNAG Continent-Ocean Transect Volume. Geological Society of America, Boulder, Colorado, pp. 129-218]. Volcanic island-arc rocks of the Avalon terrane are in contact with Central Maine terrane rocks in southern Connecticut where the latter are overthrust onto the Brompton-Cameron terrane, which is thrust over Laurentian basement. Similarities of these allochthonous island arc terranes (Brompton-Cameron, Central Maine, Avalon) in lithology, fauna and age suggest that they are faulted segments of the margin of one major late Precambrian to early Paleozoic, high latitude peri-Gondwana island arc designated as "Avalonia", which collided with Laurentia in the early to middle Paleozoic. The Brompton Cameron, Central Maine, and Avalon terranes are projected as the basement under the eastern New Jersey Coastal Plain based on drill core samples of metamorphic rocks of active margin/magmatic arc origin. A seismic reflection profile across the New York Bight traces the gentle dipping (approximately 20 degrees) Cameron's Line Taconian suture southeast beneath allochthonous Avalon and other terranes to a 4 sec TWTT depth (approximately 9 km) where the Avalonian rocks are over Laurentian crust. Gentle up-plunge (approximately 5 degrees) projections to the southwest bring the Laurentian Grenville age basement and the drift-stage early Paleozoic cover rocks to windows in Burlington Co. at approximately 1 km depth and Cape May Co. at approximately 2 km depths. The antiformal Shellburne Falls and Chester domes and Chain Lakes-Pelham dome-Bronson Hill structural trends, and the synformal Connecticut Valley-Gaspe structural trend can be traced southwest into the New Jersey Coastal Plain basement. A Mesozoic rift basin, the "Sandy Hook basin", and associated eastern boundary fault is identified, based upon gravity modeling, in the vicinity of Sandy Hook, New Jersey. The thickness of the rift-basin sedimentary rocks contained within the "Sandy Hook basin" is approximately 4.7 km, with the basin extending offshore to the east of the New Jersey coast. Gravity modeling indicates a deep rift basin and the magnetic data indicates a shallow magnetic basement caused by magnetic diabase sills and/or basalt flows contained within the rift-basin sedimentary rocks. The igneous sills and/or flows may be the eastward continuation of the Watchung and Palisades bodies. ?? 2004 Elsevier Ltd. All rights reserved.

Journal of Geodynamics

Natural resource condition assessment: Channel Islands National Park

The Natural Resource Condition Assessment (NRCA) Program aims to provide documentation about the current conditions of important Park natural resources through a spatially explicit, multidisciplinary synthesis of existing scientific data and knowledge. Findings from the NRCA will help Channel Islands National Park (CINP) managers to develop near-term management priorities. The central objectives of this assessment were to evaluate the current conditions and trends of key Park resources, identify important knowledge gaps, and determine key stressors and threats. Because a recent report provided an assessment of marine resources (Engle 2006), which constitute 50% of the park, this NRCA focused on terrestrial resources. For this NRCA, CINP resource managers identified that the greatest Park need was to know whether or not vegetation (primarily woody vegetation) is recovering on the Channel Islands and how this may be influencing recovery of the native terrestrial vertebrates. To address this, the CINP NRCA team evaluated the condition of island vegetation and vertebrate communities using both: 1) already available, summarized information from reports, theses, and journal articles and 2) new analyses of recently available, existing data collected on the islands. Descriptions of specific methods used to analyze available data by the NRCA team are provided within the assessment sections in Chapter 4 of this report. The Team used the Vital Signs indicators as identified by Fancy et al. (2009) for monitoring the condition of Natural Resources in U.S. National Parks. Using this framework, the focal study resource components were the five CINP islands, because these are the bio-geophysical resources of greatest interest to the CINP. Each of the islands has its own unique history of human impacts, ecological conditions, restoration efforts, and natural resource monitoring. The Team evaluated island scrub recovery by assessing changes over time in two metrics: 1) vegetation community structure (species composition, cover, and spatial extent), and 2) vertebrate abundance, relative to the reference condition of each island. The reference conditions for this NRCA were the particular island’s ecological state prior to European settlement or ranching, which represent the ecological state to which the NPS desires the park return. For each metric, the Team determined the Reference Conditions/Values, Current Condition, Trend, and Recovery Potential (based on a review of the literature and from the NRCA team’s analysis of newly available data), Threats and Stressors, Data Gaps, Overall Condition. A qualitative statement of overall current condition was created for each metric, after all data and literature relevant to the measures (i.e., vegetation and vertebrates) of each resource component (i.e., island) were reviewed and evaluated. All information gathered for this NRCA is available and managed by CINP, so that the compilation of resources would be available not only for the assessment team, but will be provided for future Park reference. Chapter 1 of this NRCA provides a descriptive background of NRCAs. Chapter 2 provides Park-wide natural resource background information, while Chapter 3 describes the overall methodological approach used in this assessment. Chapter 4 provides a detailed assessment of each island’s condition, which is then summarized Park-wide and for each island in Chapter 5. Based on both vegetation recovery and vertebrate population abundances, our assessment found Santa Cruz Island (SCI), Santa Rosa Island (SRI), San Miguel Island (SMI), and Anacapa Island (AI) to be, overall, in moderate condition and slowly recovering following the removal of non-native herbivores. The condition of Santa Barbara Island (SBI) is poor and remains in a state that is not likely to improve without assisted recovery by the Park. All islands have exhibited some level of native vegetation recovery, with most of the observed natural recovery occurring on SCI, SRI, SMI, and AI; very little recovery has occurred at present on SBI. Vegetation recovery is uneven across each island and largely limited to mesic areas, such as north-facing slopes. Areas that will likely require some assisted recovery efforts are those that are: at risk of erosion, isolated, lacking native seed banks, and plagued with high competition of non-native grasses and other invasives. Some plant communities, such as exposed upland vegetation on SRI, also have lower natural recovery potential and will likely require assisted recovery. The overall condition of the Park’s terrestrial island vertebrate species generally appears to be good and improving, with high potential for further recovery to reach pre-ranching era reference conditions. This is particularly true for vertebrates on SCI, SRI, and SMI. Although there are limited data for AI, the data that are available suggest that vertebrates on that island are in good condition and improving as well. However, the condition of SBI vertebrates is only moderate, with natural recovery potential being low, with many landbird species’ populations currently declining, extinct, or not improving. An overarching theme across the Park’s terrestrial island vertebrate populations is that many of them are in good condition because the particular animals are habitat generalists, so their population dynamics are not strongly tied to habitat conditions. However, those landbirds associated with riparian vegetation are not showing signs of recovery across the Park islands, and riparian areas appear to be one of the habitats with limited natural recovery potential on the islands. Overwhelmingly, the most critical data gap that emerged from our analysis was the lack of current vegetation maps for most of the islands (AI, SCI, SRI, SMI). The lack of this information made it difficult to assess spatial changes in island vegetation over time, and to relate those changes to observed changes in vertebrate species abundances. Therefore, we were limited in our spatial assessment of vegetation recovery for AI, SCI, SRI, and SMI, and were only able to relate the current vegetation with island landbird analysis for SBI. Another common theme across islands was the lack of information on best techniques for shrub reintroduction and recovery, and the success of invasive species removal and re-vegetation efforts. Abundance, distribution, and habitat affinities of reptiles and amphibians are poorly known for AI, SCI, SRI, and SMI, making it difficult to assess their population condition and trends, and their potential responses to shrub recovery. Mammal and landbird populations are reasonably well-sampled across most of the islands, but data are lacking for most vertebrates from AI, especially West and Middle AI where no vertebrate sampling occurs. The population status, cause of declines, and recovery needs are lacking for several landbird species across multiple islands, including loggerhead shrikes ( Lanius ludovicianus ssp.) on SCI and SRI, and Santa Cruz Island rufous crowned sparrows ( Aimophila ruficeps obscura ) on AI. Additionally, island song sparrows ( Melospiza melodia graminea ) remain extirpated on SBI and data are needed for assessing their recovery potential on the island. The potentially major threats facing the islands’ flora and fauna today are: areas of continued erosion, potential new non-native introductions, and climate change. The Park islands have already experienced warming over the last several decades due to human-induced climate change, and are projected to experience continued warming. These increasing stresses imposed by climate change may slow and inhibit natural recovery potentials for island fauna, and especially flora.

California

Early postcaldera rhyolite and structural resurgence at Long Valley Caldera, California

After the 767-ka caldera-forming eruption of 650 km 3 of rhyolite magma as the Bishop Tuff, 90–100 km 3 of similar rhyolite erupted in the west-central part of Long Valley caldera in as many as 40 batches spread over the 110,000-year interval from ~ 750 ka to ~ 640 ka. Centrally, this Early Rhyolite (ER) is as thick as 622 m, but it spread radially to cover much of the caldera floor, where half its area is now concealed by post-ER sediments and lavas. At least 75% of the ER is aphyric rhyolite tuff. Drillholes encountered ~ 22 (altered) ER lava flows intercalated in the pyroclastic pile, and another 11 units of (largely fresh) ER lava are exposed on the caldera's resurgent dome and at Lookout Mountain. Exposed units have been distinguished, mapped, studied petrographically and chemically, and radioisotopically dated; each is described in detail. Their phenocryst contents range from 0 to 2.5 wt%. All the phyric units have plagioclase, orthopyroxene, and ilmenite; most have biotite and rare tiny magnetite, and a few contain rare zircon. The compositional range of fresh obsidians is narrow—74.3–75.0% SiO 2 , 1.21–1.37% FeO*, and 5.12–5.26% K 2 O, but wider variations in Ti, Ba, Sr, and Zr permit distinction of individual units and eruptive groups. The limited chemical and petrographic variability shown by so many ER batches released episodically for ~ 110,000 years suggests a thermally buffered and well-stirred reservoir. The ER central area, where ER eruptions had taken place, was uplifted ~ 400 m to form a structural dome ~ 10 km in diameter. Most of the inflation is attributable to 10 sills of ER that intrude the Bishop Tuff beneath the uplift, but other processes potentially contributing to resurgence are also considered. As shown by erratics of Mesozoic rocks ice-rafted from the Sierra Nevada and dropped on ER lavas, much of the ER had erupted early enough and at low enough elevation to be inundated by the intracaldera lake and was only later lifted by the resurgence that also raised clusters of the erratics hundreds of meters higher than any shoreline. Most of the uplift was over by ~ 570 ka, but dome-crossing faults that exhibit normal throw of 10–30 m cut lavas as young as 175–125 ka. For most elements, chemical ranges of the ER lie within those of the zoned Bishop Tuff, which had erupted earlier from the same place. Only Ba, Zr, Hf, and Eu/Eu* extend to ranges outside those of the Bishop Tuff, nominally to less evolved compositions. Initial 87 Sr/ 86 Sr values of ER are likewise within the range of the Bishop Tuff, but ER ratios of 143 Nd/ 144 Nd and 206 Pb/ 204 Pb extend beyond those of the Bishop Tuff to values slightly more influenced by upper-crustal contributions. FeTi-oxide geothermometry yields 752°–844 °C for ER, compared to 700°–820 °C for the Bishop Tuff. ER fO 2 values are 0.5–1.0 log units more reduced than those of the T–fO 2 array of the Bishop Tuff. The postcaldera reduction may reflect reaction with graphite from the black lithics of Paleozoic graphitic metapelite so abundant in the Bishop Tuff. Much of the pumice emplaced during the later half of the Bishop Tuff eruption has 10–25 wt% phenocrysts, dominantly quartz and sanidine, but the 100 km 3 of ER has only 0–2.5 wt% and completely lacks quartz and sanidine. Postcaldera processes, including mixing, volatile ascent, and crystal resorption, as well as potential contaminants and magmatic inputs, are all considered.

California

Mapping cropland extent of Southeast and Northeast Asia using multi-year time-series Landsat 30-m data using Random Forest classifier on Google Earth Engine

Cropland extent maps are useful components for assessing food security. Ideally, such products are a useful addition to countrywide agricultural statistics since they are not politically biased and can be used to calculate cropland area for any spatial unit from an individual farm to various administrative unites (e.g., state, county, district) within and across nations, which in turn can be used to estimate agricultural productivity as well as degree of disturbance on food security from natural disasters and political conflict. However, existing cropland extent maps over large areas (e.g., Country, region, continent, world) are derived from coarse resolution imagery (250 m to 1 km pixels) and have many limitations such as missing fragmented and\or small farms with mixed signatures from different crop types and\or farming practices that can be, confused with other land cover. As a result, the coarse resolution maps have limited useflness in areas where fields are small (<1 ha), such as in Southeast Asia. Furthermore, coarse resolution cropland maps have known uncertainties in both geo-precision of cropland location as well as accuracies of the product. To overcome these limitations, this research was conducted using multi-date, multi-year 30-m Landsat time-series data for 3 years chosen from 2013 to 2016 for all Southeast and Northeast Asian Countries (SNACs), which included 7 refined agro-ecological zones (RAEZ) and 12 countries (Indonesia, Thailand, Myanmar, Vietnam, Malaysia, Philippines, Cambodia, Japan, North Korea, Laos, South Korea, and Brunei). The 30-m (1 pixel = 0.09 ha) data from Landsat 8 Operational Land Imager (OLI) and Landsat 7 Enhanced Thematic Mapper (ETM+) were used in the study. Ten Landsat bands were used in the analysis (blue, green, red, NIR, SWIR1, SWIR2, Thermal, NDVI, NDWI, LSWI) along with additional layers of standard deviation of these 10 bands across 1 year, and global digital elevation model (GDEM)-derived slope and elevation bands. To reduce the impact of clouds, the Landsat imagery was time-composited over four time-periods (Period 1: January- April, Period 2: May-August, and Period 3: September-December) over 3-years. Period 4 was the standard deviation of all 10 bands taken over all images acquired during the 2015 calendar year. These four period composites, totaling 42 band data-cube, were generated for each of the 7 RAEZs. The reference training data (N = 7849) generated for the 7 RAEZ using sub-meter to 5-m very high spatial resolution imagery (VHRI) helped generate the knowledge-base to separate croplands from non-croplands. This knowledge-base was used to code and run a pixel-based random forest (RF) supervised machine learning algorithm on the Google Earth Engine (GEE) cloud computing environment to separate croplands from non-croplands. The resulting cropland extent products were evaluated using an independent reference validation dataset (N = 1750) in each of the 7 RAEZs as well as for the entire SNAC area. For the entire SNAC area, the overall accuracy was 88.1% with a producer’s accuracy of 81.6% (errors of omissions = 18.4%) and user’s accuracy of 76.7% (errors of commissions = 23.3%). For each of the 7 RAEZs overall accuracies varied from 83.2 to 96.4%. Cropland areas calculated for the 12 countries were compared with country areas reported by the United Nations Food and Agriculture Organization and other national cropland statistics resulting in an R 2 value of 0.93. The cropland areas of provinces were compared with the province statistics that showed an R 2 = 0.95 for South Korea and R 2 = 0.94 for Thailand. The cropland products are made available on an interactive viewer at www.croplands.org and for download at National Aeronautics and Space Administration’s (NASA) Land Processes Distributed Active Archive Center (LP DAAC): https://lpdaac.usgs.gov/node/1281 .

International Journal of Applied Earth Observation

Oceans on Mars: An assessment of the observational evidence and possible fate

If the large Late Hesperian outflow channels were eroded by extensive floods, as appears likely, then large bodies of water must have once occupied the northern plains during that period. Previous estimates of the sizes of bodies of water in the northern lowlands range up to 3 ?? 108 km3. Several contacts have been previously mapped around the edges of the northern plains and interpreted to be shorelines remaining from these former standing bodies of water. We examine the elevations and geologic relations along these contacts in detail and find little support for their interpretation as shorelines. Some contacts are clearly of volcanic origin, and all have significant variations in elevation. Better support for the former presence of water over large parts of the northern plains is provided by the Vastitas Borealis Formation (VBF). Most of the post-Noachian fill within the northern basin is ridged plains of Lower Hesperian age, interpreted to be volcanic in origin. Overlying the ridged plains is the VBF, a thin veneer of material of Upper Hesperian age. The VBF may have been deposited from large floods. Support for this interpretation is the similarity in age between the outflow channels and the VBF, the presence of the VBF at the lower ends of the outflow channels, and identification of numerous features in the outcrop areas of the VBF that are suggestive of basal melting of an ice sheet. To cover all the area over which the VBF is exposed would require ???2.3 ?? 107 km3 of water. Spread over the entire surface of Mars, this volume is equal to a global layer (global equivalent layer, or GEL) ???156 m deep. We find no support for the larger estimates of ocean volumes that range up to 3 ?? 108 km3 and which imply comparable amounts of water per unit area as are currently present on the surface of the Earth. Under present climatic conditions on Mars an ocean would freeze in a geologically short time period (???104 years), then would sublimate away at rates strongly dependent on the presence or absence of debris on the ice surface. The present VBF is interpreted as a sublimation residue from the ponded outflow channel effluents. The fate of a volume of water thought to have been emplaced by the outflow channels (???2.3 ?? 107 km3) is largely accounted for by the presence of other existing reservoirs on the planet. An approximately 20-30 m GEL of water is estimated to be in the present polar caps, and a 50 m GEL may have escaped to space since the Hesperian, leaving ???80 m GEL unaccounted for. This amount may be partly trapped in other volatile-rich deposits on the surface, and a significant amount could have reentered the groundwater system by south polar basal melting and been progressively cold-trapped at the base of a growing cryosphere. On the basis of our assessment of the Hesperian-aged deposits, we predict that testing of the Clifford and Parker [2001] hypothesis that a Noachian-aged ocean covered up to one third of the surface of Mars will be made very difficult by the enhanced degradation rates in the Noachian and subsequent geological events in the northern lowlands.

Journal of Geophysical Research E: Planets

The geology and distribution of impact craters on Venus: What are they telling us?

Magellan has revealed an ensemble of impact craters on Venus that is unique in many important ways. We have compiled a data base describing the 842 craters on 89% of Venus' surface mapped through orbit 2578. (The craters range in diameter from 1.5 to 280 km.) We have studied the distribution, size‐density, morphology, geology, and associated surface properties of these craters both in the aggregate and, for some craters, in greater detail. We find that (1) the spatial distribution of craters is highly uniform; (2) the size‐density distribution of large craters (diameters ≥35 km) is similar to the young crater populations on other terrestrial planets but at a much lower density that indicates an average age of only about 0.5 Ga (based on the estimated population of Venus‐crossing asteroids); (3) unlike the case on other planets, the density of small craters (diameters ≤35 km) declines rapidly with decreasing diameters because of atmospheric filtering; (4) the spectrum of crater modification differs greatly from that on other planets: 62% of all craters are pristine, only 4% are embayed by lavas, and the remainder are affected by tectonism, but none are severely and progressively depleted (as extrapolated from the size‐density distribution of larger craters); (5) large craters have a progression of morphologies generally similar to those on other planets, but small craters are typically irregular or multiple rather than bowl shaped; (6) diffuse radar‐bright or ‐dark features surround some craters, and 367 similar diffuse “splotches” with no central crater are observed; and (7) other crater features unique to Venus include radar‐bright or ‐dark parabolic arcs opening westward and extensive outflows originating in crater ejecta. The first three of these observations are entirely unexpected. We interpret them as indicating that the planet's cratering record was erased by a global resurfacing event or events, the latest ending about 0.5 Ga, after which volcanic activity declined (but did not cease entirely). Since the last resurfacing event, a maximum of 10% of the planet has been resurfaced and only about 4% of the craters have been obliterated. Convective thermal evolution models support this interpretation (Arkani‐Hamed and Toksoz, 1984). Observations 3–7 confirm quantitatively the expectation that the dense atmosphere of Venus has strongly affected the production of craters. Large impactors have been relatively unaffected, intermediate‐sized ones have been fragmented and have produced overlapping or multiple craters, a narrow size range has produced shock‐induced “splotches” but no craters, and the smallest bodies have had no observable effect on the surface. The number of craters eliminated by the “atmospheric filter” is enormous, about 98% of the craters between 2 and 35 km in diameter that Magellan might have observed on a hypothetical airless Venus. Unique crater‐related features such as parabolas and outflow deposits demonstrate the roles of Venus' high atmospheric density and temperature in modifying the crater formation process. Finally, heavily fractured craters and lava‐embayed craters are found to have higher than average densities along the major fracture belts and rifted uplands connecting Aphrodite Terra and Atla, Beta, Themis, and Phoebe regiones. These craters thus provide physical evidence for recent volcanic and tectonic activity at a low level.

Journal of Geophysical Research E: Planets