A fast and global two point low storage optimization technique for tracing rays in 2D and 3D isotropic media

We present the problem of tracing rays in 2D and 3D heterogeneous isotropic media as a set of optimization problems. Each optimization problem is obtained by applying Fermat's principle to an approximation of the travel time equation from a fixed source to a fixed receiver. We assume a piecewise linear ray path that simplifies the computations of the problem, in the same way Mao and Stuart suggested in a very recent paper. Here, instead, the reflector geometry and the velocity function are computed by using nonuniformly biharmonic splines. On the other hand, to solve the optimization problem we use the Global Spectral Gradient method. This recent developed optimization scheme is a low storage optimization technique that requires very few floating point operations. It only requires the gradient of the travel time function, and it is global because it converges independently of the initial guess, that is, it does not require a close initial ray path. These three properties of the optimization method and the assumption of piecewise linear rays make this ray tracing scheme a very fast, global and effective method when estimating velocities via tomography. Moreover, in a homogeneous stratified or dipped media, any solution of the optimization problem is the best solution, i.e., it is the global minimum, no matter what numerical approach is used. We present some numerical results that show the computational advantages and the performance of this ray tracing in homogeneous and heterogeneous media.     

Propagating uncertainty from catchment experiments to estimates of streamflow reduction by invasive alien plants in southwestern South Africa

Long-term catchment experiments from South Africa have demonstrated that afforestation of grasslands and shrublands significantly reduces surface-water runoff. These results have guided the country's forestry policy and the implementation of a national Invasive Alien Plant (IAP) control programme for the past few decades. Unfortunately, woody IAP densities continue to increase, compounding existing threats to water security from population growth and climatic change. Decision makers need defensible estimates of the impacts of afforestation or invasions on runoff to weigh up alternative land use options, or guide investment of limited resources into ecosystem restoration through IAP clearing versus engineering-based water-augmentation schemes. Existing attempts to extrapolate the impacts observed in catchment afforestation experiments to broad-scale IAP impacts give no indication of uncertainty. Globally, the uncertainty inherent in the results from paired-catchment experiments is seldom propagated into subsequent analyses making use of these data. We present a fully reproducible Bayesian model that propagates uncertainty from input data to final estimates of changes in streamflow when extrapolating from catchment experiments to broader landscapes. We apply our model to South Africa's catchment experiment data, estimating streamflow losses to plantations and analogous plant invasions in the catchments of southwestern South Africa, including uncertainty. We estimate that regional streamflow is reduced by 304 million m³ or 4.14% annually as a result of IAPs, with an upper estimate of 408 million m³ (5.54%) and a lower estimate of 267 million m³ (3.63%). Our model quantifies uncertainty associated with all parameters and their contribution to overall uncertainty, helping guide future research needs. Acknowledging and quantifying inherent uncertainty enables more defensible decisions regarding water resource management.     

Investigating the micro mechanics of cemented sand using DEM

The discrete element method has been used to investigate the micro mechanics of cemented sand. High‐pressure drained triaxial tests are modelled in 3D using a flexible membrane that allows the correct deformation to develop. Simulations with up to 12 MPa confining pressure are presented, which are compared with laboratory experiments on a sand with a range of cement contents. Cementation is modelled using ‘parallel bonds’, and various parameters and strength distributions are investigated. Varying levels of cementation are successfully modelled, with the correct qualitative behaviour observed and the separate effects of cementation and confining pressures demonstrated. The triaxial behaviour is found to be highly influenced by the distribution of bond strengths. Copyright © 2014 John Wiley & Sons, Ltd.     

Airborne LiDAR and Terrestrial Laser Scanning Derived Vegetation Obstruction Factors for Visibility Models

