Understanding fracture distribution and its relation to knickpoint evolution in the Rio Icacos watershed (Luquillo Critical Zone Observatory,Puerto Rico) using landscape-scale hydrogeophysics

The Rio Icacos watershed in the Luquillo Mountains (Puerto Rico) is unique due to its extremely rapid weathering rates. The watershed is incised into a quartz diorite that has developed a large knickzone defining the river profile. Regolith thickness within the watershed generally decreases from 20 to 30 m at the ridges to several meters in the quartz diorite-dominated valley to tens of centimeters near the major river knickpoint, as determined from previous studies. Above the knickzone, we observe spheroidal corestones, but below this weathering is much less apparent. Measured erosion rates from previous studies are also high in the knickzone compared with upper elevations within the river profile. A suite of near-surface geophysical methods (i.e. ground penetrating radar and terrain conductivity) capable of fast data acquisition in rugged landscapes, was deployed at kilometer scales to characterize critical zone structure. Concentrations of chaotic ground penetrating radar (GPR) reflections and diffraction hyperbolas with low electrical conductivity were observed in vertical zones that outcrop at the land surface as areas of intense fracturing and spheroidally weathered corestones. The width of these fractured and weathered zones showed an increase with proximity to the knickpoint, and was attributed to dilation of these sub-vertical fractures near the knickpoint, as postulated theoretically by a stress model calculated for the topographic variability across the knickzone in the Rio Icacos, and that shows a release of compressive stress near the knickpoint. We hypothesize that erosion rates increase in the knickzone because of this inferred dilation of fractures. Specifically, opened fractures could enhance access of water and in turn promote spalling, erosion, and spheroidal weathering. This study shows that ground-based hydrogeophysical methods used at the landscape-scale (traditionally applied at smaller scales) can be used to explore critical zone architecture at the scales needed to explain the extreme variability in erosion rates across river profiles. © 2018 John Wiley & Sons, Ltd.     

Growing new generations of critical zone scientists

Critical zone (CZ) science is entering its second decade. A new generation of scientists is emerging trained specifically in CZ science and are contributing to advances in environmental science across disciplines. Concurrently, the global scope of CZ science is being elevated as new countries invest in CZ observatories. Global CZ science has great potential to address a diverse array of questions beyond any single discipline. In this commentary we discuss a series of CZ science grand challenges that should be targeted by early‐career researchers: understanding water availability in the CZ; expanding CZ science into new environments; communicating the societal relevance of CZ science including earthcasting to the public; seamlessly integrating biological sciences within the CZ framework; and scaling CZ processes over large spatial and temporal gradients. Targeting these grand challenges will push CZ science well into the future. We also highlight mechanisms for increased leadership within the CZ community. Copyright © 2017 John Wiley & Sons, Ltd.     

岩石的临界孔隙度反演及横波预测(英文)

临界孔隙度模型是利用岩石的临界孔隙度来计算岩石骨架的弹性模量,岩石的临界孔隙度值受到很多因素的影响,而实际应用中通常无法获得准确的临界孔隙度值,只能选取经验临界孔隙度值,就会给岩石物理建模带来误差。本文提出了一种利用纵波速度反演岩石的临界孔隙度的方法,并且把它应用于横波预测中。实验室和测井数据应用结果表明本文提出的方法可以降低以往选取经验值带来的不确定性,并且能够为横波预测提供准确的临界孔隙度值,提高了横波预测的精度。     

Geophysical constraints on deep weathering and water storage potential in the Southern Sierra Critical Zone Observatory

