Liquid water on Mars

The factors which affect fusion and evaporation of ice under a variety of Martian surface conditions are examined. It is found that a frost or ice deposit will pass through a transient liquid phase in temperate latitudes during summer, if the ice is partly or wholly dust covered. The barrier to free gaseous diffusion which the surface material presents is, under favorable (and definable) conditions, more than adequate to force the water to remain in the liquid state until its evaporation is complete. Furthermore, for a realistic range of regolith particle sizes and porosities, and depths of burial of the ice, the lifetime of the ice can be considerably longer than the duration of a single diurnal warming cycle. Current knowledge of the seosonal and diurnal behavior of the atmospheric vapor is summarized and discussed as it relates to the availability of surface ice at temperate latitudes.     

Error analysis of integrated water vapor measured by СIMEL photometer

Water vapor plays a key role in weather and climate forming, which leads to the need for continuous monitoring of its content in different parts of the Earth. Intercomparison and validation of different methods for integrated water vapor (IWV) measurements are essential for determining the real accuracies of these methods. CIMEL photometers measure IWV at hundreds of ground-based stations of the AERONET network. We analyze simultaneous IWV measurements performed by a CIMEL photometer, an RPG-HATPRO MW radiometer, and a FTIR Bruker 125-HR spectrometer at the Peterhof station of St. Petersburg State University. We show that the CIMEL photometer calibrated by the manufacturer significantly underestimates the IWV obtained by other devices. We may conclude from this intercomparison that it is necessary to perform an additional calibration of the CIMEL photometer, as well as a possible correction of the interpretation technique for CIMEL measurements at the Peterhof site.     

Heterogeneity influences on stream water–groundwater interactions in a gravel-dominated floodplain

ABSTRACT

Floodplains are composed of complex depositional patterns of ancient and recent stream sediments, and research is needed to address the manner in which coarse floodplain materials affect stream–groundwater exchange patterns. Efforts to understand the heterogeneity of aquifers have utilized numerous techniques typically focused on point-scale measurements; however, in highly heterogeneous settings, the ability to model heterogeneity is dependent on the data density and spatial distribution. The objective of this research was to investigate the correlation between broad-scale methodologies for detecting heterogeneity and the observed spatial variability in stream/groundwater interactions of gravel-dominated alluvial floodplains. More specifically, this study examined the correlation between electrical resistivity (ER) and alluvial groundwater patterns during a flood event at a site on Barren Fork Creek, in the Ozark ecoregion of Oklahoma, USA, where chert gravels were common both as streambed and as floodplain material. Water table elevations from groundwater monitoring wells for a flood event on 1–5 May 2009 were compared to ER maps at various elevations. Areas with high ER matched areas with lower water table slope at the same elevation. This research demonstrated that ER approaches were capable of indicating heterogeneity in surface water–groundwater interactions, and that these heterogeneities were present even in an aquifer matrix characterized as highly conductive. Portions of gravel-dominated floodplain vadose zones characterized by high hydraulic conductivity features can result in heterogeneous flow patterns when the vadose zone of alluvial floodplains activates during storm events.

EDITOR D. Koutsoyiannis; ASSOCIATE EDITOR X. Chen     

Fast and direct nonparametric procedures in the L-moment homogeneity test

Regional frequency analysis is an important tool to properly estimate hydrological characteristics at ungauged or partially gauged sites in order to prevent hydrological disasters. The delineation of homogeneous groups of sites is an important first step in order to transfer information and obtain accurate quantile estimates at the target site. The Hosking–Wallis homogeneity test is usually used to test the homogeneity of the selected sites. Despite its usefulness and good power, it presents some drawbacks including the subjective choice of a parametric distribution for the data and a poorly justified rejection threshold. The present paper addresses these drawbacks by integrating nonparametric procedures in the L-moment homogeneity test. To assess the rejection threshold, three resampling methods (permutation, bootstrap and Pólya resampling) are considered. Results indicate that permutation and bootstrap methods perform better than the parametric Hosking–Wallis test in terms of power as well as in time and procedure simplicity. A real-world case study shows that the nonparametric tests agree with the HW test concerning the homogeneity of the volume and the bivariate case while they disagree for the peak case, but that the assumptions of the HW test are not well respected.     

The instability of an electron beam with a finite radius in a boundless plasma

Within the framework of a linear theory, the instability of an electron beam with a finite radius in a cold magnetised boundless plasma is considered. It is shown that a finite beam dimension influences the generation of quasi-potential waves in two aspects: the perpendicular wave number is quantised so that the frequencies of the waves are subjected to strong selection; a new kind of instability appears due to wave energy losses by emission into surrounding space. Growth rate dependence of wave numbers and frequencies is investigated for typical parameters of experiments with electron beams in space.     

Evidence and Implications of Pressure Fluctuations in the ISM

A recent survey of the fine-structure excitation of neutral carbonreveals that the interstellar medium in the Galactic plane exhibits athermal pressure, nT/k, that ranges from about 10³to 10⁴cm^-3K from one location to the next, with occasional excursions in excess of about 10⁵ cm^-3K. The large excitations for small amounts of gas indicate that some regions are either subjected to shocks or must be pressurized within time scales much shorter than the time needed to reach thermal equilibrium. These rapidfluctuations probably arise from the cascade of macroscopic mechanical energyto small scales through a turbulent cascade. One consequence of thiseffect is that changes in gas temperature can arise from near adiabaticcompressions and expansions, and this may explain why investigations of21-cm emission and absorption reveal the presence of hydrogen attemperatures well below the expected values derived from the balance ofvarious known heating and cooling processes.     

An iterative scheme for use in MT 2D forward modelling — to speed up production of the Jacobian and to improve accuracy of the inverse matrix

An important stage in two-dimensional magnetotelluric modelling is the calculation of the Earth's response functions for an assumed conductivity model and the calculation of the associated Jacobian relating those response functions to the model parameters. The efficiency of the calculation of the Jacobian will affect the efficiency of the inversion modelling. Rodi (1976) produced all the Jacobian elements by inverting a single matrix and using an approximate first-order algorithm. Since only one inverse matrix required calculation the procedure speeded up the inversion. An iterative scheme to improve the approximation to the Jacobian information is presented in this paper. While this scheme takes a little longer than Rodi's algorithm, it enables a more accurate determination of the Jacobian information. It is found that the Jacobian elements can be produced in 10% of the time required to calculate an inverse matrix or to calculate a 2D starting model. A modification of the algorithm can further be used to improve the accuracy of the original inverse matrix calculated in a 2D finite difference program and hence the solution this program produces. The convergence of the iteration scheme is found to be related both to the originally calculated inverse matrix and to the change in the newly formed matrix arising from perturbation of the model parameter. A ridge regression inverse algorithm is used in conjunction with the iterative scheme for forward modelling described in this paper to produce a 2D conductivity section from field data.