Microwave Remote Sensing of the Temperature and Distribution of Sulfur Compounds in the Lower Atmosphere of Venus

A multi-wavelength radio frequency observation of Venus was performed on April 5, 1996, with the Very Large Array to investigate potential variations in the vertical and horizontal distribution of temperature and the sulfur compounds sulfur dioxide (SO₂) and sulfuric acid vapor (H₂SO₄(g)) in the atmosphere of the planet. Brightness temperature maps were produced which feature significantly darkened polar regions compared to the brighter low-latitude regions at both observed frequencies. This is the first time such polar features have been seen unambiguously in radio wavelength observations of Venus. The limb-darkening displayed in the maps helps to constrain the vertical profile of H₂SO₄(g), temperature, and to some degree SO₂. The maps were interpreted by applying a retrieval algorithm to produce vertical profiles of temperature and abundance of H₂SO₄(g) given an assumed sub-cloud abundance of SO₂. The results indicate a substantially higher abundance of H₂SO₄(g) at high latitudes (above 45°) than in the low-latitude regions. The retrieved temperature profiles are up to 25 K warmer than the profile obtained by the Pioneer Venus sounder probe at altitudes below 40 km (depending on location and assumed SO₂ abundance). For 150 ppm of SO₂, it is more consistent with the temperature profile obtained by Mariner 5, extrapolated to the surface via a dry adiabat. The profiles obtained for H₂SO₄(g) at high latitudes are consistent with those derived from the Magellan radio occultation experiments, peaking at around 8 ppm at an altitude of 46 km and decaying rapidly away from that altitude. At low latitudes, no significant H₂SO₄(g) is observed, regardless of the assumed SO₂ content. This is well below that measured by Mariner 10 (Lipa and Tyler 1979, Icarus39, 192-208), which peaked at ∼14 ppm near 47 km. Our results favor ≤100 ppm of SO₂ at low latitudes and ≤50 ppm in polar regions. The low-latitude value is statistically consistent with the results of Bézard et al. (1983, Geophs. Res. Lett.20, 1587-1590), who found that a sub-cloud SO₂ abundance of 130±40 ppm best matched their observations in the near-IR. The retrieved temperature profile and higher abundance of H₂SO₄(g) in polar regions are consistent with a strong equatorial-to-polar, cloud-level flow due to a Hadley cell in the atmosphere of Venus.     

Analytical model for aquifer response incorporating distributed stream leakage

An analytical model of stream/aquifer interaction is proposed that predicts drawdown in an aquifer with leakage from a finite-width stream induced by pumping from a well. The model is formulated based on the assumptions of stream partial penetration, a semipervious streambed, and distributed recharge across a finite-width stream. Advantages of the analytical solution include its simple structure, consisting of the Theis well function with integral modifications. The solution is derived for the semi-infinite domain between the stream and pumping well, which is of primary interest to hydrogeologists. Previous stream/aquifer analytical models are compared to the analytical solution based on dimensionless drawdown profiles. Drawdown in the aquifer near a wide stream was found to be less than that predicted by a solution that ignored stream width. Deviations between the proposed analytical solutions and previous solutions increase as stream width increases. For a hypothetical stream/aquifer system, the proposed analytical solution was equivalent to prior solutions when the ratio of the distance between the stream and aquifer to the stream width was greater than 25. This analytical solution may provide improved estimates of aquifer and streambed leakage parameters by curve fitting experimental field drawdown data.     

Influence of wave groups on SSC patterns over vortex ripples

