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.     

A greenhouse study on root dynamics of cactus pears, Opuntia ficus-indica and O. robusta

Over the last 10 years a great interest in spineless cactus pear was shown in the drier areas in terms of both fresh fruit and fodder production. However, there is a lack of knowledge on quantitative data on root dynamics of these plants needed to fully understand its potential under water limiting conditions. This study aimed at quantifying the effects of water stress on the growth of tap roots, side roots and rain roots of the species Opuntia ficus-indica (L.) Miller (cultivar Morado—green cladode) and O. robusta Wendl. (cultivar Monterey—blue cladode). One-year-old cladodes were planted in root boxes and pots (2002/2003 season) that were kept in the greenhouse at day/night temperatures of 25–30 °C/15–18 °C. Placing the cladodes flat on the soil, more areoles came in contact with the soil and therefore more roots developed in both species with an average of only 3.4% areole complexes not rooting. Each areole complex formed on average 3 roots. The highest daily tap root growth was 42 and 36 mm for O. ficus-indica and O. robusta, respectively. Tap root growth increased in the morning with water stress for both species but decreased in the afternoon. Side root growth increased with water stress, with that of O. robusta more per tap root than O. ficus-indica. O. robusta showed a finer root system than O. ficus-indica. The side roots grew as much as 8 and 5 mm per day for O. ficus-indica and O. robusta, respectively. Whitish rain roots developed on the established roots within the first hour after rewetting the soil and grew for only 3 days. Rain roots grew up to 7 and 5 mm within a day for O. ficus-indica and O. robusta, respectively. Considering all studied aspects of their roots systems, O. robusta appears to be better adapted to drought (less sensitive to water stress) than O. ficus-indica.     

Traction image method for irregular free surface boundaries in finite difference seismic wave simulation

In this study, we propose a new numerical method, named as Traction Image method, to accurately and efficiently implement the traction-free boundary conditions in finite difference simulation in the presence of surface topography. In this algorithm, the computational domain is discretized by boundary-conforming grids, in which the irregular surface is transformed into a 'flat' surface in computational space. Thus, the artefact of staircase approximation to arbitrarily irregular surface can be avoided. Such boundary-conforming gridding is equivalent to a curvilinear coordinate system, in which the first-order partial differential velocity-stress equations are numerically updated by an optimized high-order non-staggered finite difference scheme, that is, DRP/opt MacCormack scheme. To satisfy the free surface boundary conditions, we extend the Stress Image method for planar surface to Traction Image method for arbitrarily irregular surface by antisymmetrically setting the values of normal traction on the grid points above the free surface. This Traction Image method can be efficiently implemented. To validate this new method, we perform numerical tests to several complex models by comparing our results with those computed by other independent accurate methods. Although some of the testing examples have extremely sloped topography, all tested results show an excellent agreement between our results and those from the reference solutions, confirming the validity of our method for modelling seismic waves in the heterogeneous media with arbitrary shape topography. Numerical tests also demonstrate the efficiency of this method. We find about 10 grid points per shortest wavelength is enough to maintain the global accuracy of the simulation. Although the current study is for 2-D P-SV problem, it can be easily extended to 3-D problem.     

3D behaviour of contamination in landfill sites using 2D resistivity/IP imaging: case studies in Portugal

Aiming at defining a valid spatial contamination model, resistivity and induced polarization (IP) measurements were used to investigate contamination plumes in the vicinity of two municipal landfills (Ovar and ílhavo). Previous geophysical surveys and underground water samples confirmed the contamination. However 2D resistivity/IP surveys enabled in obtaining a more accurate spatial model. The Ovar survey consisted of two profiles with nine Wenner soundings each; the ílhavo survey was carried out along two individual lines using a Wenner standard pseudo-section. In both situations, negative IP values were found associated with positive IP values, which can be explained mainly by 2D or 3D geometric effects caused by the presence of the conductive plumes. The data were modelled using a 2D inversion program (RES2DINV) and the resulting resistivity and chargeability distributions were displayed as pseudo-sections. The resistivity and chargeability pseudo-sections define the contamination plumes and the sedimentary structure. These case studies illustrate the advantages of 2D resistivity/IP surveys for the mapping of shape and dimension of contamination associated with landfills.     

交互的初至波反演及静校正

一种人机交互式的初至波反演方法是在反演过程中加入人为的干预,并利用人的经验控制反演的方向,使结果更加符合实际的地质情况。再将反演的结果与其它的短波长静校正方法相结合,就能够有效地提高静校正的精度。实例证明,这种方法是有效且可行的。     

岩石复电阻率Dias模型及其参数求取方法

Dias模型是众多描述岩石复电阻率的频散特性模型中的一种,它不但参数少,且可由这些参数导出具有特定物理意义的参数。针对Dias模型参数的求取及反演解的不确定性评价,这里提出了一种模拟退火阻尼最小二乘联合反演方法。该方法首先将参数求解空间离散化,利用模拟退火找到全局最优解区间,然后利用阻尼最小二乘法在解空间进行求解。该方法有效地克服了原来两种方法的反演速度慢、反演结果依赖于初值的缺点,并实现了对反演结果多解性的评价。研究结果表明,对于岩石复电阻率数据的处理解释,Dias模型的使用比传统的Cole模型及其变形更具有潜在的优势。