The conversion of bedrock to regolith marks the inception of critical zone processes, but the factors that regulate it remain poorly understood. Although the thickness and degree of weathering of regolith are widely thought to be important regulators of the development of regolith and its water‐storage potential, the functional relationships between regolith properties and the processes that generate it remain poorly documented. This is due in part to the fact that regolith is difficult to characterize by direct observations over the broad scales needed for process‐based understanding of the critical zone. Here we use seismic refraction and resistivity imaging techniques to estimate variations in regolith thickness and porosity across a forested slope and swampy meadow in the Southern Sierra Critical Zone Observatory (SSCZO). Inferred seismic velocities and electrical resistivities image a weathering zone ranging in thickness from 10 to 35 m (average = 23 m) along one intensively studied transect. The inferred weathering zone consists of roughly equal thicknesses of saprolite (P‐velocity < 2 km s^?1) and moderately weathered bedrock (P‐velocity = 2–4 km s^?1). A minimum‐porosity model assuming dry pore space shows porosities as high as 50% near the surface, decreasing to near zero at the base of weathered rock. Physical properties of saprolite samples from hand augering and push cores are consistent with our rock physics model when variations in pore saturation are taken into account. Our results indicate that saprolite is a crucial reservoir of water, potentially storing an average of 3 m³ m^?2 of water along a forested slope in the headwaters of the SSCZO. When coupled with published erosion rates from cosmogenic nuclides, our geophysical estimates of weathering zone thickness imply regolith residence times on the order of 10⁵ years. Thus, soils at the surface today may integrate weathering over glacial–interglacial fluctuations in climate. Copyright © 2013 John Wiley & Sons, Ltd.     

CRITICAL CROSS POWER SPECTRAL DENSITY FUNCTIONS AND THE HIGHEST RESPONSE OF MULTI-SUPPORTED STRUCTURES SUBJECTED TO MULTI-COMPONENT EARTHQUAKE EXCITATIONS

The highest response of multi-supported structures subjected to partially specified multi-component earthquake support motions is considered. The seismic inputs are modelled as incompletely specified vector Gaussian random processes with known autospectral density functions but unknown cross spectral densities and these unknown functions are determined such that the steady state response variance of a given linear system is maximized. The resulting cross power spectral density functions are shown to be dependent on the system properties, autospectra of excitation and the response variable chosen for maximization. It emerges that the highest system response is associated neither with fully correlated support motions, nor with independent motions, but, instead, specific forms of cross power spectral density functions are shown to exist which produce bounds on the response of a given structure. Application of the proposed results is demonstrated by examples on a ground based extended structure, namely, a 1578 m long, three span, suspension cable bridge and a secondary system, namely, an idealized piping structure of a nuclear power plant.     

复杂结构地震动输入方向的基本原理及其影响

对于复杂结构，在动力分析中地震波的最不利输入方向对应着结构的最不利反应。本文从叠加原理出发，详细论述了用反应谱法求解结构的最不利输入方向问题，并通过工程实例来说明其影响程度。     

Seismicity simulation with a mass-spring model and a displacement hardening-softening friction law

Dieterich simulated aftershocks numerically, using a one-dimensional mass-spring model with a time-dependent friction law. But an important precursory phenomenon called quiescence cannot be produced by this model unless, as Mikumo and Miyatake showed with a three-dimensional continuum model, a somewhat arbitrary bimodal distribution of frictional strength is assumed. Here we used the friction law proposed by Stuart, which is a displacement hardening-softening model, and simulated the quiescence. By varying the parameters of the friction law in our mass-spring model, we found a variety of seismicity patterns. When we choose extremely large critical displacement we get a recurrent sequence of creep followed by mainshock without small earthquakes. But when we choose a critical displacement in the same order of magnitude as the slip-weakening critical displacement estimated by Papageorgiou and Aki from strong motion data, we get a normal seismicity pattern, including quiescence before large events. This simple model points to a promising approach for the interpretation of the rupture process during an earthquake by the same physical model.     

Turbulent transport of angular momentum and differential rotation

Abstract

With the help of simplifying approximations, we have derived expressions for the non-diffusive fluxes of the angular momentum which are brought about by the action of Coriolis forces on the convective motion. The original turbulence, which is not perturbed by the Coriolis forces, is considered given and weakly anisotropic, the anisotropy having a preferred radial direction. The eddy viscosities are evaluated. Hence, a closed equation for the angular velocity is derived, and then solved for the case of slow rotation. It is shown that the differential rotation is generated most effectively in the case of moderate rotation when the Rossby number is of order unity. At small Rossby numbers, the rotation differentiality is inhibited. A negative eddy viscosity is suggested for the case of rapid rotation. Some implications for the Sun and other astrophysical objects are discussed.     