Estimates of spatial and temporal variations in suspended sand concentrations (SSC) made with a multi-transducer Acoustic Backscatter Sensor (ABS) under a repeated wave group over a mobile rippled bed in the wave research flume at the National Hydraulics Laboratory in Ottawa, Canada, reveal an number of complex and intriguing patterns. Ensemble averages of 8 nearly identical wave groups provided much more robust estimates of SSC and allowed a detailed examination of the wave group effects. The largest SSC near the bed (< 0.10 m) occurs in phase with the largest waves in the group. Above approximately 0.10 m elevation, SSC lags behind the near bed SSC by as much as 2–3 waves; introducing significant curvature (on a semi-log plot) to the SSC profile. The log linear segments of the SSC profile grow and decay systematically on the scale of the wave group. The range in lengths of log-linear profile segments ( 0.03–0.355 m) suggest that the boundary layer thickness also fluctuates throughout the passage of the wave group. Furthermore, there are significant variations in the patterns of SSC, which occur under the largest and smallest waves in the group. Under the largest waves vertical bands of alternating high and low SSC produce an intra-wave modulation in the upper water column ( 0.075–0.30 m). The equivalent horizontal excursion of these bands scales to the ripple length. Under the smaller waves the intra-wave modulation of the SSC disappears and is replaced by temporally homogenous suspension that expands vertically through several individual wave cycles. The former pattern of homogenous suspension appears to be associated with growth of a boundary layer due to the persistent uni-directional horizontal flow during this part of the group together with the persistence of antecedent bed generated turbulence and vorticity which maintains the suspension. The latter pattern of bands of high and low SSC indicates a strong temporal and spatial constraint on the SSC (phase coupling) induced by the presence of the bedforms which may be enhanced by strong reversals in both flow and vorticity under the large waves in the group.     

MIDAS-W: a workstation-based incoherent scatter radar data acquisition system

The Millstone Hill Incoherent Scatter Data Acquisition System (MIDAS) is based on an abstract model of an incoherent scatter radar. This model is implemented in a hierarchical software system, which serves to isolate hardware and low-level software implementation details from higher levels of the system. Inherent in this is the idea that implementation details can easily be changed in response to technological advances. MIDAS is an evolutionary system, and the MIDAS hardware has, in fact, evolved while the basic software model has remained unchanged. From the earliest days of MIDAS, it was realized that some functions implemented in specialized hardware might eventually be implemented by software in a general-purpose computer. MIDAS-W is the realization of this concept. The core component of MIDAS-W is a Sun Microsystems UltraSparc 10 workstation equipped with an Ultrarad 1280 PCI bus analog to digital (A/D) converter board. In the current implementation, a 2.25 MHz intermediate frequency (IF) is bandpass sampled at 1 s intervals and these samples are multicast over a high-speed Ethernet which serves as a raw data bus. A second workstation receives the samples, converts them to filtered, decimated, complex baseband samples and computes the lag-profile matrix of the decimated samples. Overall performance is approximately ten times better than the previous MIDAS system, which utilizes a custom digital filtering module and array processor based correlator. A major advantage of MIDAS-W is its flexibility. A portable, single-workstation data acquisition system can be implemented by moving the software receiver and correlator programs to the workstation with the A/D converter. When the data samples are multicast, additional data processing systems, for example for raw data recording, can be implemented simply by adding another workstation with suitable software to the high-speed network. Testing of new data processing software is also greatly simplified, because a workstation with the new software can be added to the network without impacting the production system. MIDAS-W has been operated in parallel with the existing MIDAS-1 system to verify that incoherent scatter measurements by the two systems agree. MIDAS-W has also been used in a high-bandwidth mode to collect data on the November, 1999, Leonid meteor shower.     

Influence of uncertain boundary conditions and model structure on flood inundation predictions

In this study, the GLUE methodology is applied to establish the sensitivity of flood inundation predictions to uncertainty of the upstream boundary condition and bridges within the modelled region. An understanding of such uncertainties is essential to improve flood forecasting and floodplain mapping. The model has been evaluated on a large data set. This paper shows uncertainty of the upstream boundary can have significant impact on the model results, exceeding the importance of model parameter uncertainty in some areas. However, this depends on the hydraulic conditions in the reach e.g. internal boundary conditions and, for example, the amount of backwater within the modelled region. The type of bridge implementation can have local effects, which is strongly influenced by the bridge geometry (in this case the area of the culvert). However, the type of bridge will not merely influence the model performance within the region of the structure, but also other evaluation criteria such as the travel time. This also highlights the difficulties in establishing which parameters have to be more closely examined in order to achieve better fits. In this study no parameter set or model implementation that fulfils all evaluation criteria could be established. We propose four different approaches to this problem: closer investigation of anomalies; introduction of local parameters; increasing the size of acceptable error bounds; and resorting to local model evaluation. Moreover, we show that it can be advantageous to decouple the classification into behavioural and non-behavioural model data/parameter sets from the calculation of uncertainty bounds.     

Pore scale modeling of immiscible and miscible fluid flows using smoothed particle hydrodynamics

