Numerical study of droplet size dependent chemistry in oceanic,wintertime stratus cloud at southern mid-latitudes

Cloud droplet chemistry is modelled for the first 150 m of rise in a wintertime, mid-latitude, marine stratus cloud using observations made at and near the Cape Grim Baseline Station as a source of input parameters. The emphasis in this work was to study the variation in droplet chemistry as a function of both droplet size and nucleus composition, with a particular focus on the way in which oxidation of dissolved sulfur dioxide varied.At 150 m above the condensation level, solute concentration as a function of droplet size was found to increase by as much as 2 to 3 orders of magnitude for only a factor of 2 increase in droplet radius, primarily as a consequence of the 1/r dependence in the droplet growth equation. This type of size dependence exists at all levels in the model cloud, and has a significant influence on oxidation rate of sulfur dioxide in droplets growing on sulfate nuclei, oxidation by ozone being favoured in the smallest droplets, but oxidation by hydrogen peroxide being favoured in the larger droplets. Oxidation by ozone is favoured at all sizes in droplets formed on sea-salt nuclei as a result of the initially high alkalinity of these droplets, and in the cloud overall is calculated to be the more important oxidation pathway. Although based on a simplified chemical scheme, these results suggest that both size-dependent and nucleus-dependent chemistry of cloud droplets may need to be considered explicitly in cloud modelling work.Volume-weighted mean pH values in the range 5 to 6 were predicted from sensitivity studies in which input variables were varied over reasonable ranges, in agreement with two sets of bulk cloud-water pH data obtained by aircraft near Cape Grim.     

Increasing critical sensitivity of the Load/Unload Response Ratio before large earthquakes with identified stress accumulation pattern

The Load/Unload Response Ratio (LURR) method is proposed for short-to-intermediate-term earthquake prediction [Yin, X.C., Chen, X.Z., Song, Z.P., Yin, C., 1995. A New Approach to Earthquake Prediction — The Load/Unload Response Ratio (LURR) Theory, Pure Appl. Geophys., 145, 701–715]. This method is based on measuring the ratio between Benioff strains released during the time periods of loading and unloading, corresponding to the Coulomb Failure Stress change induced by Earth tides on optimally oriented faults. According to the method, the LURR time series usually climb to an anomalously high peak prior to occurrence of a large earthquake. Previous studies have indicated that the size of critical seismogenic region selected for LURR measurements has great influence on the evaluation of LURR. In this study, we replace the circular region usually adopted in LURR practice with an area within which the tectonic stress change would mostly affect the Coulomb stress on a potential seismogenic fault of a future event. The Coulomb stress change before a hypothetical earthquake is calculated based on a simple back-slip dislocation model of the event. This new algorithm, by combining the LURR method with our choice of identified area with increased Coulomb stress, is devised to improve the sensitivity of LURR to measure criticality of stress accumulation before a large earthquake. Retrospective tests of this algorithm on four large earthquakes occurred in California over the last two decades show remarkable enhancement of the LURR precursory anomalies. For some strong events of lesser magnitudes occurred in the same neighborhoods and during the same time periods, significant anomalies are found if circular areas are used, and are not found if increased Coulomb stress areas are used for LURR data selection. The unique feature of this algorithm may provide stronger constraints on forecasts of the size and location of future large events.     

The effects of carbon dioxide leakage on fractures, and water quality of potable aquifers during geological sequestration of CO2

