Constraints on the flux of meteoritic and cometary water on the Moon from volatile element (N–Ar) analyses of single lunar soil grains,Luna 24 core

We report new nitrogen and argon isotope and abundance results for single breccia clasts and agglutinates from four different sections of the Luna 24 drill core in order to re-evaluate the provenance of N trapped in lunar regolith, and to place limits on the flux of planetary material to the Moon’s surface. Single Luna 24 grains with ⁴⁰Ar/³⁶Ar ratios <1 show δ¹⁵N values between ?54.5‰ and +123.3‰ relative to terrestrial atmosphere. Thus, low-antiquity lunar soils record both positive and negative δ¹⁵N signatures, and the secular increase of the δ¹⁵N value previously postulated by Kerridge (Kerridge, J.F. [1975]. Science 188(4184), 162–164. doi:10.1126/science.188.4184.162) is no longer apparent when the Luna and Apollo data are combined. Instead, the N isotope signatures, corrected for cosmogenic ¹⁵N, are consistent with binary mixing between isotopically light solar wind (SW) N and a planetary N component with a δ¹⁵N value of +100‰ to +160‰. The lower δ¹⁵N values of Luna 24 grains compared to Apollo samples reflect a higher relative proportion of solar N, resulting from the higher SW fluence in the region of Mare Crisium compared to the central near side of the Moon. Carbonaceous chondrite-like micro-impactors match well the required isotope characteristics of the non-solar N component trapped in low-antiquity lunar regolith. In contrast, a possible cometary contribution to the non-solar N flux is constrained to be ?3–13%. Based on the mixing ratio of SW to planetary N obtained for recently exposed lunar soils, we estimate the flux of micro-impactors to be (2.2–5.7) × 10³ tons yr^?1 at the surface of the Moon. Our estimate for Luna 24 agrees well with that for young Apollo regolith, indicating that the supply of planetary material does not depend on lunar location. Thus, the continuous influx of water-bearing cosmic dust may have represented an important source of water for the lunar surface over the past ～1 Ga, provided that water removal rates (i.e., by meteorite impacts, photodissociation, and sputtering) do not exceed accumulation rates.     

On AGV's navigation in 3-dimensional space

This paper deals with position estimation and path control for Autonomous Guided Vehicles (AGV). To enable a vehicle or a mobile robot in following a continuous virtual path without human control, these techniques play an important role. The relationship between the vehicle's motion in 3-dimensional space and the shape of a curved surface is described. In particular, the introduction of a digital terrain model in dead reckoning is considered. Moreover, a possible nonlinear control is developed based on curvilinear path coordinates, and the proof for global stability is given. To achieve general validity, these topics are treated here independently of the cart's special mechanization (the configuration of steered wheels and driven wheels). Simulation studies are presented to illustrate the investigations.     

Westward moving dynamic substorm features observed with the IMAGE magnetometer network and other ground-based instruments

We present the ground signatures of dynamic substorm features with particular emphasis on the event interpretation capabilities provided by the IMAGE magnetometer network. This array covers the high latitudes from the sub-auroral to the cusp/cleft region. An isolated substorm on 11 Oct. 1993 during the late evening hours exhibited many of well-known features such as the Harang discontinuity, westward travelling surge and poleward leap, but also discrete auroral forms, known as auroral streamers, appeared propagating westward along the centre of the electrojet. Besides the magnetic field measurements, there were auroral observations and plasma flow and conductivity measurements obtained by EISCAT. The data of all three sets of instruments are consistent with the notion of upward field-aligned currents associated with the moving auroral patches. A detailed analysis of the electro-dynamic parameters in the ionosphere, however, reveals that they do not agree with the expectations resulting from commonly used simplifying approximations. For example, the westward moving auroral streamers which are associated with field-aligned current filaments, are not collocated with the centres of equivalent current vortices. Furthermore, there is a clear discrepancy between the measured plasma flow direction and the obtained equivalent current direction. All this suggests that steep conductivity gradients are associated with the transient auroral forms. Also self-induction effects in the ionosphere may play a role for the orientation of the plasma flows. This study stresses the importance of multi-instrument observation for a reliable interpretation of dynamic auroral processes.     

