Estimation of the 3D correlation structure of an alluvial aquifer from surface-based multi-frequency ground-penetrating radar reflection data

Knowledge about the stochastic nature of heterogeneity in subsurface hydraulic properties is critical for aquifer characterization and the corresponding prediction of groundwater flow and contaminant transport. Whereas the vertical correlation structure of the heterogeneity is often well constrained by borehole information, the lateral correlation structure is generally unknown because the spacing between boreholes is too large to allow for its meaningful inference. There is, however, evidence to suggest that information on the lateral correlation structure may be extracted from the correlation statistics of the subsurface reflectivity structure imaged by surface-based ground-penetrating radar measurements. To date, case studies involving this approach have been limited to 2D profiles acquired at a single antenna centre frequency in areas with limited complementary information. As a result, the practical reliability of this methodology has been difficult to assess. Here, we extend previous work to 3D and consider reflection ground-penetrating radar data acquired using two antenna centre frequencies at the extensively explored and well-constrained Boise Hydrogeophysical Research Site. We find that the results obtained using the two ground-penetrating radar frequencies are consistent with each other, as well as with information from a number of other studies at the Boise Hydrogeophysical Research Site. In addition, contrary to previous 2D work, our results indicate that the surface-based reflection ground-penetrating radar data are not only sensitive to the aspect ratio of the underlying heterogeneity, but also, albeit to a lesser extent, to the so-called Hurst number, which is a key parameter characterizing the local variability of the fine-scale structure.     

Atmospheric Precursors of and Response to Anomalous Arctic Sea Ice in CMIP5 Models

This study examines pre-industrial control simulations from CMIP5 climate models in an effort to better understand the complex relationships between Arctic sea ice and the stratosphere, and between Arctic sea ice and cold winter temperatures over Eurasia. We present normalized regressions of Arctic sea-ice area against several atmospheric variables at extended lead and lag times. Statistically significant regressions are found at leads and lags, suggesting both atmospheric precursors of, and responses to, low sea ice; but generally, the regressions are stronger when the atmosphere leads sea ice, including a weaker polar stratospheric vortex indicated by positive polar cap height anomalies. Significant positive midlatitude eddy heat flux anomalies are also found to precede low sea ice. We argue that low sea ice and raised polar cap height are both a response to this enhanced midlatitude eddy heat flux. The so-called "warm Arctic, cold continents" anomaly pattern is present one to two months before low sea ice, but is absent in the months following low sea ice, suggesting that the Eurasian cooling and low sea ice are driven by similar processes. Lastly, our results suggest a dependence on the geographic region of low sea ice, with low Barents–Kara Sea ice correlated with a weakened polar stratospheric vortex, whilst low Sea of Okhotsk ice is correlated with a strengthened polar vortex. Overall, the results support a notion that the sea ice, polar stratospheric vortex and Eurasian surface temperatures collectively respond to large-scale changes in tropospheric circulation.     

Laboratory kinetics of C2H radical reactions with ethane, propane, and n-butane at T = 96-296 K: implications for Titan

The kinetics of the reactions of C₂H radical with ethane (k₁), propane (k₂), and n-butane (k₃) are studied over the temperature range of T = 96-296 K with a pulsed Laval nozzle apparatus that utilizes a pulsed laser photolysis-chemiluminescence technique. The C₂H decay profiles in the presence of both the alkane reactant and O₂ are monitored by the CH(A²Δ) chemiluminescence tracer method. The results, together with available literature data, yield the following Arrhenius expressions: k₁(T) = (0.51 ± 0.06) × 10⁻¹⁰ exp[(−76 ± 30)K/T] cm³ molecule⁻¹ s⁻¹ (T = 96-800 K), k₂(T) = (0.98 ± 0.32) × 10⁻¹⁰exp[(−71 ± 60)K/T] cm³ molecule⁻¹ s⁻¹ (T = 96-361 K), and k₃(T) = (1.23 ± 0.26) × 10⁻¹⁰ cm³ molecule⁻¹ s⁻¹ (T = 96-297 K). At T = 296 K, k₁ is measured as a function of total pressure and has little or no pressure dependence. The results from this work support a direct hydrogen abstraction mechanism for the title reactions. Implications to the atmospheric chemistry of Titan are discussed.     

The use of caesium‐137 to estimate agricultural erosion on steep slopes in a tropical watershed

The estimation of erosion and sediment delivery rates in tropical mountain watersheds is difficult and most of the methods widely used for estimating soil erosion over large areas have serious limitations. The ¹³⁷Cs approach has potential for quantifying soil erosion because it can provide retrospective estimates of long‐term (since 1963) net sediment redistribution relatively quickly. Despite its great potential, ¹³⁷Cs has not yet been used in an extensive, reconnaissance level survey of erosion in complex tropical mountain environments. The objective of this study was to examine the applicability of the ¹³⁷Cs method to estimate erosion on steep tropical agricultural lands (23 to 80% slopes) in the Nizao watershed, a humid, tropical mountain area of the Dominican Republic. In this study we (i) examine the variation of ¹³⁷Cs in ten reference sites—eight coffee groves and two forested sites—and (ii) estimate erosion from 14 cultivated fields. The soil pool of ¹³⁷Cs ranged from to 150 to 192 mBq cm⁻² on reference sites with minimal erosion. Variability among reference sites was less than expected for such complex mountain terrain. The variability within coffee and forested reference (average CV=28%) sites was similar to the variability found on grassland and forested reference sites in the temperate zone. The estimated annual soil loss from 14 sampled fields ranged from 6 to 61 t ha⁻¹ year⁻¹ with an overall mean of 26 t ha⁻¹ year⁻¹. Overall, the soil erosion estimates found using the ¹³⁷Cs method were much lower than those often assumed for such steep tropical hillsides. These erosion estimates account for soil loss since 1963 only and it seems likely that soil losses may have been much higher in earlier decades immediately after initial forest clearing earlier in the 20th century. Copyright © 2000 John Wiley & Sons, Ltd.     

