Generating MODFLOW grids from boundary representation solid models

Complex stratigraphy can be difficult to simulate in MODFLOW models. MODFLOW uses a structured grid that requires that each grid layer be continuous throughout the model domain. This makes it difficult to explicitly represent common features such as pinchouts and embedded seams in a MODFLOW model. In this paper, we describe a method for automatically generating MODFLOW-compatible grids from boundary-representation solid models. Solid models are data structures developed originally for computer-aided design applications that define the geometry of three-dimensional objects. Solid models can be used to represent arbitrarily complex stratigraphy. The elevations defined by the solids are then extracted from the solids in a manner that preserves the continuous-layer requirement imposed by MODFLOW. Two basic approaches are described: The first method adjusts the MODFLOW grid dimensions (layer elevations) to fit the solid model boundaries, and the second method creates a regular MODFLOW grid and adjusts the material properties to match the changes in stratigraphy. One of the main benefits of using solid models to define stratigraphy for MODFLOW models is that it provides a grid-independent definition of the layer elevations that can be used to immediately re-create the MODFLOW grid geometry after any change to the grid resolution.     

Effect of a dolerite intrusion on triterpane stereochemistry and kerogen in Rundle oil shale,Australia

Stereochemical changes of triterpanes present in extracts from an immature oil shale sequence intruded by a 3-m dolerite sill have been studied by gas chromatography-mass spectrometry (GC-MS). The steric configuration of the hopanes was observed to change from one dominated by the thermally less stable 17β(H), 21β(H) configuration at some distance from the intrusion, to one dominated by the thermally more stable 17α(H), 21β(H) and 17β(H), 21α(H) configurating in the immediate vicinity of the intrusion. In addition, severe alteration of the kerogen appeared to have taken place as a result of the contact metamorphism, and high concentrations of extractable organic matter were observed below the intrusion. Characterization of the kerogens by Curie-point pyrolysis has enabled the effects of the intrusion on the shales to be monitored.     

Lunar occultation of Saturn: IV. Astrometric results from observations of the satellites

The method of determining local lunar limb slopes, and the consequent time scale needed for diameter studies, from accurate occultation timings at two nearby telescope is described. The results for the photoelectric observations made at Mauna Kea Observatory during the occultation of Saturn's satellites on March 30, 1974, are discussed. Analysis of all observations of occultations of Saturn's satellites during 1974 indicates possible errors in the ephemerides of Saturn and its satellites.     

Near-opposition limb darkening of solids of planetary interest

This paper presents a laboratory study of the limb darkening near opposition, of particulate materials of planetary interest and concentrates on the wavelength dependence of this limb darkening. We find that near zero phase the scattering properties of most particulate materials can be described adequately by Minnaert's law. However, there are materials for which such a representation is totally inadequate. Examples are bronzite and graphite, materials that tend to fracture into flakes having mirrorlike surfaces. In addition, there are materials, such as olivine, whose scattering properties within deep absorption bands show definite departures from Minnaert's law at large angles of incidence or emission. Our Minnaert parameters, k and B₀, were measured at a phase angle of α = 4°. For samples of comparable surface texture and roughness, k and B₀ are approximately linearly related, k increasing as B₀ increases. Very dark materials tend to have k ～ 0.5 to 0.6, while very bright materials tend to have k ～ 1. The linear relationship between k and B₀ can be explained in terms of the varying importance of multiple scattering in the surface layer. Thus for materials for which B₀ is strongly wavelength dependent, so is k. For example, for olivine, k varies from 0.73 to 0.87 between 0.4 and 1.2 μm. These variations are closely correlated with those in B₀: the value of k is relatively high outside of absorption bands and relatively low within them. For bright materials, k is very sensitive to surface roughness. For example, for BaSO₄ powder, k can be changed from ～1.0 to ～0.8 by this effect alone, a fact which has relevance to the photometry of frost-covered satellites. For dark materials, the effects of surface roughness on k are smaller and more subtle.     

The causes of bedload pulses in a gravel channel: the implications of bedload grain‐size distributions

Unsteady bedload transport was measured in two c. 5 m wide anabranches of a gravel‐bed braided stream draining the Haut Glacier d'Arolla, Switzerland, during the 1998 and 1999 melt seasons. Bedload was directly sampled using 152 mm square Helley–Smith type samplers deployed from a portable measuring bridge, and independent transport rate estimates for the coarser size fractions were obtained from the dispersion of magnetically tagged tracer pebbles. Bedload transport time series show pulsing behaviour under both marginal (1998) and partial (1999) transport regimes. There are generally weak correlations between transport rates and shear stresses determined from velocity data recorded at the measuring bridge. Characteristic parameters of the bedload grain‐size distributions (D₅₀, D₈₄) are weakly correlated with transport rates. Analysis of full bedload grain‐size distributions reveals greater structure, with a tendency for transport to become less size selective at higher transport rates. The bedload time series show autoregressive behaviour but are dif?cult to distinguish by this method. State–space plots, and associated measures of time‐series separation, reveal the structure of the time series more clearly. The measured pulses have distinctly different time‐series characteristics from those modelled using a one‐dimensional sediment routing model in which bed shear stress and grain size are varied randomly. These results suggest a mechanism of pulse generation based on irregular low‐amplitude bedforms, that may be generated in‐channel or may represent the advection of material supplied by bank erosion events. Copyright © 2003 John Wiley & Sons, Ltd.     

