Recognition of shear heating on a long‐lived major fault using Raman carbonaceous material thermometry: implications for strength and displacement history of the MTL,SW Japan

The extent to which movement on major faults causes long term shear heating is a contentious issue and an important aspect in the debate about the strength of major faults in the crust. Comparing the results of experimental work on the kinetics of crystallization of carbonaceous material with results of thermal modeling show that the Raman carbonaceous material (CM) geothermometer is well suited to studying shear heating on geological time scales in suitable lithologies exposed around exhumed major fault zones. The Median Tectonic Line (MTL), SW Japan, is the largest on‐land fault in Japan with a length of > 800 km. Application of Raman CM thermometry to pelitic schist adjacent to the fault reveals the presence of a rise in peak temperature of around 60 °C over a distance of around 150 m perpendicular to the MTL fault plane. The spatial association of this thermal anomaly with the fault implies it is due to shear heating. Thermal modeling shows the recorded thermal anomaly and steep temperature gradient is compatible with very high rates of displacement over time scales of a few thousand years. However, the implied displacement rates lie outside those generally observed. An alternative explanation is that an originally broader thermal anomaly that developed during strike slip faulting was shortened due to the effects of normal faulting. Constraints on displacement rate, width of the original anomaly, duration of heating and peak temperature imply a coefficient of friction, μ, greater than 0.4.     

A new approach to develop the Raman carbonaceous material geothermometer for low‐grade metamorphism using peak width

The Raman spectra of carbonaceous material (CM) from 19 metasediment samples collected from six widely separated areas of Southwest Japan and metamorphosed at temperatures from 165 to 655°C show systematic changes with metamorphic temperature that can be classified into four types: low‐grade CM (c. 150–280°C), medium‐grade CM (c. 280–400°C), high‐grade CM (c. 400–650°C), and well‐crystallized graphite (> c. 650°C). The Raman spectra of low‐grade CM exhibit features typical of amorphous carbon, in which several disordered bands (D‐band) appear in the first‐order region. In the Raman spectra of medium‐grade CM, the graphite band (G‐band) can be recognized and several abrupt changes occur in the trends for several band parameters. The observed changes indicate that CM starts to transform from amorphous carbon to crystallized graphite at around 280°C, and this transformation continues until 400°C. The G‐band becomes the most prominent peak at high‐grade CM suggesting that the CM structure is close to that of well‐crystallized graphite. In the highest temperature sample of 655°C, the Raman spectra of CM show a strong G‐band with almost no recognizable D‐band, implying the CM grain is well‐crystallized graphite. In the Raman spectra of low‐ to medium‐grade CM, comparisons of several band parameters with the known metamorphic temperature show inverse correlations between metamorphic temperature and the full width at half maximum (FWHM) of the D1‐ and D2‐bands. These correlations are calibrated as new Raman CM geothermometers, applicable in the range of c. 150–400°C. Details of the methodology for peak decomposition of Raman spectra from the low to medium temperature range are also discussed with the aim of establishing a robust and user‐friendly geothermometer.     

Hydrologic‐response‐driven sediment transport at a regional scale,process‐based simulation

Assessing hydrologically driven erosion at regional scales from a process‐based perspective presents a significant challenge. Most regional‐scale erosion assessments are based upon a simple steady‐state hydrology foundation. For this study, the sediment transport version of the physics‐based Integrated Hydrology Model (InHM), excited by synthetically generated rainfall, was employed to assess long‐term hydrologically driven erosion for a regional‐scale island boundary‐value problem. The spatiotemporal dynamics of runoff generation, erosion, and deposition are illustrated through saturation, water depth, velocity, and sediment concentration results. The simulations demonstrate that process‐based assessment for concept development is both feasible and tractable at regional spatial and human time scales. Copyright © 2011 John Wiley & Sons, Ltd.     

Coupled numerical–analytical approach to landscape evolution modeling

