Including spatial distribution in a data‐driven rainfall‐runoff model to improve reservoir inflow forecasting in Taiwan

Multi‐step ahead inflow forecasting has a critical role to play in reservoir operation and management in Taiwan during typhoons as statutory legislation requires a minimum of 3‐h warning to be issued before any reservoir releases are made. However, the complex spatial and temporal heterogeneity of typhoon rainfall, coupled with a remote and mountainous physiographic context, makes the development of real‐time rainfall‐runoff models that can accurately predict reservoir inflow several hours ahead of time challenging. Consequently, there is an urgent, operational requirement for models that can enhance reservoir inflow prediction at forecast horizons of more than 3 h. In this paper, we develop a novel semi‐distributed, data‐driven, rainfall‐runoff model for the Shihmen catchment, north Taiwan. A suite of Adaptive Network‐based Fuzzy Inference System solutions is created using various combinations of autoregressive, spatially lumped radar and point‐based rain gauge predictors. Different levels of spatially aggregated radar‐derived rainfall data are used to generate 4, 8 and 12 sub‐catchment input drivers. In general, the semi‐distributed radar rainfall models outperform their less complex counterparts in predictions of reservoir inflow at lead times greater than 3 h. Performance is found to be optimal when spatial aggregation is restricted to four sub‐catchments, with up to 30% improvements in the performance over lumped and point‐based models being evident at 5‐h lead times. The potential benefits of applying semi‐distributed, data‐driven models in reservoir inflow modelling specifically, and hydrological modelling more generally, are thus demonstrated. Copyright © 2012 John Wiley & Sons, Ltd.     

A model for moisture and heat flow in fractured unsaturated porous media

The present study investigates propagation of a cohesive crack in non‐isothermal unsaturated porous medium under mode I conditions. Basic points of skeleton deformation, moisture, and heat transfer for unsaturated porous medium are presented. Boundary conditions on the crack surface that consist of mechanical interaction of the crack and the porous medium, water, and heat flows through the crack are taken into consideration. For spatial discretization, the extended finite element method is used. This method uses enriched shape functions in addition to ordinary shape functions for approximation of displacement, pressure, and temperature fields. The Heaviside step function and the distance function are exploited as enrichment functions for representing the crack surfaces displacement and the discontinuous vertical gradients of the pressure and temperature fields along the crack, respectively. For temporal discretization, backward finite difference scheme is applied. Problems solved from the literature show the validity of the model as well as the dependency of structural response on the material properties and loading. Copyright © 2016 John Wiley & Sons, Ltd.     

Episodic granite accumulation and extraction from the mid‐crust

The processes of long‐range granitic magma transfer from mid‐ and lower crustal anatectic zones to upper crustal pluton emplacement sites remain controversial in the literature. This is partly because feeder networks that could have accommodated this large‐scale magma transport remain elusive in the field. Existing granite ascent models are based largely on numerical and theoretical studies that seek to demonstrate the viability of fracture‐controlled magma transport through dykes or self‐propagating hydrofractures. In most cases, the models present very little supporting field evidence, such as sufficiently voluminous near‐ or within‐source magma accumulations, to support their basic premises. We document large (deca‐ to hectometre‐scale), steeply dipping and largely homogeneous granite lenses in suprasolidus (~5 kbar, ~750 °C) mid‐crustal rocks in the Damara Belt in Namibia. The lenses are surrounded by and connected to shallowly dipping networks of stromatic leucogranites in the well‐layered gneisses of the deeply incised Husab Gorge. The outcrops define a four‐stage process from (i) the initial formation and growth of large, subvertical magma‐filled lenses as extension fractures developed at high angles to the subhorizontal regional extension in relatively competent wall‐rock layers. This stage is followed by (ii) the simultaneous lateral inflation and (iii) subcritical vertical growth of the lenses to a critical length that (iv) promotes fracture destabilization, buoyancy‐driven upward fracture mobilization and, consequently, vertical magma transport. These field observations are compared with existing numerical models and are used to constrain, by referring to the dimensions of the largest preserved inflated leucogranite lens, an estimate of the minimum fracture length (~100 m) and volume (~2.4 × 10⁵ m³) required to initiate buoyancy‐driven brittle fracture propagation in this particular mid‐crustal section. The critical values and field relationships compare favourably with theoretical models of magma ascent along vertical self‐propagating hydrofractures which close at their tails during propagation. This process leaves behind subtle wake‐like structures and thin leucogranite trails that mark the path of magma ascent. Reutilization of such conduits by repeated inflation and drainage is consistent with the episodic accumulation and removal of magma from the mid‐crust and is reflected in the sheeted nature of many upper crustal granitoid plutons.     

