Determination of physical and radiant meteor properties using PFISR interferometry measurements of head echoes

Over the past decade, High Power and Large Aperture (HPLA) radars have been widely utilized for the study of sub-millimeter extraterrestrial particles via the detection of the meteor head-echo. These observations have been a successful tool in the study of the sporadic meteor background, however, they have been limited by the lack of precise knowledge of the particle's location within the radar beam and its absolute trajectory and velocity. This limitation prevents for example the accurate determination of the meteors radiant and orbit. Interferometry measurements of the head-echo has been proven to be a detection technique that satisfies this need. Unfortunately very few radars are capable of performing them. We have developed a methodology which takes advantage of the multi-receiving capabilities of the 450 MHz Poker Flat Incoherent Scatter Radar (PFISR) enabling us to utilize the phased array of crossed-dipoles as an interferometer. This new PFISR capability allows us to determine the instantaneous position of meteors within the radar beam. This enables us to determine absolute velocities and ultimately meteor radiant and orbit around the Sun. In this work, we present initial results from 9 h of observations during which 142 particles were individually detected by the three different receiving channels simultaneously. For these meteors absolute velocities were obtained and meteor dynamical, physical and radiant properties were derived.     

A discontinuous finite element baroclinic marine model on unstructured prismatic meshes

We describe the space discretization of a three-dimensional baroclinic finite element model, based upon a discontinuous Galerkin method, while the companion paper (Comblen et al. 2010a) describes the discretization in time. We solve the hydrostatic Boussinesq equations governing marine flows on a mesh made up of triangles extruded from the surface toward the seabed to obtain prismatic three-dimensional elements. Diffusion is implemented using the symmetric interior penalty method. The tracer equation is consistent with the continuity equation. A Lax–Friedrichs flux is used to take into account internal wave propagation. By way of illustration, a flow exhibiting internal waves in the lee of an isolated seamount on the sphere is simulated. This enables us to show the advantages of using an unstructured mesh, where the resolution is higher in areas where the flow varies rapidly in space, the mesh being coarser far from the region of interest. The solution exhibits the expected wave structure. Linear and quadratic shape functions are used, and the extension to higher-order discretization is straightforward.     

Anthropogenic influences on Zuni agricultural soils

Farmers in the Zuni area of the semiarid American Southwest have successfully cultivated maize and other crops for over three millennia without using artificial fertilizers. Zuni agricultural fields are among the oldest, more or less continuously cultivated areas in the United States. Traditional Zuni agriculture is based on runoff farming, a system whereby runoff and organic‐rich sediment generated in small watersheds are captured and directed onto fields for crop use. We conducted a study to compare soil properties associated with paired and unpaired cultivated, abandoned, and uncultivated fields to evaluate the long‐term effects of cultivation on soil quality. Sampling and analytical methods of this research are especially applicable to geoarchaeological studies of anthropogenic effects on soil fertility and agricultural sustainability in ancient and traditional historical farming systems. Results of the Zuni soil study indicate that cultivation has altered some soil properties, including bulk density, organic carbon, total nitrogen, and C:N ratios in paired fields, but there is no indication that agricultural soils are degraded. This assessment supports the perception of Zuni farmers that long‐term cultivation has not caused a decline in agricultural productivity. © 2005 Wiley Periodicals, Inc.     

Effects of Stronach Dam removal on fluvial geomorphology in the Pine River, Michigan, United States

Although dam removal has been increasingly used as an option in dam management, and as a river restoration tool, few studies provide detailed quantitative assessment of the geomorphological response of rivers to dam removal. In this study, we document the response of the Pine River, Michigan, to the gradual removal of Stronach Dam. In 1996, prior to the initiation of removal, 31 permanent cross-sectional transects were established in the 10-km study area. These transects were surveyed annually during the course of the removal (1996–2003) and for the three years following removal (2004–2006). Dam removal resulted in progressive headcutting of sediments in the former impoundment, extending upstream 3.89 km of the dam. Over the course of the 10 years since dam removal was initiated, a net total of 92 000 m³ of sediment erosion occurred. The majority of sediments stored in the former reservoir remained in place, with only 12% of the estimated reservoir sediment fill being eroded. Approximately 14% of the net erosion was deposited within the stream channel 1 km downstream of the dam location, with the remainder being transported further downstream or deposited in the floodplain. Sediment fill incision resulted in a narrower and deeper channel upstream, with higher mean water velocity and somewhat coarser substrates. Downstream deposition resulted in a wider and shallower channel, with little change in substrate size composition. Counter-intuitively, water velocity also increased downstream because of the increased slope that developed. Prior to removal, bedforms in the former impoundment were dominated by runs but are showing signs of restoration toward reference conditions. Continuing changes in river geomorphology are evident even three years following removal and are likely to occur for years to come.     

