Transient response of a homogeneous half space with due regard for displacement currents

The solution for the half-space model is represented directly in the time domain as computationally stable convolution integrals. The influence of the geoelectric parameters of the earth and transmitter current waveform are then investigated for both infinitesimal and finite-dimensional transmitter loops. Simple empirical formulae are derived to account for the finite duration of the transmitter current turn off time.The whole transient process is divided into three essentially different stages: the propagation stage, the intermediate stage and the diffusion stage. The first is characterized by extremely complicated signal behavior. Apparently, interpretation of the field data using any kind of model fitting inversion algorithm is impossible in this stage. The diffusion stage virtually coincides with that used in the quasi-static case and is, therefore, unsuitable for detecting the dielectric properties of the earth. The intermediate stage is, thus, the only possible time range in which the dielectric properties can be detected using the dynamic characteristics of the signal.The duration of each stage is evaluated depending on the geoelectric parameters of the earth for different transmitter current waveforms.     

The epoch of helium reionization

We study the reionization of He  ii by quasars using a numerical approach that combines 3D radiative transfer calculations with cosmological hydrodynamical simulations. Sources producing the ionizing radiation are selected according to an empirical quasar luminosity function, and are assigned luminosities according to their intrinsic masses. The free parameters associated with this procedure are (1) a universal source lifetime, (2) a minimum mass cut-off, (3) a minimum luminosity cut-off, (4) a solid angle specifying the extent to which radiation is beamed, and (5) a tail-end spectral index for the radiative energy distribution of the sources. We present models in which these parameters are varied, and examine characteristics of the resultant reionization process that distinguish the various cases. In addition, we extract artificial spectra from the simulations and quantify statistical properties of the spectral features in each model.
We find that the most important factor affecting the evolution of He  ii reionization is the cumulative number of ionizing photons that are produced by the sources. Comparisons between He  ii opacities measured observationally and those obtained by our analysis reveal that the available ranges in plausible values for the parameters provide enough leeway to provide a satisfactory match. However, one property common to all our calculations is that the epoch of He  ii reionization must have occurred at a redshift in the range  3≲ z ≲4  . If so, future observational programmes will be able to directly trace the details of the ionization history of helium and to probe the low-density phase of the intergalactic medium during this phase of the evolution of the Universe.     

Sensitivity of Mesoscale Surface Dynamics to Surface Soil and Vegetation Contrasts over the Carolina Sandhills

A region of contrasting soils exists over the Carolinas region of the southeastern United States. Previous research has shown an increase in mesoscale summertime precipitation over this region. Numerical simulations are analyzed to investigate the relationships between mesoscale surface dynamics and the transition from clay to sandy soils over this region. Numerical modeling experiments using four different soil and vegetation patterns suggest that the presence of the clay-to-sand transition zone produces a surface heat flux gradient and enhanced convergence. The soil contrasts appear to dominate over vegetation contrasts in affecting local surface heating and convergence in the model atmosphere.     

2.5-D modeling in electromagnetic methods of geophysics

Understanding, using, or eliminating three-dimensional (3-D) effects in electromagnetic methods of geophysics are critical requirements. Numerous achievements in 3-D modeling sometimes give the impression that they are widely available today in geophysical practice. This is not necessarily true. Existing 3-D modeling packages prove that we know how to perform 3-D modeling. However, the computer resources and costs involved make the practical application of 3-D EM modeling in geophysical applications very limited.A practical compromise, or even alternative, is represented by 2.5-D modeling characterized by the use of a 3-D source in a 2-D medium. This combination allows one to mathematically describe the related boundary value problem as a sequence of independent two-dimensional problems. The typical technique leading to such a split formulation is Fourier analysis. That is why the individual terms of a split solution are often referred to as harmonics.Although each independent problem is two-dimensional, the algorithmic implementation of finite differences or integral equations for the higher harmonics has some specific features not present in the classical 2-D cases. In this paper, a hybrid scheme consisting of a combination of the finite difference technique with the integral equation approach for transient fields is described. Evaluation of algorithm accuracy is presented and a transient logging technique application is considered. The algorithm is fast and easily implemented on a personal computer     

