Co-registration and correlation of aerial photographs for ground deformation measurements

We describe and test a procedure to accurately co-register and correlate multi-temporal aerial images. We show that this procedure can be used to measure surface deformation, and explore the performance and limitations of the technique. The algorithms were implemented in a software package, COSI-Corr (available from the Caltech Tectonics Observatory website). The technique is validated on several case examples of co-seismic deformation. First, we measure co-seismic ground deformation due to the 1992, Mw 7.3, Landers, California, earthquake from 1 m resolution aerial photography of the National Aerial Photography Program (United States Geological Survey). The fault ruptures are clearly detected, including small kilometric segments with fault slip as small as a few tens of centimeters. We also obtained similar performance from images of the fault ruptures produced by the 1999 Mw 7.1 Hector Mine, California, earthquake. The measurements are shown to be biased due to the inaccuracy of the Digital Elevation Model, film distortions, scanning artifacts, and ignorance of ground displacements at the location of the tie points used to co-register the multi-temporal images. We show that some of these artifacts can be identified and corrected.     

A novel platform to evaluate the dampening of water and solute transport in an experimental channel under unsteady flow conditions

This paper presents a novel platform to study the dampening of water and solute transport in an experimental channel under unsteady flow conditions, where literature data are scarce. We address the question about what could be the smallest size of experimental platform that is useful for research, project studies, and teaching activities and that allows to do rational experiments characterized by small space occupation, short experimental duration, high measurement precision, high quality and reproducible experimental curves, low water and energy consumption, and the possibility to test a large variety of hydrograph scenarios. Whereas large scale hydraulic laboratories have focused their studies on sediment transport, our platform deals with solute transport. The objectives of our study are (a) building a platform that allows to do rational experiments, (b) enriching the lack of experimental data concerning water and solute transport under unsteady state conditions, and (c) studying the dampening of water and solute transport. We studied solute transport in a channel with lateral gain and lateral loss under different experimental configurations, and we show how the same lateral loss flow event can lead to different lateral loss mass repartitions under different configurations. In order to characterize water and solute dampening between the input and the output of the channel, we calculate dampening ratios based on peak coordinates of time flow curves and time mass curves and that express the decrease of peak amplitude and the increase of peak occurrence time between the input and output curves. Finally, we use a solute transport model coupling the diffusive wave equation for water transfer and the advection–diffusion equation for solute transport in order to simulate the experimental data. The simulations are quite good with a Nash–Sutcliffe efficiency NSE > 0.98 for water transfer and 0.84 < NSE < 0.97 for solute transport. This platform could serve hydrological modellers because it offers a variety of measured parameters (flow, water height, and solute concentration), at a fine time step under unsteady flow conditions.     

Sensitivity to domain size of mid-latitude summer simulations with a regional climate model

The issue of Regional Climate Model (RCM) domain size is studied here by using a perfect-model approach, also known as the Big-Brother experiment. It is known that the control exerted by the lateral boundary conditions (LBC) on nested simulations increases when reducing the domain size. The large-scale component of the simulation that is forced by the LBC influences the small-scale features that develop along the large-scale flow. Small-scale transient eddies need space and time to develop sufficiently however, and small domains can impede their development. Our tests performed over eastern North America in summer reveal that the small-scale features are systematically underestimated over the entire domain, even for domain as large as 140 by 140 grid points. This result differs from that obtained in winter where the small scales were mainly underestimated on the west (inflow) side of the domain. This difference is due to the circulation regime over Eastern Canada, which is characterized by weak and variable flow in summer, but strong and westerly flow in winter. For both seasons, the small-scale transient-eddy amplitudes are systematically underestimated at higher levels, but this problem is less severe in summer. Overall the model is more skilful in regenerating the small scales in summer than in winter for comparable domain sizes, which can be related to the weaker summer flow and stronger physical processes occurring in this season.     

东日本大地震启示