The Earth's topography is shaped by surface processes that operate on various scales. In particular, river processes control landscape dynamics over large length scales, whereas hillslope processes control the dynamics over smaller length scales. This scale separation challenges numerical treatments of landscape evolution that use space discretization. Large grid spacing cannot account for the dynamics of water divides that control drainage area competition, and erosion rate and slope distribution. Small grid spacing that properly accounts for divide dynamics is computationally inefficient when studying large domains. Here we propose a new approach for landscape evolution modeling that couples irregular grid‐based numerical solutions for the large‐scale fluvial dynamics and continuum‐based analytical solutions for the small‐scale fluvial and hillslope dynamics. The new approach is implemented in the landscape evolution model DAC (divide and capture). The geometrical and topological characteristics of DAC's landscapes show compatibility with those of natural landscapes. A comparative study shows that, even with large grid spacing, DAC predictions fit well an analytical solution for divide migration in the presence of horizontal advection of topography. In addition, DAC is used to study some outstanding problems in landscape evolution. (i) The time to steady‐state is investigated and simulations show that steady‐state requires much more time to achieve than predicted by fixed area calculations, due to divides migration and persistent reorganization of low‐order streams. (ii) Large‐scale stream captures in a strike‐slip environment are studied and show a distinct pattern of erosion rates that can be used to identify recent capture events. (iii) Three tectono‐climatic mechanisms that can lead to asymmetric mountains are studied. Each of the mechanisms produces a distinct morphology and erosion rate distribution. Application to the Southern Alps of New Zealand suggests that tectonic advection, precipitation gradients and non‐uniform tectonic uplift act together to shape the first‐order topography of this mountain range. Copyright © 2014 John Wiley & Sons, Ltd.     

Air‐stream temperature correlation in forested and urban headwater streams in the Southern Appalachians

Air temperature can be an effective predictor of stream temperature. However, little work has been done in studying urban impacts on air‐stream relationships in groundwater‐fed headwater streams in mountainous watersheds. We applied wavelet coherence analysis to two 13‐month continuous (1 hr interval) stream and air temperature datasets collected at Boone Creek, an urban stream, and Winkler Creek, a forest stream, in northwestern North Carolina. The main advantage of a wavelet coherence analysis approach is the ability to analyse non‐stationary dynamics for the temporal covariance between air and stream temperature over time and at multiple temporal scales (e.g. hours, days, weeks and months). The coherence is both time and scale‐dependent. Our research indicated that air temperature generally co‐varied with stream temperature at time scales greater than 0.5 day. The correlation was poor during the winter at the scales of 1 to 64 days and summer at the scales of 1.5 to 4 days, respectively. The empirical models that relate air temperature to stream temperature failed at these scales and during these periods. Finally, urbanization altered the air‐stream temperature correlation at intermediate time scales ranging from 2 to 12 days. The correlation at the urban creek increased at the 12‐day scale, whereas it decreased at scales of 2 to 7 days as compared with the forested stream, which is probably due to heated surface runoff during summer thunderstorms or leaking stormwater or wastewater collection systems. Our results provide insights into air‐stream temperature relationships over both time and scale domains. Identifying controls over time and scales are needed to predict water temperature to understand the future impacts that interacting climate and land use changes will have on aquatic ecosystem in river networks. Copyright © 2014 John Wiley & Sons, Ltd.     

Multifractal behaviour of long‐term karstic discharge fluctuations

Karstic watersheds are highly complex hydrogeological systems that are characterized by a multiscale behaviour corresponding to the different pathways of water in these systems. The main issue of karstic spring discharge fluctuations consists in the presence and the identification of characteristic time scales in the discharge time series. To identify and characterize these dynamics, we acquired, for many years at the outlet of two karstic watersheds in South of France, discharge data at 3‐mn, 30‐mn and daily sampling rate. These hydrological records constitute to our knowledge the longest uninterrupted discharge time series available at these sampling rates. The analysis of the hydrological records at different levels of detail leads to a natural scale analysis of these time series in a multifractal framework. From a universal class of multifractal models based on cascade multiplicative processes, the time series first highlights two cut‐off scales around 1 and 16 h that correspond to distinct responses of the aquifer drainage system. Then we provide estimates of the multifractal parameters α and C₁ and the moment of divergence q_D corresponding to the behaviour of karstic systems. These results constitute the first estimates of the multifractal characteristics of karstic spingflows based on 10 years of high‐resolution discharge time series and should lead to several improvements in rainfall‐karstic springflow simulation models. Copyright © 2012 John Wiley & Sons, Ltd.     