Appraisal of a new 1D weighted joint inversion of ground based and helicopter‐borne electromagnetic data

In order to couple spatial data from frequency‐domain helicopter‐borne electromagnetics with electromagnetic measurements from ground geophysics (transient electromagnetics and radiomagnetotellurics), a common 1D weighted joint inversion algorithm for helicopter‐borne electromagnetics, transient electromagnetics and radiomagnetotellurics data has been developed. The depth of investigation of helicopter‐borne electromagnetics data is rather limited compared to time‐domain electromagnetics sounding methods on the ground. In order to improve the accuracy of model parameters of shallow depth as well as of greater depth, the helicopter‐borne electromagnetics, transient electromagnetics, and radiomagnetotellurics measurements can be combined by using a joint inversion methodology. The 1D joint inversion algorithm is tested for synthetic data of helicopter‐borne electromagnetics, transient electromagnetics and radiomagnetotellurics. The proposed concept of the joint inversion takes advantage of each method, thus providing the capability to resolve near surface (radiomagnetotellurics) and deeper electrical conductivity structures (transient electromagnetics) in combination with valuable spatial information (helicopter‐borne electromagnetics). Furthermore, the joint inversion has been applied on the field data (helicopter‐borne electromagnetics and transient electromagnetics) measured in the Cuxhaven area, Germany. In order to avoid the lessening of the resolution capacities of one data type, and thus balancing the use of inherent and ideally complementary information content, a parameter reweighting scheme that is based on the exploration depth ranges of the specific methods is proposed. A comparison of the conventional joint inversion algorithm, proposed by Jupp and Vozoff ( 1975 ), and of the newly developed algorithm is presented. The new algorithm employs the weighting on different model parameters differently. It is inferred from the synthetic and field data examples that the weighted joint inversion is more successful in explaining the subsurface than the classical joint inversion approach. In addition to this, the data fittings in weighted joint inversion are also improved.     

Forward induced seismic hazard assessment: application to a synthetic seismicity catalogue from hydraulic stimulation modelling

The M _w 3.2-induced seismic event in 2006 due to fluid injection at the Basel geothermal site in Switzerland was the starting point for an ongoing discussion in Europe on the potential risk of hydraulic stimulation in general. In particular, further development of mitigation strategies of induced seismic events of economic concern became a hot topic in geosciences and geoengineering. Here, we present a workflow to assess the hazard of induced seismicity in terms of occurrence rate of induced seismic events. The workflow is called Forward Induced Seismic Hazard Assessment (FISHA) as it combines the results of forward hydromechanical-numerical models with methods of time-dependent probabilistic seismic hazard assessment. To exemplify FISHA, we use simulations of four different fluid injection types with various injection parameters, i.e. injection rate, duration and style of injection. The hydromechanical-numerical model applied in this study represents a geothermal reservoir with preexisting fractures where a routine of viscous fluid flow in porous media is implemented from which flow and pressure driven failures of rock matrix and preexisting fractures are simulated, and corresponding seismic moment magnitudes are computed. The resulting synthetic catalogues of induced seismicity, including event location, occurrence time and magnitude, are used to calibrate the magnitude completeness M _c and the parameters a and b of the frequency-magnitude relation. These are used to estimate the time-dependent occurrence rate of induced seismic events for each fluid injection scenario. In contrast to other mitigation strategies that rely on real-time data or already obtained catalogues, we can perform various synthetic experiments with the same initial conditions. Thus, the advantage of FISHA is that it can quantify hazard from numerical experiments and recommend a priori a stimulation type that lowers the occurrence rate of induced seismic events. The FISHA workflow is rather general and not limited to the hydromechanical-numerical model used in this study and can therefore be applied to other fluid injection models.     

Source parameters for the M<Subscript>w</Subscript> = 6.6, 03 February 2002, Çay Earthquake (Turkey) and aftershocks from GPS,Southwestern Turkey