Factors Influencing Adoption and Implementation of Conservation Development Ordinances in Rural United States

Conservation Development (CD) is a land-use tool to reduce the impacts of development on natural resources, yetthe use of CD is highly variable. We used a collaborative conservation model to examine factors influencing adoption and implementation to improve the future application of CD. We interviewed 2–3 individuals including planners, board members, and developers (n?=?46) from 16 rural communities in the intermountain west and northeast U.S. Motivations to preserve rural character or comply with state statutes drove adoption. Greater capacity such as outside expertise, model regulations, and iterative revisions of ordinances was associated with stronger ordinances. Adoption processes including ample dialog across diverse constituents, especially those overseeing subdivision development, facilitated higher rates of implementation. Concerns over open space management, lack of successful CD examples, and inexperience with CD posed barriers to implementation. Understanding social context revealed the strengths and shortcomings of CD and provided guidance for strategically engaging communities in private lands conservation.     

A numerical assessment of simple airblast models of impact airbursts

Asteroids and comets 10–100 m in size that collide with Earth disrupt dramatically in the atmosphere with an explosive transfer of energy, caused by extreme air drag. Such airbursts produce a strong blastwave that radiates from the meteoroid's trajectory and can cause damage on the surface. An established technique for predicting airburst blastwave damage is to treat the airburst as a static source of energy and to extrapolate empirical results of nuclear explosion tests using an energy‐based scaling approach. Here we compare this approach to two more complex models using the iSALE shock physics code. We consider a moving‐source airburst model where the meteoroid's energy is partitioned as two‐thirds internal energy and one‐third kinetic energy at the burst altitude, and a model in which energy is deposited into the atmosphere along the meteoroid's trajectory based on the pancake model of meteoroid disruption. To justify use of the pancake model, we show that it provides a good fit to the inferred energy release of the 2013 Chelyabinsk fireball. Predicted overpressures from all three models are broadly consistent at radial distances from ground zero that exceed three times the burst height. At smaller radial distances, the moving‐source model predicts overpressures two times greater than the static‐source model, whereas the cylindrical line‐source model based on the pancake model predicts overpressures two times lower than the static‐source model. Given other uncertainties associated with airblast damage predictions, the static‐source approach provides an adequate approximation of the azimuthally averaged airblast for probabilistic hazard assessment.     

Assessing the significance of wet-canopy evaporation from forests during extreme rainfall events for flood mitigation in mountainous regions of the United Kingdom

There is increased interest in the potential of tree planting to help mitigate flooding using nature-based solutions or natural flood management. However, many publications based upon catchment studies conclude that, as flood magnitude increases, benefit from forest cover declines and is insignificant for extreme flood events. These conclusions conflict with estimates of evaporation loss from forest plot observations of gross rainfall, through fall and stem flow. This study explores data from existing studies to assess the magnitudes of evaporation and attempts to identify the meteorological conditions under which they would be supported. This is achieved using rainfall event data collated from publications and data archives from studies undertaken in temperate environments around the world. The meteorological conditions required to drive the observed evaporation losses are explored theoretically using the Penman–Monteith equation. The results of this theoretical analysis are compared with the prevailing meteorological conditions during large and extreme rainfall events in mountainous regions of the United Kingdom to assess the likely significance of wet canopy evaporation loss. The collated dataset showed that event Ewc losses between approximately 2 and 38% of gross rainfall (1.5 to 39.4 mm day⁻¹) have been observed during large rainfall events (up to 118 mm day⁻¹) and that there are few data for extreme events (>150 mm day⁻¹). Event data greater than 150 mm (reported separately) included similarly high percentage evaporation losses. Theoretical estimates of wet-canopy evaporation indicated that, to reproduce the losses towards the high end of these observations, relative humidity and the aerodynamic resistance for vapour transport needed to be lower than approximately 97.5% and 0.5 to 2 s m⁻¹ respectively. Surface meteorological data during large and extreme rainfall events in the United Kingdom suggest that conditions favourable for high wet-canopy evaporation are not uncommon and indicate that significant evaporation losses during large and extreme events are possible but not for all events and not at all locations. Thus the disparity with the results from catchment studies remains.