Subsurface flow measurements using passive flux meters in variably-saturated cold-regions landscapes

To date, passive flux meters have predominantly been applied in temperate environments for tracking the movement of contaminants in groundwater. This study applies these instruments to reduce uncertainty in (typically instantaneous) flux measurements made in a low-gradient, wetland dominated, discontinuous permafrost environment. This method supports improved estimation of unsaturated and over-winter subsurface flows which are very difficult to quantify using hydraulic gradient-based approaches. Improved subsurface flow estimates can play a key role in understanding the water budget of this landscape.     

A new approach for linking event‐based upland sediment sources to downstream suspended sediment transport

In this study, we proposed a new approach for linking event sediment sources to downstream sediment transport in a watershed in central New York. This approach is based on a new concept of spatial scale, sub‐watershed area (SWA), defined as a sub‐watershed within which all eroded soils are transported out without deposition during a hydrological event. Using (rainfall) event data collected between July and November, 2007 from several SWAs of the studied watershed, we developed an empirical equation that has one independent variable, mean SWA slope. This equation was then used to determine event‐averaged unit soil erosion rate, Q_S/A, (in kg/km²/hr) for all SWAs in the studied watershed and calculate event‐averaged gross erosion E_ea (in kg/hr). The event gross erosion E_t (in kilograms) was subsequently computed as the product of E_ea and the mean event duration, T (in hours) determined using event hydrographs at the outlet of the studied watershed. Next, we developed two linear sediment rating curves (SRCs) for small and big events based on the event data obtained at the watershed outlet. These SRCs, together with T, allowed us to determine event sediment yield SY_e (in kilograms) for all events during the study period. By comparing E_t with SY_e, developing empirical equations (i) between E_t and SY_e and (ii) for event sediment delivery ratio, respectively, we revealed the event dynamic processes connecting sediment sources and downstream sediment transport. During small events, sediment transport in streams was at capacity and dominated by the deposition process, whereas during big events, it was below capacity and controlled by the erosion process. The key of applying this approach to other watersheds is establishing their empirical equations for Q_S/A and appropriately determining their numbers of SWAs. Copyright © 2011 John Wiley & Sons, Ltd.     

Evidence of Self‐Organized Criticality in riverbank mass failures: a matter of perspective?

A growing body of field, theoretical and numerical modelling studies suggests that predicting river response to even major changes in input variables is difficult. Rivers are seen to adjust rapidly and variably through time and space as well as changing independently of major driving variables. Concepts such as Self‐Organized Criticality (SOC) are considered to better reflect the complex interactions and adjustments occurring in systems than traditional approaches of cause and effect. This study tests the hypothesis that riverbank mass failures which occurred both prior to, and during, an extreme flood event in southeast Queensland (SEQ) in 2011 are a manifestation of SOC. Each wet‐flow failure is somewhat analogous to the ‘avalanche’ described in the initial sand‐pile experiments of Bak et al. (Physical Review Letters, 1987, 59(4), 381–384) and, due to the use of multitemporal LiDAR, the time period of instability can be effectively constrained to that surrounding the flood event. The data is examined with respect to the key factors thought to be significant in evaluating the existence of SOC including; non‐linear temporal dynamics in the occurrence of disturbance events within the system; an inverse power‐law relation between the magnitude and frequency of the events; the existence of a critical state to which the system readjusts after a disturbance; the existence of a cascading processes mechanism by which the same process can initiate both low‐magnitude and high‐magnitude events. While there was a significant change in the frequency of mass failures pre‐ and post‐flood, suggesting non‐linear temporal dynamics in the occurrence of disturbance events, the data did not fit an inverse power‐law within acceptable probability and other models were found to fit the data better. Likewise, determining a single ‘critical’ state is problematic when a variety of feedbacks and multiple modes of adjustment are likely to have operated throughout this high magnitude event. Overall, the extent to which the data supports a self‐organized critical state is variable and highly dependent upon inferential arguments. Investigating the existence of SOC, however, provided results and insights that are useful to the management and future prediction of these features. Copyright © 2014 John Wiley & Sons, Ltd.     

Compliance costs and the small fisher: A study of exiters from the New Zealand fishery

Over 3000 predominantly small-scale fishers have exited the New Zealand's quota management system (QMS) between its inception in 1986 and 2000. This study, based on the Ministry of Fisheries database and a questionnaire sent to the exiters, establishes that compliance costs in general, and those specifically related to the QMS, were one of the most consistent reasons for exit. Uncertainty about future QMS policy and the high cost of quota were also significant factors. It appears that the small fishers’ perception of high compliance cost can be supported by industry data.