Multiple scattering ofSH waves in 2-D media with many cavities

The full waveform synthetic seismogram of multiple scatteredSH waves by many cylindrical cavities in two-dimensional homogeneous elastic media is computed. We used the so-called single-layer potential integral representation of the scattered field and a discretization scheme with line source distribution for each cavity. The total field is the sum of the incident wave plus the field radiated from all sources, each multiplied by an unknown complex constant representing its strength. These constants are determined by imposing the appropriate boundary conditions in the least-squares sense. Here we solve scattering problems involving one, two, four, twelve and fifty cavities regularly distributed in a half-space. The seismograms computed along the free-surface show regions where the incident wave is strongly attenuated, as well as the arrivals of all multiple scattered phases. The accuray of the method is estimated from the degree of agreement of our solution for one cavity with the corresponding analytical solution, and also from the magnitude of the residual tractions along the boundaries of two cavities separated at various distances. Finally we apply the method to compute the case of fifty cylindrical cavities, each of radiusa, randomly distributed in a region 80a wide by 30a deep in a half-space. The value of scattering loss is obtained from the amplitude decay of the primary wave with distance for wavelengths in the range from 1.7a to 13.3a, using the synthetic seismogram calculated for the same distribution of 50 cavities as above, but in full-space.     

Modelling river water temperature using deterministic,empirical, and hybrid formulations in a Mediterranean stream

River water temperature is a common target of water quality models at the watershed scale, owing to its principal role in shaping biogeochemical processes and in stream ecology. Usually, models include physically‐based, deterministic formulations to calculate water temperatures from detailed meteorological information, which usually comes from meteorological stations located far from the river reaches. However, alternative empirical approaches have been proposed, that usually depend on air temperature as master variable. This study explored the performance of a semidistributed water quality application modelling river water temperature in a Mediterranean watershed, using three different approaches. First, a deterministic approach was used accounting for the different heat exchange components usually considered in water temperature models. Second, an empirical approximation was applied using the equilibrium temperature concept, assuming a linear relationship with air temperature. And third, a hybrid approach was constructed, in which the temperature equilibrium concept and the deterministic approach were combined. Results showed that the hybrid approach gave the best results, followed by the empirical approximation. The deterministic formulation gave the worst results. The hybrid approach not only fitted daily river water temperatures, but also adequately modelled the daily temperature range (maximum–minimum daily temperature). Other river water features directly dependent on water temperature, such as river intrusion depth in lentic systems (i.e. the depth at which the river inflow plunges to equilibrate density differences with lake water), were also correctly modelled even at hourly time steps. However, results for the different heat fluxes between river and atmosphere were very unrealistic. Although direct evidence of discrepancies between meteorological drivers measured at the meteorological stations and the actual river microclimate was not found, the use of models including empirical or hybrid formulations depending mainly on air temperature is recommended if only meteorological data from locations far from the river reaches are available. Copyright © 2008 John Wiley & Sons, Ltd.     

Montecarlo DLA-type simulations of wetting effects in fluid displacement in porous media

In this work, we report diffusion-limited aggregation (DLA)-type Montecarlo computations of a stochastic model of displacement of a viscous fluid by another that preferentially wets a porous medium, for the case when both fluids are immiscible in the absence of buoyancy forces. The model has the aim to simulate cooperative invasion processes found in experiments of immiscible wetting displacement. The model considers the nonlocal effects of the Laplacian pressure field and the capillary forces via hydrodynamic equations in the Darcy regime with a boundary condition for the pressure at the interface. The boundary condition contains two different types of disorder: the capillary term, which constitutes an additive random disorder, and a term containing an effective random surface tension, which couples to a curvature (it constitutes a multiplicative random term that carries nonlocal information of the whole pressure). We generate different displacement patterns for different setting of the parameters of the model. We analyze these patterns by studying the scaling properties of the interface that separate the two fluids and calculating the fractal dimension of the interface. The results show the existence of three distinct regimes of scaling. One regime at the smallest-length scales is due to the multiplicative random disorder together with the nonlocal coupling; it reveals itself in a roughness exponent α ≈ 0.80. Additionally, we find a DLA-type scaling regime with a roughness exponent α ≈ 0.60 at the largest scales and intermediate scaling regime with α ≈ 0.70 corresponding to invasion percolation with trapping. Each regime has definite scaling ranges that depend on the capillary number and the relative wetting tendency of the fluids. The behavior of the fractal dimensions of the interfaces of the aggregates constitutes a further confirmation of the existence of three scaling regimes and the multi-self-affinity of the perimeter of the interface boundaries.     

Surface temperature differences between minerals in crystalline rocks: Implications for granular disaggregation of granites through thermal fatigue

Thermal expansion differences between minerals within rocks under insolation have previously been assumed to drive breakdown by means of granular disaggregation. However, there have been no definitive demonstrations of the efficacy of this weathering mechanism. Different surface temperatures between minerals should magnify thermal expansion differences, and thus subject adjacent minerals to repeated stresses that might cause breakdown through fatigue failure. This work confirms the existence of surface temperature differences between minerals in granitic rocks under simulated short-term temperature fluctuations so as to discriminate their potential for initiating granular disaggregation. The influence of colour, as a surrogate for albedo, and crystal size, as a function of thermal mass are specifically identified because of their ease of quantification. Four rock types with a range of these properties were examined, and subjected to repeated short-term temperature cycles by radiative heating and cooling under laboratory conditions. Results show that while albedo is the main control for overall and individual maximum temperatures, crystal size is the main factor controlling higher temperature differences between minerals. Thus, stones with large differences of mineral sizes can undergo magnified stresses due to thermal expansion differences.     

