Effect of a submerged bay-mouth breakwater on tsunami behavior analyzed by 2D/3D hybrid model simulation

The hybrid numerical model had been developed to simulate a complicated 3D flow around structures generated by tsunami. In the model, the conventional 2D model is adopted for the wide region far from structures and the 3D non-hydrostatic pressure model is used in the limited region adjacent to structures. The applicability of the model is shown by comparisons of the numerical results with the experimental results of the laboratory model tests and the numerical analysis results of the conventional whole 2D simulation. In addition, the effect of a submerged structure at the opening of a breakwater is discussed from the numerical simulations by the hybrid model. The submerged structure improves the stability of the rubble mound and reduces the tsunami inflow into the bay, while it increases the water surface velocity around the opening of the breakwater. The increase of surface velocity causes the increases of impulsive forces by collision with drafts and so on.     

The 1786 earthquake-triggered landslide dam and subsequent dam-break flood on the Dadu River, southwestern China — Reply

A short reply to discussion by Chanson.     

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.     

大气CO2浓度升高对土壤中不同粒级碳的影响

不同粒级土壤中的碳有着不同的周转规律，在高CO2浓度条件下，它们含量的变化将在一定程度上反映土壤碳是累积还是减少，对明确土壤碳的变化趋势有重要意义．采用田间培养试验初步模拟研究在高CO2浓度条件下土壤不同粒级碳的分布．结果表明，加入秸秆培养1年，由于CO2浓度升高的原因导致在低氮（LN）、常规氮（NN）和高氮（HN）水平下土壤中碳分别增加0．01、1．10、1．22g／kg，表现为粒级〈53μm土壤颗粒中碳分别增加1．53、2．19、2．70g／kg．粒级〈53μmm土壤颗粒碳量的增加，主要是由于其重量分配百分数显著增加36．2％，碳浓度增加5．4％；粒级〉250μm和250～53μm土壤颗粒部分虽然其碳浓度分别增加20．8％和17．3％（P〈0．05），怛由于重量分配百分数分别显著降低22．8％和36．1％，结果碳量降低．试验表明高CO2浓度导致不同粒级土壤的分配及碳浓度的变化；高氮施肥水平下有增加土壤碳量特别是小粒级土壤碳量的趋势．     

Centennial to millennial geomagnetic secular variation

A time-varying spherical harmonic model of the palaeomagnetic field for 0–7 ka is used to investigate large-scale global geomagnetic secular variation on centennial to millennial scales. We study dipole moment evolution over the past 7 kyr, and estimate its rate of change using the Gauss coefficients of degree 1 (dipole coefficients) from the CALS7K.2 field model and by two alternative methods that confirm the robustness of the predicted variations. All methods show substantial dipole moment variation on timescales ranging from centennial to millennial. The dipole moment from CALS7K.2 has the best resolution and is able to resolve the general decrease in dipole moment seen in historical observations since about 1830. The currently observed rate of dipole decay is underestimated by CALS7K.2, but is still not extraordinarily strong in comparison to the rates of change shown by the model over the whole 7 kyr interval. Truly continuous phases of dipole decrease or increase are decadal to centennial in length rather than longer-term features. The general large-scale secular variation shows substantial changes in power in higher spherical harmonic degrees on similar timescales to the dipole. Comparisons are made between statistical variations calculated directly from CALS7K.2 and longer-term palaeosecular variation models: CALS7K.2 has lower overall variance in the dipole and quadrupole terms, but exhibits an imbalance between dispersion in   g ¹₂  and   h ¹₂  , suggestive of long-term non-zonal structure in the secular variations.     

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.