大跨斜拉桥地震动最不利输入方向分析

对于复杂结构的地震反应分析,地震波的最不利输入方向对应着结构的最不利反应。通过一座大跨斜拉桥的具体实例分析,研究了地震动输入方向不同时结构反应的变化情况,结果表明：地震动输入最不利方向和结构本身以及所采用地震波均有关,沿轴向输入并不总是最不利输入方向,同时发现,地震波的长周期成分对大跨斜拉桥结构的地震响应有很大影响。     

自组织临界性与地震预测：对目前地震预测问题争论的评述（之一）

目前国际上围绕地震预测问题的争论,很大程度上与地震的自组织临界性模型有关。根据地震ＳＯＣ模型,地震预测被认为是不可能的。事实上,这是一个误解这个误解又是关于地震的ＳＯＣ模型的诸多误解中的一个。在地震的ＳＯＣ模型中,地震的能量－频度关系具有重要的地位。     

Nonlinear baroclinic equilibration at finite supercriticality

The nonlinear equilibration of finite amplitude baroclinic waves in Phillips two-layer model is investigated at finite supercriticality. The aims are to quantify the robustness and relevance of the nonlinear theory of Warn, Gauthier and Pedlosky (WGP) for the evolution of the developing baroclinic wave, and to assess the tightness of pseudomomentum and improved pseudoenergy bounds for disturbance amplitude and energy. A high-resolution numerical model is used to perform a parameter sweep in (β,?W)-space, where β is the inverse criticality of the initial flow, and W is the ratio of the channel width to the (internal) Rossby radius. At low supercriticalities, the main predictions of WGP are found to be accurate at short times, but at long times the fully nonlinear results are found to diverge from WGP's solution. The mechanism for equilibration involves the elimination of the lower layer potential vorticity (PV) gradient, but as the supercriticality increases this is achieved by the roll-up of a train of opposite-signed vortices, rather than by coarse-grain PV homogenization as in WGP. Peak wave amplitudes are typically ≈90% of the maximum attainable under the pseudomomentum bound. New formulae are given for the pseudoenergy bound on disturbance energy which, unlike the WGP solution and the pseudomomentum bound, have non-trivial dependence on W. A detailed assessment is made of the extent to which these bounds are attained.     

地震及其前兆的观测研究与地震预测

通过对地震灾害及预测观察研究现状的分析,解析了地震预测研究中存在的问题和不足,指出了地震及其前兆观测在预测研究中的作用,同时认为,地震前兆观测既不能盲目,也不能随波逐流也进行,其进展以决于地震孕育机理的深入理解和有效观测途径或方法的选择,并简述了２０年来在震前准静电场观测研究中取得进展和体会。     

A low‐dimensional model of bedrock weathering and lateral flow coevolution in hillslopes: 2. Controls on weathering and permeability profiles,drainage hydraulics,and solute export pathways

The advance of a chemical weathering front into the bedrock of a hillslope is often limited by the rate weathering products that can be carried away, maintaining chemical disequilibrium. If the weathering front is within the saturated zone, groundwater flow downslope may affect the rate of transport and weathering—however, weathering also modifies the rock permeability and the subsurface potential gradient that drives lateral groundwater flow. This feedback may help explain why there tends to be neither “runaway weathering” to great depth nor exposed bedrock covering much of the earth and may provide a mechanism for weathering front advance to keep pace with incision of adjacent streams into bedrock. This is the second of a two‐part paper exploring the coevolution of bedrock weathering and lateral flow in hillslopes using a simple low‐dimensional model based on hydraulic groundwater theory. Here, we show how a simplified kinetic model of 1‐D rock weathering can be extended to consider lateral flow in a 2‐D hillslope. Exact and approximate analytical solutions for the location and thickness of weathering within the hillslope are obtained for a number of cases. A location for the weathering front can be found such that lateral flow is able to export weathering products at the rate required to keep pace with stream incision at steady state. Three pathways of solute export are identified: “diffusing up,” where solutes diffuse up and away from the weathering front into the laterally flowing aquifer; “draining down,” where solutes are advected primarily downward into the unweathered bedrock; and “draining along,” where solutes travel laterally within the weathering zone. For each pathway, a different subsurface topography and overall relief of unweathered bedrock within the hillslope is needed to remove solutes at steady state. The relief each pathway requires depends on the rate of stream incision raised to a different power, such that at a given incision rate, one pathway requires minimal relief and, therefore, likely determines the steady‐state hillslope profile.     

