Calibrating a Salt Water Intrusion Model with Time‐Domain Electromagnetic Data

Salt water intrusion models are commonly used to support groundwater resource management in coastal aquifers. Concentration data used for model calibration are often sparse and limited in spatial extent. With airborne and ground‐based electromagnetic surveys, electrical resistivity models can be obtained to provide high‐resolution three‐dimensional models of subsurface resistivity variations that can be related to geology and salt concentrations on a regional scale. Several previous studies have calibrated salt water intrusion models with geophysical data, but are typically limited to the use of the inverted electrical resistivity models without considering the measured geophysical data directly. This induces a number of errors related to inconsistent scales between the geophysical and hydrologic models and the applied regularization constraints in the geophysical inversion. To overcome these errors, we perform a coupled hydrogeophysical inversion (CHI) in which we use a salt water intrusion model to interpret the geophysical data and guide the geophysical inversion. We refer to this methodology as a Coupled Hydrogeophysical Inversion‐State (CHI‐S), in which simulated salt concentrations are transformed to an electrical resistivity model, after which a geophysical forward response is calculated and compared with the measured geophysical data. This approach was applied for a field site in Santa Cruz County, California, where a time‐domain electromagnetic (TDEM) dataset was collected. For this location, a simple two‐dimensional cross‐sectional salt water intrusion model was developed, for which we estimated five uniform aquifer properties, incorporating the porosity that was also part of the employed petrophysical relationship. In addition, one geophysical parameter was estimated. The six parameters could be resolved well by fitting more than 300 apparent resistivities that were comprised by the TDEM dataset. Except for three sounding locations, all the TDEM data could be fitted close to a root‐mean‐square error of 1. Possible explanations for the poor fit of these soundings are the assumption of spatial uniformity, fixed boundary conditions and the neglecting of 3D effects in the groundwater model and the TDEM forward responses.     

A comparative study of the thermal behaviour of length-fast chalcedony, length-slow chalcedony (quartzine) and moganite

The thermal behaviour of silica rocks upon heat treatment is dependent on the constituent minerals and petrographic texture types. These constituents can be shown to be mainly quartz in the form of two types of chalcedony (Length-fast (LF) chalcedony and Length-slow (LS) chalcedony, the latter also being termed quartzine) and moganite. Even though the thermal behaviour of LF-chalcedony is well understood, major uncertainties persist concerning the high-temperature behaviour of LS-chalcedony and moganite. We present here a comparative study of these three constituents of common silica rocks. Our results show that the chemical reaction is the same in all three, Si–OH + HO–Si → Si–O–Si + H₂O, but that the reaction kinetics and activation temperatures are very different. LS-chalcedony begins to react from 200 °C upwards, that is at temperatures 50 °C below the ones observed in LF-chalcedony, and shows the fastest reaction kinetics of this ‘water’ loss. Chemically bound water (SiOH) in moganite is more stable at high temperatures and no specific activation temperature is necessary for triggering the temperature-induced ‘water’ loss. Moganite is also found to act as a stabilizer in silica rocks preventing them from temperature-induced fracturing. These findings have implications for the study of potential heat treatment temperatures of silica rocks (in industry and heritage studies), but they also shed light on the different structures of SiO₂ minerals and the role of OH impurities therein.     

Hydro-micromechanical modeling of wave propagation in saturated granular crystals

Biot theory predicts wave velocities in a saturated granular medium using the pore geometry, viscosity, densities, and elastic moduli of the solid skeleton and pore fluid, neglecting the interaction between constituent particles and local flow, which becomes essential as the wavelength decreases. Here, a hydro-micromechanical model, for direct numerical simulations of wave propagation in saturated granular media, is implemented by two-way coupling the lattice Boltzmann method (LBM) and the discrete element method (DEM), which resolve the pore-scale hydrodynamics and intergranular behavior, respectively. The coupling scheme is benchmarked with the terminal velocity of a single sphere settling in a fluid. In order to mimic a small amplitude pressure wave entering a saturated granular medium, an oscillating pressure boundary on the fluid is implemented and benchmarked with the one-dimensional wave equation. The effects of input waveforms and frequencies on the dispersion relations in 3D saturated poroelastic media are investigated with granular face-centered-cubic crystals. Finally, the pressure and shear wave velocities predicted by the numerical model at various effective confining pressures are found to be in excellent agreement with Biot analytical solutions, including his prediction for slow compressional waves.     