The 03 February 2002 Çay Earthquake (M_w ～6.7) occurred on the fault segment between Eber and Ak?ehir Lakes followed by a large aftershock (M_w ～5.6) near the western end of the fault and two sequential aftershocks. We computed the coseismic surface displacements from static GPS measurements to determine the fault geometry parameters and uniform slip components. The coseismic displacements were obtained through combining the regional pre-earthquake and post-earthquake GPS data. Fault geometry and slips were acquired through the inversion of GPS data modeling the events as elastic dislocations in a half-space and assuming all four events took place on the same fault plane. Results suggest that one-segment fault of ～33 km length and dipping ～43° northward suffices to model the dislocation, assuming uniform slip distribution with 0.51 m dip slip, 0.26 m left-lateral slip extending to a depth down to ～11.5 km which is consistent with seismological evidence. The results also verify the normal faulting in the eastern flank of Isparta Angle which has long been assumed as a thrusting structure. While the available data cannot identify the four individual events on the same day, an attempted distributed slip model differentiates dip slip and left-lateral slips near the hypocenter with maximum values of ～1 and 0.6 m, respectively.     

Sea level variability at Antalya and Mente? tide gauges in Turkey: atmospheric, steric and land motion contributions

Sea level trends and interannual variability at Antalya and Menteş tide gauges are investigated during the 1985–2001 period, quantifying the roles of atmospheric, steric and local land motion contributions. Tide gauge sea level measurements, temperature/salinity climatologies and GPS data are used in the analyses and the results are compared with the output of a barotropic model forced by atmospheric pressure and wind. The overall sea level trends at two tide gauges collocated with GPS are in the range of 5.5 to 7.9 mm/yr during the study period, but showing different behaviour in the sub-periods 1985–1993 and 1993–2001 due to variations in the contributing factors both in space and time. After the removal of the atmospheric forcing and steric contribution from sea level records, the resulting trends vary between 1.9 to 4.5 mm/yr in Antalya and −1.2 to −11.6 mm/yr in Menteş depending on the period considered. Vertical land movement estimated from GPS data seems to explain the high positive residual trend in Antalya during the whole period. On the other hand, the source of the highly negative sea level trend of about −14 mm/yr in Menteş during 1985–1993 could not be resolved with the available datasets. Interannual variability of wind and atmospheric pressure appear to dominate the sea level at both tide gauges during the study period. Atmospheric and steric contributions together account for ∼50% of the total sea level variance at interannual time scales. Mass induced sea level variations which were not considered in this study may help to close the sea level trend budgets as well as to better explain the interannual sea level variance.     

Controls on sediment evacuation from glacially modified and unmodified catchments in the eastern Sierra Nevada,California

The degree of glacial modification in small catchments along the eastern Sierra Nevada, California, controls the timing and pattern of sediment flux to the adjacent fans. There is a close relationship between the depth of fan‐head incision and the pattern and degree of Late Pleistocene catchment erosion by valley glaciers; catchments with significant glacial activity are associated with deeply incised fan heads, whereas fans emerging from glacially unmodified catchments are unincised. We suggest that the depth of fan‐head incision is controlled by the potential for sediment storage during relatively dry ice‐free periods, which in turn is related to the downstream length of the glacially modified valley and creation of accommodation through valley floor slope lowering and glacial valley overdeepening and widening. Significant storage in glacially modified basins during ice‐free periods leads to sediment supply‐limited conditions at the fan head and causes deep incision. In contrast, a lack of sediment trapping allows quasi‐continuous sediment supply to the fan and prevents incision of the fan head. Sediment evacuation rates should thus show large variations in glacially modified basins, with major peaks during glacial and lows during interglacial or ice‐free periods, respectively. In contrast, sediment removal from glacially unmodified catchments in this type of setting should be free of this effect, and will be dominated instead by short‐term variations, modulated for example by changes in vegetation cover or storm frequency. This distinction may help improve our understanding of long‐term sediment yields as a measure of erosional efficiency. Copyright © 2008 John Wiley & Sons, Ltd.