提高地方地震台网地震速报水平初探

总结了西昌遥测地震台网地震速报处理技术经验，针对地方台网地震速报中常出现的地震参数误差提出了改进思路。     

重大工程设定地震动确定

概率一致反应谱通常高估设计地震动的中长周期成份,不能真实代表地震环境下的样本实现,为了合理地估计设计地震动反应谱,本文结合工程结合的动力特性和峰值加速度的经验设防标准,建议了一种重大工程设定的地震确定方法,然后依据设定地震和地震动衰减关系确定设计地震动反应谱。     

Investigations into critical earthquake load models within deterministic and probabilistic frameworks

This paper deals with the determination of critical earthquake load models for linear structures subjected to single‐point seismic inputs. The primary objective of this study is to examine the realism in critical excitations and critical responses vis a vis the framework adopted for the study and constraints that these excitations are taken to satisfy. Two alternative approaches are investigated. In the first approach, the critical earthquake is expressed in terms of a Fourier series that is modulated by an enveloping function that imparts transient nature to the inputs. The Fourier coefficients are taken to be deterministic and are constrained to satisfy specified upper and lower bounds. Estimates on these bounds, for a given site, are obtained by analysing past earthquake records from the same site or similar sites. The unknown Fourier coefficients are determined such that the response of a given structure is maximized subjected to these bounds and additional constraints on intensity, peak ground acceleration, peak ground velocity and peak ground displacement. In the second approach, the critical earthquake is modelled as a partially specified non‐stationary Gaussian random process which is defined in terms of a stationary random process of unknown power spectral density (psd) function modulated by a deterministic envelope function. The input is constrained to possess specified variance and average zero crossing rate. Additionally, a new constraint in terms of entropy rate representing the expected level of disorder in the excitation is also imposed. The unknown psd function of the stationary part of the input is determined so that the response of a given structure is maximized. The optimization problem in both these approaches is solved by using sequential quadratic programming method. The procedures developed are illustrated by considering the seismic response of a tall chimney and an earth dam. It is concluded that the imposition of lower and upper bounds on Fourier coefficients in the first approach and constraints on amount of disorder in the second approach are crucial in arriving at realistic critical excitations. Copyright © 2002 John Wiley & Sons, Ltd.     

Critical response of structures to multicomponent earthquake excitation

This paper aims to develop an improved understanding of the critical response of structures to multicomponent seismic motion characterized by three uncorrelated components that are defined along its principal axes: two horizontal and the vertical component. An explicit formula, convenient for code applications, has been derived to calculate the critical value of structural response to the two principal horizontal components acting along any incident angle with respect to the structural axes, and the vertical component of ground motion. The critical response is defined as the largest value of response for all possible incident angles. The ratio r_cr/r_srss between the critical value of response and the SRSS response—corresponding to the principal components of ground acceleration applied along the structure axes—is shown to depend on three dimensionless parameters: the spectrum intensity ratio γ between the two principal components of horizontal ground motion characterized by design spectra A(T_n) and γA(T_n); the correlation coefficient α of responses r_x and r_y due to design spectrum A(T_n) applied in the x‐ and y‐directions, respectively; and β = r_y/r_x. It is demonstrated that the ratio r_cr/r_srss is bounded by 1 and . Thus the largest value of the ratio is , 1.26, 1.13 and 1.08 for γ = 0, 0.5, 0.75 and 0.85, respectively. This implies that the critical response never exceeds times the result of the SRSS analysis, and this ratio is about 1.13 for typical values of γ, say 0.75. The correlation coefficient α depends on the structural properties but is always bounded between −1 and 1. For a fixed value of γ, the ratio r_cr/r_srss is largest if β = 1 and α = ±1. The parametric variations presented for one‐storey buildings indicate that this condition can be satisfied by axial forces in columns of symmetric‐plan buildings or can be approximated by lateral displacements in resisting elements of unsymmetrical‐plan buildings. Copyright © 2000 John Wiley & Sons, Ltd.     