Since industrial revolution, the ＂greenhouse effect＂ is one of the most important global environmental issues. Of all the greenhouse gases, CO2 is responsible for about 64% of the enhanced ＂greenhouse effect＂, making it the target for mitigation, so reducing anthropogenic discharge of carbon dioxide attracts more and more attention. Geological sequestration of CO2 in deep saline aquifers is one of the most promising options. But because unknown fractures and faults may exist in the caprock layers which can prevent the leakage of CO2, CO2 will leak upward into upper potable aquifers, and lead to adverse impacts on the shallow potable aquifers. In order to assess the potential effect of CO2 leakage from underground storage reservoirs on fractures and water quality of potable aquifers, this study used the non-isothermal reactive geochemical transport code TOUGHREACT developed by Xu et al to establish a simplified 2-D model of CO2 underground sequestration system, which includes deep saline aquifers, caprock layers, and shallow potable aquifers, and study and analyze the changes of mineral and aqueous components. The simulation results indicated that the minerals of deep saline aquifers and fractures should be mainly composed of aluminosilicate and silicate minerals, which not only enhance the mass of CO2 sequestrated by mineral trapping, but also decrease the porosity and permeability of caprock layers and fractures to prevent and reduce CO2 leakage. The results from deep saline aquifers showed that the mass of carbon dioxide trapped by minerals and solution phases is limited, the rest remained as a supercritical phase, and so once the caprock aquifers have some unknown fractures, the free carbon dioxide phase may leak from CO2 geologic sequestration reservoirs by buoyancy.     

Testing earthquake prediction methods: «The West Pacific short-term forecast of earthquakes with magnitude MwHRV ≥ 5.8»

Analyzing the tables and probability maps posted by Yan Y. Kagan and David D. Jackson in April 2002–September 2004 at http://scec.ess.ucla.edu/~ykagan/predictions_index.html and the catalog of earthquakes for the same period, the conclusion is drawn that the underlying method could be used for prediction of aftershocks, while it does not outscore random guessing when main shocks are considered.     

A laboratory experiment to monitor the contact state of a fault by transmission waves

We performed a series of laboratory experiments in which elastic waves were transmitted across a simulated fault. Two types of experiments were carried out: (1) Normal Stress Holding Test (NSHT): normal stress was kept constant for about 3 h without shear stress and transmission waves were observed. (2) Shear Stress Increasing Test (SSIT): shear stress was gradually increased until a stick-slip event occurred. Transmission waves were continuously observed throughout the process of stress accumulation. We focused on the change in transmission waves during the application of shear stress and especially during precursory slips.It was found in NSHT that the amplitude of transmission waves linearly increased with the logarithm of stationary contact time. The increase amounted to a few percent after about 3 h. Creep at asperity contacts is responsible for this phenomenon. From a theoretical consideration, it was concluded that the real contact area increased with the logarithm of stationary contact time.We observed in SSIT a significant increase in wave amplitude with shear stress application. This phenomenon cannot be attributed to the time effect observed in NSHT. Instead, it can be explained by the mechanism of “junction growth” proposed by Tabor. Junction growth yields an increase in real contact area. It is required for junction growth to occur that the material in contact is already plastic under a purely normal loading condition. A computer simulation confirmed that this requirement was satisfied in our experiments. We also found that the rate at which the amplitude increased was slightly reduced prior to a stick-slip event. The onset time of the reduction well coincides with the onset of precursory slip. The cause of the reduction is attributed to the reset of stationary contact time due to displacement. This interpretation is supported by the result of NSHT. Taking the time of stationary contact in SSIT into account, we may expect the change in wave amplitude to be, at most, only a few percent. The observed slight reduction in increasing rate is, in this sense, reasonable. The static stiffness of the fault also decreases with precursory slip. It was also found that low frequency waves are a better indicator of precursory slip than high frequency waves. This might suggest that low frequency waves with longer wavelength are a better indicator of average behavior of faults. The problem, however, merits a further investigation. The shifts in phase were also found to be a good indicator of the change in contact state of the fault. The changes in both amplitude and phase of transmission waves are unifyingly understood through the theory of transmission coefficient presented by Pyrak-Nolte et al. Rough surfaces have a tendency to give larger stick-slips than smooth surfaces. The amount of precursory slip is larger for rough surfaces than for smooth surfaces. Although it was confirmed by a computer simulation that rough surfaces have larger contact diameters than smooth surfaces, the rigorous relationship between the surface roughness (contact diameter) and the amount of precursory slips was not established.     

A way to synchronize models with seismic faults for earthquake forecasting: Insights from a simple stochastic model