Integrating the LISFLOOD‐FP 2D hydrodynamic model with the CAESAR model: implications for modelling landscape evolution

Landscape evolution models (LEMs) simulate the geomorphic development of river basins over long time periods and large space scales (100s–1000s of years, 100s of km²). Due to these scales they have been developed with simple steady flow models that enable long time steps (e.g. years) to be modelled, but not shorter term hydrodynamic effects (e.g. the passage of a flood wave). Nonsteady flow models that incorporate these hydrodynamic effects typically require far shorter time steps (seconds or less) and use more expensive numerical solutions hindering their inclusion in LEMs. The recently developed LISFLOOD‐FP simplified 2D flow model addresses this issue by solving a reduced form of the shallow water equations using a very simple numerical scheme, thus generating a significant increase in computational efficiency over previous hydrodynamic methods. This leads to potential convergence of computational cost between LEMs and hydrodynamic models, and presents an opportunity to combine such schemes. This paper outlines how two such models (the LEM CAESAR and the hydrodynamic model LISFLOOD‐FP) were merged to create the new CAESAR‐Lisflood model, and through a series of preliminary tests shows that using a hydrodynamic model to route flow in an LEM affords many advantages. The new model is fast, computationally efficient and has a stronger physical basis than a previous version of the CAESAR model. For the first time it allows hydrodynamic effects (tidal flows, lake filling, alluvial fans blocking valley floor) to be represented in an LEM, as well as producing noticeably different results to steady flow models. This suggests that the simplification of using steady flow in existing LEMs may bias their findings significantly. Copyright © 2013 John Wiley & Sons, Ltd.     

Mapping Potential Groundwater‐Dependent Ecosystems for Sustainable Management

Ecosystems which rely on either the surface expression or subsurface presence of groundwater are known as groundwater‐dependent ecosystems (GDEs). A comprehensive inventory of GDE locations at an appropriate management scale is a necessary first‐step for sustainable management of supporting aquifers; however, this information is unavailable for most areas of concern. To address this gap, this study created a two‐step algorithm which analyzed existing geospatial and remote sensing data to identify potential GDEs at both state/province and aquifer/basin scales. At the state/province scale, a geospatial information system (GIS) database was constructed for Texas, including climate, topography, hydrology, and ecology data. From these data, a GDE index was calculated, which combined vegetative and hydrological indicators. The results indicated that central Texas, particularly the Edwards Aquifer region, had highest potential to host GDEs. Next, an aquifer/basin scale remote sensing‐based algorithm was created to provide more detailed maps of GDEs in the Edwards Aquifer region. This algorithm used Landsat ETM+ and MODIS images to track the changes of NDVI for each vegetation pixel. The NDVI dynamics were used to identify the vegetation with high potential to use groundwater—such plants remain high NDVI during extended dry periods and also exhibit low seasonal and inter‐annual NDVI changes between dry and wet seasons/years. The results indicated that 8% of natural vegetation was very likely using groundwater. Of the potential GDEs identified, 75% were located on shallow soil averaging 45 cm in depth. The dominant GDE species were live oak, ashe juniper, and mesquite.     

Soil detachment and transport on field‐ and laboratory‐scale interrill areas: erosion processes and the size‐selectivity of eroded sediment