Multiple Random Walk Simulation: A Fast Method to Map Grade Uncertainty with Large Datasets

Multiple Random Walk Simulation consists of a methodology adapted to run fast simulations if close-spaced data are abundant (e.g., short-term mining models). Combining kriging with the simulation of random walks attempts to approximate traditional simulation algorithm results but at a computationally faster way when there is a large amount of conditioning samples. This paper presents this new algorithm illustrating the situations where the method can be used properly. A synthetic study case is presented in order to illustrate the Multiple Random Walk Simulation and to analyze the speed and goodness of its results against the ones from using Turning Bands Simulation and Sequential Gaussian Simulation.     

Structural characterization of deep-water deposits in a foreland basin, Silla Syncline (Chilean Patagonia), with applications to depositional processes

Based on our detailed structural characterization, we examine possible relationships between thrust faults and strike-slip faults and thrust-cored folds and depositional units in the Silla Syncline, a 4 km wide fold composed of fine-grained mudstone, coarse sandstone and conglomerate deposits of the Cerro Toro Formation in the Magallanes foreland basin, Chilean Patagonia. The syncline is bounded on its western flank by an asymmetric anticline and on its eastern flank by a broad zone of thrust faults and associated folds, which are oriented sub-parallel to the syncline axis. Deposition of the coarse-grained units of the Silla Syncline appears to have taken place in this structurally defined trough controlled primarily by thrust fault related growth structures flanking the syncline.The syncline and surrounding area have also been deformed by two sets of strike-slip faults, one right-lateral and one left-lateral. The strike-slip and thrust faulting operated contemporaneously for much of their active periods, although it appears that thrust faulting, confined within the fine-grained units, initiated slightly earlier than strike-slip faulting. In addition, younger igneous intrusions at high angle to bedding generally localize along the strike-slip faults. The cross-cutting relationships among the intrusions, strike-slip faults, and flexural slip faults show that all these structures were active during the same period, which extends beyond mid-Miocene.These conclusions support the premise that structures in deep-water sediments are important for understanding not only the deformation of a foreland basin, but also its depositional architecture.     

Biochemical alterations in juvenile carp (Cyprinus carpio) exposed to zinc: glutathione reductase as a target

The aim of this study was to investigate biochemical changes in juvenile carp (Cyprinus carpio) exposed to zinc chloride (10, 30 and 100 microM) for a period of 48 h. Zinc exposure caused a concentration-dependent reduction in glutathione reductase (GR) activity in gills, liver and brain. Gill glutathione S-transferase (GST) was reduced when animals were exposed to the highest concentration of 100 microM zinc. The phosphorylation of p38(MAPK) increased in the brain of fish exposed to zinc 100 microM, while phosphorylation of the extracellular signal-regulated protein kinase 1/2 (ERK1/2) and c-Jun N-terminal protein kinase 1/2 (JNK1/2) remained unchanged. Expression of proteins HSP60 and HSP70 were not affected by zinc exposure. Considering the significant concentration-dependent inhibition of GR in all tissues analyzed, this enzyme could be a potential biomarker of exposure to zinc, which has to be confirmed.     

Towards a comprehensive approach to parameter estimation in land surface parameterization schemes

In climate models, the land–atmosphere interactions are described numerically by land surface parameterization (LSP) schemes. The continuing improvement in realism in these schemes comes at the expense of the need to specify a large number of parameters that are either directly measured or estimated. Also, an emerging problem is whether the relationships used in LSPs are universal and globally applicable. One plausible approach to evaluate this is to first minimize uncertainty in model parameters by calibration. In this paper, we conduct a comprehensive analysis of some model diagnostics using a slightly modified version of the Simple Biosphere 3 model for a variety of biomes located mainly in the Amazon. First, the degree of influence of each individual parameter in simulating surface fluxes is identified. Next, we estimate parameters using a multi‐operator genetic algorithm applied in a multi‐objective context and evaluate simulations of energy and carbon fluxes against observations. Compared with the default parameter sets, these parameter estimates improve the partitioning of energy fluxes in forest and cropland sites and provide better simulations of daytime increases in assimilation of net carbon during the dry season at forest sites. Finally, a detailed assessment of the parameter estimation problem was performed by accounting for the decomposition of the mean squared error to the total model uncertainty. Analysis of the total prediction uncertainty reveals that the parameter adjustments significantly improve reproduction of the mean and variability of the flux time series at all sites and generally remove seasonality of the errors but do not improve dynamical properties. Our results demonstrate that error decomposition provides a meaningful and intuitive way to understand differences in model performance. To make further advancements in the knowledge of these models, we encourage the LSP community to adopt similar approaches in the future. Copyright © 2012 John Wiley & Sons, Ltd.