A Time Monte Carlo method for addressing uncertainty in land-use change models

One of the main objectives of land-use change models is to explore future land-use patterns. Therefore, the issue of addressing uncertainty in land-use forecasting has received an increasing attention in recent years. Many current models consider uncertainty by including a randomness component in their structure. In this paper, we present a novel approach for tuning uncertainty over time, which we refer to as the Time Monte Carlo (TMC) method. The TMC uses a specific range of randomness to allocate new land uses. This range is associated with the transition probabilities from one land use to another. The range of randomness is increased over time so that the degree of uncertainty increases over time. We compare the TMC to the randomness components used in previous models, through a coupled logistic regression-cellular automata model applied for Wallonia (Belgium) as a case study. Our analysis reveals that the TMC produces results comparable with existing methods over the short-term validation period (2000–2010). Furthermore, the TMC can tune uncertainty on longer-term time horizons, which is an essential feature of our method to account for greater uncertainty in the distant future.     

未来的国际天球参考架

本文简述了国际天球参考架的发展历史和现在射电参考架的现状—基准源选择的标准和参考架的稳定性。描述了地面上光学观测在依巴谷参考架的维持和加密的一系列工作。介绍由天体测量卫星GAIA和SIM给出的天球参考架可能逵到的精度。详述了在今后十年中地面天体测量的作用以及正在开展有关天球参考架的研究课题 ,同时也列出了我国正在和即将开展天体测量的几个研究课题     

Simulations of strong ground motion in SW Iberia for the 1969 February 28 (Ms= 8.0) and the 1755 November 1 (M∼ 8.5) earthquakes – II. Strong ground motion simulations

This is the second paper of a series of two concerning strong ground motion in SW Iberia due to earthquakes originating from the adjacent Atlantic area. The aim of this paper is to use the velocity model that was proposed and validated in the companion paper for seismic intensity modelling of the 1969 ( M _s= 8.0) and 1755 ( M = 8.5–8.7) earthquakes.
First, we propose a regression to convert simulated values of Peak Ground Velocity (PGV) into Modified Mercalli Intensity (MMI) in SW Iberia, and using this regression, we build synthetic isoseismal maps for a large ( M _s= 8.0) earthquake that occurred in 1969. Based on information on the seismic source provided by various authors, we show that the velocity model effectively reproduces macroseismic observations in the whole region. We also confirm that seismic intensity distribution is very sensitive to a small number of source parameters: rupture directivity, fault strike and fault dimensions. Then, we extrapolate the method to the case of the great ( M = 8.5–8.7) 1755 earthquake, for a series of hypotheses recently proposed by three authors about the location of the epicentral region. The model involving a subduction-related rupture in the Gulf of Cádiz results in excessive ground motion in northern Morocco, suggesting that the source of the 1755 earthquake should be located further west. A rupture along the western coast of Portugal, compatible with an activation of the passive western Iberian margin, would imply a relatively low average slip, which, alone, would could not account for the large tsunami observed in the whole northern Atlantic ocean. A seismic source located below the Gorringe Bank seems the most likely since it is more efficient in reproducing the distribution of high intensities in SW Iberia due to the 1755 earthquake.     

Uncertainty assessment of ephemeral gully identification,characteristics and topographic threshold when using aerial photographs in agricultural settings

Manual digitizing on aerial photographs is still commonly used for characterizing gully erosion over large areas. Even when automated detection procedures are implemented, manual digitizing is frequently being resorted to in order to constitute reference datasets used for training and validation. In both cases, manual digitizing entails some subjective decisions on behalf of the operator, which introduces uncertainty into the resulting datasets. To assess the magnitude of this uncertainty, 11 experienced operators were asked to digitize and classify ephemeral gullies (EGs) on cropland following a standardized methodology. The resulting 11 datasets were compared in terms of number, type and location of EGs. Furthermore, for EGs located on a well‐defined runoff flow concentration axis, the slope versus contributing area topographic thresholds required for initiating gully channels were assessed using four thresholding methods, and compared across the 11 datasets. The operators identified 259 different EGs. However, the number (52–139) and sum total length (8.9–23.7 km) of EGs varied widely across operators. Only 34% of the EGs were digitized by more than half of the operators, and 7% were identified by all. Identification of EGs located on a well‐defined flow concentration axis proved least subjective. The longer the EG and the more fields the EG crossed, the larger the number of operators that were able to identify it. EGs were also most easily identified when located in sugar beet fields as compared to other crops. EG classification and topographic threshold lines were also found to be strongly operator‐dependent. Quantile regression appeared to be one of the most robust thresholding methods. Operator subjectivity when digitizing EGs on orthophotographs introduces uncertainty that should be taken into account in future remote sensing‐based studies of EG erosion whenever they rely, in part or in full, on manual photograph interpretation. Copyright © 2014 John Wiley & Sons, Ltd.     

Extracting flow structures from tracer data

A method of visualizing structures in closed chaotic flows out of homogenous particle distributions is presented in the example of models of a meandering jet. To this end, the system will be leaked or opened up by defining a region of the flow, so that a particle is considered to be escaped if it leaves this region. By applying this method to an ensemble of nonescaped tracers, we are able to characterize mixing processes by visualizing the converging and stretching filamentations (stable and unstable manifolds) in the flow without using additional mathematical tools. The possibility of applying the algorithm to analyze buoy data, and a comparison with the finite time manifolds are discussed.     