Field‐ and laboratory‐scale rainfall simulation experiments were carried out in an investigation of the temporal variability of erosion processes on interrill areas, and the effects of such variation upon sediment size characteristics. Poorly aggregated sandy soils from the semi‐arid environment of Senegal, West Africa, were used on both a 40 m² field plot and a 0·25 m² laboratory plot; rainfall intensity for all experiments was 70 mm h^?1 with a duration of 1 to 2 hours. Time‐series measurements were made of the quantity and the size distribution of eroded material: these permitted an estimate of the changing temporal balance between the main erosion processes (splash and wash). Results from both spatial scales showed a similar temporal pattern of runoff generation and sediment concentration. For both spatial scales, the dominant erosional process was detachment by raindrops; this resulted in a dynamic evolution of the soil surface under raindrop impact, with the rapid formation of a sieving crust followed by an erosion crust. However, a clear difference was observed between the two scales regarding the size of particles detached by both splash and wash. While all measured values were lower than the mean weight diameter (MWD) value of the original soil (mean 0·32 mm), demonstrating the size‐selective nature of wash and splash processes, the MWD values of washed and splashed particles at the field scale ranged from 0·08 to 0·16 mm and from 0·12 to 0·30 mm respectively, whereas the MWD values of washed and splashed particles at the laboratory scale ranged from 0·13 to 0·29 mm and from 0·21 to 0·32 mm respectively. Thus only at the field scale were the soil particles detached by splash notably coarser than those transported by wash. This suggests a transport‐limited erosion process at the field scale. Differences were also observed between the dynamics of the soil loss by wash at the two scales, since results showed wider scatter in the field compared to the laboratory experiments. This scatter is probably related to the change in soil surface characteristics due to the size‐selectivity of the erosion processes at this spatial scale. Copyright © 2006 John Wiley & Sons, Ltd.     

Raindrop‐induced saltation and the enrichment of sediment discharged from sheet and interrill erosion areas

Sheet and interrill erosion areas are sources of soil material rich in nutrients and pollutants. The loss of soil, nutrients and other chemicals from these areas is a matter of concern both in terms of maintaining soil productivity and the health of offsite environments. Many experiments on rainfall erosion have shown enrichment of fine material, nutrients and other chemicals in the sediment discharged for sheet and interrill erosion areas, but often these results were obtained over short periods of time. A qualitative mechanistic model of raindrop‐induced saltation is used to illustrate how this transport mechanism influences the composition of sediment discharged by rain‐impacted flow. Initially, fast moving particles are enriched in the sediment discharge but, over time, during a rainfall event, slower moving particles become more represented. Raindrop‐induced saltation promotes the storage of material on the soil surface with a coarser composition than the original soil. Winnowing of material from this storage by the development of flow‐driven saltation during high‐intensity events can modify the composition of the sediment discharged later by raindrop‐induced saltation. Given stable soil particles, the composition of the sediment discharged at the steady state is the same as the original soil. Enrichment is a non‐steady‐state phenomenon and failure to recognize the transient nature of enrichment may lead to inappropriate interpretation of the implications of the results from short‐term experiments. Copyright © 2011 John Wiley & Sons, Ltd.     

The shape of patterns to come: from initial formation to long‐term evolution

Many studies focus on the emergence and development of rhythmic landscape patterns. In this contribution we explore the different behaviors found as patterns evolve; the trajectories that patterns exhibit as they transit from infinitesimal‐amplitude perturbation to a statistically steady state (or in some cases to continued statistical evolution). The variety of behaviors observed, either through field and laboratory experiments or numerical modeling, can be reduced to four classes: (a) simple stabilization where predictions based on the initial growth of small perturbations corresponds with the characteristics of patterns observed in nature; (b) significant pattern coarsening en route to saturated wavelength, where non‐linear interactions between finite‐amplitude pattern elements change the geometric properties of a pattern as it approaches steady‐state; (c) perpetual coarsening where the wavelength associated with the emerging pattern continues to grow over time and is only limited by physical boundaries or external constrains; (d) slow evolution toward a different attractor, a novel behavior observed in numerical modeling that involves profound temporal changes in pattern characteristics. Within these classes we also observe generalizable non‐linear behaviors: dependence on initial conditions, the emergence of pattern‐scale variables such as pattern defects, and the presence of multiple stable states. Predicting the shape of patterns to come remains a challenge – one that we suggest requires a range of modeling approaches to address both initial instabilities and the emergent properties of evolving patterns, which involve disparate forms of non‐linear interactions. Consideration of generic system behaviors at the pattern scale could enhance future pattern formation studies, facilitating appropriate pairings of analysis approaches and pattern‐evolution modes. Copyright © 2013 John Wiley & Sons, Ltd.     

New insights on the origin of troctolites from the breakaway area of the Godzilla Megamullion (Parece Vela back‐arc basin): The role of melt‐mantle interaction on the composition of the lower crust