Chemo‐mechanics of cemented granular solids subjected to precipitation and dissolution of mineral species

This paper studies the chemo‐mechanics of cemented granular solids in the context of continuum thermodynamics for fluid‐saturated porous media. For this purpose, an existing constitutive model formulated in the frame of the Breakage Mechanics theory is augmented to cope with reactive processes. Chemical state variables accounting for the reactions between the solid constituents and the solutes in the pore fluid are introduced to enrich the interactions among the microstructural units simulated by the model (i.e., grains and cement bonds). Two different reactive processes are studied (i.e., grain dissolution and cement precipitation), using the chemical variables to describe the progression of the reactions and track changes in the size of grains and bonds. Finally, a homogenization strategy is used to derive the energy potentials of the solid mixture, adopting probability density functions that depend on both mechanical and chemical indices. It is shown that the connection between the statistics of the micro‐scale attributes and the continuum properties of the solid enables the mathematical capture of numerous mechanical effects of lithification and chemical deterioration, such as changes in stiffness, expansion/contraction of the elastic domain, and development of inelastic strains during reaction. In particular, the model offers an interpretation of the plastic strains generated by aggressive environments, which are here interpreted as an outcome of chemically driven debonding and comminution. As a result, the model explains widely observed macroscopic signatures of geomaterial degradation by reconciling the energetics of the deformation/reaction processes with the evolving geometry of the microstructural attributes. Copyright © 2015 John Wiley & Sons, Ltd.     

Late Quaternary sea‐level change and evolution of Belfast Lough,Northern Ireland: new offshore evidence and implications for sea‐level reconstruction

The interplay of eustatic and isostatic factors causes complex relative sea‐level (RSL) histories, particularly in paraglacial settings. In this context the past record of RSL is important in understanding ice‐sheet history, earth rheology and resulting glacio‐isostatic adjustment. Field data to develop sea‐level reconstructions are often limited to shallow depths and uncertainty exists as to the veracity of modelled sea‐level curves. We use seismic stratigraphy, 39 vibrocores and 26 radiocarbon dates to investigate the deglacial history of Belfast Lough, Northern Ireland, and reconstruct past RSL. A typical sequence of till, glacimarine and Holocene sediments is preserved. Two sea‐level lowstands (both max. ?40 m) are recorded at c. 13.5 and 11.5k cal a bp . Each is followed by a rapid transgression and subsequent periods of RSL stability. The first transgression coincides temporally with a late stage of Meltwater Pulse 1a and the RSL stability occurred between c. 13.0 and c. 12.2k cal a bp (Younger Dryas). The second still/slowstand occurred between c. 10.3 and c. 11.5k cal a bp . Our data provide constraints on the direction and timing of RSL change during deglaciation. Application of the Depth of Closure concept adds an error term to sea‐level reconstructions based on seismic stratigraphic reconstructions.     

Evaluation of NLDAS‐2 evapotranspiration against tower flux site observations