Research presented here explores the feasibility of leveraging vegetation data derived from airborne light detection and ranging (LiDAR) and terrestrial laser scanning (TLS) for visibility modeling. Using LiDAR and TLS datasets of a lodgepole pine (Pinus contorta) dominant ecosystem, tree canopy and trunk obstructions were isolated relevant to a discrete visibility beam in a short‐range line‐of‐sight model. Cumulative obstruction factors from vegetation were compared with reference visibility values from digital photographs along sightline paths. LiDAR‐derived tree factors were augmented with single‐scan TLS data for obstruction prediction. Good correlation between datasets was found up to 10 m from the terrestrial scanner, but fine scale visibility modeling was problematic at longer distances. Analysis of correlation and regression results reveal the influence of obstruction shadowing inherent to discrete LiDAR and TLS, potentially limiting the feasibility of modeling visibility over large areas with similar technology. However, the results support the potential for TLS‐derived subcanopy metrics for augmenting large amounts of aerial LiDAR data to significantly improve models of forest structure. Subtle LiDAR processing improvements, including more accurate tree delineation through higher point density aerial data, combined with better vegetation quantification processes for TLS data, will advance the feasibility and accuracy of data integration.     

Assessment of carbon in woody plants and soil across a vineyard-woodland landscape

Background

Quantification of ecosystem services, such as carbon (C) storage, can demonstrate the benefits of managing for both production and habitat conservation in agricultural landscapes. In this study, we evaluated C stocks and woody plant diversity across vineyard blocks and adjoining woodland ecosystems (wildlands) for an organic vineyard in northern California. Carbon was measured in soil from 44 one m deep pits, and in aboveground woody biomass from 93 vegetation plots. These data were combined with physical landscape variables to model C stocks using a geographic information system and multivariate linear regression.

Results

Field data showed wildlands to be heterogeneous in both C stocks and woody tree diversity, reflecting the mosaic of several different vegetation types, and storing on average 36.8 Mg C/ha in aboveground woody biomass and 89.3 Mg C/ha in soil. Not surprisingly, vineyard blocks showed less variation in above- and belowground C, with an average of 3.0 and 84.1 Mg C/ha, respectively.

Conclusions

This research demonstrates that vineyards managed with practices that conserve some fraction of adjoining wildlands yield benefits for increasing overall C stocks and species and habitat diversity in integrated agricultural landscapes. For such complex landscapes, high resolution spatial modeling is challenging and requires accurate characterization of the landscape by vegetation type, physical structure, sufficient sampling, and allometric equations that relate tree species to each landscape. Geographic information systems and remote sensing techniques are useful for integrating the above variables into an analysis platform to estimate C stocks in these working landscapes, thereby helping land managers qualify for greenhouse gas mitigation credits. Carbon policy in California, however, shows a lack of focus on C stocks compared to emissions, and on agriculture compared to other sectors. Correcting these policy shortcomings could create incentives for ecosystem service provision, including C storage, as well as encourage better farm stewardship and habitat conservation.

     

Surmounting luminescence age overestimation in Alaska-margin Arctic Ocean sediments by use of ‘micro-hole’ quartz dating

Because of several difficulties with the application of radiocarbon (¹⁴C) dating to Arctic Ocean sediments, numeric dating techniques are needed that can complement, supplant and reach beyond the ¹⁴C method. However, large age overestimates (often >7 kyr) for near-sediment-water-interface horizons from Arctic Ocean cores have been almost universal when luminescence dating has been applied to multigrain aliquots of fine silt (4–11 μm) quartz and feldspar grains. Here micro-hole quartz-grain photon-stimulated-luminescence (PSL) dating is applied to the 0.5–2.0 cm horizons of multicores from high-sedimentation-rate sites spanning depths from 87 m to 1140 m at the Alaska margin of the Arctic Ocean. Expected near zero ages (0–200 a) result when grains larger than ～11 μm are used, demonstrating that fine-silt age overestimations here and perhaps elsewhere in the Arctic Ocean are a function of grain and aliquot size. At the 87 m site, the micro-hole PSL approach revealed no significant gradient in age estimates over the 1–26 cm horizon range, implying that bioturbation reached to at least 26 cm. Micro-hole PSL dating of 25–62 μm quartz grains from trans-ocean sea-ice sediment also produced expected near zero ages, in contrast to earlier reported long-bleach multigrain PSL results from 4–11 μm fractions of the same samples. The micro-hole PSL approach thus surmounts the age overestimation problem associated with the use of multigrain silt fractions, and overcomes limitations of the ¹⁴C method in this region. Finally, results unexpectedly suggest the potential of micro-hole quartz PSL for use in provenance studies of Arctic Ocean sea-ice sediment.     