一组高质量的数据记录了2011年3月东日本大地震的震前、地震发生时以及震后的地表形变。但目前还没有适当的模型可以利用类似数据来对地震和海啸危险性进行评估。随着基于空间的测量板块构造边界和大型断层表面形变技术的发展,地震科学已经跨入一个全新的时代。日本在实施这些技术方面一直走在前列,特别是随着大约15年     

Chironomid assemblages in cores from multiple water depths reflect oxygen-driven changes in a deep French lake over the last 150 years

We sampled modern chironomids at multiple water depths in Lake Annecy, France, before reconstructing changes in chironomid assemblages at sub-decadal resolution in sediment cores spanning the last 150 years. The lake is a large, deep (z_max = 65 m), subalpine waterbody that has recently returned to an oligotrophic state. Comparison between the water-depth distributions of living chironomid larvae and subfossil head capsules (HC) along three surface-sediment transects indicated spatial differences in the influence of external forcings on HC deposition (e.g. tributary effects). The transect with the lowest littoral influence and the best-preserved, depth-specific chironomid community characteristics was used for paleolimnological reconstructions at various water depths. At the beginning of the twentieth century, oxygen-rich conditions prevailed in the lake, as inferred from M. contracta-type and Procladius sp. at deep-water sites (i.e. cores from 56 to 65 m) and Paracladius sp. and H. grimshawi-type in the core from 30 m depth. Over time, chironomid assemblages in cores from all three water depths converged toward the dominance of S. coracina-type, indicating enhanced hypoxia. The initial change in chironomid assemblages from the deep-water cores occurred in the 1930s, at the same time that an increase in lake trophic state is inferred from an increase in total organic carbon (TOC) concentration in the sediment. In the 1950s, an assemblage change in the core from 30 m water depth reflects the rapid expansion of the hypoxic layer into the shallower region of the lake. Lake Annecy recovered its oligotrophic state in the 1990s. Chironomid assemblages, however, still indicate hypoxic conditions, suggesting that modern chironomid assemblages in Lake Annecy are decoupled from the lake trophic state. Recent increases in both TOC and the hydrogen index indicate that changes in pelagic functioning have had a strong indirect influence on the composition of the chironomid assemblage. Finally, the dramatic decrease in HC accumulation rate over time suggests that hypoxic conditions are maintained through a feedback loop, wherein the accumulation of (un-consumed) organic matter and subsequent bacterial respiration prevent chironomid re-colonization. We recommend study of sediment cores from multiple water depths, as opposed to investigation of only a single core from the deepest part of the lake, to assess the details of past ecological changes in large deep lakes.     

Teleseismic waveform modelling with a one-way wave equation

It is now widely accepted that elastic properties of the continental lithosphere and the underlying sublithospheric mantle are both anisotropic and laterally heterogeneous at a range of scales. To fully exploit modern three-component broad-band array data sets requires the use of comprehensive modelling tools. In this work, we investigate the use of a wide-angle, one-way wave equation to model variations in teleseismic 3-D waveforms due to 2-D elastic heterogeneity and anisotropy. The one-way operators are derived based on a high-frequency approximation of the square-root operator and include the effects of wave propagation as well as multiple scattering. Computational cost is reduced through a number of physically motivated approximations. We present synthetic results from simple 1-D (layer over a half-space) and 2-D (subduction zone) models that are compared with reference solutions. The algorithm is then used to model data from an array of broad-band seismograph stations deployed in northwestern Canada as part of the IRIS-PASSCAL/LITHOPROBE CANOE experiment. In this region radial-component receiver functions show a clear continental Moho and the presence of crustal material dipping into the mantle at the suture of two Palaeo-Proterozoic terranes. The geometry of the suture is better defined on the transverse component where subduction is associated with a ∼10 km thick layer exhibiting strong elastic anisotropy. The modelling reproduces the main features of the receiver functions, including the effects of anisotropy, heterogeneity and finite-frequency scattering.     

High-end aerial digital cameras and their impact on the automation and quality of the production workflow

