High‐accuracy practical spline‐based 3D and 2D integral transformations in potential‐field geophysics

Potential, potential field and potential‐field gradient data are supplemental to each other for resolving sources of interest in both exploration and solid Earth studies. We propose flexible high‐accuracy practical techniques to perform 3D and 2D integral transformations from potential field components to potential and from potential‐field gradient components to potential field components in the space domain using cubic B‐splines. The spline techniques are applicable to either uniform or non‐uniform rectangular grids for the 3D case, and applicable to either regular or irregular grids for the 2D case. The spline‐based indefinite integrations can be computed at any point in the computational domain. In our synthetic 3D gravity and magnetic transformation examples, we show that the spline techniques are substantially more accurate than the Fourier transform techniques, and demonstrate that harmonicity is confirmed substantially better for the spline method than the Fourier transform method and that spline‐based integration and differentiation are invertible. The cost of the increase in accuracy is an increase in computing time. Our real data examples of 3D transformations show that the spline‐based results agree substantially better or better with the observed data than do the Fourier‐based results. The spline techniques would therefore be very useful for data quality control through comparisons of the computed and observed components. If certain desired components of the potential field or gradient data are not measured, they can be obtained using the spline‐based transformations as alternatives to the Fourier transform techniques.     

In situ tensile fracture toughness of surficial cohesive marine sediments

This study reports the first in situ measurements of tensile fracture toughness, K _IC, of soft, surficial, cohesive marine sediments. A newly developed probe continuously measures the stress required to cause tensile failure in sediments to depths of up to 1 m. Probe measurements are in agreement with standard laboratory methods of K _IC measurements in both potter’s clay and natural sediments. The data comprise in situ depth profiles from three field sites in Nova Scotia, Canada. Measured K _IC at two muddy sites (median grain size of 23–50 μm) range from near zero at the sediment surface to >1,800 Pa m^1/2 at 0.2 m depth. These profiles also appear to identify the bioturbated/mixed depth. K _IC for a sandy site (>90% sand) is an order of magnitude lower than for the muddy sediments, and reflects the lack of cohesion/adhesion. A comparison of K _IC, median grain size, and porosity in muddy sediments indicates that consolidation increases fracture strength, whereas inclusion of sand causes weakening; thus, sand-bearing layers can be easily identified in K _IC profiles. K _IC and vane-measured shear strength correlate strongly, which suggests that the vane measurements should perhaps be interpreted as shear fracture toughness, rather than shear strength. Comparison of in situ probe-measured values with K _IC of soils and gelatin shows that sediments have a K _IC range intermediate between denser compacted soils and softer, elastic gelatin.     

Contrasting rainfall generated debris flows from adjacent watersheds at Forest Falls, southern California, USA

Debris flows are widespread and common in many steeply sloping areas of southern California. The San Bernardino Mountains community of Forest Falls is probably subject to the most frequently documented debris flows in southern California. Debris flows at Forest Falls are generated during short-duration high-intensity rains that mobilize surface material. Except for debris flows on two consecutive days in November 1965, all the documented historic debris flows have occurred during high-intensity summer rainfall, locally referred to as ‘monsoon’ or ‘cloudburst’ rains. Velocities of the moving debris range from about 5 km/h to about 90 km/h. Velocity of a moving flow appears to be essentially a function of the water content of the flow. Low velocity debris flows are characterized by steep snouts that, when stopped, have only small amounts of water draining from the flow. In marked contrast are high-velocity debris flows whose deposits more resemble fluvial deposits. In the Forest Falls area two adjacent drainage basins, Snow Creek and Rattlesnake Creek, have considerably different histories of debris flows. Snow Creek basin, with an area about three times as large as Rattlesnake Creek basin, has a well developed debris flow channel with broad levees. Most of the debris flows in Snow Creek have greater water content and attain higher velocities than those of Rattlesnake Creek. Most debris flows are in relative equilibrium with the geometry of the channel morphology. Exceptionally high-velocity flows, however, overshoot the channel walls at particularly tight channel curves. After overshooting the channel, the flows degrade the adjacent levee surface and remove trees and structures in the immediate path, before spreading out with decreasing velocity. As the velocity decreases the clasts in the debris flows pulverize the up-slope side of the trees and often imbed clasts in them. Debris flows in Rattlesnake Creek are relatively slow moving and commonly stop in the channel. After the channel is blocked, subsequent debris flows cut a new channel upstream from the blockage that results in the deposition of new debris-flow deposits on the lower part of the fan. Shifting the location of debris flows on the Rattlesnake Creek fan tends to prevent trees from becoming mature. Dense growths of conifer seedlings sprout in the spring on the late summer debris flow deposits. This repeated process results in stands of even-aged trees whose age records the age of the debris flows.     