Dendrochronology and lichenometry: colonization, growth rates and dating of geomorphological events on the east side of the North Patagonian Icefield, Chile

This paper highlights the importance for dating accuracy of initial studies of delay before colonization for both trees and lichens and tree age below core height, particularly in recently deglaciated terrain where colonization and growth rates may vary widely due to differences in micro-environment. It demonstrates, for the first time, how dendrochronology and lichenometry can be used together in an assessment of each other's colonization and growth rates, and then cross-correlated to provide a supportive dating framework. The method described for estimating tree age below core height is also new. The results show that on the east side of the North Patagonian Icefield in the Arco and Colonia valleys, Nothofagus age below a core height of 112 cm can vary from 5 to 41 years and delay before colonization may range from a maximum of 22 years near water to a minimum of 93 years on the exposed flanks of the Arenales and Colonia Glaciers. Tree age plus colonization delay supplied a maximum growth rate of 4.7 mm/year for the lichen Placopsis perrugosa and lichen colonization is estimated to take from 2.5 to approximately 13 years. A minimum lichenometric date of 1883 was estimated for an ice-formed trimline at the junction of the Arenales and Colonia glaciers and a maximum dendrochronological date of 1881 for a water-formed trimline in the Arco valley. Tree and lichen ages around the valley suggest that a glacial outburst drained the 1881 high level lake releasing approximately 265 million cubic metres of water. Repeated flooding, with a minimum of 38 high lake levels, is suggested by horizontal sediment lines on the Arco valley walls and moraine flanks. Dating confirmed diminishing flood levels with a last minor flood in 1963. The wider significance of the work is that it should produce more accurate dating of recent glacier fluctuations around the North Patagonia Icefield, an area where dated reference surfaces are extremely scarce.     

Understanding how environmental heterogeneity and elevation drives the distribution of woody communities across vegetation types within the campo rupestre in South America

The campo rupestre sensu lato is a vegetation type that occurs in South American mountains, supports a distinctive flora characterized by high rates of endemism, high herbaceous species richness and often-neglected but also species-rich of the arboreal stratum. We aimed to investigate how environmental factors and elevation are associated with the distribution and diversity of woody species in different rupestrian vegetation types across South America. Using a database of 2,049 woody species from 185 sites across four vegetation types within the campo rupestre, we assessed how the vegetation types were grouped according to their floristic composition and number of shared indicator species, as well as by using different beta diversity indices. The most important variables from a set of 27 variables(e.g. altitude, geo-edaphic and climatic) explaining species distribution were identified using redundancy analysis(RDA) and variation partitioning methods. The distribution of vegetation types was related to both environmental and spatial fractions, with a set of 17 variables retained(e.g. rockiness, grass cover and temperature seasonality as the most important variables). There was an association between the floristic composition of each vegetation type and the elevation range. Although the identified vegetation types are floristically related, they are distinguished by exclusive and habitat-specialist woody species. This uniqueness of vegetation types should be considered in terms of complementarity for the conservation of campos rupestres.     

The Relationship between Geomorphic River Structure and Coarse Particulate Organic Matter (CPOM) Storage along the Kangaroo River,New South Wales,Australia

Fluvial landforms provide a physical template upon which to appraise biophysical relationships along river courses. In this study, the spatial pattern of organic matter storage along the Kangaroo River, NSW, is related to geomorphic controls that operate at a range of scales within a nested hierarchy. This snapshot study of CPOM storage found that at the catchment scale the longitudinal pattern of coarse particulate organic matter (CPOM) storage is dependent on the type and downstream pattern of River Styles. At the reach scale, CPOM storage is dependent on the geomorphic unit structure and physical heterogeneity of the river, and associated energy conditions along the reach. At the geomorphic unit scale, CPOM storage capacity is related to the position of geomorphic units relative to the thalweg (i.e. flow characteristics) and associated roughness attributes. At the hydraulic unit scale, CPOM storage capacity is related to local flow velocity and substrate characteristics (clast size and distribution).     