Petrology of lavas from the Puu Oo eruption of Kilauea Volcano: III. The Kupaianaha episode (1986–1992)

 The Puu Oo eruption has been remarkable in the historical record of Kilauea Volcano for its duration (over 13 years), volume (>1 km³) and compositional variation (5.7–10 wt.% MgO). During the summer of 1986, the main vent for lava production moved 3 km down the east rift zone and the eruption style changed from episodic geyser-like fountaining at Puu Oo to virtually continuous, relatively quiescent effusion at the Kupaianaha vent. This paper examines this next chapter in the Puu Oo eruption, episodes 48 and 49, and presents new ICP-MS trace element and Pb-, Sr-, and Nd-isotope data for the entire eruption (1983–1994). Nearly aphyric to weakly olivine-phyric lavas were erupted during episodes 48 and 49. The variation in MgO content of Kupaianaha lavas erupted before 1990 correlates with changes in tilt at the summit of Kilauea, both of which probably were controlled by variations in Kilauea's magma supply rate. These lavas contain euhedral olivines which generally are in equilibrium with whole-rock compositions, although some of the more mafic lavas which erupted during 1990, a period of frequent pauses in the eruption, accumulated 2–4 vol.% olivine. The highest forsterite content of olivines (∼85%) in Kupaianaha lavas indicates that the parental magmas for these lavas had MgO contents of ∼10 wt.%, which equals the highest observed value for lavas during this eruption. The composition of the Puu Oo lavas has progressively changed during the eruption. Since early 1985 (episode 30), when mixing between an evolved rift zone magma and a more mafic summit reservoir-derived magma ended, the normalized (to 10 wt.% MgO) abundances of highly incompatible elements and CaO have systematically decreased with time, whereas ratios of these trace elements and Pb, Sr, and Nd isotopes, and the abundances of Y and Yb, have remained relatively unchanged. These results indicate that the Hawaiian plume source for Puu Oo magmas must be relatively homogeneous on a scale of 10–20 km³ (assuming 5–10% partial melting), and that localized melting within the plume has apparently progressively depleted its incompatible elements and clinopyroxene component as the eruption continued. The rate of variation of highly incompatible elements in Puu Oo lavas is much greater than that observed for Kilauea historical summit lavas (e.g., Ba/Y 0.09 a^–1 vs ∼0.03 a^–1). This rapid change indicates that Puu Oo magmas did not mix thoroughly with magma in the summit reservoir. Thus, except for variable amounts of olivine fractionation, the geochemical variation in these lavas is predominantly controlled by mantle processes. Received: 8 March 1996 / Accepted: 30 April 1996     

Factors influencing the calculation of periodic secondary circulation in a tidal river: numerical modelling of the lower Sacramento River,USA

Secondary circulation is the component of three‐dimensional (3D) flow in river channels perpendicular to the primary flow direction. Secondary circulation calculated from acoustic Doppler current profiler (ADCP) transects is sensitive to the calculation method and is affected by the transect angle relative to the mean flow direction and variations in the flow direction along a transect. To quantify bounds on transect alignment relative to river flow for field data collection and examine tidal time‐scale variability in secondary circulation, the 3D hydrodynamic model UnTRIM was applied to simulate the hydrodynamics in the lower reach of the Sacramento River (CA, USA). Secondary circulation was calculated using the Rozovskii and the zero net discharge methods on repeated transects extracted from the model results in regions of both relatively uniform and complex flows. When the depth‐averaged flow direction along a transect varied by more than about 5 °, occurring when the transect was as little as 10 to 20 ° out of normal to the mean flow direction, the Rozovskii method produced more realistic secondary circulation than the zero net discharge method. Analysis indicated that ADCP transects should be within 20 ° of perpendicular to the mean flow direction when calculating secondary circulation. Secondary circulation strength around two tidally influenced bends generally increased with increasing flow and broke down near slack water. However, the strength of the secondary circulation was not only a function of the flow magnitude, but also depended on the direction of the water flow and the transect location relative to the river curvature, which varied with the tidal flow direction. Copyright © 2015 John Wiley & Sons, Ltd.     

