On different types of adjustment usable to calculate the parameters of the stream power law

Model parameterization through adjustment to field data is a crucial step in the modeling and the understanding of the drainage network response to tectonic or climatic perturbations. Using as a test case a data set of 18 knickpoints that materialize the migration of a 0.7-Ma-old erosion wave in the Ourthe catchment of northern Ardennes (western Europe), we explore the impact of various data fitting on the calibration of the stream power model of river incision, from which a simple knickpoint celerity equation is derived. Our results show that statistical least squares adjustments (or misfit functions) based either on the stream-wise distances between observed and modeled knickpoint positions at time t or on differences between observed and modeled time at the actual knickpoint locations yield significantly different values for the m and K parameters of the model. As there is no physical reason to prefer one of these approaches, an intermediate least-rectangles adjustment might at first glance appear as the best compromise. However, the statistics of the analysis of 200 sets of synthetic knickpoints generated in the Ourthe catchment indicate that the time-based adjustment is the most capable of getting close to the true parameter values. Moreover, this fitting method leads in all cases to an m value lower than that obtained from the classical distance adjustment (for example, 0.75 against 0.86 for the real case of the Ourthe catchment), corresponding to an increase in the non-linear character of the dependence of knickpoint celerity on discharge.     

Twin structure in coesite studied by high resolution electron microscopy

A synthetic crystal of coesite, SiO₂, has been found to contain microtwins. The (100) reflection twin has been made evident on incoherent (010) planes. The twin structure is directly imaged by high resolution electron microscopy at 2,5 Å resolution, allowing an atomic model to be proposed. This twin could be incoherent on large surfaces and the dislocation density is predicted. The X-ray crystallographic data are analyzed critically in view of the presence of such defects in single crystals. The mutual orientation between a silicon matrix and coesite precipitates in Czochralski-silicon annealed specimens is proposed to be non unique as a consequence of the easy twinning of coesite.     

Reconciling pyroclastic flow and surge: the multiphase physics of pyroclastic density currents

Two end-member types of pyroclastic density current are commonly recognized: pyroclastic surges are dilute currents in which particles are carried in turbulent suspension and pyroclastic flows are highly concentrated flows. We provide scaling relations that unify these end-members and derive a segregation mechanism into basal concentrated flow and overriding dilute cloud based on the Stokes number (S_T), the stability factor (Σ_T) and the dense-dilute condition (D_D). We recognize five types of particle behaviors within a fluid eddy as a function of S_T and Σ_T: (1) particles sediment from the eddy, (2) particles are preferentially settled out during the downward motion of the eddy, but can be carried during its upward motion, (3) particles concentrate on the periphery of the eddy, (4) particles settling can be delayed or ‘fast-tracked’ as a function of the eddy spatial distribution, and (5) particles remain homogeneously distributed within the eddy. We extend these concepts to a fully turbulent flow by using a prototype of kinetic energy distribution within a full eddy spectrum and demonstrate that the presence of different particle sizes leads to the density stratification of the current. This stratification may favor particle interactions in the basal part of the flow and D_D determines whether the flow is dense or dilute. Using only intrinsic characteristics of the current, our model explains the discontinuous features between pyroclastic flows and surges while conserving the concept of a continuous spectrum of density currents.     

Subsurface DC resistivity mapping: approximate 1-D interpretation

Resistivity prospecting is the main tool used to investigate the shallow structure of the ground. A series of new techniques for determining the 2-D and 3-D geometry of the ground is now finding increasing use, but the light and simple Wenner prospecting technique remains a practical and efficient tool for rapidly mapping lateral variations in resistivity. When the resistivity changes are smooth, 1-D modelling can be used to interpret the data, and the criteria governing this approximation can be defined from synthetic data generated by a 3-D slab-model. For a Wenner array, two quadripole configurations can be used, Normal and Dipole-Dipole. For these two configurations the width of the transition zone, the apparent anisotropy effect and the precision of the resistivity values recovered from 1-D inversion differ. However the simultaneous inversion of both sets of data gives better results than for either configuration by itself. Two examples illustrate that in geological contexts where the thickness of the weathered zone causes the changes in the apparent resistivity value, this parameter can be recovered from 1-D inversion.     

Paleoclimate Changes during the Early Oiigocene in the Hoh Xil Region, Northern Tibetan Plateau

Sedimentological, cyclic-stratigraphic, paleomagnetic, and clay-mineralogical studies on the early Oligocene Yaxicuo Group in the Hoh Xil Basin, the largest Cenozoic sedimentary basin in the hinterland of the Tibetan Plateau, provide abundant information of paleoclimate changes. A 350-m thick section in the middle-lower Yaxicuo Group was analyzed to reveal the climatic history that occurred in the Hoh Xil region during the early Oligocene interval 31.30-30.35 Ma, dated with the paleomagnetic chronostratigraphy. The results indicate that add and cold climate dominated the Hoh Xil region during the early Oligocene in general, being related to the global cooling and drying events that occurred in the earliest Oligocene. Within this period, relatively warm and wet climate accompanied by strong tectonic activity occurred in the 31.05-30.75 Ma interval; while add and cold climate and relatively inactive tectonics occurred in the 31.30-31.05 and 30.75-30.35 Ma intervals. Furthermore, spectral analyses of high-temporal resolution paleoclimatic records show orbital periods including eccentricity, obliquity, and precession. It is concluded that paleoclimate changes during the early Oligocene in the Hoh Xil region were forced by both tectonic activity and orbital periods.     

