Nonuniform cratering of the terrestrial planets

We estimate the impact flux and cratering rate as a function of latitude on the terrestrial planets using a model distribution of planet crossing asteroids and comets [Bottke, W.F., Morbidelli, A., Jedicke, R., Petit, J.-M., Levison, H.F., Michel, P., Metcalfe, T.S., 2002. Icarus 156, 399-433]. After determining the planetary impact probabilities as a function of the relative encounter velocity and encounter inclination, the impact positions are calculated analytically, assuming the projectiles follow hyperbolic paths during the encounter phase. As the source of projectiles is not isotropic, latitudinal variations of the impact flux are predicted: the calculated ratio between the pole and equator is 1.05 for Mercury, 1.00 for Venus, 0.96 for the Earth, 0.90 for the Moon, and 1.14 for Mars over its long-term obliquity variation history. By taking into account the latitudinal dependence of the impact velocity and impact angle, and by using a crater scaling law that depends on the vertical component of the impact velocity, the latitudinal variations of the cratering rate (the number of craters with a given size formed per unit time and unit area) is in general enhanced. With respect to the equator, the polar cratering rate is about 30% larger on Mars and 10% on Mercury, whereas it is 10% less on the Earth and 20% less on the Moon. The cratering rate is found to be uniform on Venus. The relative global impact fluxes on Mercury, Venus, the Earth and Mars are calculated with respect to the Moon, and we find values of 1.9, 1.8, 1.6, and 2.8, respectively. Our results show that the relative shape of the crater size-frequency distribution does not noticeably depend upon latitude for any of the terrestrial bodies in this study. Nevertheless, by neglecting the expected latitudinal variations of the cratering rate, systematic errors of 20-30% in the age of planetary surfaces could exist between equatorial and polar regions when using the crater chronology method.     

226Ra and222Rn contents of Galapagos Rift hydrothermal waters —the importance of low-temperature interactions with crustal rocks

Hydrothermal waters collected by “Alvin” from the Galapagos Spreading Center are enriched in²²²Rn by factors of 50–200 over bottom waters. The²²⁶Ra in the same samples, however, is enriched by less than a factor of four over bottom waters. Enrichments of²²²Rn result primarily from α-recoil from rock surfaces while²²⁶Ra enrichments are dominantly produced by high-temperature alteration of cooling ridge volcanics. The abundances of both nuclides exhibit positive correlations with temperature. The data extrapolate to bottom water temperatures and compositions, demonstrating the importance of seawater mixing. Different vents, however, have different mixing lines, and vents with high²²²Rn have low²²⁶Ra. We propose these patterns result from variations in the extent of low-temperature crustal interaction with the hydrothermal fluids. Low-temperature crustal waters can maintain high steady state²²²Rn contents due to the α-recoil additions to the fluids. The²²⁶Ra, however, is strongly adsorbed at low-temperatures resulting in low concentrations of this nuclide in low-temperature crustal waters. Thus, physical mixing of a crustal water component with hydrothermal waters or variable crustal path lengths of the hydrothermal fluids can account for the variable mixing lines and²²²Rn/²²⁶Ra values of the hot springs.The²²²Rn/²²⁶Ra value appears to be a sensitive indicator of low-temperature crustal interaction. Values > 100 have experienced extensive crustal interaction and are indicative of diffuse hydrothermal flow. Values between 1 and 10 are indicative of primary hydrothermal fluids which have not experienced significant interaction with the crust. Values of²²²Rn/²²⁶Ra between 10³ and 10⁴ are indicative of interaction of the hydrothermal fluids with sediments. Such values are observed in water samples from the Galapagos hydrothermal mounds.     

Long-term increase in diffuse groundwater recharge following expansion of rainfed cultivation in the Sahel,West Africa

Rapid population growth in sub-Saharan West Africa and related cropland expansion were shown in some places to have increased focused recharge through ponds, raising the water table. To estimate changes in diffuse recharge, the water content and matric potential were monitored during 2009 and 2010, and modeling was performed using the Hydrus-1D code for two field sites in southwest Niger: (1) fallow land and (2) rainfed millet cropland. Monitoring results of the upper 10 m showed increased water content and matric potential to greater depth under rainfed cropland (>2.5 m) than under fallow land (≤1.0 m). Model simulations indicate that conversion from fallow land to rainfed cropland (1) increases vadose-zone water storage and (2) should increase drainage flux (～25 mm year^?1) at 10-m depth after a 30–60 year lag. Therefore, observed regional increases in groundwater storage may increasingly result from diffuse recharge, which could compensate, at least in part, groundwater withdrawal due to observed expansion in irrigated surfaces; and hence, contribute to mitigate food crises in the Sahel.     

