Temporal and spatial pattern of thermokarst lake area changes at Yukon Flats,Alaska

To better understand the linkage between lake area change, permafrost conditions and intra‐annual and inter‐annual variability in climate, we explored the temporal and spatial patterns of lake area changes for a 422 382‐ha study area within Yukon Flats, Alaska using Landsat images of 17 dates between 1984 and 2009. Only closed basin lakes were used in this study. Among the 3529 lakes greater than 1 ha, closed basin lakes accounted for 65% by number and 50% by area. A multiple linear regression model was built to quantify the temporal change in total lake area with consideration of its intra‐annual and inter‐annual variability. The results showed that 80.7% of lake area variability was attributed to intra‐annual and inter‐annual variability in local water balance and mean temperature since snowmelt (interpreted as a proxy for seasonal thaw depth). Another 14.3% was associated with long‐term change. Among 2280 lakes, 350 lakes shrank, and 103 lakes expanded. The lakes with similar change trends formed distinct clusters, so did the lakes with similar short term intra‐annual and inter‐annual variability. By analysing potential factors driving lake area changes including evaporation, precipitation, indicators for regional permafrost change, and flooding, we found that ice‐jam flooding events were the most likely explanation for the observed temporal pattern. In addition to changes in the frequency of ice jam flooding events, the observed changes of individual lakes may be influenced by local variability in permafrost distributions and/or degradation. Copyright © 2012 John Wiley & Sons, Ltd.     