The North American Land Data Assimilation System project phase 2 (NLDAS‐2) has run four land surface models for a 30‐year (1979–2008) retrospective period. Land surface evapotranspiration (ET) is one of the most important model outputs from NLDAS‐2 for investigating land–atmosphere interaction or to monitor agricultural drought. Here, we evaluate hourly ET using in situ observations over the Southern Great Plains (Atmospheric Radiation Measurement/Cloud and Radiation Testbed network) for 1 January 1997–30 September 1999 and daily ET u‐sing in situ observations at the AmeriFlux network over the conterminous USA for an 8‐year period (2000–2007). The NLDAS‐2 models compare well against observations, with the National Centers for Environmental Prediction's Noah land surface model performing best, followed, in order, by the Variable Infiltration Capacity, Sacramento Soil Moisture Accounting, and Mosaic models. Daily evaluation across the AmeriFlux network shows that for all models, performance depends on season and vegetation type; they do better in spring and fall than in winter or summer and better for deciduous broadleaf forest and grasslands than for croplands or evergreen needleleaf forest. Copyright © 2014 John Wiley & Sons, Ltd.     

Late Holocene environmental conditions in Coronation Gulf,southwestern Canadian Arctic Archipelago: evidence from dinoflagellate cysts,other non‐pollen palynomorphs,and pollen

Boxcore 99LSSL‐001 (68.095° N, 114.186° W; 211 m water depth) from Coronation Gulf represents the first decadal‐scale marine palynology and late Holocene sediment record for the southwestern part of the Northwest Passage. The record was studied for organic‐walled microfossils (dinoflagellate cysts, non‐pollen palynomorphs), pollen, terrestrial spores, and sediment characteristics. ²¹⁰Pb, ¹³⁷Cs, and three accelerator mass spectrometry ¹⁴C dates constrain the chronology. Three prominent palaeoenvironmental zones were identified. During the interval AD 1470–1680 (Zone I), the climate was warmer and wetter than at present, and environmental conditions were more favourable to biological activity and northward boreal forest migration, with reduced sea‐ice and a longer open‐water (growing) season. The interval AD 1680–1940 (Zone II) records sea‐ice increase, and generally cool, polar conditions during the Little Ice Age. During AD 1940–2000 (Zone III), organic microfossils indicate an extended open‐water season and decreased sea‐ice, with suggested amelioration surpassing that of Zone I. Although more marine studies are needed to place this record into an appropriate context, the succession from ameliorated (Zone I) to cooler, sea‐ice influenced conditions (Zone II) and finally to 20th‐century warming (Zone III) corresponds well with several terrestrial climatic records from the neighbouring mainland and Victoria Island, and with lower‐resolution marine records to the west. Copyright © 2011 John Wiley & Sons, Ltd.     

Vegetation history of the penultimate glacial period (Marine isotope stage 6) at Ioannina,north‐west Greece

The vegetational history of the penultimate glacial period, Marine Isotope Stage (MIS) 6 (c. 185–135 ka), has remained relatively unexplored. Here we present a new record from the Ioannina basin, north‐west Greece, which constitutes the highest‐resolution terrestrial pollen record for this interval produced to date. It shows that the vegetation history of MIS 6 in this region can be divided into two parts: an early period (185–155 ka) with pronounced oscillations in tree population extent, and a later period (155–135 ka) with much smaller tree populations and subdued oscillations. This pattern is analogous to the MIS 3/MIS 2 division during the last glacial in the same sequence, although the early part of MIS 6 had larger Pinus populations and fewer temperate trees relative to the equivalent interval in MIS 3. This implies cooler and wetter conditions, which is somewhat counterintuitive given the high summer insolation during MIS 6e, but is in line with other palaeoclimatic evidence from the Mediterranean. Comparison with North Atlantic records suggests that despite the absence of pronounced iceberg discharges during MIS 6, North Atlantic millennial‐scale variability had a significant downstream impact on tree populations in north‐west Greece. Copyright © 2011 John Wiley & Sons, Ltd.