Coastal conglomerates around the Hadjer el Khamis inselbergs (western Chad,central Africa): new evidence for Lake Mega‐Chad episodes

This paper reports on a study dealing with the rhyolitic inselbergs of Hadjer el Khamis that formed palaeoislands during Lake Mega‐Chad events. Field observations have shown that: (1) conglomeratic patches of immature to mature clasts are preserved at the feet of the Hadjer el Khamis inselbergs; (2) in cross‐section, their pro?le reveals a well de?ned cliff–platform junction at a constant elevation (325 m). The monolithological clasts show all degrees of roundness, from angular cobbles to well rounded pebbles. This wide range of maturity suggests a coastal origin for these cobbles. The system was permanently fed with angular clasts, which were progressively worn by waves. Cobbles that were wave‐worked for the longest time are the best rounded. The cliff–platform junction is the result of erosion by waves, which attacked and undercut the inselberg cliffs during Lake Mega‐Chad events. Asymmetrical erosion pro?les moreover suggest a wind regime dominated by SW to NE oriented winds. These interpretations have two implications. The elevation of the cliff–platform junction is an indication of the highest water level of Lake Mega‐Chad at 320–325 m, which is in agreement with other observations elsewhere in the basin. The SW to NE oriented winds show that monsoon‐related winds were prevalent during Lake Mega‐Chad events, suggesting the Inter‐Tropical Convergence Zone was located higher in latitude than today. Copyright © 2003 John Wiley & Sons, Ltd.     

Integrity monitoring of fixed ambiguity Precise Point Positioning (PPP) solutions

Traditional positioning methods, such as conventional Real Time Kinematic (cRTK) rely upon local reference networks to enable users to achieve high-accuracy positioning. The need for such relatively dense networks has significant cost implications. Precise Point Positioning (PPP) on the other hand is a positioning method capable of centimeter-level positioning without the need for such local networks, hence providing significant cost benefits especially in remote areas. This paper presents the state-of-the-art PPP method using both GPS and GLONASS measurements to estimate the float position solution before attempting to resolve GPS integer ambiguities. Integrity monitoring is carried out using the Imperial College Carrier-phase Receiver Autonomous Integrity Monitoring method. A new method to detect and exclude GPS base-satellite failures is developed. A base-satellite is a satellite whose measurements are differenced from other satellite’s measurements when using between-satellite-differenced measurements to estimate position. The failure detection and exclusion methods are tested using static GNSS data recorded by International GNSS Service stations both in static and dynamic processing modes. The results show that failure detection can be achieved in all cases tested and failure exclusion can be achieved for static cases. In the kinematic processing cases, failure exclusion is more difficult because the higher noise in the measurement residuals increases the difficulty to distinguish between failures associated with the base-satellite and other satellites.     

Effects of sample size on the accuracy of geomorphological models

Commonly, the most costly part of geomorphological distribution modelling studies is gathering the data. Thus, guidance for researchers concerning the quantity of field data needed would be extremely practical. This paper scrutinises the relationship between the sample size (the number of observations varied from 20 to 600) and the predictive ability of the generalized linear model (GLM), generalized additive model (GAM), generalized boosting method (GBM) and artificial neural network (ANN) in two data settings, i.e., independent and split-sample approaches. The study was performed using empirical data of periglacial processes from an area of 600 km² in northernmost Finland at grid resolutions of 1 ha (100 × 100 m) and 25 ha (500 × 500 m). A rather sharp increase in the predictive ability of the models was observed when the number of observations increased from 20 to 100, and the level of robust predictions was reached with 200 observations. The result indicates that no more than a few hundred observations are needed in geomorphological distribution modelling at a medium scale resolution (ca. 0.01–1 km²).     

A rapid method for converting medical Computed Tomography scanner topogram attenuation scale to Hounsfield Unit scale and to obtain relative density values

Medical Computed Tomography scanners permit to rapidly obtain qualitative and quantitative information on objects in a non-destructive manner in both longitudinal and transversal directions. CT-scan scales used to express attenuation on the images are different for longitudinal (topograms) and transversal (tomograms) views, restraining the complementary use of the information extracted from both views. Moreover, topograms are of full use in core analysis as they reveal fine-scale information on stratigraphical, sedimentological and physical properties of the material. This paper presents a method to convert CT-scan topogram intensity scale to Hounsfield Unit (HU; scale used on tomograms) and to later extract a relative density from it. The method is based on relationships between topogram and tomogram values, absolute density and effective atomic numbers obtained on a series of minerals. Results show excellent correlations between converted topogram values and HU values obtained on the tomograms, as well as between relative densities derived from the converted topogram values, and the absolute density of the minerals.