The IGN digital camera project was established in the early 1990s. The first research surveys were carried out in 1996 and the digital camera was first used in production in 2000. In 2004 approximately 10 French departments (accounting for 10% of the territory) were covered with a four-head camera system and since summer 2005 all IGN imagery has been acquired digitally. Nevertheless the camera system is still evolving, with tests on new geometric configurations being continuously undertaken. The progressive integration of the system in IGN production workflows has allowed IGN to keep the system evolving in accordance with production needs. Remaining problems are due to specific camera characteristics such as CCD format, the optical quality of off-the-shelf lenses, and because some production tools are ill-adapted to digital images with a large dynamic range. However, when considering the pros and cons of integrating these images into production lines, the disadvantages are largely balanced by the numerous benefits this technology offers.     

The coupling of optimal economic growth and climate dynamics

In this paper, we study optimal economic growth programs coupled with climate change dynamics. The study is based on models derived from MERGE, a well established integrated assessment model (IAM). We discuss first the introduction in MERGE of a set of “tolerable window” constraints which limit both the temperature change and the rate of temperature change. These constraints, obtained from ensemble simulations performed with the Bern 2.5-D climate model, allow us to identity a domain intended to preserve the Atlantic thermohaline circulation. Next, we report on experiments where a two-way coupling is realized between the economic module of MERGE and an intermediate complexity “3-D-” climate model (C-GOLDSTEIN) which computes the changes in climate and mean temperature. The coupling is achieved through the implementation of an advanced “oracle based optimization technique” which permits the integration of information coming from the climate model during the search for the optimal economic growth path. Both cost-effectiveness and cost-benefit analysis modes are explored with this combined “meta-model” which we refer to as GOLDMERGE. Some perspectives on future implementations of these approaches in the context of “collaborative” or “community” integrated assessment modules are derived from the comparison of the different approaches.     

Projection of future climate change conditions using IPCC simulations, neural networks and Bayesian statistics. Part 1: Temperature mean state and seasonal cycle in South America

Projections for South America of future climate change conditions in mean state and seasonal cycle for temperature during the twenty-first century are discussed. Our analysis includes one simulation of seven Atmospheric-Ocean Global Circulation Models, which participated in the Intergovernmental Panel on Climate Change Project and provided at least one simulation for the twentieth century (20c3m) and one simulation for each of three Special Report on Emissions Scenarios (SRES) A2, A1B, and B1. We developed a statistical method based on neural networks and Bayesian statistics to evaluate the models’ skills in simulating late twentieth century temperature over continental areas. Some criteria [model weight indices (MWIs)] are computed allowing comparing over such large regions how each model captures the temperature large scale structures and contributes to the multi-model combination. As the study demonstrates, the use of neural networks, optimized by Bayesian statistics, leads to two major results. First, the MWIs can be interpreted as optimal weights for a linear combination of the climate models. Second, the comparison between the neural network projection of twenty-first century conditions and a linear combination of such conditions allows the identification of the regions, which will most probably change, according to model biases and model ensemble variance. Model simulations in the southern tip of South America and along the Chilean and Peruvian coasts or in the northern coasts of South America (Venezuela, Guiana) are particularly poor. Overall, our results present an upper bound of potential temperature warming for each scenario. Spatially, in SRES A2, our major findings are that Tropical South America could warm up by about 4°C, while southern South America (SSA) would also undergo a near 2–3°C average warming. Interestingly, this annual mean temperature trend is modulated by the seasonal cycle in a contrasted way according to the regions. In SSA, the amplitude of the seasonal cycle tends to increase, while in northern South America, the amplitude of the seasonal cycle would be reduced leading to much milder winters. We show that all the scenarios have similar patterns and only differ in amplitude. SRES A1B differ from SRES A2 mainly for the late twenty-first century, reaching more or less an 80–90% amplitude compared to SRES A2. SRES B1, however, diverges from the other scenarios as soon as 2025. For the late twenty-first century, SRES B1 displays amplitudes, which are about half those of SRES A2.