Storage of sediment‐associated heavy metals along the channelized Odra River,Poland

The importance of long‐term storage of heavy metals in groyne fields, functioning over 150 years, is investigated for the River Odra (Oder), western Poland. Construction of groynes along the Odra preceded rapid development of heavy industrialization in the largest coal mine districts in Poland and the Czech Republic that resulted in persistent riverine pollution. The 187 km long Middle Odra reach was repeatedly channelized from the first half of the eighteenth century to the turn of the twentieth century, during which time partially filled groyne fields were dissected by new bank lines and groyne systems, with older groyne fields partially keyed into the floodplain. Consequently, concentrations of zinc, lead, cadmium, and copper within historically deposited groyne field sediments exceed local geochemical background levels by more than 60, 40, 15 and 10 times, respectively. Sediments contaminated with heavy metals occur within three distinctive geomorphic zones: zone I is up to 250 m wide and furthest from the present channel, comprising decimeter‐thick polluted sediments, overlying eighteenth century sand and gravel bars; zone II represents the former nineteenth century groyne fields, with widths between 10 and 100 m, filled with as much as 3 m of polluted sediments; zone III represents the twentieth century groyne fields, which are several to a dozen metres wide and filled with polluted sediments averaging depths of more than 2 m. This investigation indicates that large and extensive sediment quantities of moderately polluted sediments are stored immediately along the banks of the River Odra. These sediments could be a significant secondary pollution source and therefore careful maintenance of contemporary bank protection structures is required. Copyright © 2009 John Wiley & Sons, Ltd.     

Impacts of stormwater runoff in the Southern California Bight: Relationships among plume constituents

The effects from two winter rain storms on the coastal ocean of the Southern California Bight were examined as part of the Bight ‘03 program during February 2004 and February–March 2005. The impacts of stormwater from fecal indicator bacteria, water column toxicity, and nutrients were evaluated for five major river discharges: the Santa Clara River, Ballona Creek, the San Pedro Shelf (including the Los Angeles, San Gabriel, and Santa Ana Rivers), the San Diego River, and the Tijuana River. Exceedances of bacterial standards were observed in most of the systems. However, the areas of impact were generally spatially limited, and contaminant concentrations decreased below California Ocean Plan standards typically within 2–3 days. The largest bacterial concentrations occurred in the Tijuana River system where exceedances of fecal indicator bacteria were noted well away from the river mouth. Maximum nitrate concentrations (~40 μM) occurred in the San Pedro Shelf region near the mouth of the Los Angeles River. Based on the results of general linear models, individual sources of stormwater differ in both nutrient concentrations and the concentration and composition of fecal indicator bacteria. While nutrients appeared to decrease in plume waters due to simple mixing and dilution, the concentration of fecal indicator bacteria in plumes depends on more than loading and dilution rates. The relationships between contaminants (nutrients and fecal indicator bacteria) and plume indicators (salinity and total suspended solids) were not strong indicating the presence of other potentially important sources and/or sinks of both nutrients and fecal indicator bacteria. California Ocean Plan standards were often exceeded in waters containing greater than 10% stormwater (<28–30 salinity range). The median concentration dropped below the standard in the 32–33 salinity range (1–4% stormwater) for total coliforms and Enterococcus spp. and in the 28–30 salinity range (10–16% stormwater) for fecal coliforms. Nutrients showed a similar pattern with the highest median concentrations in water with greater than 10% stormwater. Relationships between colored dissolved organic matter (CDOM) and salinity and between total suspended solids and beam attenuation indicate that readily measurable, optically active variables can be used as proxies to provide at least a qualitative, if not quantitative, evaluation of the distribution of the dissolved, as well as the particulate, components of stormwater plumes. In this context, both CDOM absorption and the beam attenuation coefficient can be derived from satellite ocean color measurements of inherent optical properties suggesting that remote sensing of ocean color should be useful in mapping the spatial areas and durations of impacts from these contaminants.     

A fluid dynamics approach to modelling the 18th March 2007 lahar at Mt. Ruapehu,New Zealand

