Data assimilation: From photon counts to Earth System forecasts

Data assimilation – the set of techniques whereby information from observing systems and models is combined optimally – is rapidly becoming prominent in endeavours to exploit Earth Observation for Earth sciences, including climate prediction. This paper explains the broad principles of data assimilation, outlining different approaches (optimal interpolation, three-dimensional and four-dimensional variational methods, the Kalman Filter), together with the approximations that are often necessary to make them practicable. After pointing out a variety of benefits of data assimilation, the paper then outlines some practical applications of the exploitation of Earth Observation by data assimilation in the areas of operational oceanography, chemical weather forecasting and carbon cycle modelling. Finally, some challenges for the future are noted.

Subseasonal extremes of precipitation and active-break cycles of the Indian summer monsoon in a climate-change scenario

Changes to the behaviour of subseasonal precipitation extremes and active-break cycles of the Indian summer monsoon are assessed in this study using pre-industrial and 2 × CO2 integrations of the Hadley Centre coupled model HadCM3, which is able to simulate the monsoon seasonal cycle reasonably. At 2 × CO2, mean summer rainfall increases slightly, especially over central and northern India. The mean intensity of daily precipitation during the monsoon is found to increase, consistent with fewer wet days, and there are increases to heavy rain events beyond changes in the mean alone. The chance of reaching particular thresholds of heavy rainfall is found to approximately double over northern India, increasing the likelihood of damaging floods on a seasonal basis. The local distribution of such projections is uncertain, however, given the large spread in mean monsoon rainfall change and associated extremes amongst even the most recent coupled climate models. The measured increase of the heaviest precipitation events over India is found to be broadly in line with the degree of atmospheric warming and associated increases in specific humidity, lending a degree of predictability to changes in rainfall extremes. Active-break cycles of the Indian summer monsoon, important particularly due to their effect on agricultural output, are shown to be reasonably represented in HadCM3, in particular with some degree of northward propagation. We note an intensification of both active and break events, particularly when measured against the annual cycle, although there is no suggestion of any change to the duration or likelihood of monsoon breaks. Copyright © 2009 Royal Meteorological Society

Fast regridding of large, complex geospatial datasets

In the earth sciences, data are commonly cast on complex grids in order to model irregular domains such as coastlines, or to evenly distribute grid points over the globe. It is common for a scientist to wish to re-cast such data onto a grid that is more amenable to manipulation, visualization, or comparison with other data sources. The complexity of the grids presents a significant technical difficulty to the regridding process. In particular, the regridding of complex grids may suffer from severe performance issues, in the worst case scaling with the product of the sizes of the source and destination grids. We present a mechanism for the fast regridding of such datasets, based upon the construction of a spatial index that allows fast searching of the source grid. We discover that the most efficient spatial index under test (in terms of memory usage and query time) is a simple look-up table. A kd-tree implementation was found to be faster to build and to give similar query performance at the expense of a larger memory footprint. Using our approach, we demonstrate that regridding of complex data may proceed at speeds sufficient to permit regridding on-the-fly in an interactive visualization application, or in a Web Map Service implementation. For large datasets with complex grids the new mechanism is shown to significantly outperform algorithms used in many scientific visualization packages.

Basement depths in the Parecis Basin (Amazon) with receiver functions from small local earthquakes in the Porto dos Gauchos seismic zone

Receiver functions from small local earthquakes were used to determine sediment thicknesses in Porto dos Gauchos seismic zone (PGSZ), Parecis basin, Amazonian craton, Brazil. The high velocity contrast between basement and sediments (P-wave velocities of 6.1 and 3.0 km/s, respectively) favors the generation of clear P-to-S converted phases (Ps) seen in the radial component, and also S-to-P conversions (Sp) seen in the vertical component. A reference 10 velocity model determined with shallow refraction experiment in PGSZ helped to convert Ps P time differences to basement depths at 15 stations deployed for aftershocks studies. The results of receiver function integrated with the shallow refraction reveal that the basement depths in the PGSZ increases from the basin border in the north up to about 600 m depth in the south. The basement topography, however, does not vary smoothly and a basement high with a steep topography was detected near the epicentral area. A 400 m elevation difference within 1.7 km distance suggests a possible border fault of a buried graben. This feature seems to be oriented roughly WSW-ENE and could indicate basement structures related to the seismicity of the Porto dos Gauchos Seismic Zone. (C) 2011 Elsevier Ltd. All rights reserved.

Major shear zones of southern Brazil and Uruguay: escape tectonics in the eastern border of Rio de La plata and Paranapanema cratons during the Western Gondwana amalgamation

