Liquefaction potential evaluations by energy-based method and stress-based method for various ground motions

An energy-based liquefaction potential evaluation method (EBM) previously developed was applied to a uniform sand model shaken by seismic motions recorded at different sites during different magnitude earthquakes. It was also applied to actual liquefaction case histories in Urayasu city during the 2011 M9.0 Tohoku earthquake and in Tanno-cho during the 2003 M8.0 Tokachi-oki earthquake. In all these evaluations, the results were compared with those by the currently used stress-based method (SBM) under exactly the same seismic and geotechnical conditions. It was found that EBM yields similar results with SBM for several ground motions of recent earthquakes but has easier applicability without considering associated parameters. In Urayasu city, the two methods yielded nearly consistent results by using an appropriate coefficient in SBM for the M9.0 earthquake, though both overestimated the actual liquefaction performance, probably because effects of plasticity and aging on in situ liquefaction strength were not taken into account. In Tanno-cho, EBM could evaluate actual liquefaction performance due to a small-acceleration motion during a far-field large magnitude earthquake while SBM could not.     

Computational model coupling mode II discrete fracture propagation with continuum damage zone evolution

We propose a numerical method that couples a cohesive zone model (CZM) and a finite element‐based continuum damage mechanics (CDM) model. The CZM represents a mode II macro‐fracture, and CDM finite elements (FE) represent the damage zone of the CZM. The coupled CZM/CDM model can capture the flow of energy that takes place between the bulk material that forms the matrix and the macroscopic fracture surfaces. The CDM model, which does not account for micro‐crack interaction, is calibrated against triaxial compression tests performed on Bakken shale, so as to reproduce the stress/strain curve before the failure peak. Based on a comparison with Kachanov's micro‐mechanical model, we confirm that the critical micro‐crack density value equal to 0.3 reflects the point at which crack interaction cannot be neglected. The CZM is assigned a pure mode II cohesive law that accounts for the dependence of the shear strength and energy release rate on confining pressure. The cohesive shear strength of the CZM is calibrated by calculating the shear stress necessary to reach a CDM damage of 0.3 during a direct shear test. We find that the shear cohesive strength of the CZM depends linearly on the confining pressure. Triaxial compression tests are simulated, in which the shale sample is modeled as an FE CDM continuum that contains a predefined thin cohesive zone representing the idealized shear fracture plane. The shear energy release rate of the CZM is fitted in order to match to the post‐peak stress/strain curves obtained during experimental tests performed on Bakken shale. We find that the energy release rate depends linearly on the shear cohesive strength. We then use the calibrated shale rheology to simulate the propagation of a meter‐scale mode II fracture. Under low confining pressure, the macroscopic crack (CZM) and its damaged zone (CDM) propagate simultaneously (i.e., during the same loading increments). Under high confining pressure, the fracture propagates in slip‐friction, that is, the debonding of the cohesive zone alternates with the propagation of continuum damage. The computational method is applicable to a range of geological injection problems including hydraulic fracturing and fluid storage and should be further enhanced by the addition of mode I and mixed mode (I+II+III) propagation. Copyright © 2016 John Wiley & Sons, Ltd.     

Pixel‐oriented land use classification in energy balance modelling

Mass and energy transfer between soil, vegetation and atmosphere is the process that allows to maintain an adequate energy and water balance in the earth–atmosphere system. However, the evaluation of the energy balance components, such as the net radiation and the sensible and latent heat fluxes, is characterized by significant uncertainties related to both the dynamic nature of heat transfer processes and surfaces heterogeneity. Therefore, a detailed land use classification and an accurate evaluation of vegetation spatial distribution are required for an accurate estimation of these variables. For this purpose, in the present article, a pixel‐oriented supervised classification was applied to obtain land use maps of the Basilicata region in Southern Italy by processing three Landsat TM and ETM+ satellite images. An accuracy analysis based on the overall accuracy index and the agreement Khat of Cohen coefficient showed a good performance of the applied classification methodology and a good quality of the obtained maps. Subsequently, these maps were used in the application of a simplified two‐source energy balance model for estimating the actual evapotranspiration at a regional scale. The comparison between the simulations made by applying the simplified two‐source energy balance model and the measurements of evapotranspiration at a lysimetric station located in the study area showed the applicability and the validity of the proposed methodology. Copyright © 2012 John Wiley & Sons, Ltd.     

Energy partitioning over a semi‐arid shrubland in northern China