The troctolites and olivine‐gabbros from the Dive 6 K‐1147 represent the most primitive gabbroic rocks collected at the Godzilla Megamullion, a giant oceanic core complex formed at an extinct spreading segment of the Parece Vela back‐arc basin (Philippine Sea). Previous investigations have shown that these rocks have textural and major elements mineral compositions consistent with a formation through multistage interaction between mantle‐derived melts and a pre‐existing ultramafic matrix. New investigations on trace element mineral compositions basically agree with this hypothesis. Clinopyroxenes and plagioclase have incompatible element signatures similar to that of typical‐MORB. However, the clinopyroxenes show very high Cr contents (similar to those of mantle clinopyroxene) and rim having sharply higher Zr/REE ratios with respect to the core. These features are in contrast with an evolution constrained by fractional crystallization processes, and suggest that the clinopyroxene compositions are controlled by melt‐rock interaction processes. The plagioclase anorthite versus clinopyroxene Mg#[Mg/(Mg + Fe^Tot)] correlation of the Dive 6 K‐1147 rocks shows a trend much steeper than those depicted by other oceanic gabbroic sections. Using a thermodynamic model, we show that this trend is reproducible by fractionation of melts assimilating 1 g of mantle peridotite per 1 °C of cooling. This model predicts the early crystallization of high Mg# clinopyroxene, consistent with our petrological observation. The melt‐peridotite interaction process produces Na‐rich melts causing the crystallization of plagioclase with low anorthite component, typically characterizing the evolved gabbros from Godzilla Megamullion.     

Refinement of Bouwer‐Rice Baildown Test Analysis

The recoverability of light nonaqueous phase liquids (LNAPL) in the subsurface can be estimated using LNAPL transmissivity. LNAPL transmissivity is analogous to aquifer transmissivity in that it represents the volume of LNAPL that flows through a unit width of a porous medium for a unit gradient in a unit time. Methods for estimating LNAPL transmissivity from baildown test data have been modified from the Bouwer and Rice (1976) slug test method by Lundy and Zimmerman (1996) and Huntley (2000). The primary assumptions when estimating LNAPL transmissivity with the Bouwer‐Rice method include, a quasi‐steady‐state model for recharge to the well (the model assumes steady‐state radial flow to the well with rate dependent well drawdown and no storage effects) and that the ratio of change in LNAPL drawdown to change in LNAPL thickness at the well (ds/db) is constant. This ratio will be referred to as the j‐ratio. Rather than having to meet boundary conditions for a predetermined j‐ratio value, each baildown test provides the data to estimate the j‐ratio value that is unique to that data set. This calculation methodology in turn results in the Bouwer‐Rice method being applicable to a wider range of baildown tests where it is not required for the potentiometric surface or LNAPL/water interface to remain constant.     

The dynamics of subsurface water divides—watersheds of Lake Stechlin and neighbouring lakes

Undoubtedly watersheds and their limiting groundwater divides are not steady but change. Nevertheless the temporal and spatial scales of the changes are seldom explored and rarely can be quantified. The Lake Stechlin region offers the exceptional chance to study the transient dynamics of catchments, the movement of subsurface water divides and their disappearance and reappearance. Based on long‐term data of surface water and groundwater levels a computer model delivers illustrative pictures of the dynamics of natural water bodies. A simple analytically derived criterion for the occurrence of groundwater divides shows remarkable coincidence with observations and with the output from the numerical model. Copyright © 2001 John Wiley & Sons, Ltd.     

Scale invariance and self‐similarity in kinematic wave overland flow in space and time