An ecohydrological approach to predicting regional woody species distribution patterns in dryland ecosystems

This paper presents a quantitative ecohydrological framework for predicting regional distribution patterns of woody species in dryland ecosystems. The framework is based on an existing stochastic model for the daily mass balance of water that represents the interactions between soils, climate, and vegetation. Individual species selection is based on an optimality trade-off hypothesis, which states that dryland vegetation patterns are constrained by maximization of water use and simultaneous minimization of water stress. The relative importance of water use and stress avoidance to the overall fitness of three Acacia species is determined from the heterogeneous basin, the Upper Ewaso Ng’iro river basin, of the central Kenya highlands. The model results indicate that overall fitness is more strongly influenced by water use than stress avoidance but that consideration of both stress avoidance and water use is critical to predicting basin-scale patterns of species distribution. We identify a linear trend in the frequency and intensity of storms with the same annual total using a basin-wide gauge precipitation dataset. After calibration, we apply the basin average linear trends in time for average rain per storm and storm arrival rates. The model results indicate the upslope migration of two species, Acacia tortilis and Acacia xanthophloea to areas with higher total rainfall. Lastly, we explore the modeled changes of species cover in the basin influenced by changes in rainfall total holding growing season rainfall variability constant and changes in growing season rainfall variability holding total rainfall constant. We find that changes in dryland species distribution patterns and relative abundance may be as sensitive to growing season rainfall variability as they are to changes in total rainfall amounts.     

Velocity and Temperature Dissimilarity in the Surface Layer Uncovered by the Telegraph Approximation

We consider the Janjic (NCEP Office Note 437:61, 2001) boundary-layer model, which is one of the most widely used in numerical weather prediction models. This boundary-layer model is based on a number of length scales that are, in turn, obtained from a master length multiplied by constants. We analyze the simulation results obtained using different sets of constants with respect to measurements using sonic anemometers, and interpret these results in terms of the turbulence processes in the atmosphere and of the role played by the different length scales. The simulations are run on a virtual machine on the Chameleon cloud for low-wind-speed, unstable, and stable conditions.

     

Paleowind velocity and paleocurrents of pluvial Lake Manly,Death Valley,USA

Pluvial lake deposits are found throughout western North America and are frequently used to reconstruct regional paleoclimate. In Death Valley, California, USA, we apply the beach particle technique (BPT) of Adams (2003), Sedimentology, 50, 565–577 and Adams (2004), Sedimentology, 51, 671–673 to Lake Manly deposits at the Beatty Junction Bar Complex (BJBC), Desolation Canyon, and Manly Terraces and calculate paleowind velocities of 14–27 m/s. These wind velocities are within the range of present-day wind velocities recorded in the surrounding area. Sedimentary structures and clast provenance at Desolation Canyon and the Manly Terraces indicate sediment transport from north to south. Lake level, based on the elevation of constructional features, indicates that the hill west of the BJBC was an island and that the BJBC spits formed during simple lake regression. The data are consistent with the hypothesis that the present wind regime (velocity and direction) formed the pluvial Lake Manly features.     

Understanding groundwater processes by representing aquifer heterogeneity in the Maules Creek Catchment, Namoi Valley (New South Wales, Australia)