An improved earthquake catalogue (<Emphasis Type="Italic">M</Emphasis> ≥ 4.0) for Turkey and its near vicinity (1900–2012)

Many catalogues, agency reports and research articles have been published on seismicity of Turkey and its surrounding since 1950s. Given existing magnitude heterogeneity, erroneous information on epicentral location, event date and time, this past published data however is far from fulfilling the required standards. Paucity of a standardized format in the available catalogues have reinforced the need for a refined and updated catalogue for earthquake related hazard and risk studies. During this study, ~37,000 earthquakes and related parametric data were evaluated by utilizing more than 41 published studies and, an integrated database was prepared in order to analyse all parameters acquired from the catalogues and references for each event. Within the scope of this study, the epicentral locations of M ≥ 5.0 events were firstly reappraised based on the updated Active Fault Map of Turkey. An improved catalogue of 12.674 events for the period 1900–2012 was as a result recompiled for the region between 32–45N° and 23–48E° by analyzing in detail accuracy of all seismological parameters available for each event. The events consist of M ≥ 4.0 are reported in several magnitude scales (e.g. moment magnitude, M_w; surface wave magnitude, M_S; body-wave magnitude m_b; local magnitude M_L and duration magnitude M_d) whereas the maximum focal depth reaches up to 225-km. In order to provide homogenous data, the improved catalogue is unified in terms of M_w. Fore-and aftershocks were also removed from the catalogue and completeness analyses were performed both separately for various tectonic sources and as a whole for the study region of interest. Thus, the prepared homogenous and declustered catalogue consisting of 6573 events provides the basis for a reliable input to the seismic hazard assessment studies for Turkey and its surrounding areas.     

A Note on Extension Variances in IR2

Matheron (1971) proposed an approximation of the extension variance in IR. We propose in this note an extension of this formula in IR ² , based on a MacLaurin formula. Its application is shown in an example, the estimation of the maximum depressional storage of a soil surface.     

Recent Climate Warming: Surface Air Temperature Series and Geothermal Evidence

Long-term (1961 – 1996) meteorological air temperature series together with the reconstructed ground surface temperature histories, obtained by inverting borehole temperature-depth profiles, were used to project regional patterns of the recent (climate) warming rate on the territory of the Czech Republic. The characteristic magnitude of the warming rate of 0.02 –0.03 K/yr was confirmed by the results of several years of monitoring the temperature in two experimental boreholes. The monitoring of shallow temperatures at depths of about 30 –40 m, i.e. below the reach of the seasonal surface temperature variations, can serve as an alternative tool of direct quantitative assessment of the present warming rate. The data also seem to sustain a potential man-made component contributing to the more pronounced recent warming rate observed in the areas of large agglomeration.     

Role of reservoir regulation and groundwater feedback in a simulated ground-soil-vegetation continuum: A long-term regional scale analysis

The regional terrestrial water cycle is strongly altered by human activities. Among them, reservoir regulation is a way to spatially and temporally allocate water resources in a basin for multi-purposes. However, it is still not sufficiently understood how reservoir regulation modifies the regional terrestrial- and subsequently, the atmospheric water cycle. To address this question, the representation of reservoir regulation into the terrestrial component of fully coupled regional Earth system models is required. In this study, an existing process-based reservoir network module is implemented into NOAH-HMS, that is, the terrestrial component of an atmospheric–hydrologic modelling system, namely, the WRF-HMS. It allows to quantitatively differentiate role of reservoir regulation and of groundwater feedback in a simulated ground-soil-vegetation continuum. Our study focuses on the Poyang Lake basin, where the largest freshwater lake of China and reservoirs of different sizes are located. As compared to streamflow observations, the newly extended NOAH-HMS slightly improves the streamflow and streamflow duration curves simulation for the Poyang Lake basin for the period 1979–1986. The inclusion of reservoir regulation leads to major changes in the simulated groundwater recharges and evaporation from reservoirs at local scale, but has minor effects on the simulated soil moisture and surface runoff at basin scale. The performed groundwater feedback sensitivity analysis shows that the strength of the groundwater feedback is not altered by the consideration of reservoir regulation. Furthermore, both reservoir regulation and groundwater feedback modify the partitioning of the simulated evapotranspiration, thus affecting the atmospheric water cycle in the Poyang Lake region. This finding motivates future research with our extended fully coupled atmospheric–hydrologic modelling system by the community.     

生物质燃烧颗粒物有机示踪化合物的测定和应用

使用气相色谱-质谱法对2002年7月-2003年7月的北京市大气粗、细颗粒物样品中的左旋葡聚糖等糖类化合物进行了测定。结果表明，左旋葡聚糖主要存在于细颗粒中，可以作为示踪化合物来研究生物质燃烧现象。生物质燃烧对北京市大气颗粒物有较重要的贡献。对应于12％-40％的PM2.5有机碳和10％-33％的PM10有机碳。北京市在2002年10月和11月受到明显的生物质燃烧的影响，可能由于农田秸杆焚烧和秋季落叶焚烧。2003年5月7日颗粒物样品受到直线距离约为1000km以外的内蒙古自治区呼伦贝尔市大兴安岭林区森林大火烟雾的影响。生物质燃烧事件具有突发性，可以长距离传输；生物质作为农村生物燃料的使用其燃烧排放具有经常性和持久性的特征。