Probing the deep critical zone beneath the Luquillo Experimental Forest,Puerto Rico

Recent work has suggested that weathering processes occurring in the subsurface produce the majority of silicate weathering products discharged to the world's oceans, thereby exerting a primary control on global temperature via the well‐known positive feedback between silicate weathering and CO₂. In addition, chemical and physical weathering processes deep within the critical zone create aquifers and control groundwater chemistry, watershed geometry and regolith formation rates. Despite this, most weathering studies are restricted to the shallow critical zone (e.g. soils, outcrops). Here we investigate the chemical weathering, fracturing and geomorphology of the deep critical zone in the Bisley watershed in the Luquillo Critical Zone Observatory, Puerto Rico, from two boreholes drilled to 37.2 and 27.0 m depth, from which continuous core samples were taken. Corestones exposed aboveground were also sampled. Weathered rinds developed on exposed corestones and along fracture surfaces on subsurface rocks slough off of exposed corestones once rinds attain a thickness up to ~1 cm, preventing the corestones from rounding due to diffusion limitation. Such corestones at the land surface are assumed to be what remains after exhumation of similar, fractured bedrock pieces that were observed in the drilled cores between thick layers of regolith. Some of these subsurface corestones are massive and others are highly fractured, whereas aboveground corestones are generally massive with little to no apparent fracturing. Subsurface corestones are larger and less fractured in the borehole drilled on a road where it crosses a ridge compared with the borehole drilled where the road crosses the stream channel. Both borehole profiles indicate that the weathering zone extends to well below the stream channel in this upland catchment; hence weathering depth is not controlled by the stream level within the catchment and not all of the water in the watershed is discharged to the stream. Copyright © 2013 John Wiley & Sons, Ltd.     

Subduction zone seismicity and the thermo-mechanical evolution of downgoing lithosphere

In this paper we discuss characteristic features of subduction zone seismicity at depths between about 100 km and 700 km, with emphasis on the role of temperature and rheology in controlling the deformation of, and the seismic energy release in downgoing lithosphere. This is done in two steps. After a brief review of earlier developments, we first show that the depth distribution of hypocentres at depths between 100 km and 700 km in subducted lithosphere can be explained by a model in which seismic activity is confined to those parts of the slab which have temperatures below a depth-dependent critical valueT _cr.Second, the variation of seismic energy release (frequency of events, magnitude) with depth is addressed by inferring a rheological evolution from the slab's thermal evolution and by combining this with models for the system of forces acting on the subducting lithosphere. It is found that considerable stress concentration occurs in a reheating slab in the depth range of 400 to 650–700 km: the slab weakens, but the stress level strongly increases. On the basis of this stress concentration a model is formulated for earthquake generation within subducting slabs. The model predicts a maximum depth of seismic activity in the depth range of 635 to 760 km and, for deep earthquake zones, a relative maximum in seismic energy release near the maximum depth of earthquakes. From our modelling it follows that, whereas such a maximum is indeed likely to develop in deep earthquake zones, zones with a maximum depth around 300 km (such as the Aleutians) are expected to exhibit a smooth decay in seismic energy release with depth. This is in excellent agreement with observational data. In conclusion, the incoroporation of both depth-dependent forces and depth-dependent rheology provides new insight into the generation of intermediate and deep earthquakes and into the variation of seismic activity with depth.Our results imply that no barrier to slab penetration at a depth of 650–700 km is required to explain the maximum depth of seismic activity and the pattern of seismic energy release in deep earthquake zones.