Double-difference Relocation of the Aftershocks of the Tecom��n, Colima, Mexico Earthquake of 22 January 2003

On 22 January 2003, the M _w?=?7.6 Tecomán earthquake struck offshore of the state of Colima, Mexico, near the diffuse triple junction between the Cocos, Rivera, and North American plates. Three-hundred and fifty aftershocks of the Tecomán earthquake with magnitudes between 2.6 and 5.8, each recorded by at least 7 stations, are relocated using the double difference method. Initial locations are determined using P and S readings from the Red Sismológica Telemétrica del Estado de Colima (RESCO) and a 1-D velocity model. Because only eight RESCO stations were operating immediately following the Tecomán earthquake, uncertainties in the initial locations and depths are fairly large, with average uncertainties of 8.0?km in depth and 1.4?km in the north?Csouth and east?Cwest directions. Events occurring between 24 January and 31 January were located using not only RESCO phase readings but also additional P and S readings from 11 temporary stations. Average uncertainties decrease to 0.8?km in depth, 0.3?km in the east?Cwest direction, and 0.7?km in the north?Csouth direction for events occurring while the temporary stations were deployed. While some preliminary studies of the early aftershocks suggested that they were dominated by shallow events above the plate interface, our results place the majority of aftershocks along the plate interface, for a slab dipping between approximately 20° and 30°. This is consistent with the slab positions inferred from geodetic studies. We do see some upper plate aftershocks that may correspond to forearc fault zones, and faults inland in the upper plate, particularly among events occurring more than 3?months after the mainshock.     

Biogeochemical Behavior of Arsenic Species at Paranaguá Estuarine Complex, Southern Brazil

The behavior and dynamics of arsenic at Ilha do Mel, Guaraque?aba, and Paranaguá, located in the Paranaguá Estuarine Complex, Brazil, were investigated in this work. Samples were collected in March, September, and December of 2005. With respect to arsenic behavior, the total concentrations were higher for Paranaguá (22.5?±?2.5?μg?L^?1) on average, followed by the samples of Guaraque?aba (14.4?±?3.2?μg?L^?1) and Ilha do Mel (8.7?±?1.1?μg?L^?1). The concentrations found in this work were consistently greater when compared to other estuaries. The results can be attributed to geological factors. Nevertheless, it was also possible to realize that higher concentrations were found in Paranaguá, the more impacted environment. These results represent a strong indication that human activities in the area contribute to the enrichment of the estuary with this metalloid. Arsenic was mainly present in the dissolved phase, which corresponded to more than 70% of the total concentration. This behavior can represent a greater residence time of this metalloid in the water column. The marine phytoplankton presence was decisive in dictating the distribution of As in the waters of Paranaguá estuary. Organic species, which are considered the less toxic, were found in higher concentrations in the eutrophic areas such as Guaraque?aba and Paranaguá. This was a typical behavior during the summer probably due to higher biological activity.     

Storage of carbon in U.S. forests predicted from satellite data,ecosystem modeling,and inventory summaries

A plant and soil simulation model based on satellite observations of vegetation and climate data was used to estimate the potential carbon pools in standing wood biomass across all forest ecosystems of the conterminous United States up to the year 1997. These modeled estimates of vegetative carbon potential were compared to aggregated measurements of standing wood biomass from the U. S. Forest Service’s national Forest Inventory and Analysis (FIA) data set and the Carbon Online Estimator (COLE) to understand: 1) predominant geographic variations in tree growth rate and 2) local land cover and land use history including the time since the last stand-replacing disturbance (e.g., from wildfire or harvest). Results suggest that although wood appears to be accumulating at high rates in many areas of the U.S. (Northwest and Southeast), there are still extensive areas of relatively low biomass forest in the late 1990s according to FIA records. We attribute these low biomass accumulation levels to the high frequency of disturbances, which can be observed even in high production areas such as the Southeast due to frequent forest harvests. Ecosystem models like the one presented in this study have been coupled with satellite observations of land cover and green plant density to uniquely differentiate areas with a high potential for vegetative carbon storage at relatively fine spatial resolution.