Lahars are water-sediment mass flows from a volcanic source. They can be triggered by a variety of mechanisms and span a continuum of flow rheology and hydraulic properties, even within the same event. Lahars are extremely powerful landscaping agents and represent a considerable hazard potential. However, this highly dynamic character and a lack of direct measurements has made modelling lahars difficult. This study therefore applies a fluid dynamics model; Delft3D, to analyse the 18th March 2007 dam break lahar at Mount Ruapehu, New Zealand. The modelled lahar routed through the Whangaehu gorge in ~30 min, crossed the Whangaehu fan in ~60 min, and then over a further 3 h travelled an additional ~22 km distance along the Whangaehu River to the Tangiwai bridge. The modelled mean frontal velocity was 6.5 m s⁻¹ along the gorge although peak velocity reached up to 19.6 m s⁻¹. The modelled lahar flow front progressively slowed across the fan but along the River it accelerated from 2.1–3.3 m s⁻¹. Calculated peak velocity along the River was <4.5 m s⁻¹. These results generally compare well with gauged records, with historical records, and with other modelling approaches. However, discrepancies in frontal velocity and time to peak stage arise due to (1) specifying roughness, which arises from slope variations between adjacent computational nodes, and which is stage-dependant, and (2) due to rapid topographic changes that produce frequent hydraulic jumps, which are inadequately accommodated in the numerical scheme. The overall pattern of discharge attenuation, and of relationships between topographic and hydraulic variables, is similar to that calculated for lahars on other volcanoes. This modelling method could be applied at other similar sites where a likely source hydrograph and high-resolution topographic data are available. These results have important implications for hazard management at Ruapehu and for examining geomorphic and sedimentary impacts of this lahar.     

Thermal characterization of the Vulcano fumarole field

Ground-based thermal infrared surveys can contribute to complete heat budget inventories for fumarole fields. However, variations in atmospheric conditions, plume condensation and mixed-pixel effects can complicate vent area and temperature measurements. Analysis of vent temperature frequency distributions can be used, however, to characterise and quantify thermal regions within a field. We examine this using four thermal infrared thermometer and thermal image surveys of the Vulcano Fossa fumarole field (Italy) during June 2004 and July 2005. These surveys show that regions occupied by low temperature vents are characterised by distributions that are tightly clustered around the mean (i.e., the standard deviation is low), highly peaked (positive kurtosis) and skewed in the low temperature direction (negative skewness). This population is associated with wet fumaroles, where boiling controls maximum temperature to cause a narrow distribution with a mode at 90–100°C. In contrast, high temperature vent regions have distributions that are widely spread about the mean (i.e., the standard deviation is high), relatively flat (negative kurtosis) and skewed in the high temperature direction (positive skewness). In this dry case, fumaroles are water-free so that maximum temperatures are not fixed by boiling. As a result greater temperature variation is possible. We use these results to define two vent types at Vulcano on the basis of their thermal characteristics: (1) concentrated (localized) regions of high temperature vents, and (2) dispersed low temperature vents. These occur within a much larger region of diffuse heat emission across which surfaces are heated by steam condensation, the heat from which causes elevated surface temperatures. For Vulcano's lower fumarole zone, high and low temperature vents occupied total areas of 3 and 6 m², respectively, and occurred within a larger (430 m²) vent-free zone of diffuse heat emission. For this lower zone, we estimate that 21–43 × 10³ W of heat was lost by diffuse heat emission. A further 4.5 × 10³ W was lost by radiation from high temperature vents, and 6.5 × 10³ W from low temperature vents. Thus, radiative heat losses from high and low temperature vents within Vulcano's lower fumarole zone respectively account for 10% and 15% of the total heat lost from this zone. This shows that radiation from open vents can account for a non-trivial portion of the total fumarole field heat budget.     

An empirical function to estimate the depths of linear hot sources: Laboratory modeling and field measurements of lava tubes

Estimating depths of buried lava tubes is important for determining the thermal budgets and effusion rates of basaltic volcanic systems. This research used a laboratory experiment scaled to a lava tube system to measure the 3D temperature field surrounding a hot viscous fluid flowing through a buried glass tube while varying conditions such as flow rate and temperature. The depth of the glass tube was changed for different experimental runs. Numerical techniques were applied to model the laboratory experiment. The surface thermal distributions from 166 thermal traverses, constrained to a depth to width ratio of 0.6 to 1.6, were analyzed to empirically derive a depth estimation function using regression techniques. This “Linear Anomaly Surface Transect (LAST)” depth function is a scaleable depth estimation technique which can be solved with thermal imaging data alone. The minimum temperature, maximum temperature and width of a Lorentzian distribution fit to a surface thermal transect, are the only inputs required for the LAST function to estimate the depths of the hot source. The input parameters were then applied to non-laboratory situations including the Kuhio lava tube system in Hawai’i. The LAST function produced depth estimates of ∼ 0.3 m for the Kuhio lava tube in Hawai’i, which did not agree with observations on the ground. This is the result of the complex composition and geometry of an actual lava tube where heat transfer is controlled by more than a simple fluid filling a tube, but also by convection of gasses and fluids in a partially filled passage. Though not effective for lava tubes at this time, the model provides promising results for simple cases applied to engineering and underground coal fires.