Fractals are famous for their self‐similar nature at different spatial scales. Similar to fractals, solutions of scale invariant processes are self‐similar at different space–time scales. This unique property of scale‐invariant processes can be utilized to translate the solution of the processes at a much larger or smaller space–time scale (domain) based on the solution calculated on the original space–time scale. This study investigates scale invariance conditions of kinematic wave overland flow process in one‐parameter Lie group of point transformations framework. Scaling (stretching) transformation is one of the one‐parameter Lie group of point transformations and it has a unique importance among the other transformations, as it leads to the scale invariance or scale dependence of a process. Scale invariance of a process yields a self‐similar solution at different space–time scales. However, the conditions for the process to be scale invariant usually dictate various relationships between the scaling coefficients of the dependent and independent variables of the process. Therefore, the scale invariance of a process does not assure a self‐similar solution at any arbitrary space and time scale. The kinematic wave overland flow process is modelled mathematically as initial‐boundary value problem. The conditions to be satisfied by the system of governing equations as well as the initial and boundary conditions of the kinematic wave overland flow process are established in order for the process to be scale invariant. Also, self‐similarity of the solution of the kinematic wave overland flow under the established invariance conditions is demonstrated by various numerical example problems. Copyright © 2011 John Wiley & Sons, Ltd.     

Analytical studies on transient groundwater flow induced by land reclamation using different fill materials

Land reclamation may have a significant influence on groundwater regimes. Analytical solutions have been developed in the past to study the impact of land reclamation on a steady‐state groundwater flow and transient flow in fill materials, assuming that the reclamation site consists of a single zone of uniform hydraulic parameters. In this paper, we derive analytical solutions to describe the transient water table change in response to multi‐stage land reclamation where the fill material is uniform in each stage but the hydraulic conductivity of the fill material varies from stage to stage. By introducing the method of separation of variables, we develop a transient analytical solution to study the impact of land reclamation consisting of fill material with different hydraulic properties on groundwater dynamics. The results show that the water table first increases significantly into the reclaimed zone following the fill material deposition, and then the increase gradually propagates into the original aquifer. The change of water table in the original aquifer mainly depends on the value of hydraulic conductivity of the fill materials. Examples in this paper illustrate how the aquifer system experiences a long time unsteady‐state flow as a result of the reclamation, and it takes at least tens of years for the system to approach a new equilibrium. It is suggested that for a large‐scale reclamation project, the response of the groundwater regime to reclamation should be carefully studied. Copyright © 2013 John Wiley & Sons, Ltd.     

Form and growth of bars in a wandering gravel‐bed river

The changing form of developing alluvial river bars has rarely been studied in the field, especially in the context of the fixed, compound, mainly alternate gravel bars that are the major morphological feature of the wandering style. Century scale patterns of three‐dimensional growth and development, and the consequent scaling relations of such bars, are examined along the gravel‐bed reach of lower Fraser River, British Columbia, Canada. A retrospective view based on maps and aerial photographs obtained through the twentieth century shows that individual bars have a life history of about 100 years, except in certain, protected positions. A newly formed gravel bar quickly assumes its ultimate thickness and relatively quickly approaches its equilibrium length. Growth continues mainly by lateral accretion of unit bars, consistent with the lateral style of instability of the river. Bar growth is therefore allometric. Mature bars approach equilibrium dimensions and volume that scale with the overall size of the channel. Accordingly, the bars conform with several published criteria for the ultimate dimensions of alternate barforms. Sand bars, observed farther downstream, have notably different morphology. Fraser River presents a typical wandering channel planform, exhibiting elements of both meandered and low‐order braided channels. Hydraulic criteria to which the Fraser bars conform illustrate why this planform develops and persists. The modest rate of bed material transfer along the channel – typical of the wandering type – determines a century‐length time scale for bar development. This time scale is consistent with estimates that have been made for change of the macroform elements that determine the overall geometry of alluvial channels. Copyright © 2009 John Wiley & Sons, Ltd.     

A multi‐geochronological study of the Hakusan volcano,central Japan