During the last decade, the widely distributed shrublands in northern China have shown significant signs of recovery from desertification, the result of widespread conservation practices. However, to support the current efforts in conservation, more knowledge is needed on surface energy partitioning and its biophysical controls. Using eddy‐covariance measurements made over a semi‐arid shrubland in northwest China in 2012, we examined how surface energy‐balance components vary on diurnal and seasonal scales, and how biophysical factors control bulk surface parameters and energy exchange. Sensible heat flux (H) exceeded latent heat flux (λE) during most of the year, resulting in an annual Bowen ratio (β, i.e. H/λE) of 2.0. λE exceeded H only in mid‐summer when frequent rainfall co‐occurred with the seasonal peak in leaf area index (LAI). Evapotranspiration reached a daily maximum of 3.3 mm day^?1, and summed to 283 mm yr^?1. The evaporative fraction (EF, i.e. λE/R_n), Priestley–Taylor coefficient (α), surface conductance (g_s) and decoupling coefficient (Ω) were all positively correlated with soil water content (SWC) and LAI. The direct enhancement of λE by high vapour pressure deficit (VPD) was buffered by a concurrent suppression of g_s. The g_s played a direct role in controlling EF and α by mediating the effects of LAI, SWC and VPD. Our results highlight the importance of adaptive plant responses to water scarcity in regulating ecosystem energy partitioning, and suggest an important role for revegetation in the reversal of desertification in semi‐arid areas. Copyright © 2015 John Wiley & Sons, Ltd.     

Exploring hydroclimatic change disparity via the Budyko framework

The Budyko framework characterizes landscape water cycles as a function of climate. We used this framework to identify regions with contrasting hydroclimatic change during the past 50 years in Sweden. This analysis revealed three distinct regions: the mountains, the forests, and the areas with agriculture. Each region responded markedly different to recent climate and anthropogenic changes, and within each region, we identified the most sensitive subregions. These results highlight the need for regional differentiation in climate change adaptation strategies to protect vulnerable ecosystems and freshwater resources. Further, the Budyko curve moved systematically towards its water and energy limits, indicating augmentation of the water cycle driven by changing vegetation, climate and human interactions. This finding challenges the steady state assumption of the Budyko curve and therefore its ability to predict future water cycles. Copyright © 2013 John Wiley & Sons, Ltd.     

Long-Term Variability of the Wind Field over the Indian Ocean Based on ERA-Interim Reanalysis

A detailed study of long-term variability of winds using 30 years of data from the European Centre for Medium-range Weather Forecasts global reanalysis (ERA-Interim) over the Indian Ocean has been carried out by partitioning the Indian Ocean into six zones based on local wind extrema. The trend of mean annual wind speed averaged over each zone shows a significant increase in the equatorial region, the Southern Ocean, and the southern part of the trade winds. This indicates that the Southern Ocean winds and the southeast trade winds are becoming stronger. However, the trend for the Bay of Bengal is negative, which might be caused by a weakening of the monsoon winds and northeast trade winds. Maximum interannual variability occurs in the Arabian Sea due to monsoon activity; a minimum is observed in the subtropical region because of the divergence of winds. Wind speed variations in all zones are weakly correlated with the Dipole Mode Index (DMI). However, the equatorial Indian Ocean, the southern part of the trade winds, and subtropical zones show a relatively strong positive correlation with the Southern Oscillation Index (SOI), indicating that the SOI has a zonal influence on wind speed in the Indian Ocean. Monsoon winds have a decreasing trend in the northern Indian Ocean, indicating monsoon weakening, and an increasing trend in the equatorial region because of enhancement of the westerlies. The negative trend observed during the non-monsoon period could be a result of weakening of the northeast trade winds over the past few decades. The mean flux of kinetic energy of wind (FKEW) reaches a minimum of about 100?W?m^?2 in the equatorial region and a maximum of about 1500?W?m^?2 in the Southern Ocean. The seasonal variability of FKEW is large, about 1600?W?m^?2, along the coast of Somalia in the northern Indian Ocean. The maximum monthly variability of the FKEW field averaged over each zone occurs during boreal summer. During the onset and withdrawal of monsoon, FKEW is as low as 50?W?m^?2. The Southern Ocean has a large variation of about 1280?W?m^?2 because of strong westerlies throughout the year.     

Modelling long-term HFC emissions from India's residential air-conditioning sector: exploring implications of alternative refrigerants,best practices,and a sustainable lifestyle within an integrated assessment modelling framework