New strategies to detect life on Mars

If we are to find unequivocal evidence for life on Mars, we will need new ways to search for it. Jeff L Bada and the MOD team describe the innovative strategy developed for the ExoMars mission.     

Grazing impacts on gully dynamics indicate approaches for gully erosion control in northeast Australia

Drainage network extension in semi‐arid rangelands has contributed to a large increase in the amount of fine sediment delivered to the coastal lagoon of the Great Barrier Reef, but gully erosion rates and dynamics are poorly understood. This study monitored annual erosion, deposition and vegetation cover in six gullies for 13 years, in granite‐derived soils of the tropical Burdekin River basin. We also monitored a further 11 gullies in three nearby catchments for 4 years to investigate the effects of grazing intensity. Under livestock grazing, the long‐term fine sediment yield from the planform area of gullies was 6.1 t ha^‐1 yr^‐1. This was 7.3 times the catchment sediment yield, indicating that gullies were erosion hotspots within the catchment. It was estimated that gully erosion supplied between 29 and 44% of catchment sediment yield from 4.5% of catchment area, of which 85% was derived from gully wall erosion. Under long‐term livestock exclusion gully sediment yields were 77% lower than those of grazed gullies due to smaller gully extent, and lower erosion rates especially on gully walls. Gully wall erosion will continue to be a major landscape sediment source that is sensitive to grazing pressure, long after gully length and depth have stabilised. Wall erosion was generally lower at higher levels of wall vegetation cover, suggesting that yield could be reduced by increasing cover. Annual variations in gully head erosion and net sediment yield were strongly dependent on annual rainfall and runoff, suggesting that sediment yield would also be reduced if surface runoff could be reduced. Deposition occurred in the downstream valley segments of most gullies. This study concludes that reducing livestock grazing pressure within and around gullies in hillslope drainage lines is a primary method of gully erosion control, which could deliver substantial reductions in sediment yield. Copyright © 2018 John Wiley & Sons, Ltd.     

Using stable isotopes to estimate travel times in a data‐sparse Arctic catchment: Challenges and possible solutions

Use of isotopes to quantify the temporal dynamics of the transformation of precipitation into run‐off has revealed fundamental new insights into catchment flow paths and mixing processes that influence biogeochemical transport. However, catchments underlain by permafrost have received little attention in isotope‐based studies, despite their global importance in terms of rapid environmental change. These high‐latitude regions offer limited access for data collection during critical periods (e.g., early phases of snowmelt). Additionally, spatio‐temporal variable freeze–thaw cycles, together with the development of an active layer, have a time variant influence on catchment hydrology. All of these characteristics make the application of traditional transit time estimation approaches challenging. We describe an isotope‐based study undertaken to provide a preliminary assessment of travel times at Siksik Creek in the western Canadian Arctic. We adopted a model–data fusion approach to estimate the volumes and isotopic characteristics of snowpack and meltwater. Using samples collected in the spring/summer, we characterize the isotopic composition of summer rainfall, melt from snow, soil water, and stream water. In addition, soil moisture dynamics and the temporal evolution of the active layer profile were monitored. First approximations of transit times were estimated for soil and streamwater compositions using lumped convolution integral models and temporally variable inputs including snowmelt, ice thaw, and summer rainfall. Comparing transit time estimates using a variety of inputs revealed that transit time was best estimated using all available inflows (i.e., snowmelt, soil ice thaw, and rainfall). Early spring transit times were short, dominated by snowmelt and soil ice thaw and limited catchment storage when soils are predominantly frozen. However, significant and increasing mixing with water in the active layer during the summer resulted in more damped steam water variation and longer mean travel times (~1.5 years). The study has also highlighted key data needs to better constrain travel time estimates in permafrost catchments.