A numerical model based on smoothed particle hydrodynamics (SPH) was developed and used to simulate immiscible and miscible fluid flows in porous media and to study effects of pore scale heterogeneity and anisotropy on such flows.     

Gullying and erosion control at archaeological sites in Grand Canyon,Arizona

Gully erosion of cultural sites in Grand Canyon National Park is an urgent management problem that has intensified in recent decades, potentially related to the effects of Glen Canyon Dam. We studied 25 gullies at nine sites in Grand Canyon over the 2002 monsoon–erosion season to better understand the geomorphology of the gully erosion and the effectiveness of erosion‐control structures (ECS) installed by the park under the direction of the Zuni Conservation Program. Field results indicate that Hortonian overland flow leads to concentrated flow in gullies and erosion focused at knickpoints along channels as well as at gully heads. Though groundcover type, soil shear strength and permeability vary systemat‐ically across catchments, gradient and, to a lesser degree, contributing drainage area seem to be the first‐order controls on gully extent, location of new knickpoints, and ECS damage. The installed ECS do reduce erosion relative to reaches without them and initial data suggest woody checkdams are preferable to rock linings, but maintenance is essential because damaged structures can exacerbate erosion. Topographic data from intensive field surveys and detailed photogrammetry provide slope–contributing area data for gully heads that have a trend consistent with previous empirical and theoretical formulations from a variety of landscapes. The same scaling holds below gully heads for knickpoint and ECS topographic data, with threshold coefficients the lowest for gully heads, slightly higher for knickpoints, and notably higher for damaged ECS. These topographic thresholds were used with 10‐cm digital elevation models to create simple predictive models for gully extent and structure damage. The model predictions accounted for the observed gullies but there are also many false‐positives. Purely topographical models are probably inadequate at this scale and application, but models that also parameterize the variable soil properties across sites would be useful for predicting erosion problems and ECS failure. Copyright © 2006 John Wiley & Sons, Ltd.     

青藏高原典型地区的地貌量化分析——兼对高原“夷平面”的讨论

利用新近公布的SRTM数字高程模型(DEM),选取青藏高原北部及高原内、外流区域进行大尺度定量地貌分析。分析表明,青藏高原不同地区的地貌差异反映了它们在地貌演化上的阶段性。在高原北部的祁连山地区,局部地区绝对高程增加的同时,地势起伏反而变缓。这些地区水系的发育程度普遍较低,剥蚀物质往往只经过近距离的搬运仍滞留在逆冲褶皱带内,山间盆地和平地成为山间小河的侵蚀堆积基准,因此“削高填低”的过程有效地降低了局部地形起伏。高原平坦地势是伴随着造山过程及之后的高海拔侵蚀基准和内流型水系条件下“削高填低”剥蚀过程的结果。我们认为高原内部为现今仍在承受剥蚀的地势平坦面。它的形成具穿时性,是内流型水系河流下切侵蚀能力弱化的结果,不直接反映海拔的高低。如果平坦侵蚀面的形成与海拔高程无必然的关联,或侵蚀面可以在任何海拔高度形成,而不一定代表以海平面为基准的准平原,那么它就不能作为一个可靠的参照面用于直接示踪和约束高原的抬升量和抬升时间。     

Ground viscosity and stiffness measurements for near surface seismic velocity

In land surveys, the weathering layer can often distort the seismic signal due to it passing through rapid velocity and density changes, dispersion, scattering and inelastic absorption. In a simple spring‐dashpot model for the earth response, an equivalent medium groups these complex phenomena into two parameters only; these are called ground viscosity and ground stiffness. The most recent controllers for vibrators can estimate both parameters. To validate these measurements, Saudi Aramco conducted an experiment measuring ground viscosity and stiffness from two different vibrator control systems over an area of varying terrain conditions, including unconsolidated sand and limestone outcrop. The two systems measured different values, but detected similar trends that correlated well with weathering conditions and surface geology, e.g. lower viscosity values on the outcrop than on the sand. The ratio of ground viscosity to ground stiffness can approximate the shallow S‐wave velocity, which we converted into P‐wave velocity through calibration with sparse uphole data. Static corrections incorporating this velocity information somewhat improved the focusing of seismic time sections. This new approach does not require additional acquisition efforts, and can model shallow complex formations in arid areas where classical static methods often fail.