Age, geology, geophysics, and geochemistry of Mahukona Volcano, Hawai`i

The size, shape, and magmatic history of the most recently discovered shield volcano in the Hawaiian Islands, Mahukona, have been controversial. Mahukona corresponds to what was thought to be a gap in the paired sequence (Loa and Kea trends) of younger Hawaiian volcanoes (<4?Ma). Here, we present the results of marine expeditions to Mahukona where new bathymetry, sidescan sonar, gravity data, and lava samples were collected to address these controversies. Modeling of bathymetric and gravity data indicate that Mahukona is one of the smallest Hawaiian volcanoes (～6,000?km³) and that its magmatic system was not focused in a long-lived central reservoir like most other Hawaiian volcanoes. This lack of a long-lived magmatic reservoir is reflected by the absence of a central residual gravity high and the random distribution of cones on Mahukona Volcano. Our reconstructed subsidence history for Mahukona suggests it grew to at least ～270?m below sea level but probably did not form an island. New ⁴⁰Ar–³⁹Ar plateau ages range from 350 to 654?ka providing temporal constraints for Mahukona’s post-shield and shield stages of volcanism, which ended prematurely. Mahukona post-shield lavas have high ³He/⁴He ratios (16–21?Ra), which have not been observed in post-shield lavas from other Hawaiian volcanoes. Lava compositions range widely at Mahukona, including Pb isotopic values that straddle the boundary between Kea and Loa sequences of volcanoes. The compositional diversity of Mahukona lavas may be related to its relatively small size (less extensive melting) and the absence of a central magma reservoir where magmas would have been homogenized.     

Use of the SHMAK clarity tube for measuring water clarity: Comparison with the black disk method

A metre‐long transparent tube for measuring visual water clarity was developed as part of the New Zealand Stream Health Monitoring and Assessment Kit (SHMAK) for use by non‐scientists. Water clarity measured using the tube (the horizontal sighting range of a black target) was compared with that measured using the standard “black disk” method (in which a black target is viewed horizontally underwater) to determine the feasibility of estimating black disk clarities from clarity tube readings. Readings using the two methods showed good agreement for waters with clarity <c. 50 cm. In clearer waters the slope of the relationship depends on the reflectivity of the background material used in the clarity tube. With a reflective background (white, reflectance c. 75%; or grey, reflectance c. 55%), the clarity tube readings are directly proportional and approximately equal to the black disk readings. With a black background (reflectance c. 5.5%), clarity tube readings >c. 50 cm are lower than the black disk visibility and follow a logarithmic relationship with visibility, which enables useful readings to be made in clearer water. Independent testing of the relationship showed that black disk visibilities may be predicted with reasonable confidence in waters that contain very low concentrations of dissolved humic materials (yellow substance). In waters that are noticeably brown‐coloured, the relationship is less consistent especially at black disk visibilities of more than 1.5 m. Recommendations are given for the use of the clarity tube.     

Retrieval of local surface wave velocities from traffic noise – an example from the La Barge basin (Wyoming)

In regions where active source seismic exploration is constrained by limitations of energy penetration and recovery, cost and logistical concerns, or regulatory restrictions, analysis of natural source seismic data may provide an alternative. In this study, we investigate the feasibility of using locally‐generated seismic noise in the 2–6 Hz band to obtain a subsurface model via interferometric analysis. We apply this technique to three‐component data recorded during the La Barge Passive Seismic Experiment, a local deployment in south‐western Wyoming that recorded continuous seismic data between November 2008 and June 2009. We find traffic noise from a nearby state road to be the dominant source of surface waves recorded on the array and observe surface wave arrivals associated with this source up to distances of 5 kms. The orientation of the road with respect to the deployment ensures a large number of stationary points, leading to clear observations on both in‐line and cross‐line virtual source‐receiver pairs. This results in a large number of usable interferograms, which in turn enables the application of standard active source processing methods like signal processing, common offset stacking and traveltime inversion. We investigate the dependency of the interferograms on the amount of data, on a range of processing parameters and on the choice of the interferometry algorithm. The obtained interferograms exhibit a high signal‐to‐noise ratio on all three components. Rotation of the horizontal components to the radial/transverse direction facilitates the separation of Rayleigh and Love waves. Though the narrow frequency spectrum of the surface waves prevents the inversion for depth‐dependent shear‐wave velocities, we are able to map the arrival times of the surface waves to laterally varying group and phase velocities for both Rayleigh and Love waves. Our results correlate well with the known geological structure. We outline a scheme for obtaining localized surface wave velocities from local noise sources and show how the processing of passive data benefits from a combination with well‐established exploration seismology methods. We highlight the differences with interferometry applied to crustal scale data and conclude with recommendations for similar deployments.     

Continental-scale patterns in soil geochemistry and mineralogy: Results from two transects across the United States and Canada

In 2004, the US Geological Survey (USGS) and the Geological Survey of Canada (GSC) initiated a pilot study that involved collection of more than 1500 soil samples from 221 sites along two continental transects across Canada and the United States. The pilot study was designed to test and refine protocols for a soil geochemical survey of North America. The two transects crossed a wide array of soil parent materials, soil ages, climatic conditions, landforms, land covers and land uses. Sample sites were selected randomly at approximately 40-km intervals from a population defined as all soils of the continent. At each site, soils representing 0 to 5 cm depth, and the O, A, and C horizons, if present, were collected and analyzed for their near-total content of over 40 major and trace elements. Soils from 0–5 cm depth were also collected for analysis of organic compounds. Results from the transects confirm that soil samples collected at a 40-km spacing reveal coherent, continental- to subcontinental-scale geochemical and mineralogical patterns that can be correlated to aspects of underlying soil parent material, soil age and climate influence. The geochemical data also demonstrate that at the continental-scale the dominance of any of these major factors that control soil geochemistry can change across the landscape. Along both transects, soil mineralogy and geochemistry change abruptly with changes in soil parent materials. However, the chemical influence of a soil’s parent material can be obscured by changing climatic conditions. For the transects, increasing precipitation from west to east and increasing temperature from north to south affect both soil mineralogy and geochemistry because of climate effects on soil weathering and leaching, and plant productivity. Regional anomalous metal concentrations can be linked to natural variations in soil parent materials, such as high Ni and Cr in soils developed on ultramafic rocks in California or high P in soils formed on weathered Ordovician limestones in central Kentucky. On local scales, anomalous metal concentrations recognized in soil profiles, such as high P in soils from animal confinement sites, are consistent with local anthropogenic disturbances. At a larger scale, the distribution of Hg across the west to east transect demonstrates that it can be difficult to distinguish between natural or anthropogenic contributions and that many factors can contribute to an element’s spatial distribution.Only three samples in a subset of seventy-three 0–5 cm depth soil samples from the north to south transect had organochlorine pesticides values above the method detection limit, apparently related to historic usage of the pesticides DDT and dieldrin.     

A coupled oscillator model of shelf and ocean tides

The resonances of tides in the coupled open ocean and shelf are modeled by a mechanical analogue consisting of a damped driven larger mass and spring (the open-ocean) connected to a damped smaller mass and spring (the shelf). When both masses are near resonance, the addition of even a very small mass can significantly affect the oscillations of the larger mass. The influence of the shelf is largest if the shelf is resonant with weak friction. In particular, an increase of friction on a near-resonant shelf can, perhaps surprisingly, lead to an increase in ocean tides. On the other hand, a shelf with large friction has little effect on ocean tides. Comparison of the model predictions with results from numerical models of tides during the ice ages, when lower sea levels led to a much reduced areal extent of shelves, suggests that the predicted larger tidal dissipation then is related to the ocean basins being close to resonance. New numerical simulations with a forward global tide model are used to test expectations from the mechanical analogue. Setting friction to unrealistically large values in Hudson Strait yields larger North Atlantic _M2$M_2$ amplitudes, very similar to those seen in a simulation with the Hudson Strait blocked off. Thus, as anticipated, a shelf with very large friction is nearly equivalent in its effect on the open ocean to the removal of the shelf altogether. Setting friction in shallow waters throughout the globe to unrealistically large values yields even larger open ocean tidal amplitudes, similar to those found in simulations of ice-age tides. It thus appears that larger modeled tides during the ice ages can be a consequence of enhanced friction in shallower water on the shelf in glacial times as well as a reduced shelf area then. Single oscillator and coupled oscillator models for global tides show that the maximum extractable power for human use is a fraction of the present dissipation rate, which is itself a fraction of global human power consumption.     

Density variations in the lower thermosphere from analysis of the AE-C accelerometer measurements

Density measurements at 140, 160, 180 and 200km from the AE-C accelerometer experiment during 1974 are analyzed by least-squares multiple linear regression. The resulting empirical model is compared to the widely-used Jacchia (1971) model which is primarily based on satellite orbital decay data above 200 km. The semiannual and geomagnetic (K_p) variations derived from the AE-C data set are in good agreement with Jacchia's predictions. However, individual magnetic storm responses often exceed such average representations of the geomagnetic effect by more than a factor of two. The density is more variable with respect to the daily 10.7cm solar flux than Jacchia's model predicts. The regression analysis reveals a positive correlation with the absolute value of the geographic latitude. The statistical weakness of this latter effect reflects the transient, aperiodic nature of high-latitude heating events related to magnetic activity, which can only be properly analyzed on an individual basis. Finally, the sub-solar bulge effect enters the regression equation with about half the amplitude predicted by the Jacchia model at these heights.     

Estimating groundwater recharge in lowland watersheds

Little attention has been given to the role of groundwater in the hydrological cycle of lowland watersheds. Our objective in this study was to estimate total recharge to groundwater by analysing water table response to storm events and the rate at which water was transferred into the shallow aquifer. This was conducted at three sites in a rural watershed in the lower Atlantic coastal plain near Charleston, South Carolina, USA. A novel version of the water table fluctuation method was used to estimate total recharge to the shallow aquifer by comparing hourly data of water table position following storm events and measuring water table recession behavior, rather than subjective graphical analysis methods. Also, shallow aquifer recharge rates (vertical fluxes) were estimated using Darcy's Law by comparing static water levels in a water table well and in a shallow piezometer during dry periods. The total annual recharge estimated ranged from 107 ± 39 mm·yr^–1 (5–10% of annual precipitation) at a poorly drained topographic low area to 1140 ± 230 mm·yr^–1 (62–94% of annual precipitation) for a moderately well‐drained upland site. The average aquifer recharge rate was 114 ± 60 mm·yr^–1, which is similar to previous estimations of base flow for the ephemeral third‐order streams in this watershed. The difference in the two methods may have been caused by processes not accounted for in the Darcy flux method, soil moisture deficits, and average evapotranspiration demand, which is about 1100 mm·yr^–1 for this region. Although other factors also can affect partitioning of recharge, an integrated approach to inspecting easily gathered groundwater data can provide information on an often neglected aspect of water budget estimation. We also discuss the effects of land use change on recharge reduction, given a typical development scenario for the region. Copyright © 2011 John Wiley & Sons, Ltd.     

Statistics of the variations of the high-energy electron population between 7 and 28 jovian radii as measured by the Galileo spacecraft

Recent measurements of the high-energy, omni-directional electron environment by the Galileo spacecraft Energetic Particle Detector (EPD) have been analyzed in the range from 7 to 28 Jupiter radii. 10-min averages of these data between Jupiter orbit insertion in 1995 to the end of the mission have been analyzed to provide estimates of the electron differential fluxes at 1.5, 2, and 11 MeV in the jovian equatorial plane as a function of radial distance. These data provide a long term picture of the variations in the high-energy electron environment over the ∼8 years of the Galileo mission. This paper reviews those measurements and the statistics associated with them for the 8 year period. In general, the data variations are well behaved with variations being within a factor of ∼2 of a median value at a given distance from Jupiter. These results are analyzed in detail and the orbit variations discussed in the context of the overall data set. The results of this analysis of the long-term statistical variations in high-energy electron fluxes are directly applicable to models that estimate the effects of the radiation environment on Jupiter's moons and their atmospheres as they permit estimates of the possible range of radiation effects that might be expected.