Inverse reconstructions of ecosystem flows in investigating regime shifts: impact of the choice of objective function

Inverse analysis is increasingly used in ecosystem modelling to objectively reconstruct a large number of unknown flows or interactions from a small number of observations. This type of analysis may be useful in relating observed regime shifts in ecosystem structure to underlying processes. Inversions of ecosystem flow networks currently use a constrained least-squares solution which at the same time minimizes the squared norm (the sum of squares) of the reconstructed flows. This minimum norm (MN) inversion is thought to be a parsimonious solution to the ecosystem flow inverse problem, but it may well not reflect how ecosystems are organised. It has been proposed instead that ecosystems evolve to maximize energy/mass flows or that they maximize the information content of the network weighted by ecosystem flows (ascendancy). We used simulated inverse experiments, where inverse analyses are applied to simulations of flow networks, to explore objective functions different than the MN generally used. We could not compute inverse solutions that maximize ascendancy because the objective function is unbounded. We could calculate inversions that maximize flows; however, these generally overestimated the simulated flows, even though the simulations were designed to maximize flows. It appears that the ecosystem flow inverse problem is too under-determined (too few data relative to the number of unknowns) to allow the use of these maximizing goal functions. We introduce a new minimization that simultaneously minimizes the squared flows and the squared differences between flows. This smoothing minimization makes the inverse flows as even as possible and it helps with some technical issues with MN inversions. The simulated inverse experiments indicated that this smoothed norm (SM) is the most robust in comparative analyses of contrasting ecosystem states, such as those that can be associated with regime shifts. Like the MN inversion, the SM inversion has no ecological basis. However, it is a conservative norm that is less likely to produce false differences between the dynamics of regimes.     

Risk assessment of groundwater pollution using sensitivity analysis and a worst-case scenario analysis

This paper illustrates how sensitivity analysis and a worst-case scenario analysis can be useful tools in risk assessment of groundwater pollution. The approach is applied to a study area in Hungary with several known groundwater pollution sources and nearby drinking water production wells. The main concern is whether the contamination sources threaten the drinking water wells of the area. A groundwater flow and transport model is set up to answer this question. Due to limited data availability, the results of this model are associated with large uncertainty. Sensitivity analysis and a worst-case scenario analysis are applied to estimate this uncertainty and build confidence in the model results.     

Data assimilation (4D-VAR) to forecast flood in shallow-waters with sediment erosion

In this paper, the four-dimensional variational data assimilation technique (4D-VAR) is presented as a tool to forecast floods. Our study is limited to purely hydrological flows and supposes that the weather, here a big rain, has been already forecasted by meteorological services. The technique consists in minimizing, in the sense of Lagrange, the cost function: a measure of the difference between calculated data and available observations, here the water level. This is done under constraints that are the equations of the physical model. In our case, we modified the shallow-water equations to include a simplified sediment transport model. The steepest descent algorithm is then used to find the minimum. This is made possible because we can compute analytically the gradient of the cost function by using the adjoint equations of the model. As an application of the 4D-VAR technique, the overflowing of the Chicoutimi River at the Chute-Garneau dam, during the 1996 flood, is investigated. It is found that the 4D-VAR method reduces the error in the water height forecast even when the erosion model is not activated. In terms of Lyapunov exponents, we estimate the predictability horizon of such an event to be about half-an-hour after a big rain. However, this limit of predictability can be increased by using more observations or by using a finer computational grid.     

Changes in stream water quality due to logging of the boreal forest in the Montmorency Forest,Québec

Summer stream water quality was monitored before and following the logging of 50% of the boreal forest within three small watersheds (<50 ha) nested in the ‘Ruisseau des Eaux‐Volées’ Experimental Watershed, Montmorency Forest (Québec, Canada). Logging was conducted in winter, on snow cover according to recommended best management practices (BMPs) to minimize soil disturbance and protect advance growth. A 20‐m forest buffer was maintained along perennial streams. In watershed 7·2, cut‐blocks were located near the stream network and logging was partially allowed within the riparian buffer zone. In watersheds 7·5 and 7·7, logging occurred farther away from the stream network. Observations were also made for watershed 7·3 that collected the runoff from watersheds 7·2 and 7·5, and watershed 7·6, the uproad portion of watershed 7·7. The control watershed 0·2 was contiguous to the impacted watersheds and remained undisturbed. Following clearcutting, changes in summer daily maximum and minimum stream temperatures remained within ± 1 °C while changes in diurnal variation did not decrease by more than 0·5 °C. Concentrations of NO₃^? greatly increased by up to 6000% and concentrations of K⁺ increased by up to 300% during the second summer after logging. Smaller increases were observed for Fe_total (up to 71%), specific conductance (up to 26%), and Mg²⁺ (up to 19%). Post‐logging pH decreased slightly by no more than 7% while PO₄^3? concentration remained relatively constant. Suspended sediment concentrations appeared to increase during post‐logging, but there was not enough pre‐logging data to statistically confirm this result. Logging of moderate intensity and respecting established BMPs may account for the limited changes of water quality parameters and the low exceedances of the criteria for the protection of aquatic life. The proximity of the cutover to the stream network and logging within the riparian zone did not appear to affect water quality. Copyright © 2008 John Wiley & Sons, Ltd.