Numerical models are starting to be used for determining the future behaviour of seismic faults and fault networks. Their final goal would be to forecast future large earthquakes. In order to use them for this task, it is necessary to synchronize each model with the current status of the actual fault or fault network it simulates (just as, for example, meteorologists synchronize their models with the atmosphere by incorporating current atmospheric data in them). However, lithospheric dynamics is largely unobservable: important parameters cannot (or can rarely) be measured in Nature. Earthquakes, though, provide indirect but measurable clues of the stress and strain status in the lithosphere, which should be helpful for the synchronization of the models.The rupture area is one of the measurable parameters of earthquakes. Here we explore how it can be used to at least synchronize fault models between themselves and forecast synthetic earthquakes. Our purpose here is to forecast synthetic earthquakes in a simple but stochastic (random) fault model. By imposing the rupture area of the synthetic earthquakes of this model on other models, the latter become partially synchronized with the first one. We use these partially synchronized models to successfully forecast most of the largest earthquakes generated by the first model. This forecasting strategy outperforms others that only take into account the earthquake series. Our results suggest that probably a good way to synchronize more detailed models with real faults is to force them to reproduce the sequence of previous earthquake ruptures on the faults. This hypothesis could be tested in the future with more detailed models and actual seismic data.     

大变形有限元法分析边坡稳定性

边坡稳定性分析的重要任务是确定边坡潜在滑面的位置及相应的安全系数,鉴于坡体位移可在现场监测中直接得到,提出边坡潜在滑面确定的位移判据,采用大变形有限元首先计算出边坡的位移场,然后在开挖面的不同高度作塑性区内的水平位移参考线,依次连接参考线上的位移突变点即得到边坡失稳的潜在滑面位置,分析滑面上的应力以及岩体抗剪强度即得到边坡的整体安全系数。     

Towards Formulation of a Space-borne System for Early Warning of Floods: Can Cost-Effectiveness Outweigh Prediction Uncertainty?

The three most important components necessary for functioning of an operational flood warning system are: (1) a rainfall measuring system; (2) a soil moisture updating system; and, (3) a surface discharge measuring system. Although surface based networks for these systems can be largely inadequate in many parts of the world, this inadequacy particularly affects the tropics, which are most vulnerable to flooding hazards. Furthermore, the tropical regions comprise developing countries lacking the financial resources for such surface-based monitoring. The heritage of research conducted on evaluating the potential for measuring discharge from space has now morphed into an agenda for a mission dedicated to space-based surface discharge measurements. This mission juxtaposed with two other upcoming space-based missions: (1) for rainfall measurement (Global Precipitation Measurement, GPM), and (2) soil moisture measurement (Hydrosphere State, HYDROS), bears promise for designing a fully space-borne system for early warning of floods. Such a system, if operational, stands to offer tremendous socio-economic benefit to many flood-prone developing nations of the tropical world. However, there are two competing aspects that need careful assessment to justify the viability of such a system: (1) cost-effectiveness due to surface data scarcity; and (2) flood prediction uncertainty due to uncertainty in the remote sensing measurements. This paper presents the flood hazard mitigation opportunities offered by the assimilation of the three proposed space missions within the context of these two competing aspects. The discussion is cast from the perspective of current understanding of the prediction uncertainties associated with space-based flood prediction. A conceptual framework for a fully space-borne system for early-warning of floods is proposed. The need for retrospective validation of such a system on historical data comprising floods and its associated socio-economic impact is stressed. This proposal for a fully space-borne system, if pursued through wide interdisciplinary effort as recommended herein, promises to enhance the utility of the three space missions more than what their individual agenda can be expected to offer.     

A Bay–River Coupled Model for Storm Surge Prediction Along the Andhra Coast of India

The paper describes a two-dimensional bay–river coupled numerical model for storm surges along the Andhra coast of India. The effect of the Krishna and Godavari rivers on the surge development is analysed. A comparative study of the surge generated by a tropical cyclone with and without the inclusion of rivers is done in detail. Three cyclones that struck the Andhra coast in November 1977, May 1990 and November 1996 were used for the simulation studies. It is found that the idealized model without a river overestimates the sea-level elevation as compared to a more realistic bay–river coupled model. The temporal variation of surge values at the mouth of the rivers is also studied for all three cyclone cases. It is found that the effect of the presence of rivers depends on the strength of the cyclone, its point of landfall and the location of the rivers with respect to the landfall point.