The effect of ammonia ice on the outgoing thermal radiance from the atmosphere of Jupiter

We examine the effects of NH₃ ice particle clouds in the atmosphere of Jupiter on outgoing thermal radiances. The cloud models are characterized by a number density at the cloud base, by the ratio of the scale height of the vertical distribution of particles (H_p) to the gas scale height (H_g), and by an effective particle radius. NH₃ ice particle-scattering properties are scaled from laboratory measurements. The number density for the various particle radius and scale height models is inferred from the observed disk average radiance at 246 cm^?1, and preliminary lower limits on particle sizes are inferred from the lack of apparent NH₃ absorption features in the observed spectral radiances as well as the observed minimum flux near 2100 cm^?1. We find lower limits on the particle size of 3 μm if H_p/H_g = 0.15, or 10μmif H_p/H_g = 0.50 or 0.05. NH₃ ice particles are relatively dark near the far-infrared and 8.5-μm atmospheric windows, and the outgoing thermal radiances are not very sensitive to various assumptions about the particle-scattering function as opposed to radiances at 5 μm, where particles are relatively brighter. We examined observations in these three different spectral window regions which provide, in principle, complementary constraints on cloud parameters. Characterization of the cloud scale height is difficult, but a promising approach is the examination of radiances and their center-to-limb variation in spectral regions where there is significant opacity provided by gases of known vertical distribution. A blackbody cloud top model can reduce systematic errors due to clouds in temperature sounding to the level of 1K or less. The NH₃ clouds provide a substantial influence on the internal infrared flux field near the 600-mbar level.     

Oxygen isotope and 26Al‐26Mg systematics of aluminum‐rich chondrules from unequilibrated enstatite chondrites

Abstract— Correlated in situ analyses of the oxygen and magnesium isotopic compositions of aluminum‐rich chondrules from unequilibrated enstatite chondrites were obtained using an ion microprobe. Among eleven aluminum‐rich chondrules and two plagioclase fragments measured for ²⁶Al‐²⁶Mg systematics, only one aluminum‐rich chondrule contains excess ²⁶Mg from the in situ decay of ²⁶Al; the inferred initial ratio (²⁶Al/²⁷Al)_o = (6.8 ± 2.4) × 10^?6 is consistent with ratios observed in chondrules from carbonaceous chondrites and unequilibrated ordinary chondrites. The oxygen isotopic compositions of five aluminum‐rich chondrules and one plagioclase fragment define a line of slope ?0.6 ± 0.1 on a three‐oxygen‐isotope diagram, overlapping the field defined by ferromagnesian chondrules in enstatite chondrites but extending to more ¹⁶O‐rich compositions with a range in δ¹⁸O of about ?12‰. Based on their oxygen isotopic compositions, aluminum‐rich chondrules in unequilibrated enstatite chondrites are probably genetically related to ferromagnesian chondrules and are not simple mixtures of materials from ferromagnesian chondrules and calcium‐aluminum‐rich inclusions (CAIs). Relative to their counterparts from unequilibrated ordinary chondrites, aluminum‐rich chondrules from unequilibrated enstatite chondrites show a narrower oxygen isotopic range and much less resolvable excess ²⁶Mg from the in situ decay of ²⁶Al, probably resulting from higher degrees of equilibration and isotopic exchange during post‐crystallization metamorphism. However, the presence of ²⁶Al‐bearing chondrules within the primitive ordinary, carbonaceous, and now enstatite chondrites suggests that ²⁶Al was at least approximately homogeneously distributed across the chondrite‐forming region.