The Mantiqueira Province represents a series of supracrustal segments of the South-American counterpart formed during the Gondwana Supercontinent agglutination. In this crustal domain, the process of escape tectonics played a conspicuous role, generating important NE-N-S-trending lineaments. The oblique component of the motions of the colliding tectonic blocks defined the transpressional character of the main suture zones: Lancinha-Itariri, Cubato-Arcadia-Areal, Serrinha-Rio Palmital in the Ribeira Belt and Sierra Ballena-Major Gercino in the Dom Feliciano Belt. The process as a whole lasted for ca. 60 Ma, since the initial collision phase until the lateral escape phase predominantly marked by dextral and subordinate sinistral transpressional shear zones. In the Dom Feliciano Belt, southern Brazil and Uruguay, transpressional event at 630-600 Ma is recognized and in the Ribeira Belt, despite less coevally, the transpressional event occurred between 590 and 560 Ma in its northern-central portion and between ca. 625 and 595 Ma in its central-southern portion. The kinematics of several shear zones with simultaneous movement in opposite directions at their terminations is explained by the sinuosity of these lineaments in relation to a predominantly continuous westward compression.

Evaluation of petrogenetic models for intermediate and silicic plutonic rocks from the Sierra de Valle Fertil-La Huerta, Argentina: Petrologic constraints on the origin of igneous rocks in the Ordovician Famatinian-Puna paleoarc

The whole Valle Fertil-La Huerta section appears as a calc-alkaline plutonic suite typical of a destructive plate margin. New Sr and Nd isotopic whole-rock data and published whole-rock geochemistry suggest that the less-evolved intermediate (dioritic) rocks can be derived by magmatic differentiation, mainly by hornblende + plagioclase +/- Fe-Ti oxide fractional crystallization, from mafic (gabbroic) igneous precursors. Closed-system differentiation, however, cannot produce the typical intermediate (tonalitic) and silicic (granodioritic) plutonic rocks, which requires a preponderant contribution of crustal components. Intermediate and silicic plutonic rocks from Valle Fertil-La Huerta section have formed in a plate subduction setting where the thermal and material input of mantle-derived magmas promoted fusion of fertile metasedimentary rocks and favored mixing of gabbroic or dioritic magmas with crustal granitic melts. Magma mixing is observable in the field and evident in variations of chemical elemental parameters and isotopic ratios, revealing that hybridization coupled with fractionation of magmas took place in the crust. Consideration of the whole-rock geochemical and isotopic data in the context of the Famatinian-Puna magmatic belt as a whole demonstrates that the petrologic model postulated for the Sierra Valle Fertil-La Huerta section has the potential to explain the generation of plutonic and volcanic rocks across the Early Ordovician paleoarc from central and northwestern Argentina. As the petrologic model does not require the intervention of old Precambrian continental crust, the nature of the basement on which thick accretionary turbiditic sequences were deposited remains a puzzling aspect. Discussion in this paper provides insights into the nature of magmatic source rocks and mechanisms of magma generation in Cordilleran-type volcano-plutonic arcs of destructive plate margins. (C) 2010 Elsevier Ltd. All rights reserved.

The Guarguaraz Complex and the Neoproterozoic-Cambrian evolution of southwestern Gondwana: Geochemical signatures and geochronological constraints

The Guarguaraz Complex, in western Argentina, comprises a metasedimentary assemblage, associated with mafic sills and ultramafic bodies intruded by basaltic dikes, which are interpreted as Ordovician dismembered ophiolites. Two kinds of dikes are recognized, a group associated with the metasediments and the other ophiolite-related. Both have N-MORB signatures, with epsilon(Nd) between +3.5 and +8.2, indicating a depleted source, and Grenville model ages between 0.99 and 1.62 Ga. A whole-rock Sm-Nd isochron yielded an age of 655 +/- 76 Ma for these mafic rocks, which is compatible with cianobacteria and acritarchae recognized in the clastic metasedimentary platform sequences, that indicate a Neoproterozoic (Vendian)-Cambrian age of deposition. The Guarguaraz metasedimentary-ophiolitic complex represents, therefore, a remnant of an oceanic basin developed to the west of the Grenville-aged Cuyania terrane during the Neoproterozoic. The southernmost extension of these metasedimentary sequences in Cordon del Portillo might represent part of this platform and not fragments of the Chilenia terrane. An extensional event related to the fragmentation of Rodinia is represented by the mafic and ultramafic rocks. The Devonian docking of Chilenia emplaced remnants of ocean floor and slices of the Cuyania terrane (Las Yaretas Gneisses) in tectonic contact with the Neoproterozoic metasediments, marking the Devonian western border of Gondwana. (C) 2009 Elsevier Ltd. All rights reserved.