India's growing role in the global climate debate makes it imperative to analyse emission reduction policies and strategies across a range of GHGs, especially for under-researched non-CO₂ gases. Hydrofluorocarbons' (HFCs) usage in cooling equipment and subsequent emissions are expected to increase dramatically in India with the phase-out of hydrochlorofluorocarbons (HCFCs) as coolants in air-conditioning equipment. We focus on the residential air-conditioning sector in India and analyse a suite of HFC and alternative coolant gas scenarios for understanding the implications for GHG emissions from this sector within an integrated assessment modelling framework. We find that, if unabated, HFC410A emissions will contribute to 36% of the total global warming impact from the residential air-conditioner sector in India in 2050, irrespective of the future economic growth trajectory, and the remaining 64% is from energy to power residential air-conditioners. A move towards more efficient, low global warming potential (GWP) alternative refrigerants will significantly reduce the cumulative global warming footprint of this sector by 37% during the period 2010–2050, due to gains both from energy efficiency as well as low GWP alternatives. Best practices for reducing direct emissions are important, but only of limited utility, and if a sustainable lifestyle is adopted by consumers with lower floorspace, low GWP refrigerants, and higher building envelope efficiencies, cumulative emissions during 2010–2050 can be reduced by 46% compared to the Reference scenario.

Policy relevance

Our analysis has important implications for Indian climate policy. We highlight that the Indian government's amendment proposal to the Montreal Protocol is a strong signal to the Indian market that the transition away from high GWP refrigerants towards low/zero GWP alternatives will happen sooner or later. The Bureau of Energy Efficiency should extend building energy conservation code policy to residential buildings immediately, and the government should mandate it. Government authorities should set guidelines and mandate reporting of data related to air-conditioner coolant recharge frequency and recovery of scrapped air-conditioner units. For contentious issues like flammability where there is no consensus within the industry, the government needs to undertake an independent technical assessment that can provide unbiased and reliable information to the market.     

Climate change mitigation scenarios and policies and measures: the case of Kazakhstan

This article illustrates the main difficulties encountered in the preparation of GHG emission projections and climate change mitigation policies and measures (P&M) for Kazakhstan. Difficulties in representing the system with an economic model have been overcome by representing the energy system with a technical-economic growth model (MARKAL-TIMES) based on the stock of existing plants, transformation processes, and end-use devices. GHG emission scenarios depend mainly on the pace of transition in Kazakhstan from a planned economy to a market economy. Three scenarios are portrayed: an incomplete transition, a fast and successful one, and even more advanced participation in global climate change mitigation, including participation in some emission trading schemes. If the transition to a market economy is completed by 2020, P&M already adopted may reduce emissions of CO₂ from combustion by about 85 MtCO₂ by 2030 – 17% of the emissions in the baseline (WOM) scenario. One-third of these reductions are likely to be obtained from the demand sectors, and two-thirds from the supply sectors. If every tonne of CO₂ not emitted is valued up to US$10 in 2020 and $20 in 2030, additional P&M may further reduce emissions by 110 MtCO₂ by 2030.     

地壳运动驱动力的探讨——核能与地球演化

本文综合核物理、天文和地质学的最新研究成果,推导出在45亿年前太阳系的全部原始类地行星及其某些卫星(包括地球及月球)上,发生过大规模的、持续的铀、钚链式核裂变,释出了巨大的热能,熔化了整个原始星球。新星球的物质发生了重力分异作用、形成了按密度划分的圈层结构。地球中心是内地核,在高压下压成了固态。铀、钚下沉到内地核顶部停留,继续发生链式核裂变,但改变不了内地核的固态性质。从此处产生的核裂变能,主要以热对流的形式外导。直到前45亿年以后,铀的链式核裂变停止了,只剩下钚的链式核裂变,生成的热能减少了,地球开始冷凝。当熔浆表面温度下降到700~800℃时,较轻的花岗岩质岩浆首先凝固成薄层。它很脆弱、经不起大风浪和潮汐力的冲击,破碎成小块,随着下面对流的玄武岩质熔浆汇集起来,形成大片飘浮物,后来就成为大陸;小片的便成为岛屿。温度逐渐下降,地幔冷凝成固体。温度下降至100℃以下,大气中的水蒸气冷凝成水,下落汇集在地表低凹处,便有了海洋和湖泊。所以海洋底主要分布着玄武岩。从内地核顶部不断产生的核能以"地幔柱"的形式穿过固体地幔,上升至地壳。受坚硬地壳的隔挡、便集中在地壳下构成软流层。当地壳岩石受热软化和可以流动之后,便开始向压力小的方向流动。又受日、月引力和地球自身自转力的作用,软流层的动能又有了增强。这种巨大的能量支配了地壳的升降、褶曲、断裂、变质、岩浆、火山和地震等活动。地球是天体之一,无时不受自然界各种作用力的影响,但支配它演化的主要作用力是核能。本文还简述了核能、核素与太阳系星体演化的关系。