A FEFLOW three-dimensional (3D) groundwater model is developed to enhance the understanding of groundwater processes in the complex alluvial stratigraphy of Maules Creek Catchment (New South Wales, Australia). The aquifer vertical heterogeneity is replicated by indexing 204 lithological logs into units of high or low hydraulic conductivity, and by developing a 3D geological conceptual model with a vertical resolution based on the average lithological unit thickness for the region. The model mesh is populated with the indexed geology using nearest neighbour gridding. The calibrated model is successful in simulating the observed flow dynamics and in quantifying the important water-budget components. This indicates that the lateral groundwater flow from the mountainous region is the main inflow component of the system. Under natural conditions, the Namoi River acts as a sink of water, but groundwater abstraction increasingly removes a large amount of water each year causing dewatering of the system. The pumping condition affects the river–aquifer interaction by reversing the flow, from gaining to losing river conditions during the simulation period. The procedure is relevant for the development of groundwater models of heterogeneous systems in order to improve the understanding of the interplay between aquifer architecture and groundwater processes.     

Observations of Thin Layers in Coastal Hawaiian Waters

Thin layers of plankton have been documented in a wide variety of environments. The growing body of observations indicates that these features are a critical component of marine ecosystem dynamics and functioning. In the past two decades, much of the research on thin layers was undertaken in temperate coastal waters. Here, we report the first known observations of thin layers of phytoplankton in tropical Hawaiian waters. We conducted an overnight shipboard study during which time we made high-resolution observations of physical and optical structure in the water column. During the overnight cruise, we observed the greatest number of thin layers in the early evening hours when thermal stratification was strongest and most persistent due to a combination of warm air and surface water, as well as light winds. A comparison of these observations with those from temperate regions leads us to hypothesize that the nature and persistence of the physical structure is very important in determining the persistence of thin layered structures. Because plankton biomass is generally lower in tropical regions, the heterogeneous aggregation of food in thin subsurface layers may be more critical to the marine ecosystem than it is in temperate regions where plankton are generally more abundant.     

Retrospective and prospective model simulations of sea level rise impacts on Gulf of Mexico coastal marshes and forests in Waccasassa Bay, Florida

The State of Florida (USA) is especially threatened by sea level rise due to extensive low elevation coastal habitats (approximately 8,000?km²?<?1?m above sea level) where the majority of the human population resides. We used the Sea Level Affecting Marshes Model (SLAMM) simulation to improve understanding of the magnitude and location of these changes for 58,000?ha of the Waccasassa Bay region of Florida??s central Gulf of Mexico coast. To assess how well SLAMM portrays changes in coastal wetland systems resulting from sea level rise, we conducted a hindcast in which we compared model results to 30?years of field plot data. Overall, the model showed the same pattern of coastal forest loss as observed. Prospective runs of SLAMM using 0.64?m, 1?m and 2?m sea level rise scenarios predict substantial changes over this century in the area covered by coastal wetland systems including net losses of coastal forests (69%, 83%, and 99%, respectively) and inland forests (33%, 50%, and 88%), but net gains of tidal flats (17%, 142%, and 3,837%). One implication of these findings at the site level is that undeveloped, unprotected lands inland from the coastal forest should be protected to accommodate upslope migration of this natural community in response to rising seas. At a broader scale, our results suggest that coastal wetland systems will be unevenly affected across the Gulf of Mexico as sea level rises. Species vulnerable to these anticipated changes will experience a net loss or even elimination.     

Exploring weathering and regolith transport controls on Critical Zone development with models and natural experiments

The architecture of the Critical Zone, including mobile regolith thickness and depth to the weathering front, is first order controlled by advance of a weathering front at depth and transport of sediment at the surface. Differences in conditions imposed by slope aspect in the Gordon Gulch catchment of the Boulder Creek Critical Zone Observatory present a natural experiment to explore these interactions. The weathering front is deeper and saprolite more decayed on north-facing than on south-facing slopes. Simple numerical models of weathering front advance, mobile regolith production, and regolith transport are used to test how weathering and erosion rates interact in the evolution of weathered profiles. As the processes which attempt are being made to mimic are directly tied to climate variables such as mean annual temperature, the role of Quaternary climate variation in governing the evolution of Critical Zone architecture can be explored with greater confidence.