The Hakusan volcano, central Japan, is located in a region where two subducting plates (the Pacific Plate and the Philippine Sea Plate) overlap near the junction of four plates adjacent to the Japanese Islands (the Pacific Plate, the Philippine Sea Plate, the Eurasia Plate, and the North American Plate). The Hakusan volcano consists of products from four major volcanic episodes: Kagamuro, Ko‐hakusan, and Shin‐Hakusan I and II. To date the eruption events of the Hakusan volcano we applied thermoluminescence and fission track methods. ²³⁸U(²³⁴U)–²³⁰Th disequilibrium and ²⁰⁶Pb/²³⁸U methods were applied to date the zircon crystallization ages for estimating the magma residence time before the eruptions. The eruption ages we obtained are ca 250 ka for Kagamuro, ca 100 ka and ca 60 ka for Ko‐Hakusan, ca 50 ka for Shin‐Hakusan I, and <10 ka for Shin‐Hakusan II. They are concordant with previous reports based on K–Ar dating. Some of the pyroclastic rocks, possibly originating from Shin‐Hakusan II activities, are dated to be ca 36 ka or 50 ka, and belong to the Shin‐Hakusan I activity. The zircon crystallization ages show several clusters prior to eruption. The magma residence time was estimated for each volcanic activity by comparing the major crystallization events and eruption ages, and we found a gradual decrease from ca. 500 ky for the Kagamuro activity to ca. 5 ky for the Shin‐Hakusan II activity. This decrease in residence time may be responsible for the decrease in volume of erupted material estimated from the current topography of the region. The scale of volcanic activity, which was deduced from the number of crystallized zircons, is more or less constant throughout the Hakusan volcanic activity. Therefore, the decrease in magma residence time is most likely the result of stress field change.     

A Practical Modeling Framework for Non‐Fickian Transport and Multi‐Species Sequential First‐Order Reaction

Many studies indicate that small‐scale heterogeneity and/or mobile–immobile mass exchange produce transient non‐Fickian plume behavior that is not well captured by the use of the standard, deterministic advection‐dispersion equation (ADE). An extended ADE modeling framework is presented here that is based on continuous time random walk theory. It can be used to characterize non‐Fickian transport coupled with simultaneous sequential first‐order reactions (e.g., biodegradation or radioactive decay) for multiple degrading contaminants such as chlorinated solvents, royal demolition explosive, pesticides, and radionuclides. To demonstrate this modeling framework, new transient analytical solutions are derived and are inverted in Laplace space. Closed‐form, steady‐state, multi‐species analytical solutions are also derived for non‐Fickian transport in highly heterogeneous aquifers with linear sorption–desorption and matrix diffusion for use in spreadsheets. The solutions are general enough to allow different degradation rates for the mobile and immobile zones. The transient solutions for multi‐species transport are applied to examine the effects of source remediation on the natural attenuation of downgradient plumes of both parent and degradation products in highly heterogeneous aquifers. Results for representative settings show that the use of the standard, deterministic ADE can over‐estimate cleanup rates and under‐predict the cleanup timeframe in comparison to the extended ADE analytical model. The modeling framework and calculations introduced here are also applied for a 30 year groundwater cleanup program at a site in Palm Bay, Florida. The simulated plume concentrations using the extended ADE exhibited agreement with observed long concentration tails of trichloroethene, cis 1,2 DCE, and VC that remained above cleanup goals.     

Wellbore flow‐rate solution for a constant‐head test in two‐zone finite confined aquifers

The solution describing the wellbore flow rate in a constant‐head test integrated with an optimization approach is commonly used to analyze observed wellbore flow‐rate data for estimating the hydrogeological parameters of low‐permeability aquifers. To our knowledge, the wellbore flow‐rate solution for the constant‐head test in a two‐zone finite‐extent confined aquifer has never been reported so far in the literature. This article is first to develop a mathematical model for describing the head distribution in the two‐zone aquifer. The Laplace domain solutions for the head distributions and wellbore flow rate in a two‐zone finite confined aquifer are derived using the Laplace transform, and their corresponding time domain solutions are then obtained using the Bromwich integral method and residue theorem. These new solutions are expressed in terms of an infinite series with Bessel functions and not straightforward to calculate numerically. A large‐time solution for the wellbore flow rate is therefore developed by employing the relationship of small Laplace variable versus large time variable and L'Hospital's rule. The result shows that the large‐time solution is identical to the steady‐state solution obtained after applying the Tauberian theorem into the Laplace domain solution. This large‐time solution can reduce to the Thiem equation in the case of no skin. Finally, the newly developed solution is used to investigate the effects of outer boundary distance and conductivity ratio on the wellbore flow rate. Copyright © 2011 John Wiley & Sons, Ltd.