ALBERT

Description: 〈b〉Paleomagnetic determination of paleolatitude and rotation of Bering Island (Komandorsky Islands) Russia: comparison with rotations in the Aleutian Islands and Kamchatka〈/b〉〈br〉 P. S. Minyuk and D. B. Stone〈br〉 Stephan Mueller Spec. Publ. Ser., 4, 329-348, https://doi.org/10.5194/smsps-4-329-2009, 2009〈br〉 A paleomagnetic study was carried out on Paleogene sedimentary rocks from Bering Island, Komandorsky islands, located at the far western end of the Aleutian Island Arc. The age of these sediments has been debated at length, but the combination of magnetostratigraphy with the fossil record indicates that the base of the section is of early Eocene (approximately 55 Ma) and the top latest Eocene age. Paleomagnetic data were obtained from 260 samples from 60 individual bedding units. The combined data show a clockwise rotation 〈i〉R〈/i〉=26.3°±8.5°, 〈i〉F〈/i〉=8.1°±2.5° with respect to the North American Plate and 〈i〉R〈/i〉=38°±8.8°, 〈i〉F〈/i〉=8.7°±2.7° with respect to the Eurasian Plate. They also show a shallowing of the inclination which yields a paleolatitude of 53°, 12° south of its expected latitude. The shallowing may have a component due to compaction, but the wide variation in sampled lithologies, combined with internal consistency of the data set, would argue against the shallowing being significant. To compare these data with other Aleutian Arc data we compiled a comprehensive survey of all available data sets. Out of these we selected four islands for which the data passed basic reliability criteria, namely Umnak, Amlia, Amchitka and Medny islands. All four showed significant clockwise rotation with respect to both North American and Eurasian polar wander paths. Several mechanisms can generate the observed rotation, ranging from block rotation driven by oblique relative motion of the Pacific plate, through lateral transport along the curve of the arc, to whole-arc rotation about its eastern end. The distribution and age spread of the rotation data are insufficient to discriminate between mechanisms, but it seems likely that different mechanism may have operated at different times and in different locations.

Description: 〈b〉Paleomagnetism of the Cretaceous rocks from Cape Kronotskiy, East Kamchatka and reconstruction of terrane trajectories in the NE Pacific area〈/b〉〈br〉 W. Harbert, N. V. Tsukanov, D. V. Alexeiev, C. Gaedicke, R. Freitag, B. V. Baranov, S. G. Skolotnev, W. Kramer, and W. Seifert〈br〉 Stephan Mueller Spec. Publ. Ser., 4, 313-327, https://doi.org/10.5194/smsps-4-313-2009, 2009〈br〉 The Kamchatka Peninsula of northeastern Russia is located along the northwestern margin of the Bering Sea and consists of zones of complexly deformed accreted terranes. Paleomagnetic samples were collected for study from a Late Cretaceous aged locality at Cape Kronotskiy (λ=54°44.8´ N, φ=162°1.29´ E). Two components of magnetization were observed. During stepwise thermal demagnetization, the B-magnetic component was observed up to 600°C having a direction and associated uncertainty in stratigraphic coordinates of 〈i〉D〈sub〉s〈/sub〉〈/i〉=300.7°, 〈i〉I〈sub〉s〈/sub〉〈/i〉=48.7°, α〈sub〉95〈/sub〉=10.9°, k-value=11.8, n=17. The B component paleolatitude calculated from the Fisher mean in stratigraphic coordinates and associated statistics are λ〈sub〉obs〈/sub〉=30.4° N or S, λ〈sub〉95〈/sub〉=8.9°, n=17 (sites), k-value=11. Our overall study paleolatitude result is similar to a previously reported paleomagnetic study completed within this unit. Terrane trajectories calculated using the finite rotation poles of Engebretson et al. (1985), which are corrected for either Pacific-hotspot drift or True Polar Wander hotspot-spin axis relative motion, show that the sampled unit represents a far traveled tectonostratigraphic terrane and support a model in which accretion (docking) events of this composite or superterrane with the North America plate occur at approximately 40 Ma.

Description: 〈b〉Comparison of Cretaceous granitoids of the Chaun tectonic zone to those of the Taigonos Peninsula, NE Asia: rock chemistry, composition of rock forming minerals, and conditions of formation〈/b〉〈br〉 P. L. Tikhomirov, M. V. Luchitskaya, and I. R. Kravchenko-Berezhnoy〈br〉 Stephan Mueller Spec. Publ. Ser., 4, 289-311, https://doi.org/10.5194/smsps-4-289-2009, 2009〈br〉 The Cretaceous granitoid complexes of the Eastern Taigonos and the Prybrezhny Taigonos belts (southern part of the Taigonos Peninsula), Tanyurer pluton of the Okhotsk-Chukotka volcanic belt, and the Peekiney, Moltykan, and Telekay plutons of the Chaun tectonic zone are discussed in relation to their structural position, petrography, rock and mineral chemistry and physicochemical conditions of melt crystallization. These granitoid plutons were generated through melting of a compositionally heterogeneous crustal source, with direct contribution from mafic melts produced in the mantle wedge above active or extinct Benioff zones. Variations of the trace-element composition of granitoids are controlled to a greater extent by local compositional peculiarities of the source regions than by the geodynamic regime as such. The final crystallization of these plutons occurred at comparatively shallow depths, between 1–2 and 6–7 km, in a temperature interval of 700–770°C. The depth of emplacement of the bodies decreases with increasing distance from the areas with oceanic and transitional type crust, as does the degree of incompatible element enrichment of the mantle and crustal sources of melts. Variations in fo〈sub〉2〈/sub〉 values at the late stages of crystallization of the plutons reach 3–4 orders of magnitude, exceeding the limits of the quartz-fayalite-magnetite (QFM) and nickel-nickel oxide (NNO) buffer equilibria, which likely results from local variations of the source composition.

Description: 〈b〉Thermogeodynamic manifestations in the Caucasus and their genesis〈/b〉〈br〉 G. E. Gugunava, J. K. Kiria, and T. B. Bochorishvili〈br〉 eEarth Discuss., 4, 77-89, doi:10.5194/eed-4-77-2009, 2009〈br〉 〈b〉Publication in eE not foreseen〈/b〉 (discussion: closed, 2 comments)〈br〉 In the work two aspects of thermal character are considered: first of all this is the connection of subduction phenomena with thermal life of the Caucasus on the basis of over interpreted data of magnetotelluric sounding, and secondly, origin of thermostressed condition of the Caucasus and its geological aspects which is manifested in the following: 〈br〉〈br〉 1. in the zones of anomalous thermodisplacements thermofaults should occur (Le Pishon et al., 1977). These thermofaults are in good correlation with deep faults which are distinguished by geological and seismic methods, these thermofaults may be earthquake sources (Spitak, Racha, etc. earthquakes), also may be channels through which magma derivates (giving mineral deposits) may penetrate on surface (Gugunava and Gijeishvili, 1989); 〈br〉〈br〉 2. in the body of sedimentary complex thermostressed seals and seal failures occur, which are apparently traps for oil-gas fluids. Good correlation of thermodense anomalies with oil deposits of the Caucasus is shown (Alexidze et al., 1985; Gugunava, 1980). 〈br〉〈br〉 Everything above mentioned was carried out within frames of stationary thermal model which did not allow us to reveal time characteristics of interconnection of geological medium and thermal field. 〈br〉〈br〉 Now investigations are being carried out within the frames of stationary thermal model and its interconnection with geological environment.

Description: 〈b〉Thermogeodynamic manifestations in the Caucasus and their genesis〈/b〉〈br〉 G. E. Gugunava, J. K. Kiria, and T. B. Bochorishvili〈br〉 eEarth Discuss., 4, 77-89, doi:10.5194/eed-4-77-2009, 2009〈br〉 〈b〉Publication in eE not foreseen〈/b〉 (discussion: closed, 2 comments)〈br〉 In the work two aspects of thermal character are considered: first of all this is the connection of subduction phenomena with thermal life of the Caucasus on the basis of over interpreted data of magnetotelluric sounding, and secondly, origin of thermostressed condition of the Caucasus and its geological aspects which is manifested in the following: 〈br〉〈br〉 1. in the zones of anomalous thermodisplacements thermofaults should occur (Le Pishon et al., 1977). These thermofaults are in good correlation with deep faults which are distinguished by geological and seismic methods, these thermofaults may be earthquake sources (Spitak, Racha, etc. earthquakes), also may be channels through which magma derivates (giving mineral deposits) may penetrate on surface (Gugunava and Gijeishvili, 1989); 〈br〉〈br〉 2. in the body of sedimentary complex thermostressed seals and seal failures occur, which are apparently traps for oil-gas fluids. Good correlation of thermodense anomalies with oil deposits of the Caucasus is shown (Alexidze et al., 1985; Gugunava, 1980). 〈br〉〈br〉 Everything above mentioned was carried out within frames of stationary thermal model which did not allow us to reveal time characteristics of interconnection of geological medium and thermal field. 〈br〉〈br〉 Now investigations are being carried out within the frames of stationary thermal model and its interconnection with geological environment.

Description: 〈b〉Holocene evolution and sedimentation rate of Alikes Lagoon, Zakynthos island, Western Greece: preliminary results〈/b〉〈br〉 P. Avramidis and N. Kontopoulos〈br〉 eEarth, 4, 23-29, doi:10.5194/ee-4-23-2009, 2009〈br〉 In the present study we present preliminary results from Alikes lagoon in Zakynthos island, an area that is one of the most seismically active regions of Greece. In order to estimate – interpret the Holocene evolution of the area and to reconstruct the palaeoenvironmental changes, we based on data derived from a 21 m sediment core. Sediment types, structure, colour, as well as contact depths and bed characteristics were recorded in the field. Standarised sedimentological analysis was carried out, on 46 samples including grain size analysis, calculation of moment measures, and micro- and molluscan fossils of 17 selected samples. Moreover, radiocarbon age determinations have been made on individual 〈i〉Cardium〈/i〉 shells from two horizons and whole – core Magnetic Susceptibility (MS) measurements were taken. The interpretation of depositional environments suggests a coastal environment (restricted-shallow) with reduced salinity such as a lagoon margin and in a tidal flat and/or marsh particularly. The maximum age of the studied sediments is about 8500 BP. The rate of sedimentation between 8280 BP while 5590 BP was 5.3 mm/yr and between 5590 BP and modern times is on the order of 1.03 mm/yr. These sedimentation rates results are similar to other coastal areas of western Greece.

Description: 〈b〉Morphology of the pore space in claystones – evidence from BIB/FIB ion beam sectioning and cryo-SEM observations〈/b〉〈br〉 G. Desbois, J. L. Urai, and P. A. Kukla〈br〉 eEarth, 4, 15-22, doi:10.5194/ee-4-15-2009, 2009〈br〉 The morphology of pore space has a strong effect on mechanical and transport properties of mudrocks and clay-rich fault gouge, but its characterization has been mostly indirect. We report on a study of Boom clay from a proposed disposal site of radioactive waste (Mol site, Belgium) using high resolution SEM at cryogenic temperature, with ion beam cross-sectioning to prepare smooth, damage free surfaces. Pores commonly have crack-like tips, preferred orientation parallel to bedding and power law size distribution. We define a number of pore types depending on shape and location in the microstructure: large jagged pores in strain shadows of clastic grains, high aspect ratio pores between similarly oriented phyllosilicate grains and crescent-shaped pores in saddle reefs of folded phyllosilicates. 3-D reconstruction by serial cross-sectioning shows 3-D connectivity of the pore space. These findings offer a new insight into the morphology of pores down to nano-scale in comparison to traditional pore size distributions calculated from mercury Injection experiments, explain slaking of clays by successive wetting and drying and provide the basis for microstructure-based models of transport in clays.

Description: 〈b〉Holocene evolution and sedimentation rate of Alikes Lagoon, Zakynthos island, Western Greece: preliminary results〈/b〉〈br〉 P. Avramidis and N. Kontopoulos〈br〉 eEarth, 4, 23-29, doi:10.5194/ee-4-23-2009, 2009〈br〉 In the present study we present preliminary results from Alikes lagoon in Zakynthos island, an area that is one of the most seismically active regions of Greece. In order to estimate – interpret the Holocene evolution of the area and to reconstruct the palaeoenvironmental changes, we based on data derived from a 21 m sediment core. Sediment types, structure, colour, as well as contact depths and bed characteristics were recorded in the field. Standarised sedimentological analysis was carried out, on 46 samples including grain size analysis, calculation of moment measures, and micro- and molluscan fossils of 17 selected samples. Moreover, radiocarbon age determinations have been made on individual 〈i〉Cardium〈/i〉 shells from two horizons and whole – core Magnetic Susceptibility (MS) measurements were taken. The interpretation of depositional environments suggests a coastal environment (restricted-shallow) with reduced salinity such as a lagoon margin and in a tidal flat and/or marsh particularly. The maximum age of the studied sediments is about 8500 BP. The rate of sedimentation between 8280 BP while 5590 BP was 5.3 mm/yr and between 5590 BP and modern times is on the order of 1.03 mm/yr. These sedimentation rates results are similar to other coastal areas of western Greece.

Description: 〈b〉Map – geographic locations related to the manuscript sequence〈/b〉〈br〉 〈br〉 Stephan Mueller Spec. Publ. Ser., 4, 3-3, https://doi.org/10.5194/smsps-4-3-2009, 2009〈br〉

Description: 〈b〉Age and paleomagnetism of the Okhotsk-Chukotka Volcanic Belt (OCVB) near Lake El'gygytgyn, Chukotka, Russia〈/b〉〈br〉 D. B. Stone, P. W. Layer, and M. I. Raikevich〈br〉 Stephan Mueller Spec. Publ. Ser., 4, 243-260, https://doi.org/10.5194/smsps-4-243-2009, 2009〈br〉 Paleomagnetic results from the upper two thirds of the whole section of the Okhotsk-Chukotka Volcanic Belt (OCVB) volcanics exposed in the area around Lake El'gygytgyn, Chukotka yield stable, consistent magnetic vectors and well-preserved reversed directions. The magnetostratigraphy and 〈sup〉40〈/sup〉Ar/〈sup〉39〈/sup〉Ar geochronologic data reported here indicate that the sampled OCVB volcanics were erupted between about 90 and 67 Ma, and show no significant change in the apparent pole position over that time. The OCVB extends from northeast China, across Northeast Russia to the Bering Straight. This belt is made up of both extrusive and intrusive rocks, with the extrusive rocks and their associated sediments being dominant. The whole belt important in interpreting the paleogeography of the region because it overlies many of the accreted terranes of Northeast Russia. Most importantly, it overlies parts of the Chukotka-Alaska block which is thought to have moved out of the Arctic Ocean region, as well as terranes accreted from the south. These latter terranes have been rafted northwards on the paleo-plates of the Pacific, implying that the present relative paleogeography of all of the terranes overlain by the OCVB were essentially in place by 67 Ma, and possibly as early as 90 Ma. However, comparing our paleomagnetic pole position for the OCVB with those for North America and Eurasia (a proxy for Siberia) shows a statistically significant displacement of the OCVB pole to the south west. This implies that not only the OCVB, but the underlying terranes of northeast Russia, experienced southerly displacement with respect to the Siberian and North American platforms since the Late Cretaceous.

Description: 〈b〉Tectonic reconstruction of Uda-Murgal arc and the Late Jurassic and Early Cretaceous convergent margin of Northeast Asia–Northwest Pacific〈/b〉〈br〉 S. D. Sokolov, G. Ye. Bondarenko, A. K. Khudoley, O. L. Morozov, M. V. Luchitskaya, M. I. Tuchkova, and P. W. Layer〈br〉 Stephan Mueller Spec. Publ. Ser., 4, 273-288, https://doi.org/10.5194/smsps-4-273-2009, 2009〈br〉 A long tectonic zone composed of Upper Jurassic to Lower Cretaceous volcanic and sedimentary rocks is recognized along the Asian continent margin from the Mongol-Okhotsk fold and thrust belt on the south to the Chukotka Peninsula on the north. This belt represents the Uda-Murgal arc, which was developed along the convergent margin between Northeast Asia and Northwest Meso-Pacific. Several segments are identified in this arc based upon the volcanic and sedimentary rock assemblages, their respective compositions and basement structures. The southern and central parts of the Uda-Murgal arc were a continental margin belt with heterogeneous basement represented by metamorphic rocks of the Siberian craton, the Verkhoyansk terrigenous complex of Siberian passive margin and the Koni-Taigonos Late Paleozoic to Early Mesozoic island arc with accreted oceanic terranes. At the present day latitude of the Pekulney and Chukotka segments there was an ensimatic island arc with relicts of the South Anyui oceanic basin in a backarc basin. Accretionary prisms of the Uda-Murgal arc and accreted terranes contain fragments of Permian, Triassic to Jurassic and Jurassic to Cretaceous (Tithonian–Valanginian) oceanic crust and Jurassic ensimatic island arcs. Paleomagnetic and faunal data show significant displacement of these oceanic complexes and the terranes of the Taigonos Peninsula were originally parts of the Izanagi oceanic plate.

Description: 〈b〉Provenance analysis and tectonic setting of the Triassic clastic deposits in Western Chukotka, Northeast Russia〈/b〉〈br〉 M. I. Tuchkova, S. Sokolov, and I. R. Kravchenko-Berezhnoy〈br〉 Stephan Mueller Spec. Publ. Ser., 4, 177-200, https://doi.org/10.5194/smsps-4-177-2009, 2009〈br〉 The study area is part of the Anyui subterrane of the Chukotka microplate, a key element in the evolution of the Amerasia Basin, located in Western Chukotka, Northeast Russia. The subterrane contains variably deformed, folded and cleaved rhythmic Triassic terrigenous deposits which represent the youngest stage of widespread marine deposition which form three different complexes: Lower-Middle Triassic, Upper Triassic (Carnian) and Upper Triassic (Norian). All of the complexes are represented by rhythmic interbeds of sandstone, siltstone and mudstone. Macrofaunas are not numerous, and in some cases deposits are dated by analogy to, or by their relationship with, other units dated with macrofaunas. The deposits are composed of pelagic sediments, low-density flows, high-density flows, and shelf facies associations suggesting that sedimentation was controlled by deltaic progradation on a continental shelf and subsequent submarine fan sedimentation at the base of the continental slope. Petrographic study of the mineral composition indicates that the sandstones are lithic arenites. Although the Triassic sandstones appear similar in outcrop and by classification, the constituent rock fragments are of diverse lithologies, and change in composition from lower grade metamorphic rocks in the Lower-Middle Triassic to higher grade metamorphic rocks in the Upper Triassic. This change suggests that the Triassic deposits represent an unroofing sequence as the source of the clastic material came from more deeply buried rocks with time.

Description: 〈b〉Multi-phase tectonic structures in the collision zone of the Kolyma-Omolon microcontinent and the eastern margin of the North Asian craton, Northeastern Russia〈/b〉〈br〉 A. V. Prokopiev and V. S. Oxman〈br〉 Stephan Mueller Spec. Publ. Ser., 4, 65-70, https://doi.org/10.5194/smsps-4-65-2009, 2009〈br〉 The sequence of formation of structures is established in the zone of junction of the eastern margin of the North Asian craton and the northeastern flank of the Kolyma-Omolon microcontinent, in the area of bend of the Kolyma structural loop. Detailed structural studies revealed two phases in the formation of Mesozoic structures – an early thrust phase and a late strike-slip phase. Structures formed during each of the phases are described. Thrust structures are represented by the Setakchan nappe on which the minimum amount of horizontal displacement is estimated at 13–15 km. Later superposed left-lateral strike-slip faults have a north strike. Formation of these latter structures occurred during the second phase of collision between the Kolyma-Omolon microcontinent and the eastern margin of the North Asian craton.

Description: 〈b〉Geology of the Shelves surrounding the New Siberian Islands, Russian Arctic〈/b〉〈br〉 D. Franke and K. Hinz〈br〉 Stephan Mueller Spec. Publ. Ser., 4, 35-44, https://doi.org/10.5194/smsps-4-35-2009, 2009〈br〉 A total of 11 700 km of multichannel seismic reflection data were acquired during recent reconnaissance surveys of the wide, shallow shelves of the Laptev and western East Siberian Seas around the New Siberian Islands. To the north of the Laptev Sea, the Gakkel Ridge, an active mid-ocean ridge which separates the North American and Eurasian Plates, meets abruptly the steep slope of the continental shelf. Extension has affected the Laptev Shelf since at least the Early Tertiary and has resulted in the formation of three major, generally north-south trending rift basins: the Ust' Lena Rift, the Anisin Basin, and the New Siberian Basin. Our data indicate that the rift basins on the Laptev Shelf are not continuous with those on the East Siberian Shelf. The latter shelf can best be described as an epicontinental platform which has undergone continuous subsidence since the Late Cretaceous. The greatest subsidence occurred in the NE, manifested by a major depocentre filled with inferred (?)Late Cretaceous to Tertiary sediments up to 5 s (twt) thick. On the basis of deep reflection data we revise and adjust Mesozoic domain boundaries around the New Siberian Islands.

Description: 〈b〉A seismic swarm near Neshkan, Chukotka, northeastern Russia, and implications for the boundary of the Bering plate〈/b〉〈br〉 K. G. Mackey, K. Fujita, B. M. Sedov, L. V. Gounbina, and S. Kurtkin〈br〉 Stephan Mueller Spec. Publ. Ser., 4, 261-271, https://doi.org/10.5194/smsps-4-261-2009, 2009〈br〉 A seismic swarm lasting over two years occurred near the village of Neshkan, Chukotka, far northeastern Russia, beginning with a 〈i〉M〈sub〉L〈/sub〉〈/i〉, 4.2 (4.1 〈i〉m〈sub〉b〈/sub〉〈/i〉) earthquake on 9 December 2002. The swarm generated considerable anxiety among the local populace and authorities. Two temporary seismic stations were deployed during the latter part of September 2003, and recorded over 150 events with magnitudes up to 3.0. Eighteen locatable events appear to form a northeast striking linear trend, parallel to other seismicity trends in Chukotka, extending 20 km to the southwest from the village. We interpret this trend as a previously unknown fault. A small pond located ~1 km west of the village drained and some apparent surface deformation was observed over the course of the earthquake sequence. Relocation of historic seismicity in the region shows that a magnitude 6.0 in 1996 may have ruptured an adjacent fault segment. Other, less well located but larger, teleseismic events earlier in the 20th century may also have occurred on or near this fault. The seismicity is consistent a proposed region of transtension along the northern boundary of a Bering plate.

Description: 〈b〉South Anyui suture: tectono-stratigraphy, deformations, and principal tectonic events〈/b〉〈br〉 S. D. Sokolov, G. Ye. Bondarenko, P. W. Layer, and I. R. Kravchenko-Berezhnoy〈br〉 Stephan Mueller Spec. Publ. Ser., 4, 201-221, https://doi.org/10.5194/smsps-4-201-2009, 2009〈br〉 Geochronologic and structural data from the terranes of the South Anyui suture zone record a protracted deformational history before, during and after an Early Cretaceous collision of the passive margin of the Chukotka-Arctic Alaska continental block with the active continental margin of the North Asian continent. Preceding this collision, the island arc complexes of the Yarakvaam terrane on the northern margin of the North Asian craton record Early Carboniferous to Neocomian ages in ophiolite, sedimentary, and volcanic rocks. Triassic to Jurassic amphibolites constrain the timing of subduction and intraoceanic deformation along this margin. The protracted (Neocomian to Aptian) collision of the Chukotka passive margin with the North Asian continent is preserved in a range of structural styles including first north verging folding, then south verging folding, and finally late collisional dextral strike slip motions which likely record a change from orthogonal collision to oblique collision. Due to this collision, the southern passive margin of Chukotka was overthrust by tectonic nappes composed of tectono-stratigraphic complexes of the South Anyui terrane. Greenschists with ages of 115–119 Ma are related to the last stages of this collision. The postcollisional orogenic stage (Albian to Cenomanian) is characterized by sinistral strike slip faults and an extensional environment.

Description: 〈b〉Deformation of the Northwestern Okhotsk Plate: How is it happening?〈/b〉〈br〉 D. Hindle, K. Fujita, and K. Mackey〈br〉 Stephan Mueller Spec. Publ. Ser., 4, 147-156, https://doi.org/10.5194/smsps-4-147-2009, 2009〈br〉 The Eurasia (EU) – North America (NA) plate boundary zone across Northeast Asia still presents many open questions within the plate tectonic paradigm. Constraining the geometry and number of plates or microplates present in the plate boundary zone is especially difficult because of the location of the EU-NA euler pole close to or even upon the EU-NA boundary. One of the major challenges remains the geometry of the Okhotsk plate (OK). whose northwestern portion terminates on the EU-OK-NA triple junction and is thus caught and compressed between converging EU and NA. We suggest that this leads to a coherent and understandable large scale deformation pattern of mostly northwest-southeast trending strike-slip faults which split Northwest OK into several extruding slivers. When the fault geometry is analysed together with space geodetic and focal mechanism data it suggests a central block which is extruding faster bordered east and west by progressively slower extruding blocks until the OK plate boundary faults are encountered. Taking into account elastic loading from both the intra-OK faults and the OK-Pacific (PA) boundary reconciles geodetic motions with geologic slip rates on at least the OK-NA boundary which corresponds to the Ulakhan fault.

Description: 〈b〉The Tommot pluton: a Middle Paleozoic rift-related alkaline gabbro and syenite complex, Yakutia, northeast Russia〈/b〉〈br〉 V. A. Trunilina, P. W. Layer, L. M. Parfenov, A. I. Zaitsev, and Y. S. Orlov〈br〉 Stephan Mueller Spec. Publ. Ser., 4, 97-109, https://doi.org/10.5194/smsps-4-97-2009, 2009〈br〉 The Tommot pluton is located within the continental Omulevka terrane of the inner zone of the Verkhoyansk-Kolyma Mesozoic orogen. It is a small complex (~12 km〈sup〉2〈/sup〉) composed of alkaline-ultramafic rocks, alkaline and subalkaline gabbroids, and alkaline and quartz syenites. The pluton is unique both in the composition and age of its constituent rocks. Mineralogical-petrographical and geochemical studies of the rocks indicate that the alkaline rocks resulted from the melting of depleted mantle horizons. K-Ar, Rb-Sr, and 〈sup〉40〈/sup〉Ar/〈sup〉39〈/sup〉Ar age determinations confirm a Paleozoic age of the rocks. Formation of the alkaline rocks is related to Middle Paleozoic rifting which occurred as two discrete events: a Late Devonian event, which affected the marginal part of the Siberian continent, and a Late Carboniferous event that reflects internal deformation of the Omulevka terrane or late-stage extension. A spatially associated alkali granite, the Somnitel'nyy pluton, is Late Jurassic–Early Cretaceous in age and is synchronous with accretion of the Kolyma-Omolon Superterrane to Siberia in the Mesozoic.

Description: 〈b〉Biography – Leonid M. Parfenov (1937–2002)〈/b〉〈br〉 K. Fujita, A. V. Prokopiev, and D. B. Stone〈br〉 Stephan Mueller Spec. Publ. Ser., 4, 5-9, https://doi.org/10.5194/smsps-4-5-2009, 2009〈br〉

Description: 〈b〉Geochronology and thermochronology of Cretaceous plutons and metamorphic country rocks, Anyui-Chukotka fold belt, North East Arctic Russia〈/b〉〈br〉 E. L. Miller, S. M. Katkov, A. Strickland, J. Toro, V. V. Akinin, and T. A. Dumitru〈br〉 Stephan Mueller Spec. Publ. Ser., 4, 157-175, https://doi.org/10.5194/smsps-4-157-2009, 2009〈br〉 U-Pb isotopic dating of seven granitoid plutons and associated intrusions from the Bilibino region (Arctic Chukotka, Russia) was carried out using the SHRIMP-RG. The crystallization ages of these granitoids, which range from approximately 116.9±2.5 to 108.5±2.7 Ma, bracket two regionally significant deformational events. The plutons cut folds, steep foliations and thrust-related structures related to sub-horizontal shortening at lower greenschist facies conditions (〈i〉D〈/i〉〈sub〉1〈/sub〉), believed to be the result of the collision of the Arctic Alaska-Chukotka microplate with Eurasia along the South Anyui Zone (SAZ). Deformation began in the Late Jurassic, based on fossil ages of syn-orogenic clastic strata, and involves strata as young as early Cretaceous (Valanginian) north of Bilibino and as young as Hauterivian-Barremian, in the SAZ. The second phase of deformation (〈i〉D〈/i〉〈sub〉2〈/sub〉) is developed across a broad region around and to the east of the Lupveem batholith of the Alarmaut massif and is interpreted to be coeval with magmatism. 〈i〉D〈/i〉〈sub〉2〈/sub〉 formed gently-dipping, high-strain foliations (〈i〉S〈/i〉〈sub〉2〈/sub〉). Growth of biotite, muscovite and actinolite define 〈i〉S〈/i〉〈sub〉2〈/sub〉 adjacent to the batholith, while chlorite and white mica define 〈i〉S〈/i〉〈sub〉2〈/sub〉 away from the batholith. Sillimanite (± andalusite) at the southeastern edge the Lupveem batholith represent the highest grade metamorphic minerals associated with 〈i〉D〈/i〉〈sub〉2〈/sub〉. 〈i〉D〈/i〉〈sub〉2〈/sub〉 is interpreted to have developed during regional extension and crustal thinning. Extension directions as measured by stretching lineations, quartz veins, boudinaged quartz veins is NE-SW to NW-SE. Mapped dikes associated with the plutons trend mostly NW-SE and indicate NE-SW directed extension. 〈sup〉40〈/sup〉Ar/〈sup〉39〈/sup〉Ar ages from 〈i〉S〈/i〉〈sub〉2〈/sub〉 micas range from 109.3±1.2 to 103.0±1.8 Ma and are interpreted as post-crystallization cooling ages following a protracted period of magmatism and high heat flow. Regional uplift and erosion of many kilometers of cover produced a subdued erosional surface prior to the eruption of volcanic rocks of the Okhotsk-Chukotka volcanic belt (OCVB) whose basal units (~87 Ma) overlie this profound regional unconformity. A single fission track age on apatite from granite in the Alarmaut massif yielded an age of 90±11 Ma, in good agreement with this inference.

Description: 〈b〉Summary of Northeast Asia geodynamics and tectonics*〈/b〉〈br〉 L. M. Parfenov, G. Badarch, N. A. Berzin, A. I. Khanchuk, M. I. Kuzmin, W. J. Nokleberg, A. V. Prokopiev, M. Ogasawara, and H. Yan〈br〉 Stephan Mueller Spec. Publ. Ser., 4, 11-33, https://doi.org/10.5194/smsps-4-11-2009, 2009〈br〉 The compilation, synthesis, description, and interpretation of regional geology and tectonics of major regions, such as Northeast Asia (Eastern Russia, Mongolia, Northern China, South Korea, and Japan) and the Circum-North Pacific (the Russian Far East, Alaska, and Canadian Cordillera), requires a complex methodology that includes: (1) definitions of key terms; (2) compilation of a regional geodynamics map that can be constructed according to modern tectonic concepts and definitions; and (3) formulation of a comprehensive tectonic model. This article presents a summary of the regional dynamics and tectonics of Northeast Asia as developed for a new summary geodynamics map of the region. This article also describes how a high-quality summary geodynamics map and companion tectonic analysis substantially aids in the understanding of the origin of major rock units, major structures, and contained mineral and fuel resources, and provides important guidelines for new research. 〈br〉〈br〉 * Prepared in memory of Leonid M. Parfenov, the leader of the geodynamics map team for the International collaborative project on NE Asia tectonics and metallogenesis.

Description: 〈b〉Detrital zircon geochronologic tests of the SE Siberia-SW Laurentia paleocontinental connection〈/b〉〈br〉 J. S. MacLean, J. W. Sears, K. R. Chamberlain, A. K. Khudoley, A. V. Prokopiev, A. P. Kropachev, and G. G. Serkina〈br〉 Stephan Mueller Spec. Publ. Ser., 4, 111-116, https://doi.org/10.5194/smsps-4-111-2009, 2009〈br〉 Strikingly similar Late Mesoproterozoic stratigraphic sequences and correlative U-Pb detrital-zircon ages may indicate that the Sette Daban region of southeastern Siberia and the Death Valley region of southwestern North America were formerly contiguous parts of a Grenville foreland basin. The Siberian section contains large numbers of detrital zircons that correlate with Grenville, Granite-Rhyolite, and Yavapai basement provinces of North America. The sections in both Siberia and Death Valley exhibit west-directed thrust faults that may represent remnants of a Grenville foreland thrust belt. North American detrital-zircon components do not occur in Siberian samples above a ~600 Ma breakup unconformity, suggesting that rifting and continental separation blocked transfer of clastic sediment between the cratons by 600 Ma. Faunal similarities suggest, however, that the two cratons remained within the breeding ranges of Early Cambrian trilobites and archeocyathans.

Description: 〈b〉Seismotectonics of the Chersky Seismic Belt, eastern Sakha Republic (Yakutia) and Magadan District, Russia〈/b〉〈br〉 K. Fujita, B. M. Koz'min, K. G. Mackey, S. A. Riegel, M. S. McLean, and V. S. Imaev〈br〉 Stephan Mueller Spec. Publ. Ser., 4, 117-145, https://doi.org/10.5194/smsps-4-117-2009, 2009〈br〉 The Chersky seismic belt represents a zone of deformation between the North American and Eurasian plates in northeast Russia. The belt extends from the Laptev Sea into the Chersky Range where it splits into two branches. One branch extends to Kamchatka and the Aleutian-Kurile Junction, while the other branch extends south towards Sakhalin Island. Focal mechanisms indicate a change from extension to transpression in the northern Verkhoyansk Range and generally left-lateral transpression in the Chersky Range extending to northern Kamchatka. The few focal mechanisms on the second branch suggest right-lateral transpression. A large number of faults, sub-parallel to the seismicity and presumed to be strike-slip, are visible in satellite imagery and topographic maps and are also associated with seismically generated landslides. 〈br〉〈br〉 These data support a model in which the Sea of Okhotsk forms the core of a separate Okhotsk microplate surrounded by diffuse boundaries on the north and west. Microseismicity in continental northeast Russia is most heavily concentrated within and between the fault systems along the northern boundary of the proposed Okhotsk plate and indicates a high level of deformation. The sense of slip on the faults (both from focal mechanisms and geology) are also generally consistent with the extrusion of the Okhotsk plate to the southeast as it is compressed between its larger neighbors. The northernmost part of the Okhotsk plate may be decoupled to some degree from the more stable central Sea of Okhotsk.

Description: 〈b〉Review of geology of the New Siberian Islands between the Laptev and the East Siberian Seas, North East Russia〈/b〉〈br〉 M. Kos'ko and E. Korago〈br〉 Stephan Mueller Spec. Publ. Ser., 4, 45-64, https://doi.org/10.5194/smsps-4-45-2009, 2009〈br〉 The New Siberian Islands comprise De Long Islands, Anjou Islands, and Lyakhov Islands. Early Paleozoic, Mesozoic and Cenozoic sediments and igneous rocks are known on the De Long Islands. Cambrian slate, siltstone, mudstone and silicified limestone occur on Bennett Island. Ordovician volcanogenic turbidites, lavas, and small intrusions of andesite-basalt, basalt, dolerite, and porphyritic diorite were mapped on Henrietta Island. The igneous rocks are of calc-alkaline island arc series. The Ordovician age of the sequence was defined radiometrically. Early Paleozoic strata were faulted and folded presumably in the Caledonian time. Early Cretaceous sandstone and mudstone are known on Bennett Island. They are overlain by a 106–124 Ma basalt unit. Cenozoic volcanics are widespread on the De Long Islands. Zhokhov Island is an eroded stratovolcano. The volcanics are mostly of picrite-olivine type and limburgite. Radiometric dating indicates Miocene to Recent ages for Cenozoic volcanism. 〈br〉〈br〉 On the Anjou islands Lower-Middle Paleozoic strata consist of carbonates, siliciclastics, and clay. A Northwest-southeast syn-sedimentary facies zonation has been reconstructed. Upper Paleozoic strata are marine carbonate, clay and siliciclastic facies. Mudstone and clay predominate in the Triassic to Upper Jurassic section. Aptian-Albian coal bearing deposits uconformably overlap lower strata indicating Early Cretaceous tectonism. Upper Cretaceous units are mostly clay and siltstone with brown coal strata resting on Early Cretaceous weathered rhyolite. Cenozoic marine and nonmarine silisiclastics and clay rest upon the older units with a transgressive unconformity including a weathering profile in the older rocks. 〈br〉〈br〉 Manifestations of Paleozoic and Triassic mafic and Cretaceous acidic magmatism are also found on these islands. The pre-Cretaceous structure of the Anjou islands is of a block and fold type Late Cimmerian in age followed by faulting in Cenozoic time. 〈br〉〈br〉 The Lyakhov islands are located at the western end of the Late Cimmerian South Anyui suture. Sequences of variable age, composition, and structural styles are known on the Lyakhov Islands. These include an ancient metamorphic sequence, Late Paleozoic ophiolitic sequence, Late Mesozoic turbidite sequence, Cretaceous granites, and Cenozoic sediments. Fold and thrust imbricate structures have been mapped on southern Bol'shoi Lyakhov Island. North-northwestern vergent thrusts transect the Island and project offshore. Open folds of Jurassic–Early Cretaceous strata are characteristic of Stolbovoi and Malyi Lyakhov islands. 〈br〉〈br〉 Geology of the New Siberian Islands supports the concept of a circum Arctic Phanerozoic fold belt. The belt is comprised of Caledonian, Ellesmerian, Early Cimmerian and Late Cimmerian fold systems, manifested in many places on the mainland and on islands around the Arctic Ocean. Knowledge of the geology of the New Siberian Islands has been used to interpret anomalous gravity and magnetic field maps and Multi Channel Seismic (MCS) lines. Two distinguishing structural stages are universally recognized within the offshore sedimentary cover which correlate with the onshore geology of the New Siberian Islands. Dating of the upper structural stage and constituent seismic units is based on structural and stratigraphic relationships between Late Mesozoic and Cenozoic units in the archipelago. The Laptev Sea–western East Siberian Sea seismostratigraphic model for the upper structural stage has much in common with the seismostratigraphic model in the American Chukchi Sea.

Description: 〈b〉Editors and Acknowledgement〈/b〉〈br〉 〈br〉 Stephan Mueller Spec. Publ. Ser., 4, 1-1, https://doi.org/10.5194/smsps-4-1-2009, 2009〈br〉

Description: 〈b〉Morphology of the pore space in claystones – evidence from BIB/FIB ion beam sectioning and cryo-SEM observations〈/b〉〈br〉 G. Desbois, J. L. Urai, and P. A. Kukla〈br〉 eEarth, 4, 15-22, doi:10.5194/ee-4-15-2009, 2009〈br〉 The morphology of pore space has a strong effect on mechanical and transport properties of mudrocks and clay-rich fault gouge, but its characterization has been mostly indirect. We report on a study of Boom clay from a proposed disposal site of radioactive waste (Mol site, Belgium) using high resolution SEM at cryogenic temperature, with ion beam cross-sectioning to prepare smooth, damage free surfaces. Pores commonly have crack-like tips, preferred orientation parallel to bedding and power law size distribution. We define a number of pore types depending on shape and location in the microstructure: large jagged pores in strain shadows of clastic grains, high aspect ratio pores between similarly oriented phyllosilicate grains and crescent-shaped pores in saddle reefs of folded phyllosilicates. 3-D reconstruction by serial cross-sectioning shows 3-D connectivity of the pore space. These findings offer a new insight into the morphology of pores down to nano-scale in comparison to traditional pore size distributions calculated from mercury Injection experiments, explain slaking of clays by successive wetting and drying and provide the basis for microstructure-based models of transport in clays.

Description: 〈b〉Late Pleistocene palaeoproductivity patterns during the last climatic cycle in the Guyana Basin as revealed by calcareous nannoplankton〈/b〉〈br〉 G.-E. López-Otálvaro, J. A. Flores, F. J. Sierro, I. Cacho, J.-O. Grimalt, E. Michel, E. Cortijo, and L. Labeyrie〈br〉 eEarth, 4, 1-13, doi:10.5194/ee-4-1-2009, 2009〈br〉 Variations in the assemblage and abundance of coccoliths reveal changes in oceanic and atmospheric dynamics in the Guyana Basin over the last climatic cycle, mainly linked to latitudinal variations in the ITCZ (Intertropical Convergence Zone). Records of the N ratio (a palaeoproductivity index of coccolithophores) allowed us to monitor nutri-thermocline fluctuations. Nannofossil accumulation rates (NAR) vary closely with the N ratio, indicating a strong correlation between these two palaeoproductivity proxies. Decreases in the N ratio and NAR values suggest lower palaeoproductivity during glacial substages, indicating a deep nutri-thermocline (deep stratification of the mixed layer) as a consequence of the piling up of warm water dragged by the NEC. This setting was favoured by the southern shift of the ITCZ and Trade winds which blew perpendicular to the Guyana coast. By contrast, increases in the N ratio and NAR values revealed higher palaeoproductivity during interglacial substages, suggesting a shoaling of the nutri-thermocline. This scenario is favoured by a northward displacement of the ITCZ with the southeast Trade winds blowing alongshore. Additionally, palaeoproductivity changes during substages of Marine Isotope Stage (MIS) 6-5 are of much higher amplitude than those recorded in substages of MIS 4-2 and the early Holocene. Similarities between the palaeoproductivity and the 65° N summer insolation records, suggest a link between the depth of nutri-thermocline, the latitudinal migration of the ITCZ and ice-sheet changes in the Northern Hemisphere.

Description: 〈b〉Late Pleistocene palaeoproductivity patterns during the last climatic cycle in the Guyana Basin as revealed by calcareous nannoplankton〈/b〉〈br〉 G.-E. López-Otálvaro, J. A. Flores, F. J. Sierro, I. Cacho, J.-O. Grimalt, E. Michel, E. Cortijo, and L. Labeyrie〈br〉 eEarth, 4, 1-13, doi:10.5194/ee-4-1-2009, 2009〈br〉 Variations in the assemblage and abundance of coccoliths reveal changes in oceanic and atmospheric dynamics in the Guyana Basin over the last climatic cycle, mainly linked to latitudinal variations in the ITCZ (Intertropical Convergence Zone). Records of the N ratio (a palaeoproductivity index of coccolithophores) allowed us to monitor nutri-thermocline fluctuations. Nannofossil accumulation rates (NAR) vary closely with the N ratio, indicating a strong correlation between these two palaeoproductivity proxies. Decreases in the N ratio and NAR values suggest lower palaeoproductivity during glacial substages, indicating a deep nutri-thermocline (deep stratification of the mixed layer) as a consequence of the piling up of warm water dragged by the NEC. This setting was favoured by the southern shift of the ITCZ and Trade winds which blew perpendicular to the Guyana coast. By contrast, increases in the N ratio and NAR values revealed higher palaeoproductivity during interglacial substages, suggesting a shoaling of the nutri-thermocline. This scenario is favoured by a northward displacement of the ITCZ with the southeast Trade winds blowing alongshore. Additionally, palaeoproductivity changes during substages of Marine Isotope Stage (MIS) 6-5 are of much higher amplitude than those recorded in substages of MIS 4-2 and the early Holocene. Similarities between the palaeoproductivity and the 65° N summer insolation records, suggest a link between the depth of nutri-thermocline, the latitudinal migration of the ITCZ and ice-sheet changes in the Northern Hemisphere.

Description: 〈b〉Geotectonic setting of the Tertiary Uyandina and Indigirka-Zyryanka basins, Republic Sakha (Yakutia), Northeast Russia, using coal rank data〈/b〉〈br〉 H.-J. Paech〈br〉 Stephan Mueller Spec. Publ. Ser., 4, 85-96, https://doi.org/10.5194/smsps-4-85-2009, 2009〈br〉 Outcrops along the Inach River in the Uyandina basin and those along the Myatis' River in the Indigirka-Zyryanka basin were studied in detail and sampled for coal rank determinations. The Uyandina basin is an intramontane pull-apart basin characterized by extensional structures within the Moma rift system. The coal rank is below 0.3% vitrinite reflectance (〈i〉R〈sub〉r〈/sub〉〈/i〉), which indicates shallow, immature conditions of basin formation and very low subsidence. The Myatis' River coal-bearing outcrops in the Indigirka-Zyryanka basin reveal compression induced by continent collision. The compressive deformation includes also lowermost Pliocene strata. Due to the position in the Verkhoyansk-Chersky fold belt adjacent to the Kolyma-Omolon microcontinent the Indigirka-Zyryanka basin has much in common with a foredeep, i.e. the asymmetry in thickness and tectonic structure. The vitrinite reflectance data (〈i〉R〈sub〉r〈/sub〉〈/i〉) which range from 0.25% to more than 5% reinforce the accepted models that describe basin subsidence and geothermal history and the tectonic deformation.

Description: 〈b〉Structural studies near Pevek, Russia: implications for formation of the East Siberian Shelf and Makarov Basin of the Arctic Ocean〈/b〉〈br〉 E. L. Miller and V. E. Verzhbitsky〈br〉 Stephan Mueller Spec. Publ. Ser., 4, 223-241, https://doi.org/10.5194/smsps-4-223-2009, 2009〈br〉 The Pevek region of Arctic Russia provides excellent beach cliff exposure of sedimentary and igneous rocks that yield detailed information on the nature, progression and timing of structural events in this region. Regional folding and thrust faulting, with the development of a south-dipping axial plane cleavage/foliation developed during N-S to NE-SW directed shortening and formation of the Chukotka-Anyui fold belt. This deformation involves strata as young as Valanginian (136–140 Ma, Gradstein et al., 2004). Fold-related structures are cut by intermediate to silicic batholiths, plutons and dikes of Cretaceous age. Reported K-Ar whole rock and mineral ages on the granitoids range from 144 to 85 Ma, but to the south, more reliable U-Pb zircon ages on compositionally similar plutons yield a much narrower age range of ~120–105 Ma (Miller et al., this volume) and a pluton in Pevek yields a U-Pb age on zircon of 108.1±1.1 Ma with evidence for inheritance of slightly older 115 Ma zircons. Magmas were intruded during an episode of E-W to ENE-WSW directed regional extension based on the consistent N-S to NNW-SSE orientation of over 800 mapped dikes and quartz veins. Analysis of small-offset faults and slickensides yield results compatible with those inferred from the dikes. Younger tectonic activity across this region is minor and the locus of magmatic activity moved southward towards the Pacific margin as represented by the 〈90 Ma Okhotsk-Chukotsk volcanic belt (OCVB). A lengthy period of uplift and erosion occurred after emplacement of Cretaceous plutons and produced the peneplain beneath the younger OCVB. 〈br〉〈br〉 Based on our studies, we speculate that ~120–105 Ma magmatism, which heralds a change in tectonic regime from compression to extension, could represent one of the consequences of the inception of rifting in the Amerasian Basin of the Arctic, forming the Makarov Basin north of the Siberian shelf at this longitude. A synthesis of available seismic reflection, gravity and magnetic data for the offshore Siberian Shelf reveals a widespread, seismically mappable basement-sedimentary cover contact that deepens northward towards the edge of the shelf with few other significant basins. Various ages have been assigned to the oldest strata above the unconformity, ranging from Cretaceous (Albian – 112–100 Ma) to Tertiary (Paleocene–Eocene – ~60–50 Ma). The period of uplift and erosion documented along the Arctic coast of Russia at this longitude could represent the landward equivalent of the (yet undrilled) offshore basement-sedimentary cover contact, thus overlying sedimentary sequences could be as old as early Late Cretaceous. Although quite speculative, these conclusions suggest that land-based geologic, structural, petrologic and geochronologic studies could provide useful constraints to help resolve the plate tectonic history of the Arctic Ocean.

Description: 〈b〉Middle Paleozoic-Mesozoic boundary of the North Asian craton and the Okhotsk terrane: new geochemical and geochronological data and their geodynamic interpretation〈/b〉〈br〉 A. V. Prokopiev, J. Toro, J. K. Hourigan, A. G. Bakharev, and E. L. Miller〈br〉 Stephan Mueller Spec. Publ. Ser., 4, 71-84, https://doi.org/10.5194/smsps-4-71-2009, 2009〈br〉 The Okhotsk terrane, located east of the South Verkhoyansk sector of the Verkhoyansk fold-and-thrust belt, has Archean crystalline basement and Riphean to Early Paleozoic sedimentary cover similar to that of the adjacent the North Asian craton. However, 2.6 Ga biotite orthogneisses of the Upper Maya uplift of the Okhotsk terrane yielded Early Devonian 〈sup〉40〈/sup〉Ar/〈sup〉39〈/sup〉Ar cooling ages, evidence of a Mid-Paleozoic metamorphic event not previously known. These gneisses are also intruded by 375±2 Ma (Late Devonian) calc-alkaline granodiorite plutons that we interpret as part of a continental margin volcanic arc. Therefore, Late Devonian rifting, which affected the entire eastern margin of North Asia separating the Okhotsk terrane from the North Asian craton, was probably a back-arc event. 〈br〉〈br〉 Our limited 〈sup〉40〈/sup〉Ar/〈sup〉39〈/sup〉Ar data from the South Verkhoyansk metamorphic belt suggests that low grade metamorphism and deformation started in the Late Jurassic due to accretion of the Okhotsk terrane to the North Asia margin along the Bilyakchan fault. Shortening and ductile strain continued in the core of the South Verkhoyansk metamorphic belt until about 120 Ma due to paleo-Pacific subduction along the Uda-Murgal continental margin arc.

Description: Hyrrokkin sarcophaga is a parasitic foraminifera that is commonly found in cold-water coral reefs where it infests the file clam Acesta excavata and the scleractinian coral Desmophyllum pertusum (formerly known as Lophelia pertusa). Here, we present measurements of the trace element and isotopic composition of these parasitic foraminifera, analyzed by inductively coupled optical emission spectrometry (ICP-OES), electron probe microanalysis (EPMA) and mass spectrometry (gas-source MS and inductively-coupled-plasma MS). Our results reveal that the geochemical signature of H. sarcophaga depends on the host organism it infests. Sr / Ca ratios are 1.1 mmol mol−1 higher in H. sarcophaga that infest D. pertusum, which could be an indication that dissolved host carbonate material is utilized in shell calcification, given that the aragonite of D. pertusum has a naturally higher Sr concentration compared to the calcite of A. excavata. Similarly, we measure 3.1 ‰ lower δ13C and 0.25 ‰ lower δ18O values in H. sarcophaga that lived on D. pertusum, which might be caused by the direct uptake of the host's carbonate material with a more negative isotopic composition or different pH regimes in these foraminifera (pH can exert a control on the extent of CO2 hydration/hydroxylation) due to the uptake of body fluids of the host. We also observe higher Mn / Ca ratios in foraminifera that lived on A. excavata but did not penetrate the host shell compared to specimen that penetrated the shell, which could be interpreted as a change in food source, changes in the calcification rate, Rayleigh fractionation or changing oxygen conditions. While our measurements provide an interesting insight into the calcification process of this unusual foraminifera, these data also indicate that the geochemistry of this parasitic foraminifera is unlikely to be a reliable indicator of paleoenvironmental conditions using Sr / Ca, Mn / Ca, δ18O or δ13C unless the host organism is known and its geochemical composition can be accounted for.

Description: In this work we have performed a detailed study of vectors to ore to a representative volcanic-rock-hosted replacive volcanogenic massive sulfide (VMS) deposit located in the northern Iberian Pyrite Belt (Spain), the Aguas Teñidas deposit. The investigated vectors include the following: (1) mineralogical zoning, (2) host sequence characterization and mineralized unit identification based on whole rock geochemistry discrimination diagrams, (3) study of the characteristics and behaviour of whole rock geochemical anomalies around the ore (e.g. alteration-related compositional changes, characteristics and extent of geochemical halos of indicative elements such as Cu, Zn, Pb, Sb, Tl, and Ba around the deposit), and (4) application of portable X-ray fluorescence (p-XRF) analysis to the detection of the previous vectors. In the footwall, a concentric cone-shaped hydrothermal alteration zone bearing the stockwork passes laterally, from core to edge, from quartz (only local) to chlorite–quartz, sericite–chlorite–quartz, and sericite–quartz alteration zones. The hydrothermal alteration is also found in the hanging wall despite being tectonically allochthonous to the orebody: a proximal sericite alteration zone is followed by a more distal albite-rich one. Whole rock major elements show an increase in alteration indexes (e.g. AI, CCPI) towards the mineralization, a general SiO2 enrichment, and FeO enrichment as well as K2O and Na2O depletion towards the centre of the hydrothermal system, with MgO showing a less systematic behaviour. K2O and Na2O leached from the centre of the system are transported and deposited in more external areas. Copper, Pb, and Zn produce proximal anomalies around mineralized areas, with the more mobile Sb, Tl, and Ba generating wider halos. Whereas Sb and Tl halos form around all mineralized areas, Ba anomalies are restricted to areas around the massive sulfide body. Our results show that proposed vectors, or adaptations designed to overcome p-XRF limitations, can be confidently used by analysing unprepared hand specimens, including the external rough curved surface of drill cores. The data presented in this work are not only applicable to VMS exploration in the Iberian Pyrite Belt, but on a broader scale they will also contribute to improving our general understanding of vectors to ore in replacive-type VMS deposits.

Description: We combine satellite data products to provide a first and general overview of the physical sea ice conditions along the drift of the international Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition and a comparison with previous years (2005–2006 to 2018–2019). We find that the MOSAiC drift was around 20 % faster than the climatological mean drift, as a consequence of large-scale low-pressure anomalies prevailing around the Barents–Kara–Laptev sea region between January and March. In winter (October–April), satellite observations show that the sea ice in the vicinity of the Central Observatory (CO; 50 km radius) was rather thin compared to the previous years along the same trajectory. Unlike ice thickness, satellite-derived sea ice concentration, lead frequency and snow thickness during winter months were close to the long-term mean with little variability. With the onset of spring and decreasing distance to the Fram Strait, variability in ice concentration and lead activity increased. In addition, the frequency and strength of deformation events (divergence, convergence and shear) were higher during summer than during winter. Overall, we find that sea ice conditions observed within 5 km distance of the CO are representative for the wider (50 and 100 km) surroundings. An exception is the ice thickness; here we find that sea ice within 50 km radius of the CO was thinner than sea ice within a 100 km radius by a small but consistent factor (4 %) for successive monthly averages. Moreover, satellite acquisitions indicate that the formation of large melt ponds began earlier on the MOSAiC floe than on neighbouring floes.

Description: Across the Qinghai–Tibet Plateau (QTP) there is a narrow engineering corridor with widely distributed slopes called the Qinghai–Tibet Engineering Corridor (QTEC), where a variety of important infrastructures are concentrated. These facilities are transportation routes for people, materials, energy, etc. from inland China to the Tibet Autonomous Region. From Golmud to Lhasa, the engineering corridor covers 632 km of permafrost containing the densely developed Qinghai–Tibet Railway and Qinghai–Tibet Highway, as well as power and communication towers. Slope failure in permafrost regions, caused by permafrost degradation, ground ice melting, etc., affects the engineering construction and permafrost environments in the QTEC. We implement a variety of sensors to monitor the hydrological and thermal deformation between permafrost slopes and permafrost engineering projects in the corridor. In addition to soil temperature and moisture sensors, the global navigation satellite system (GNSS), terrestrial laser scanning (TLS), and unmanned aerial vehicles (UAVs) were adopted to monitor the spatial distribution and changes in thermal deformation. An integrated dataset of hydrological and thermal deformation in permafrost engineering and slopes in the QTEC from the 1950s to 2020, including meteorological and ground observations, TLS point cloud data, and RGB and thermal infrared (TIR) images, can be of great value for estimating the hydrological and thermal impact and stability between engineering and slopes under the influence of climate change and engineering disturbance. The dataset and code were uploaded to the Zenodo repository and can be accessed through https://zenodo.org/communities/qtec (last access: 23 June 2021), including meteorological and ground observations at https://doi.org/10.5281/zenodo.5009871 (Luo et al., 2020d), TLS measurements at https://doi.org/10.5281/zenodo.5009558 (Luo et al., 2020a), UAV RGB and TIR images at https://doi.org/10.5281/zenodo.5016192 (Luo et al., 2020b), and R code for permafrost indices and visualisation at https://doi.org/10.5281/zenodo.5002981 (Luo et al., 2020c).

Description: We present the first incorporation of the Common Representative Intermediates version 2.2 tropospheric chemistry mechanism, CRI v2.2, combined with stratospheric chemistry, into the global chemistry–climate United Kingdom Chemistry and Aerosols (UKCA) model to give the CRI-Strat 2 mechanism. A rigorous comparison of CRI-Strat 2 with the earlier version, CRI-Strat, is performed in UKCA in addition to an evaluation of three mechanisms, CRI-Strat 2, CRI-Strat and the standard UKCA chemical mechanism, StratTrop v1.0, against a wide array of surface and airborne chemical data. CRI-Strat 2 comprises a state-of-the-art isoprene scheme, optimized against the Master Chemical Mechanism v3.3.1, which includes isoprene peroxy radical isomerization, HOx recycling through the addition of photolabile hydroperoxy aldehydes (HPALDs), and isoprene epoxy diol (IEPOX) formation. CRI-Strat 2 also features updates to several rate constants for the inorganic chemistry, including the reactions of inorganic nitrogen and O(1D). The update to the isoprene chemistry in CRI-Strat 2 increases OH over the lowest 500 m in tropical forested regions by 30 %–50 % relative to CRI-Strat, leading to an improvement in model–observation comparisons for surface OH and isoprene relative to CRI-Strat and StratTrop. Enhanced oxidants also cause a 25 % reduction in isoprene burden and an increase in oxidation fluxes of isoprene and other biogenic volatile organic compounds (BVOCs) at low altitudes with likely impacts on subsequent aerosol formation, atmospheric lifetime, and climate. By contrast, updates to the rate constants of O(1D) with its main reactants relative to CRI-Strat reduces OH in much of the free troposphere, producing a 2 % increase in the methane lifetime, and increases the tropospheric ozone burden by 8 %, primarily from reduced loss via O(1D)+H2O. The changes to inorganic nitrogen reaction rate constants increase the NOx burden by 4 % and shift the distribution of nitrated species closer to that simulated by StratTrop. CRI-Strat 2 is suitable for multi-decadal model integrations and the improved representation of isoprene chemistry provides an opportunity to explore the consequences of HOx recycling in the United Kingdom Earth System Model (UKESM1). This new mechanism will enable a re-evaluation of the impact of BVOCs on the chemical composition of the atmosphere and further probe the feedback between the biosphere and the climate.

Description: The use of mass flow simulations in volcanic hazard zonation and mapping is often limited by model complexity (i.e. uncertainty in correct values of model parameters), a lack of model uncertainty quantification, and limited approaches to incorporate this uncertainty into hazard maps. When quantified, mass flow simulation errors are typically evaluated on a pixel-pair basis, using the difference between simulated and observed (“actual”) map-cell values to evaluate the performance of a model. However, these comparisons conflate location and quantification errors, neglecting possible spatial autocorrelation of evaluated errors. As a result, model performance assessments typically yield moderate accuracy values. In this paper, similarly moderate accuracy values were found in a performance assessment of three depth-averaged numerical models using the 2012 debris avalanche from the Upper Te Maari crater, Tongariro Volcano, as a benchmark. To provide a fairer assessment of performance and evaluate spatial covariance of errors, we use a fuzzy set approach to indicate the proximity of similarly valued map cells. This “fuzzification” of simulated results yields improvements in targeted performance metrics relative to a length scale parameter at the expense of decreases in opposing metrics (e.g. fewer false negatives result in more false positives) and a reduction in resolution. The use of this approach to generate hazard zones incorporating the identified uncertainty and associated trade-offs is demonstrated and indicates a potential use for informed stakeholders by reducing the complexity of uncertainty estimation and supporting decision-making from simulated data.

Description: In this study, we analyse the climatology of ionosphere over Nepal based on GPS-derived vertical total electron content (VTEC) observed from four stations as defined in Table 1: KKN4 (27.80∘ N, 85.27∘ E), GRHI (27.95∘ N, 82.49∘ E), JMSM (28.80∘ N, 83.74∘ E) and DLPA (28.98∘ N, 82.81∘ E) during the years 2008 to 2018. The study illustrates the diurnal, monthly, annual, seasonal and solar cycle variations in VTEC during all times of solar cycle 24. The results clearly reveal the presence of equinoctial asymmetry in TEC, which is more pronounced in maximum phases of solar cycle in the year 2014 at KKN4 station, followed by descending, ascending and minimum phases. Diurnal variations in VTEC showed the short-lived day minimum which occurs between 05:00 to 06:00 LT (local time) at all the stations considered, with diurnal peaks between 12:00 and 15:00 LT. The maximum value of TEC is observed more often during the spring equinox than the autumn equinox, with a few asymmetries. Seasonal variation in TEC is observed to be a manifestation of variations in solar flux, particularly regarding the level of solar flux in consecutive solstices.

Description: We report on applications of the ultraviolet-light-emitting-diode-based incoherent broadband cavity-enhanced absorption spectroscopy (UV-LED-IBBCEAS) technique for optical monitoring of HONO, NO2 and CH2O in a simulation chamber. Performance intercomparison of UV-LED-IBBCEAS with a wet chemistry-based NitroMAC sensor and a Fourier transform infrared (FTIR) spectrometer has been carried out on real-time simultaneous measurement of HONO, NO2 and CH2O concentrations during the reaction of NO2 with H2O vapour in CESAM (French acronym for Experimental Multiphasic Atmospheric Simulation Chamber). The 1σ (signal-to-noise ratio (SNR) = 1) detection limits of 112 pptv for NO2, 56 pptv for HONO and 41 ppbv for CH2O over 120 s were found for the UV-LED-IBBCEAS measurement. On the contrary to many set-ups where cavities are installed outside the simulation chamber, we describe here an original in situ permanent installation. The intercomparison results demonstrate that IBBCEAS is a very well suitable technique for in situ simultaneous measurements of multiple chemically reactive species with high sensitivity and high precision even if the absorption bands of these species are overlapped. It offers excellent capacity for non-invasive optical monitoring of chemical reactions without any perturbation. For the application to simulation chambers, it has the advantage to provide a spatially integrated measurement across the reactor and hence to avoid point-sampling-related artefacts.

Description: Assimilation of weather radar measurements including radar reflectivity and radial wind data has been operational at the Deutscher Wetterdienst, with a diagonal observation error (OE) covariance matrix. For an implementation of a full OE covariance matrix, the statistics of the OE have to be a priori estimated, for which the Desroziers method has been often used. However, the resulted statistics consists of contributions from different error sources and are difficult to interpret. In this work, we use an approach that is based on samples for truncation error in radar observation space to approximate the representation error due to unresolved scales and processes (RE) and compare its statistics with the OE statistics estimated by the Desroziers method. It is found that the statistics of the RE help the understanding of several important features in the variances and correlation length scales of the OE for both reflectivity and radial wind data and the other error sources from the microphysical scheme, radar observation operator and the superobbing technique may also contribute, for instance, to differences among different elevations and observation types. The statistics presented here can serve as a guideline for selecting which observations are assimilated and for assignment of the OE covariance matrix that can be diagonal or full and correlated.

Description: Determining controls on the temperature sensitivity of heterotrophic soil respiration remains critical to incorporating soil–climate feedbacks into climate models. Most information on soil respiratory responses to temperature comes from laboratory incubations of isolated soils and typically subsamples of individual horizons. Inconsistencies between field and laboratory results may be explained by microbial priming supported by cross-horizon exchange of labile C or N. Such exchange is feasible in intact soil profiles but is absent when soils are isolated from surrounding depths. Here we assess the role of soil horizon connectivity, by which we mean the degree to which horizons remain layered and associated with each other as they are in situ, on microbial C and N substrate use and its relationship to the temperature sensitivity of respiration. We accomplished this by exploring changes in C : N, soil organic matter composition (via C : N, amino acid composition and concentration, and nuclear magnetic resonance spectroscopy), and the δ13C of respiratory CO2 during incubations of organic horizons collected across boreal forests in different climate regions where soil C and N compositions differ. The experiments consisted of two treatments: soil incubated (1) with each organic horizon separately and (2) as a whole organic profile, permitting cross-horizon exchange of substrates during the incubation. The soils were incubated at 5 and 15 ∘C for over 430 d. Enhanced microbial use of labile C-rich, but not N-rich, substrates were responsible for enhanced, whole-horizon respiratory responses to temperature relative to individual soil horizons. This impact of a labile C priming mechanism was most emergent in soils from the warmer region, consistent with these soils' lower C bioreactivity relative to soils from the colder region. Specifically, cross-horizon exchange within whole soil profiles prompted increases in mineralization of carbohydrates and more 13C-enriched substrates and increased soil respiratory responses to warming relative to soil horizons incubated in isolation. These findings highlight that soil horizon connectivity can impact microbial substrate use in ways that affect how soil effluxes of CO2 are controlled by temperature. The degree to which this mechanism exerts itself in other soils remains unknown, but these results highlight the importance of understanding mechanisms that operate in intact soil profiles – only rarely studied – in regulating a key soil–climate feedback.

Description: Plumes from the boreal spring biomass burning (BB) in northern peninsular Southeast Asia (nPSEA) are lifted into the subtropical jet stream and transported and deposited across nPSEA, South China, Taiwan and even the western North Pacific Ocean. This paper as part of the Seven SouthEast Asian Studies (7-SEAS) project effort attempts to improve the chemical weather prediction capability of the Weather Research and Forecasting coupled with the Community Multiscale for Air Quality (WRF–CMAQ) model over a vast region, from the mountainous near-source burning sites at nPSEA to its downwind region. Several sensitivity analyses of plume rise are compared in the paper, and it is discovered that the initial vertical allocation profile of BB plumes and the plume rise module (PLMRIM) are the main reasons causing the inaccuracies of the WRF–CMAQ simulations. The smoldering emission from the Western Regional Air Partnership (WRAP) empirical algorithm included has improved the accuracies of PM10, O3 and CO at the source. The best performance at the downwind sites is achieved with the inline PLMRIM, which accounts for the atmospheric stratification at the mountainous source region with the FINN burning emission dataset. Such a setup greatly improves not only the BB aerosol concentration prediction over near-source and receptor ground-based measurement sites but also the aerosol vertical distribution and column aerosol optical depth of the BB aerosol along the transport route. The BB aerosols from nPSEA are carried by the subtropical westerlies in the free troposphere to the western North Pacific, while BB aerosol has been found to interact with the local pollutants in the Taiwan region through three conditions: (a) overpassing western Taiwan and entering the central mountain area, (b) mixing down to western Taiwan, (c) transport of local pollutants upwards and mixing with a BB plume on higher ground. The second condition, which involves the prevailing high-pressure system from Asian cold surge, is able to impact most of the population in Taiwan.

Description: The German Bight located within the central North Sea is a hydro- and morphodynamically highly complex system of estuaries, barrier islands, and part of the world's largest coherent tidal flats, the Wadden Sea. To identify and understand challenges faced by coastal stakeholders, such as harbor operators or governmental agencies, to maintain waterways and employ numerical models for further analyses, it is imperative to have a consistent database for both bathymetry and surface sedimentology. Current commercial and public data products are insufficient in spatial and temporal resolution and coverage for recent analysis methods. Thus, this first part of a two-part publication series of the German joint project EasyGSH-DB describes annual bathymetric digital terrain models at a 10 m gridded resolution for the German North Sea coast and German Bight from 1996 to 2016 (Sievers et al., 2020a, https://doi.org/10.48437/02.2020.K2.7000.0001), as well as surface sedimentological models of discretized cumulative grain size distribution functions for 1996, 2006, and 2016 on 100 m grids (Sievers et al., 2020b, https://doi.org/10.48437/02.2020.K2.7000.0005). Furthermore, basic morphodynamic and sedimentological processing analyses, such as the estimation of, for example, bathymetric stability or surface maps of sedimentological parameters, are provided (Sievers et al., 2020a, b, see respective download links).

Description: The Phanerozoic tectonothermal evolution of the SW slope of the East European Platform (EEP) in Poland is reconstructed by means of thermal maturity, low-temperature thermochronometry, and thermal modelling. We provide a set of new thermochronometric data and integrate stratigraphic and thermal maturity information to constrain the burial and thermal history of sediments. Apatite fission track (AFT) analysis and zircon (U-Th)/He (ZHe) thermochronology have been carried out on samples of sandstones, bentonites, diabase, and crystalline basement rocks collected from 17 boreholes located in central and NE Poland. They penetrated sedimentary cover of the EEP subdivided from the north to south into the Baltic, Podlasie, and Lublin basins. The average ZHe ages from Proterozoic basement rocks as well as Ordovician to Silurian bentonites and Cambrian to lower Carboniferous sandstones range from 848 ± 81 to 255 ± 22 Ma with a single early Permian age of 288 Ma, corresponding to cooling after a thermal event. The remaining ZHe ages represent partial reset or source ages. The AFT ages of samples are dispersed in the range of 235.8 ± 17.3 Ma (Middle Triassic) to 42.1 ± 11.1 Ma (Paleogene) providing a record of Mesozoic and Cenozoic cooling. The highest frequency of the AFT ages is in the Jurassic and Early Cretaceous prior to Alpine basin inversion. Thermal maturity results are consistent with the SW-ward increase of the Paleozoic and Mesozoic sediments thickness. An important break in a thermal maturity profile exists across the base Permian–Mesozoic unconformity. Thermal modelling showed that significant heating of Ediacaran to Carboniferous sedimentary successions occurred before the Permian with maximum paleotemperatures in the earliest and latest Carboniferous for Baltic–Podlasie and Lublin basins, respectively. The results obtained suggest an important role of early Carboniferous uplift and exhumation at the SW margin of the EEP. The SW slope of the latter was afterward overridden in the Lublin Basin by the Variscan orogenic wedge. Its tectonic loading interrupted Carboniferous uplift and caused resumption of sedimentation in the late Viséan. Consequently, a thermal history of the Lublin Basin is different from that in the Podlasie and Baltic basins but similar to other sections of the Variscan foreland, characterized by maximum burial at the end of Carboniferous. The Mesozoic thermal history was characterized by gradual cooling from peak temperatures at the transition from Triassic to Jurassic due to decreasing heat flow. Burial caused maximum paleotemperatures in the SW part of the study area, where the EEP was covered by an extensive sedimentary pile. However, further NE, due to low temperatures caused by shallow burial, the impact of fluids can be detected by vitrinite reflectance, illite/smectite, and thermochronological data. Our new results emphasize the importance of using multiple low-temperature thermochronometers and thermal modelling in connection with thermal maturity analysis to elucidate the near-surface evolution of platform margins.

Description: Snow stands out from materials at the Earth’s surface owing to its unique optical properties. Snow optical properties are sensitive to the snow microstructure, triggering potent climate feedbacks. The impacts of snow microstructure on its optical properties such as reflectance are, to date, only partially understood. However, precise modelling of snow reflectance, particularly bidirectional reflectance, are required in many problems, e.g. to correctly process satellite data over snow-covered areas. This study presents a dataset that combines bidirectional reflectance measurements over 500–2500 nm and the X-ray tomography of the snow microstructure for three snow samples of two different morphological types. The dataset is used to evaluate the stereological approach from Malinka (2014) that relates snow optical properties to the chord length distribution in the snow microstructure. The mean chord length and specific surface area (SSA) retrieved with this approach from the albedo spectrum and those measured by the X-ray tomography are in excellent agreement. The analysis of the 3D images has shown that the random chords of the ice phase obey the gamma distribution with the shape parameter m taking the value approximately equal to or a little greater than 2. For weak and intermediate absorption (high and medium albedo), the simulated bidirectional reflectances reproduce the measured ones accurately but tend to slightly overestimate the anisotropy of the radiation. For such absorptions the use of the exponential law for the ice chord length distribution instead of the one measured with the X-ray tomography does not affect the simulated reflectance. In contrast, under high absorption (albedo of a few percent), snow microstructure and especially facet orientation at the surface play a significant role in the reflectance, particularly at oblique viewing and incidence.

Description: Attribution in its general definition aims to quantify drivers of change in a system. According to IPCC Working Group II (WGII) a change in a natural, human or managed system is attributed to climate change by quantifying the difference between the observed state of the system and a counterfactual baseline that characterizes the system's behavior in the absence of climate change, where “climate change refers to any long-term trend in climate, irrespective of its cause” (IPCC, 2014). Impact attribution following this definition remains a challenge because the counterfactual baseline, which characterizes the system behavior in the hypothetical absence of climate change, cannot be observed. Process-based and empirical impact models can fill this gap as they allow us to simulate the counterfactual climate impact baseline. In those simulations, the models are forced by observed direct (human) drivers such as land use changes, changes in water or agricultural management but a counterfactual climate without long-term changes. We here present ATTRICI (ATTRIbuting Climate Impacts), an approach to construct the required counterfactual stationary climate data from observational (factual) climate data. Our method identifies the long-term shifts in the considered daily climate variables that are correlated to global mean temperature change assuming a smooth annual cycle of the associated scaling coefficients for each day of the year. The produced counterfactual climate datasets are used as forcing data within the impact attribution setup of the Inter-Sectoral Impact Model Intercomparison Project (ISIMIP3a). Our method preserves the internal variability of the observed data in the sense that factual and counterfactual data for a given day have the same rank in their respective statistical distributions. The associated impact model simulations allow for quantifying the contribution of climate change to observed long-term changes in impact indicators and for quantifying the contribution of the observed trend in climate to the magnitude of individual impact events. Attribution of climate impacts to anthropogenic forcing would need an additional step separating anthropogenic climate forcing from other sources of climate trends, which is not covered by our method.

Description: Biocides used in film protection products leaching from facades are known to be a potential threat to the environment. This study identifies individual sources and entry pathways in a small-scale urban area. We investigate emissions of commonly used biocides (terbutryn, diuron, and octylisothiazolinone – OIT) and some of their transformation products (TPs; diuron-desmethyl, terbumeton, terbuthylazine-2-hydroxy, and terbutryn-desethyl) from a 2 ha residential area 13 years after construction has ended. Sampling utilizes existing urban water infrastructure representative for decentralized storm water management in central and northern Europe and applies a two-step approach to (a) determine the occurrence of biocides above water quality limits (i.e., predicted no-effect concentration, PNEC) and (b) identify source areas and characterize entry pathways into surface and groundwater. Monitoring focuses on the analysis of selected biocides and TPs by liquid chromatography–mass spectrometry/mass spectrometry (LC-MS/MS) in water samples taken from facades, rainwater pipes, drainage, and storm water infiltration systems. In standing water in a swale, we found high concentrations of diuron (174 ng L−1) and terbutryn (40 ng L−1) above PNEC for surface water. We confirmed expected sources, i.e., facades. Sampling of rain downpipes from flat roofs identified additional sources of all biocides and two TPs of terbutryn and one TP of diuron. Diuron and terbutryn were found in three drainage pipes representing different entry pathways of biocides. In one drainage pipe collecting road runoff, only diuron-desmethyl and terbutryn-desethyl were detected. In two other drainage pipes collecting infiltrated water through soil, terbuthylazine-2-hydroxy was additionally detected. One of the pipes collecting infiltrated water through soil concentration showed the highest concentrations of terbutryn and two of its TPs (terbutryn-desethyl and terbuthylazine-2-hydroxy). This suggests a high leaching potential of terbutryn. The applied two-step approach determined sources and pathways of biocide and their TPs. This study contributes to expanding knowledge on their entry and distribution and, thus, eventually towards reducing emissions.

Description: Regional anomalies of steric sea level are either due to redistribution of heat and freshwater anomalies or due to ocean–atmosphere buoyancy fluxes. Interannual to decadal variability in sea level across the tropical Pacific is mainly due to steric variations driven by wind stress anomalies. The importance of air–sea buoyancy fluxes is less clear. We use a global, eddy-permitting ocean model and a series of sensitivity experiments with quasi-climatological momentum and buoyancy fluxes to identify the contribution of buoyancy fluxes for interannual to decadal sea level variability in the tropical Pacific. We find their contribution on interannual timescales to be strongest in the central tropical Pacific at around a 10∘ latitude in both hemispheres and also relevant in the very east of the tropical domain. Buoyancy-flux-forced anomalies are correlated with variations driven by wind stress changes, but their effect on the prevailing anomalies and the importance of heat and freshwater fluxes vary locally. In the eastern tropical basin, interannual sea level variability is amplified by anomalous heat fluxes, while the importance of freshwater fluxes is small, and neither has any impact on decadal timescales. In the western tropical Pacific, the variability on interannual and decadal timescales is dampened by both heat and freshwater fluxes. The mechanism involves westward-propagating Rossby waves that are triggered during El Niño–Southern Oscillation (ENSO) events by anomalous buoyancy fluxes in the central tropical Pacific and counteract the prevailing sea level anomalies once they reach the western part of the basin.

Description: The Nioghalvfjerdsfjorden glacier (also known as the 79∘ North Glacier) drains approximately 8 % of the Greenland Ice Sheet. Supraglacial lakes (SGLs), or surface melt ponds, are a persistent summertime feature and are thought to drain rapidly to the base of the glacier and influence seasonal ice velocity. However, seasonal development and spatial distribution of SGLs in the north-east of Greenland are poorly understood, leaving a substantial error in the estimate of meltwater and its impacts on ice velocity. Using results from an automated detection of melt ponds, atmospheric and surface mass balance modelling, and reanalysis products, we investigate the role of specific climatic conditions in melt onset, extent, and duration from 2016 to 2019. The summers of 2016 and 2019 were characterised by above-average air temperatures, particularly in June, as well as a number of rainfall events, which led to extensive melt ponds to elevations up to 1600 m. Conversely, 2018 was particularly cold, with a large accumulated snowpack, which limited the development of lakes to altitudes less than 800 m. There is evidence of inland expansion and increases in the total area of lakes compared to the early 2000s, as projected by future global warming scenarios.

Description: Coincident wide-angle and multi-channel seismic data acquired within the scope of the PAMELA Moz3-5 project allow us to reconsider the formation mechanism of East African margins offshore of southern Mozambique. This study specifically focuses on the sedimentary and deep-crustal architecture of the Limpopo margin (LM) that fringes the eastern edge of the Mozambique’s Coastal Plain (MCP) and its offshore southern prolongation the North Natal Valley (NNV). It relies primarily on the MZ3 profile that runs obliquely from the northeastern NNV towards the Mozambique basin (MB) with additional inputs from a tectonostratigraphy analysis of industrial onshore–offshore seismic lines and nearby or crossing velocity models from companion studies. Over its entire N–S extension the LM appears segmented into (1) a western domain that shows the progressive eastward crustal thinning and termination of the MCP/NNV continental crust and its overlying pre-Neocomian volcano-sedimentary basement and (2) a central corridor of anomalous crust bounded to the east by the Mozambique fracture zone (MFZ) and the oceanic crust of the MB. A prominent basement high marks the boundary between these two domains. Its development was most probably controlled by a steep and deeply rooted fault, i.e., the Limpopo fault. We infer that strike-slip or slightly transtensional rifting occurred along the LM and was accommodated along this Limpopo fault. At depth we propose that ductile shearing was responsible for the thinning of the continental crust and an oceanward flow of lower crustal material. This process was accompanied by intense magmatism that extruded to form the volcanic basement and gave the corridor its peculiar structure and mixed nature. The whole region remained at a relative high level during the rifting period and a shallow marine environment dominated the pre-Neocomian period during the early phase of continent–ocean interaction. It is only some time after break-up in the MB and the initiation of the MFZ that decoupling occurred between the MCP/NNV and the corridor, allowing for the latter to subside and become covered by deep marine sediments. A scenario for the early evolution and formation of the LM is proposed taking into account both recent kinematic and geological constraints. It implies that no or little change in extensional direction occurred between the intra-continental rifting and subsequent phase of continent–ocean interaction.

Description: Models are an important tool to predict Earth system dynamics. An accurate prediction of future states of ecosystems depends on not only model structures but also parameterizations. Model parameters can be constrained by data assimilation. However, applications of data assimilation to ecology are restricted by highly technical requirements such as model-dependent coding. To alleviate this technical burden, we developed a model-independent data assimilation (MIDA) module. MIDA works in three steps including data preparation, execution of data assimilation, and visualization. The first step prepares prior ranges of parameter values, a defined number of iterations, and directory paths to access files of observations and models. The execution step calibrates parameter values to best fit the observations and estimates the parameter posterior distributions. The final step automatically visualizes the calibration performance and posterior distributions. MIDA is model independent, and modelers can use MIDA for an accurate and efficient data assimilation in a simple and interactive way without modification of their original models. We applied MIDA to four types of ecological models: the data assimilation linked ecosystem carbon (DALEC) model, a surrogate-based energy exascale earth system model: the land component (ELM), nine phenological models and a stand-alone biome ecological strategy simulator (BiomeE). The applications indicate that MIDA can effectively solve data assimilation problems for different ecological models. Additionally, the easy implementation and model-independent feature of MIDA breaks the technical barrier of applications of data–model fusion in ecology. MIDA facilitates the assimilation of various observations into models for uncertainty reduction in ecological modeling and forecasting.

Description: In this study we present results from the factor analysis of 43 aerosol mass spectrometer (AMS) datasets and provide an overview of worldwide organic aerosol (OA) components and their evolution in the atmosphere. At most sites, the OA can be separated into oxygenated OA (OOA), hydrocarbon-like OA (HOA), and sometimes other components such as biomass burning OA (BBOA). In many analyses, the OOA can be further deconvolved into low-volatility OOA (LV-OOA) and semi-volatile OOA (SV-OOA). A wide range of f44 (ratio of m/z 44 to total signal in the component mass spectrum) and O:C ratios are observed for both LV-OOA (0.17±0.04, 0.73±0.14) and SV-OOA (0.07±0.04, 0.35±0.14) components, reflecting the fact that there is a continuum of OOA properties in ambient aerosol. Differences in the mass spectra of these components are characterized in terms of the two main ions m/z 44 (CO2+) and m/z 43 (mostly C2H3O+). The LV-OOA component spectra have higher f44 and lower f43 than SV-OOA. The OOA components (OOA, LV-OOA, and SV-OOA) from all sites cluster within a well defined triangular region in the f44 vs. f43 space, which can be used as a standardized means of comparing and characterizing any OOA components (laboratory or ambient) observed with the AMS. Examination of the OOA components in this triangular space indicates that OOA component spectra become increasingly similar to each other and to fulvic acid and HULIS sample spectra as f44 (a surrogate for O:C and an indicator of photochemical aging) increases. This indicates that ambient OA converges towards highly aged LV-OOA with atmospheric oxidation. The common features of the transformation between SV-OOA and LV-OOA at multiple sites potentially enables a simplified description of the oxidation of OA in the atmosphere. Comparison of laboratory SOA data with ambient OOA indicates that laboratory SOA are more similar to SV-OOA, and rarely become as oxidized as ambient LV-OOA, likely due to the higher loadings employed in the experiments and/or limited oxidant exposure in most chamber experiments.

Description: Airborne lidar and in-situ measurements of the aerosol properties were conducted between Svalbard Island and Scandinavia in April 2008. Evidence of aerosol transport from Europe and Asia is given. The analysis of the aerosol optical properties based on a multiwavelength lidar (355, 532, 1064 nm) including depolarization at 355 nm aims at distinguishing the role of the different aerosol sources (Siberian wild fires, Eastern Asia and European anthropogenic emissions). Combining, first aircraft measurements, second FLEXPART simulations with a calculation of the PBL air fraction originating from the three different mid-latitude source regions, and third level-2 CALIPSO data products (i.e. backscatter coefficient, depolarisation and color ratio in aerosol layers) along the transport pathways, appears a valuable approach to identify the role of the different aerosol sources even after a transport time larger than 4 days. Above Asia, CALIPSO data indicate more depolarisation (up to 15%) and largest color ratio (〉0.5) for the northeastern Asia emissions (i.e. an expected mixture of Asian pollution and dust), while low depolarisation together with smaller and quasi constant color ratio (≈0.3) are observed for the Siberian biomass burning emissions. A similar difference is visible between two layers observed by the aircraft above Scandinavia. The analysis of the time evolution of the aerosol optical properties revealed by CALIPSO between Asia and Scandinavia shows a gradual decrease of the aerosol backscatter, depolarisation ratio and color ratio which suggests the removal of the largest particles in the accumulation mode. A similar study conducted for a European plume has shown aerosol optical properties intermediate between the two Asian sources with color ratio never exceeding 0.4 and moderate depolarisation ratio being always less than 8%, i.e. less aerosol from the accumulation mode.

Description: Secondary organic aerosol (SOA) affects atmospheric composition, air quality and radiative transfer. However major difficulties are encountered in the development of reliable models for SOA formation. Constraints on processes involved in SOA formation can be obtained by interpreting the speciation and evolution of organics in the gaseous and condensed phase simultaneously. In this study we investigate SOA formation from dark α-pinene ozonolysis with particular emphasis upon the mass distribution of gaseous and particulate organic species. A detailed model for SOA formation is compared with the results from experiments performed in the EUropean PHOtoREactor (EUPHORE) simulation chamber, including on-line gas-phase composition obtained from Chemical-Ionization-Reaction Time-Of-Flight Mass-Spectrometry measurements, and off-line analysis of SOA samples performed by Electrospray Ionisation Ion Trap Mass Spectrometry. The temporal profile of SOA mass concentration is relatively well reproduced by the model. Sensitivity analysis highlights the importance of the choice of vapour pressure estimation method. Comparisons of the simulated gaseous- and condensed-phase mass distributions with those observed show a generally good agreement. The simulated speciation has been used to (i) propose a chemical structure for the principal gaseous semi-volatile organic compounds and condensed monomer organic species and (ii) explore the possible contribution of a range of accretion reactions occurring in the condensed phase. We find that oligomer formation through esterification reactions gives the best agreement between the observed and simulated mass spectra.

Description: The first aerosol indirect effect over a clean, northern high-latitude site was investigated by determining the aerosol cloud interaction (ACI) using three different approaches; ground-based in situ measurements, combined ground-based in situ measurements and satellite retrievals and using only satellite retrievals. The obtained values of ACI were highest for in situ ground-based data, clearly lower for combined ground-based and satellite data, and lowest for data relying solely on satellite retrievals. One of the key findings of this study was the high sensitivity of ACI to the definition of the aerosol burden. We showed that at least a part of the variability in ACI can be explained by how different investigators have related different cloud properties to "aerosol burden".

Description: Methyl chloride (CH3Cl) is a chlorine-containing trace gas in the atmosphere contributing significantly to stratospheric ozone depletion. Large uncertainties in estimates of its source and sink magnitudes and temporal and spatial variations currently exist. GEIA inventories and other bottom-up emission estimates are used to construct a priori maps of the surface fluxes of CH3Cl. The Model of Atmospheric Transport and Chemistry (MATCH), driven by NCEP interannually varying meteorological data, is then used to simulate CH3Cl mole fractions and quantify the time series of sensitivities of the mole fractions at each measurement site to the surface fluxes of various regional and global sources and sinks. We then implement the Kalman filter (with the unit pulse response method) to estimate the surface fluxes on regional/global scales with monthly resolution from January 2000 to December 2004. High frequency observations from the AGAGE, SOGE, NIES, and NOAA/ESRL HATS in situ networks and low frequency observations from the NOAA/ESRL HATS flask network are used to constrain the source and sink magnitudes. The inversion results indicate global total emissions around 4100±470 Gg yr−1 with very large emissions of 2200±390 Gg yr−1 from tropical plants, which turn out to be the largest single source in the CH3Cl budget. Relative to their a priori annual estimates, the inversion increases global annual fungal and tropical emissions, and reduces the global oceanic source. The inversion implies greater seasonal and interannual oscillations of the natural sources and sink of CH3Cl compared to the a priori. The inversion also reflects the strong effects of the 2002/2003 globally widespread heat waves and droughts on global emissions from tropical plants, biomass burning and salt marshes, and on the soil sink.

Description: The goal of this research is to quantify diesel- and gasoline-powered motor vehicle emissions within the Mexico City Metropolitan Area (MCMA) using on-road measurements captured by a mobile laboratory combined with positive matrix factorization (PMF) receptor modeling. During the MCMA-2006 ground-based component of the MILAGRO field campaign, the Aerodyne Mobile Laboratory (AML) measured many gaseous and particulate pollutants, including carbon dioxide, carbon monoxide (CO), nitrogen oxides (NOx), benzene, toluene, alkylated aromatics, formaldehyde, acetaldehyde, acetone, ammonia, particle number, fine particulate mass (PM2.5), and black carbon (BC). These serve as inputs to the receptor model, which is able to resolve three factors corresponding to gasoline engine exhaust, diesel engine exhaust, and the urban background. Using the source profiles, we calculate fuel-based emission factors for each type of exhaust. The MCMA's gasoline-powered vehicles are considerably dirtier, on average, than those in the US with respect to CO and aldehydes. Its diesel-powered vehicles have similar emission factors of NOx and higher emission factors of aldehydes, particle number, and BC. In the fleet sampled during AML driving, gasoline-powered vehicles are responsible for 97% of mobile source emissions of CO, 22% of NOx, 95–97% of aromatics, 72–85% of carbonyls, 74% of ammonia, negligible amounts of particle number, 26% of PM2.5, and 2% of BC; diesel-powered vehicles account for the balance. Because the mobile lab spent 17% of its time waiting at stoplights, the results may overemphasize idling conditions, possibly resulting in an underestimate of NOx and overestimate of CO emissions. On the other hand, estimates of the inventory that do not correctly account for emissions during idling are likely to produce bias in the opposite direction. Nevertheless, the fuel-based inventory suggests that mobile source emissions of CO and NOx are overstated in the official inventory while emissions of VOCs may be understated. For NOx, the fuel-based inventory is lower for gasoline-powered vehicles but higher for diesel-powered ones compared to the official inventory.

Description: A cross-platform field campaign, OP3, was conducted in the state of Sabah in Malaysian Borneo between April and July of 2008. Among the suite of observations recorded, the campaign included measurements of NOx and O3–crucial outputs of any model chemistry mechanism. We describe the measurements of these species made from both the ground site and aircraft. We examine the output from the global model p-TOMCAT at two resolutions for this location during the April campaign period. The models exhibit reasonable ability in capturing the NOx diurnal cycle, but ozone is overestimated. We use a box model containing the same chemical mechanism to explore the weaknesses in the global model and the ability of the simplified global model chemical mechanism to capture the chemistry at the rainforest site. We achieve a good fit to the data for all three species (NO, NO2, and O3), though the model is much more sensitive to changes in the treatment of physical processes than to changes in the chemical mechanism. Indeed, without some parameterization of the nighttime boundary layer-free troposphere mixing, a time dependent box model will not reproduce the observations. The final simulation uses this mixing parameterization for NO and NO2 but not O3, as determined by the vertical structure of each species, and matches the measurements well.

Description: The elemental composition of laboratory chamber secondary organic aerosol (SOA) from glyoxal uptake, α-pinene ozonolysis, isoprene photooxidation, single-ring aromatic photooxidation, and naphthalene photooxidation is evaluated using Aerodyne high-resolution time-of-flight mass spectrometer data. SOA O/C ratios ranged from 1.13 for glyoxal to 0.30–0.43 for α-pinene ozonolysis. The elemental composition of α-pinene and naphthalene SOA was also confirmed by offline mass spectrometry. The fraction of organic signal at m/z 44 is generally a good measure of SOA oxygenation for all systems except for glyoxal uptake, in which m/z 44 substantially underpredicts O/C. Although chamber SOA has generally been considered less oxygenated than ambient SOA, single-ring aromatic- and naphthalene-derived SOA can reach O/C ratios upward of 0.7, well within the range of OOA, though still not as high as some ambient measurements. The spectra of aromatic- and isoprene-high-NOx SOA resemble that of OOA, but the spectrum of glyoxal uptake does not resemble that of any PMF component.

Description: This paper presents estimates of the spectral solar absorption due to atmospheric aerosols during the 2006 MILAGRO/INTEX-B (Megacity Initiative-Local And Global Research Observations/Phase B of the Intercontinental Chemical Transport Experiment) field campaign. The aerosol absorption was derived from measurements of the spectral solar radiation and the spectral aerosol optical depth made on the J31 aircraft flying over the Gulf of Mexico and over Mexico City. We present the spectral single scattering albedo (SSA) and aerosol absorption optical depth (AAOD) for two flights over the Gulf of Mexico and three flights over Mexico City for wavelengths from 350 to approximately 1650 nm. The spectral aerosol optical properties of each case are different and illustrate the variability of the aerosol optical properties in the Mexico City area. The results can be described in terms of three different wavelength region: The 350–500 nm region where the aerosol absorption often falls off sharply presumably due to organic carbonaceous particles and windblown dust; the 500–1000 nm region where the decrease with wavelength is slower presumably due to black carbon; and the near infrared spectral region (1000 nm to 1650 nm) where it is difficult to obtain reliable results since the aerosol absorption is relatively small and the gas absorption dominates. However, there is an indication of a small and somewhat wavelength independent absorption in the region beyond 1000 nm. For one of the flights over the Gulf of Mexico near the coastline it appears that a cloud/fog formation and evaporation led to an increase of absorption possibly due to a water shell remaining on the particles after the cloud/fog had dissipated. For two of the Mexico City cases, the single scattering albedo is roughly constant between 350–500 nm consistent with other Mexico City results. In three of the cases a single absorption Angstrom exponent (AAE) fits the aerosol absorption optical depth over the entire wavelength range of 350 to 1650 nm relatively well (r2〉0.86).

Description: To understand and quantify the impact of local, regional, and distant atmospheric mercury sources to rainfall mercury deposition in the Pensacola, Florida watershed, a program of event-based rainfall sampling was started in late 2004. Modified Aerochem-Metrics wet/dry rainfall samplers were deployed at three sites in the region around the Crist coal-fired power plant and event-based samples were collected continuously for three years. Samples were analyzed for total mercury, volatile elements As, Se, and Sn, and a suite of trace elements including Al, Bi, Cd, Ce, Co, Cr, Cs, Cu, Dy, Er, Eu, Fe, Ho, Ga, Gd, La, Li, Lu, Mg, Mn, Na, Nb, Nd, Ni, P, Pb, Pr, Rb, Sb, Sc, Si, Sm, Sr, Tb, Th, Ti, Tm, U, V, Y, Yb, and Zn. Nutrients and major ions were also measured on each sample. Multivariate statistical methods are used to sort these tracers into factors that represent potential source components that contribute to the rainfall chemistry. Hg, As, Se, Sn, Sb, and non sea-salt sulfate were all significantly correlated with one anthropogenic factor. Using various Hg/element ratios, we can estimate that 22–33% of the rainfall mercury in the region results from coal combustion.

Description: Estimation of a trend of an atmospheric state variable is often performed by fitting a linear regression line to a set of data of this variable sampled at different times. Often these data are irregularly sampled in space and time and clustered in a sense that error correlations among data points cause a similar error of data points sampled at similar times. Since this can affect the estimated trend, we suggest to take the full error covariance matrix of the data into account. Superimposed periodic variations can be jointly fitted in a straight forward manner, even if the shape of the periodic function is not known. Global data sets, particularly satellite data, can form the basis to estimate the error correlations.

Description: We compare Tropospheric Emission Spectrometre (TES) versions 3 and 4, V003 and V004, respectively, nadir-stare ozone profiles with ozonesonde profiles from the Arctic Intensive Ozonesonde Network Study (ARCIONS, http://http://croc.gsfc.nasa.gov/arcions/) during the Arctic Research on the Composition of the Troposphere from Aircraft and Satellites (ARCTAS) field mission. The ozonesonde data are from launches timed to match the Aura's overpass, where 11 coincidences spanned 44° N to 71° N from April to July 2008. Using the TES "stare" observation mode, 32 observations are taken over each coincident ozonesonde launch. By effectively sampling the same air mass 32 times, comparisons are made between the empirically-calculated random errors to the expected random errors from measurement noise, temperature and interfering species, such as water. This study represents the first validation of high latitude (〉60°) TES ozone. We find that the calculated errors are consistent with the actual errors with a similar vertical distribution that varies between 5% and 20% for V003 and V004 TES data. In general, TES ozone profiles are positively biased (by less than 15%) from the surface to the upper troposphere (~1000 to 100 hPa) and negatively biased (by less than 20%) from the upper troposphere to the lower stratosphere (100 to 30 hPa) when compared to the ozonesonde data. Lastly, for V003 and V004 TES data between 44° N and 71° N there is a small variability in the mean biases (from −14 to +15%), mean theoretical errors (from 6 to 13%), and mean random errors (from 9 to 19%).

Description: The Antarctic near-coastal sub-micrometre aerosol particle features in summer were characterised based on measured data on aerosol hygroscopicity, size distributions, volatility and chemical ion and organic carbon mass concentrations. Hysplit model was used to calculate the history of the air masses to predict the particle origin. Additional measurements of meteorological parameters were utilised. The hygroscopic properties of particles mostly resembled those of marine aerosols. The measurements took place at 130 km from the Southern Ocean, which was the most significant factor affecting the particle properties. This is explained by the lack of additional sources on the continent of Antarctica. The Southern Ocean was thus a likely source of the particles and nucleating and condensing vapours. The particles were very hygroscopic (HGF 1.75 at 90 nm) and very volatile. Most of the sub-100 nm particle volume volatilised below 100 °C. Based on chemical data, particle hygroscopic and volatile properties were explained by a large fraction of non-neutralised sulphuric acid together with organic material. The hygroscopic growth factors assessed from chemical data were similar to measured. Hygroscopicity was higher in dry continental air masses compared with the moist marine air masses. This was explained by the aging of the marine organic species and lower methanesulphonic acid volume fraction together with the changes in the inorganic aerosol chemistry as the aerosol had travelled long time over the continental Antarctica. Special focus was directed in detailed examination of the observed new particle formation events. Indications of the preference of negative over positive ions in nucleation could be detected. However, in a detailed case study, the neutral particles dominated the particle formation process. Freshly nucleated particles had the smallest hygroscopic growth factors, which increased subsequent to particle aging.

Description: We present carbon dioxide (CO2) estimates from the Tropospheric Emission Spectrometer (TES) on the EOS-Aura satellite launched in 2004. For observations between 40° S and 45° N, we find about 1 degree of freedom with peak sensitivity at 511 hPa. The estimated error is ~10 ppm for a single target and about 1.3 ppm for monthly averages on spatial scales of 20°×30°. Monthly spatially-averaged TES results from 2005–2008 processed with a uniform initial guess and prior are compared to CONTRAIL aircraft data over the Pacific ocean, aircraft data at the Southern Great Plains (SGP) ARM site in the southern US, and the Mauna Loa and Samoa surface stations. Comparisons to Mauna Loa observatory show a correlation of 0.92, a standard deviation of 1.3 ppm, a predicted error of 1.2 ppm, and a ~2% low bias, which is subsequently corrected, and comparisons to SGP aircraft data over land show a correlation of 0.67 and a standard deviation of 2.3 ppm. TES data between 40° S and 45° N for 2006–2007 are compared to surface flask data, GLOBALVIEW, the Atmospheric Infrared Sounder (AIRS), and CarbonTracker. Comparison to GLOBALVIEW-CO2 ocean surface sites shows a correlation of 0.60 which drops when TES is offset in latitude, longitude, or time. At these same locations, TES shows a 0.62 and 0.67 correlation to CarbonTracker with TES observation operator at the surface and 5 km, respectively. We also conducted an observing system simulation experiment to assess the potential utility of the TES data for inverse modeling of CO2 fluxes. We find that if biases in the data and model are well characterized, the averaged data have the potential to provide sufficient information to significantly reduce uncertainty on annual estimates of regional CO2 sources and sinks. Averaged pseudo-data at 10°×10° reduced uncertainty in flux estimates by as much as 70% for some tropical regions.

Description: The spatial distribution of aerosol chemical composition and the evolution of the Organic Aerosol (OA) fraction is investigated based upon airborne measurements of aerosol chemical composition in the planetary boundary layer across Europe. Sub-micron aerosol chemical composition was measured using a compact Time-of-Flight Aerosol Mass Spectrometer (cToF-AMS). A range of sampling conditions were evaluated, including relatively clean background conditions, polluted conditions in North-Western Europe and the near-field to far-field outflow from such conditions. Ammonium nitrate and OA were found to be the dominant chemical components of the sub-micron aerosol burden, with mass fractions ranging from 20–50% each. Ammonium nitrate was found to dominate in North-Western Europe during episodes of high pollution, reflecting the enhanced NOx and ammonia sources in this region. OA was ubiquitous across Europe and concentrations generally exceeded sulphate by 50–100%. A factor analysis of the OA burden was performed in order to probe the evolution across this large range of spatial and temporal scales. Two separate Oxygenated Organic Aerosol (OOA) components were identified; one representing an aged-OOA, termed Low Volatility-OOA and another representing fresher-OOA, termed Semi Volatile-OOA on the basis of their mass spectral similarity to previous studies. The factors derived from different flights were not chemically the same but rather reflect the range of OA composition sampled during a particular flight. Significant chemical processing of the OA was observed downwind of major sources in North-Western Europe, with the LV-OOA component becoming increasingly dominant as the distance from source and photochemical processing increased. The measurements suggest that the aging of OA can be viewed as a continuum, with a progression from a less oxidised, semi-volatile component to a highly oxidised, less-volatile component. Substantial amounts of pollution were observed far downwind of continental Europe, with OA and ammonium nitrate being the major constituents of the sub-micron aerosol burden. Such anthropogenically perturbed air masses can significantly perturb regional climate far downwind of major source regions.

Description: Results from a regional air quality forecast model, AIRPACT-3, are compared to OMI tropospheric NO2 integrated column densities for an 18 month period over the Pacific Northwest. AIRPACT column densities were well correlated with cloud-free monthly averages of tropospheric NO2 (R=0.75) to NASA retrievals for months without wildfires, but were poorly correlated with significant model overpredictions (R=0.21) for months with wildfires when OMI and AIRPACT were compared over the entire domain. AIRPACT forecasted higher NO2 in some US urban areas, and lower NO2 in many Canadian urban areas, when compared to OMI. There are significant changes in results after spatially averaging model results to the daily OMI swath. Also, it is shown that applying the averaging kernel to model results in cloudy conditions has a large effect, but applying the averaging kernel in cloud free conditions has little effect. The KNMI and NASA retrievals of tropospheric NO2 from OMI (collection 3) are compared. The NASA product is shown to be significantly different than the KNMI tropospheric NO2 product, i.e. July 2007 (R=0.60) and January 2008 (R=0.69).

Description: We present a kinetic double layer model coupling aerosol surface and bulk chemistry (K2-SUB) based on the PRA framework of gas-particle interactions (Pöschl et al., 2007). K2-SUB is applied to a popular model system of atmospheric heterogeneous chemistry: the interaction of ozone with oleic acid. We show that our modelling approach allows de-convoluting surface and bulk processes, which has been a controversial topic and remains an important challenge for the understanding and description of atmospheric aerosol transformation. In particular, we demonstrate how a detailed treatment of adsorption and reaction at the surface can be coupled to a description of bulk reaction and transport that is consistent with traditional resistor model formulations. From literature data we have derived a consistent set of kinetic parameters that characterise mass transport and chemical reaction of ozone at the surface and in the bulk of oleic acid droplets. Due to the wide range of rate coefficients reported from different experimental studies, the exact proportions between surface and bulk reaction rates remain uncertain. Nevertheless, the model results suggest an important role of chemical reaction in the bulk and an upper limit of ~1×10−11 cm2 s−1 for the surface reaction rate coefficient. Sensitivity studies show that the surface accommodation coefficient of the gas-phase reactant has a strong non-linear influence on both surface and bulk chemical reactions. We suggest that K2-SUB may be used to design, interpret and analyse future experiments for better discrimination between surface and bulk processes in the oleic acid-ozone system as well as in other heterogeneous reaction systems of atmospheric relevance.

Description: Vertical number fluxes of aerosol particles and vertical fluxes of CO2 were measured with the eddy covariance method at the top of a 53 m high tower in the Amazon rain forest as part of the LBA (The Large Scale Biosphere Atmosphere Experiment in Amazonia) experiment. The observed aerosol number fluxes included particles with sizes down to 10 nm in diameter. The measurements were carried out during the wet and dry season in 2008. In this study focus is on the dry season aerosol fluxes, with significant influence from biomass burning, and these are compared with aerosol fluxes measured during the wet season. The primary goal is to quantify the dry deposition sink and to investigate whether particle deposition velocities change when going from the clean wet season into the more polluted dry season. Furthermore, it is tested whether the rain forest is always a net sink of particles in terms of number concentrations, or if particle emission from the surface under certain circumstances may dominate over the dry deposition sink. The particle deposition velocity vd increased linearly with increasing friction velocity in both seasons and the relations are described by vdd=(2.7 u* −0.2)×10−3 (dry season) and vdw=2.5 u*×10−3 (wet season), where u* is the friction velocity. The fact that the two relations are very similar to each other indicates that the seasonal change in aerosol number size distribution is not enough for causing any significant change in deposition velocity. In general, particle deposition velocities in this study are low compared to studies over boreal forests. The reason is probably domination of accumulation mode particles in the Amazon boundary layer, both in the dry and wet season, and low wind speeds in the tropics compared to the midlatitudes. Net particle deposition fluxes prevailed in daytime in both seasons and the deposition flux was considerably larger in the dry season due to the much higher dry season particle concentration. In the dry season, nocturnal particle fluxes behaved very similar to the nocturnal CO2 fluxes. Throughout the night, the measured particle flux at the top of the tower was close to zero, but early in the morning there was an upward particle flux peak that is not likely a result of entrainment or local pollution. It is possible that these morning upward particle fluxes are associated with emission of natural biogenic particles from the rain forest. Emitted particles may be stored within the canopy during stable conditions at nighttime, similarly to CO2, and being released from the canopy when conditions become more turbulent in the morning.

Description: Different approaches to estimate organic sampling artifacts have been adopted by US long-term air quality monitoring networks. This study documents field blank (bQF) and backup filter (quartz-fiber behind quartz-fiber filter; QBQ) carbon levels for the: 1) Interagency Monitoring of Protected Visual Environments (IMPROVE); 2) Speciation Trends Network (STN; part of the Chemical Speciation Network [CSN]; and 3) Southeastern Aerosol Research and Characterization (SEARCH) networks. Filter pack sample handling procedures and blank correction methods are examined. Due to a relatively short (1–15 min) passive exposure period, STN/CSN and SEARCH network bQF organic carbon (OC; 0.8–1 μg/cm2) may underestimate positive and negative OC artifacts, respectively, resulting from passive adsorption or volatilization of volatile or semi-volatile organic compounds on quartz-fiber filters while they are in the sampler. This is evidenced by a lack of temporal or spatial variability and low bQF levels. With ~7 d of ambient passive exposure, average IMPROVE bQF and QBQ OC are comparable (2.4±0.5 and 3.1±0.8 μg/cm2) and more than twice those found in the STN/CSN and SEARCH networks. Lower STN/CSN flow rates and larger filter deposit areas result in 9–20% of the areal density (μg/cm2) compared to IMPROVE areal deposits. STN/CSN bQF values are 11–34% lower than linear regression intercepts derived from collocated IMPROVE-STN/CSN data pairs. Using a preceding organic denuder in the SEARCH network reduces the organic vapor adsorption on QBQ, though the longer QBQ period in the sampler may result from passive organic vapor adsorption as well as evaporated OC from the front filter deposits.

Description: HUmic-LIke Substances (HULIS) have been identified as major contributors to the organic carbon in atmospheric aerosol. The term HULIS is used to describe the organic material found in aerosol particles which resembles the humic organic material in river and sea water and in soils. In this study two sets of filter samples from atmospheric aerosols were collected at different sites. One sample was collected at the K-puszta rural site in Hungary, about 80 km SE of Budapest, and a second set of samples was collected at a site in Rondônia, Amazonia, Brazil, during the LBA-SMOCC biomass burning season experiment. HULIS were extracted from the samples, and their hygroscopic properties were studied using a Hygroscopicity Tandem Differential Mobility Analyzer (H-TDMA) at relative humidity (RH) 100%. The H-TDMA measurements were carried out at a dry diameter of 100 nm and for RH ranging from 30 to 98%. At 90% RH the HULIS samples showed diameter growth factors between 1.04 and 1.07, reaching values of 1.4 at 98% RH. The cloud nucleating properties of the two sets of aerosol samples were analyzed using two types of thermal static cloud condensation nucleus counters (CCNC). Two different parameterization models were used to investigate the potential effect of HULIS surface activity, both yielding similar results. For the K-puszta winter HULIS sample, the surface tension at the point of activation was estimated to be lowered by between 34% (47.7 mN/m) and 31% (50.3 mN/m) for dry sizes between 50 and 120 nm in comparison to pure water. A moderate lowering was also observed for the entire water soluble aerosol sample, including both organic and inorganic compounds, where the surface tension was decreased by between 2% (71.2 mN/m) and 13% (63.3 mN/m).

Description: Air quality in the outflow from urban centers affects millions of people, as well as, natural and managed ecosystems downwind. In locations where there are large sources of biogenic VOCs downwind of urban centers, the outflow is characterized by a high VOC reactivity due to biogenic emissions and low NOx. However most field and chamber studies have focused on limiting cases of high NOx or of near zero NOx. Recent measurements of a wide suite of VOCs, O3 and meteorological parameters at several locations within the Sacramento urban plume have provided a detailed benchmark for testing our understanding of chemistry in a plume transitioning from high NOx to low NOx and high VOC reactivity. As an additional simplification, the strong mountain valley circulation in the region makes this urban plume a physical realization of a nearly idealized Lagrangian plume. Here, we describe a model of this plume. We use a Lagrangian model representing chemistry based on the Master Chemical Mechanism (MCM) v3.1 along with mixing and deposition. We discuss the effects of entrainment of background air, the branching ratio for the production of isoprene nitrates and the effects of soil NOx emissions on the composition of the evolving plume. The model predicts that after 2–3 h of chemical processing only 45% of the peroxynitrates (ΣPNs) are PAN and that most (69%) RONO2 are secondary alkyl nitrate products of the reaction of OH with RONO2. We find the model is more consistent with the observations if: a) the yield of ΣPNs from large and multi-functional aldehydes is close to zero; and b) the reaction between OH and RONO2 produces multifunctional nitrates as opposed to either HNO3 or NO2 as is typical in most currently adopted reaction mechanisms. Model results also show that adding NOx emissions throughout the transect increases the available NOx in the downwind regions, but modeled ozone concentrations were little affected by the increased NOx.

Description: Data collected at four ground-level sites are analyzed (1) to determine the surface cloud radiative effect (CRE) induced by cirrus clouds at regional scale for shortwave (CRESW) and longwave (CRELW) fluxes and (2) to derive the sensitivity of surface CRESW to the cloud optical thickness (COT) modulated by the solar zenith angle and the atmospheric turbidity (noted CRESW*) and the sensitivity of surface CRELW to the infrared emissive power of cirrus cloud modulated by the water vapor content (noted CRELW*). The average CRESW* is −120 W m−2 COT−1 but it ranges from −80 to −140 m−2 COT−1 depending on the solar illumination with a residual variability ranges from +40 and −40 W m−2 COT−1 from pristine to turbid conditions, respectively. The CRELW*, that corresponds to the infrared transmissivity of the atmosphere, ranges from 3% to 40% from dry to wet atmospheric conditions, respectively. The subvisible cirrus class (COT

Description: In this study, the Community Multiscale Air Quality (CMAQ) modeling system is used to simulate the ozone (O3) episodes during the Program of Regional Integrated Experiments of Air Quality over the Pearl River Delta, China, in October 2004 (PRIDE-PRD2004). The simulation suggests that O3 pollution is a regional phenomenon in the PRD. Elevated O3 levels often occurred in the southwestern inland PRD, Pearl River estuary (PRE), and southern coastal areas during the 1-month field campaign. Three evolution patterns of simulated surface O3 are summarized based on different near-ground flow conditions. More than 75% of days featured interaction between weak synoptic forcing and local sea-land circulations. Integrated process rate (IPR) analysis shows that photochemical production is the dominant contributor to O3 enhancement from 09:00 to 15:00 LST (local standard time) in the atmospheric boundary layer over most areas with elevated O3 occurrence in the mid-afternoon. The simulated ozone production efficiency is 2–8 O3 molecules per NOx molecule oxidized in areas with high O3 chemical production. Precursors of O3 originating from different source regions in the central PRD are mixed during transport to downwind rural areas during nighttime and early morning, where they then contribute to the daytime O3 photochemical production. Such close interactions among precursor emissions, transports, and O3 photochemical production result in the regional O3 pollution over the PRD. Sensitivity studies suggest that O3 formation is volatile organic compound-limited in the central inland PRD, PRE, and surrounding coastal areas with less chemical aging (NOx/NOy〉0.6), but is NOx-limited in the rural southwestern PRD with photochemically aged air (NOx/NOy

Description: Dust storms from major Asian sources are usually carried by northwesterly or westerly winds over northern and southeastern China to the Pacific Ocean. These pathways leave central and southwestern China nearly free of incursions. But a strong dust event on 5–6 May 2005 was captured in a 15-month series of weekly filter samples of PM2.5 at three sites in Chongqing. It illustrated that desert dust can be transported to this region, and sometimes strongly. Annual PM2.5 and dust were similar at the three sites, but higher than in simultaneous samples in Beijing. High correlations of dust concentration were found between the cites during spring, indicating that Asian dust affects a broader swath of China than is often realized. During the event, the concentrations of mineral dust were high at all sites (20–30 μg m−3; 15%–20% of PM2.5 in Chongqing, and 15 μg m−3; 20%–30% of PM2.5 in Beijing), and were part of a broader spring maximum. The proportions of crustal elements and pollution-derived components such as Pb, SO42−, and organic carbon indicated that the sources for this dust differed from Beijing. The dust was considerably enriched in Ca and Mg, characteristic of western deserts, whereas Beijing's dust had the lower Ca and Mg of eastern deserts. This observation agrees with synoptic patterns and back-trajectories. Driven by a cold air outbreak from the northwest, dust from the western Gobi Desert was transported at lower altitudes (

Description: Anthropogenic sources of mercury emissions have a significant impact on global pollution. Therefore, finding uncharacterised sources and assessing the emissions from these sources are important. However, limited data are available worldwide on mercury emissions from crematories. In Japan, 99.9% of dead bodies are cremated, which is the highest percentage in the world, and more than 1600 crematories are in operation. We thus focused on emissions from crematories in Japan. The number of targeted facilities was seven, and we used continuous emission monitoring to measure the mercury concentrations and investigate mercury behaviour. The total mercury concentrations in stack gases were a few μg/Nm3 (normal cubic meters). Considering the time profile of mercury and its species in cremations, the findings confirmed that the mercury in stack gas originated from dental amalgam. The amount of mercury emissions was calculated using the total concentration and gas flow rate. Furthermore, the annual amount of mercury emission from crematories in Japan was estimated by using the total number of corpses. The emission amount was considerably lower than that estimated in the UK. From statistical analyses on population demographics and measurements, future total emissions from crematories were also predicted. As a result, the amount of mercury emitted by crematories will likely increase by 2.6-fold from 2007 to 2037.

Description: Owing to recent industrialization, Central East China has become a significant source of air pollutants. To examine the processes controlling the chemistry and transport of tropospheric ozone, we continuously measured non-methane volatile organic compounds (NMVOCs) as part of an intensive field campaign at Mount Tai, China, in June 2006 (MTX2006), using proton transfer reaction mass spectrometry (PTR-MS). Temporal variations of NMVOCs were recorded in mass-scan mode from m/z 17 to m/z 300 during 12–30 June 2006. More than thirty kinds of NMVOCs were detected up to m/z 160, including alkenes, aromatics, alcohols, aldehydes, and ketones. Oxygenated VOCs were the predominant NMVOCs. During the night of 12 June, we observed an episode of high NMVOCs concentrations attributed to the burning of agricultural biomass. The ΔNMVOCs/ΔCO ratios derived by PTR-MS measurements for this episode are compared to emission ratios from various types of biomass burning as reviewed by Andreae and Merlet (2001) and to ratios recently measured by PTR-MS in tropical forests (Karl et al., 2007) and at urban sites (Warneke et al., 2007).

Description: While Antarctica is often described as a pristine environment, the potential threats from local pollution sources including tourist ships and emissions associated with scientific activities have recently been raised. However, to date there has been no systematic attempt to model the impacts of such pollutants at the continental scale. Indeed, until very recently there was not even a sulphur emission budget available for Antarctica. Here we present the first comprehensive study of atmospheric pollution in Antarctica using a limited area chemistry climate model, and a monthly emissions inventory for sulphur from maintenance of research stations, ground and air traffic, shipping and the active volcano Mt. Erebus. We find that ship emissions, both sulphurous and black carbon, dominate anthropogenic pollution near the ground. These are likely to rise considerably if recent trends in tourism continue.

Description: Two 3-dimensional global atmospheric transport models for persistent organic pollutants (POPs) have been employed to investigate the association between the large-scale atmospheric motions and poleward transports of persistent semi-volatile organic chemicals (SVOCs). We examine the modeled daily air concentration of α- and γ-hexachlorocyclohexane (HCH) over a period from 1997–1999 during which a number of episodic atmospheric transport events were detected in this modeling study. These events provide modeling evidence for improving the interpretation on the cold condensation effect and poleward atmospheric transport of SVOCs at the mid-troposphere. Two episodic transport events of γ-HCH (lindane) to the high Arctic (80–90° N), one from Asian and another from Eurasian sources, are reported in this paper. The both events suggest that the episodic atmospheric transports occurring at the mid-troposphere (e.g. from 3000–5500 m height) are driven by atmospheric horizontal and vertical motions. The association of the transport events with atmospheric circulation is briefly discussed. Strong southerly winds, forced by the evolution of two semi-permanent high pressure systems over mid-high latitudes in the Northern Hemisphere, play an important role in the long-range transport (LRT) of HCHs to the high latitudes from its sources. Being consistent with the cold condensation effect and poleward atmospheric transport in a mean meridional atmospheric circulation simulated by a 2-D atmospheric transport model, as reported by the first part of this study, this modeling study indicates that cold condensation is likely occurring more intensively at the mid-troposphere where rapid declining air temperature results in condensed phase of the chemical over and near its source regions and where stronger winds convey the chemical more rapidly to the polar region during the episodic poleward atmospheric transport events.

Description: Time profiles of molecular iodine emissions from seven species of seaweed have been measured at high time resolution (7.5 s) by direct spectroscopic quantification of the gas phase I2 using broadband cavity enhanced absorption spectroscopy. Substantial differences were found between species, both in the amounts of I2 emitted when the plants were exposed to air and in the shapes of their emission time profiles. Two species of kelp, Laminaria digitata and Laminaria hyperborea, were found to be the most potent emitters, producing an intense burst of I2 when first exposed to air. I2 was also observed from Saccharina latissima and Ascophyllum nodosum but in lower amounts and with broader time profiles. I2 mixing ratios from two Fucus species and Dictyopteris membranacea were at or below the detection limit of the present instrument (25 pptv). A further set of experiments investigated the time dependence of I2 emissions and aerosol particle formation when fragments of L. digitata were exposed to desiccation in air, to ozone and to oligoguluronate stress factors. Particle formation occurred in all L. digitata stress experiments where ozone and light were present, subject to the I2 mixing ratios being above certain threshold amounts. Moreover, the particle number concentrations closely tracked variations in the I2 mixing ratios, confirming the results of previous studies that the condensable particle-forming gases derive from the photochemical oxidation of the plant's I2 emissions. This work also supports the theory that particle nucleation in the coastal atmosphere occurs in "hot-spot" regions of locally elevated concentrations of condensable gases: the greatest atmospheric concentrations of I2 and hence of condensable iodine oxides are likely to be above plants of the most efficiently emitting kelp species and localised in time to shortly after these seaweeds are uncovered by a receding tide.

Description: Gaseous elemental mercury (GEM) measurements at Alert, Canada, from 1995 to 2007 were analyzed for statistical time trends and for correlations with meteorological and climate data. A significant decreasing trend in annual GEM concentration is reported at Alert, with an estimated slope of −0.0086 ng m−3 yr−1 (−0.6% yr−1) over this 13-year period. It is shown that there has been a shift in the month of minimum mean GEM concentration from May to April due to a change in the timing of springtime atmospheric mercury depletion events (AMDEs). These AMDEs are found to decrease with increasing local temperature within each month, both at Alert and at Amderma, Russia. These results agree with the temperature dependence suggested by previous experimental results and theoretical kinetic calculations and highlight the potential for changes in Arctic mercury chemistry with climate. A correlation between total monthly AMDEs at Alert and the Polar/Eurasian Teleconnection Index was observed only in March, perhaps due to higher GEM inputs in early spring in those years with a weak polar vortex. A correlation of AMDEs at Alert with wind direction supports the origin of mercury depletion events over the Arctic Ocean, in agreement with a previous trajectory study of ozone depletion events. Interannual variability in total monthly depletion event frequency at Alert does not appear to correlate significantly with total or first-year northern hemispheric sea ice area or with other major teleconnection patterns. Nor do AMDEs at either Alert or Amderma correlate with local wind speed, as might be expected if depletion events are sustained by stable, low-turbulence atmospheric conditions. The data presented here – both the change in timing of depletion events and their relationship with temperature – can be used as additional constraints to improve the ability of global models to predict the cycling and deposition of mercury in the Arctic.

Description: Polar ecosystems are considered to be the last pristine environments of the Earth relatively uninfluenced by human activities. Antarctica in particular, compared to the Arctic is considered to be even less affected by any kind of anthropogenic influences. Once contaminants reach the polar regions, their lifetime in the troposphere depends on local removal processes. Atmospheric mercury, in particular, has unique characteristics that include long-range transport to polar regions and the transformation to more toxic and water-soluble compounds that may potentially become bioavailable. These chemical-physical properties have given mercury on the priority list of an increasing number of international, European and national conventions and agreements aimed at the protection of the ecosystems including human health (i.e., GEO, UNEP, AMAP, UN-ECE, HELCOM, OSPAR) thus stimulating a significant amount of research including measurements of Hg0 reaction rate constant with atmospheric oxidants, experimental and modelling studies in order to understand the cycling of Hg in polar regions and its impact to these ecosystems. Special attention in terms of contamination of polar regions, is paid to the consequences of the springtime phenomena, referred to as ''atmospheric mercury depletion event'' (AMDE), during which elemental gaseous mercury (GEM or Hg0) through a series of photochemically-initiated reactions involving halogens, may be converted to a reactive form that may accumulate in polar ecosystems. The discovery of the AMDE, first noted in the Arctic, has also been observed at both poles and was initially considered to result in an important net input of atmospheric Hg into the polar surfaces. However, recent studies point out that complex processes take place after deposition that may result in less significant net-inputs from the atmosphere since a fraction, sometimes significant of deposited Hg may be recycled. Therefore, the contribution of this unique reactivity occurring in polar atmospheres to the global budget of atmospheric Hg and the role played by snow and ice surfaces of these regions are important issues. This paper presents a review of atmospheric mercury studies conducted in the Antarctic troposphere, both at coastal locations and on the Antarctic Plateau since 1985. Our current understanding of atmospheric reactivity in this region is also presented.

Description: Recent thinning of glaciers over the Himalayas (sometimes referred to as the third polar region) have raised concern on future water supplies since these glaciers supply water to large river systems that support millions of people inhabiting the surrounding areas. Black carbon (BC) aerosols, released from incomplete combustion, have been increasingly implicated as causing large changes in the hydrology and radiative forcing over Asia and its deposition on snow is thought to increase snow melt. In India BC from biofuel combustion is highly prevalent and compared to other regions, BC aerosol amounts are high. Here, we quantify the impact of BC aerosols on snow cover and precipitation from 1990 to 2010 over the Indian subcontinental region using two different BC emission inventories. New estimates indicate that Indian BC from coal and biofuel are large and transport is expected to expand rapidly in coming years. We show that over the Himalayas, from 1990 to 2000, simulated snow/ice cover decreases by ~0.9% due to aerosols. The contribution of the enhanced Indian BC to this decline is ~30%, similar to that simulated for 2000 to 2010. Spatial patterns of modeled changes in snow cover and precipitation are similar to observations (from 1990 to 2000), and are mainly obtained with the newer BC estimates.

Description: The detailed mechanism of secondary new particle formation in the atmosphere is still under debate. It is proposed that particle formation happens via activation of 1–2 nm atmospheric neutral molecular clusters and/or large molecules. Since traditional instrumentation does not reach these sizes, the hypothesis has not yet been verified. By directly measuring particle size distributions down to mobility diameters of about 1.3 nm with a pulse-height CPC we provide evidence of the nucleation mechanism in coastal environment (Mace Head, Ireland) and in boreal forest (Hyytiälä, Finland). In both places neutral sub-3 nm condensation nuclei (nano-CN) were continuously present, even when no new particle formation was detected. In Mace Head, however, the concentration of the nano-CN was far too low to account for the particle formation rates during particle bursts. Thus the results imply that on coastal sites new particle formation initiates, as proposed earlier, via homogenous nucleation from biogenic iodine vapors. In contrary, activation of pre-existing nano-CN remains a possible explanation in the boreal forest, but the observed concentrations are not the limiting factor for the particle formation events.

Description: This paper presents a summary of the measurements that were made during the heavily-instrumented Reactive Halogens in the Marine Boundary Layer (RHaMBLe) coastal study in Roscoff on the North West coast of France. It was clearly demonstrated that iodine-mediated coastal particle formation occurs, driven by daytime low tide emission of molecular iodine, I2, by macroalgal species fully or partially exposed by the receding waterline. Ultrafine particle concentrations strongly correlate with the rapidly recycled reactive iodine species, IO, produced at high concentrations following photolysis of I2. The heterogeneous macroalgal I2 sources lead to variable relative concentrations of iodine species observed by path-integrated and in situ measurement techniques. Apparent particle emission fluxes were associated with an enhanced apparent depositional flux of ozone, consistent with both a direct O3 deposition to macroalgae and involvement of O3 in iodine photochemistry and subsequent particle formation below the measurement height. The magnitude of the particle formation events was observed to be greatest at the lowest tides with higher concentrations of ultrafine particles growing to much larger sizes, probably by the condensation of anthropogenically-formed condensable material. At such sizes the particles should be able to act as cloud condensation nuclei at reasonable atmospheric supersaturations.

Description: In general, regional and global chemistry transport models apply instantaneous mixing of emissions into the model's finest resolved scale. In case of a concentrated source, this could result in erroneous calculation of the evolution of both primary and secondary chemical species. Several studies discussed this issue in connection with emissions from ships and aircrafts. In this study, we present an approach to deal with the non-linear effects during dispersion of NOx emissions from ships. It represents an adaptation of the original approach developed for aircraft NOx emissions, which uses an exhaust tracer to trace the amount of the emitted species in the plume and applies an effective reaction rate for the ozone production/destruction during the plume's dilution into the background air. In accordance with previous studies examining the impact of international shipping on the composition of the troposphere, we found that the contribution of ship induced surface NOx to the total reaches 90% over remote ocean and makes 10–30% near coastal regions. Due to ship emissions, surface ozone increases by up to 4–6 ppbv making 10% contribution to the surface ozone budget. When applying the ship plume parameterization, we showed that the large scale NOx decreases and the ship NOx contribution is reduced by up to 20–25%. Similar decrease was found in case of O3. The plume parameterization suppressed the ship induced ozone production by 15–30% over large areas of the focused region. To evaluate the presented parameterization, nitrogen oxide measurements over the English Channel were compared with modeled values and it was found that after activating the parameterization the model accuracy increases.

Description: Chlorine dioxide, OClO, column amounts retrieved from measurements of the SCIAMACHY satellite instrument are presented and validated by comparison with simultaneous ground-based DOAS observations. In addition, the measurements are compared to model calculations taking into account the photochemical change along the light path. Although OClO does not participate directly in the destruction of ozone, its accurate measurement as well as modelling is crucial to understand the highly perturbed chlorine chemistry in the polar vortices. SCIAMACHY OClO slant columns retrieved during spring 2005 have been quantitatively validated by comparison with slant columns retrieved from measurements made in Ny-Ålesund (79° N, 12° E), Summit (73° N, 38° W) and Bremen (53° N, 9° E). Overall, good agreement is found. OClO slant column densities modelled with a set of stacked box models and considering the light path through the atmosphere are also included in this comparison. The model predictions differ significantly from the measured quantities. OClO amounts are underestimated for conditions of strong chlorine activation and at large solar zenith angles. Sensitivity studies for several parameters in the stacked box model have been performed and it is inferred that using the chemistry known to date, the observed OClO cannot be adequately reproduced within the range of uncertainties given for the various model parameters.

Description: Recent work has hypothesized that tropical cyclones in the deep Atlantic and eastern Pacific basins develop from the cyclonic Kelvin cat's eye of a tropical easterly wave critical layer located equatorward of the easterly jet axis that typifies the trade wind belt. The cyclonic critical layer is thought to be important to tropical cyclogenesis because its cat's eye provides (i) a region of cyclonic vorticity and weak deformation by the resolved flow, (ii) containment of moisture entrained by the developing flow and/or lofted by deep convection therein, (iii) confinement of mesoscale vortex aggregation, (iv) a predominantly convective type of heating profile, and (v) maintenance or enhancement of the parent wave until the developing proto-vortex becomes a self-sustaining entity and emerges from the wave as a tropical depression. This genesis sequence and the overarching framework for describing how such hybrid wave-vortex structures become tropical depressions/storms is likened to the development of a marsupial infant in its mother's pouch, and for this reason has been dubbed the "marsupial paradigm". Here we conduct the first multi-scale test of the marsupial paradigm in an idealized setting by revisiting the problem of the transformation of an easterly wave-like disturbance into a tropical storm vortex using the WRF model. An analysis of the evolving winds, equivalent potential temperature, and relative vertical vorticity is presented from coarse (28 km) and high resolution (3.1 km) simulations. The results are found to support key elements of the marsupial paradigm by demonstrating the existence of a vorticity dominant region with minimal strain/shear deformation within the critical layer pouch that contains strong cyclonic vorticity and high saturation fraction. This localized region within the pouch serves as the "attractor" for an upscale "bottom up" development process while the wave pouch and proto-vortex move together. Implications of these findings are discussed in relation to an upcoming field experiment for the most active period of the Atlantic hurricane season in 2010 that is to be conducted collaboratively between the National Oceanic and Atmospheric Administration (NOAA), the National Science Foundation (NSF), and the National Aeronautics and Space Adminstration (NASA).

Description: We use observations of total particle number concentration at 36 worldwide sites and a global aerosol model to quantify the primary and secondary sources of particle number. We show that emissions of primary particles can reasonably reproduce the spatial pattern of observed condensation nuclei (CN) (R2=0.51) but fail to explain the observed seasonal cycle at many sites (R2=0.1). The modeled CN concentration in the free troposphere is biased low (normalised mean bias, NMB=−88%) unless a secondary source of particles is included, for example from binary homogeneous nucleation of sulfuric acid and water (NMB=−25%). Simulated CN concentrations in the continental boundary layer (BL) are also biased low (NMB=−74%) unless the number emission of anthropogenic primary particles is increased or an empirical BL particle formation mechanism based on sulfuric acid is used. We find that the seasonal CN cycle observed at continental BL sites is better simulated by including a BL particle formation mechanism (R2=0.3) than by increasing the number emission from primary anthropogenic sources (R2=0.18). Using sensitivity tests we derive optimum rate coefficients for this nucleation mechanism, which agree with values derived from detailed case studies at individual sites.

Description: Most parameterizations for precipitating convection in use today are bulk schemes, in which an ensemble of cumulus elements with different properties is modelled as a single, representative entraining-detraining plume. We review the underpinning mathematical model for such parameterizations, in particular by comparing it with spectral models in which elements are not combined into the representative plume. The chief merit of a bulk model is that the representative plume can be described by an equation set with the same structure as that which describes each element in a spectral model. The equivalence relies on an ansatz for detrained condensate introduced by Yanai et al. (1973) and on a simplified microphysics. There are also conceptual differences in the closure of bulk and spectral parameterizations. In particular, we show that the convective quasi-equilibrium closure of Arakawa and Schubert (1974) for spectral parameterizations cannot be carried over to a bulk parameterization in a straightforward way. Quasi-equilibrium of the cloud work function assumes a timescale separation between a slow forcing process and a rapid convective response. But, for the natural bulk analogue to the cloud-work function (the dilute CAPE), the relevant forcing is characterised by a different timescale, and so its quasi-equilibrium entails a different physical constraint. Closures of bulk parameterization that use the non-entraining parcel value of CAPE do not suffer from this timescale issue. However, the Yanai et al. (1973) ansatz must be invoked as a necessary ingredient of those closures.

Description: We validate OMI ozone profiles between 0.22–215 hPa and stratospheric ozone columns down to 215 hPa (SOC215) against v2.2 MLS data from 2006. The validation demonstrates convincingly that SOC can be derived accurately from OMI data alone, with errors comparable to or smaller than those from current MLS retrievals, and it demonstrates implicitly that tropospheric ozone column can be retrieved accurately from OMI or similar nadir-viewing ultraviolet measurements alone. The global mean biases are within 2.5% above 100 hPa and 5–10% below 100 hPa; the standard deviations (1σ) are 3.5–5% between 1–50 hPa, 6–9% above 1 hPa and 8–15% below; 50 hPa. OMI shows some latitude and solar zenith angle dependent biases, but the mean biases are mostly within 5% and the standard deviations are mostly within 2–5% except for low altitudes and high latitudes. The excellent agreement with MLS data shows that OMI retrievals can be used to augment the validation of MLS and other stratospheric ozone measurements made with even higher vertical resolution than that for OMI. OMI SOC215 shows a small bias of −0.6% with a standard deviation of 2.8%. When compared as a function of latitude and solar zenith angle, the mean biases are within 2% and the standard deviations range from 2.1 to 3.4%. Assuming 2% precision for MLS SOC215, we deduce that the upper limits of random-noise and smoothing errors for OMI SOC215 range from 0.6% in the southern tropics to 2.8% at northern middle latitudes.

Description: Atmospheric mercury concentrations were measured during two oceanographic cruise campaigns covering the Adriatic Sea, the first during the autumn in 2004 and the second in the summer of 2005. The inclement weather during the autumn campaign meant that no clear in-situ production of oxidised gas phase mercury was seen, and that events where high values of HgII(g) and/or Hg associated with particulates (HgP) were observed, were the result of plumes from anthropogenic emission sources. During the summer campaign however, the by now rather familiar diurnal variation of HgII(g) concentration, with maxima around midday, was observed. Again there were events when high HgII(g) and particulates (HgP) concentrations were seen which did not fit with the pattern of daily in-situ HgII(g) production, which were traceable, with the help of back trajectory calculations, to anthropogenic emission sources. All the emission plumes encountered could be traced back to ports, not all of which are associated with major industrial installations. It therefore seems likely in theses cases that the emissions are either due to shipping or to port activities. Box modelling studies of the summer 2005 campaign show that although the in-situ production of HgII(g) occurs in the MBL, the exact chemical mechanism responsible is difficult to determine. However given the high O3 concentrations encountered during this campaign it seems clear that if Hg0 does react with O3, it does not produce gas phase HgII, and the reaction between Hg0 and OH if it occurs, does not contribute appreciably to HgII(g) production.

Description: Windsor (Ontario) – the automotive capital of Canada does not have any significant mercury (Hg) sources. However, Windsor experiences trans-boundary air pollution as it is located immediately downwind of industrialized regions of the United States of America. A study was conducted in 2007 aimed to identify the potential regional sources of total gaseous mercury (TGM) and investigate the effects of regional sources and other factors on seasonal variability of TGM concentrations in Windsor. TGM concentration was measured at the University of Windsor campus using a Tekran® 2537A Hg vapour analyzer. An annual mean of 2.02±1.63 ng/m3 was observed in 2007. The average TGM concentration was high in the summer (2.48 ng/m3) and winter (2.17 ng/m3), compared to spring (1.88 ng/m3) and fall (1.76 ng/m3). Hybrid receptor modeling potential source contribution function (PSCF) was used by incorporating 72-h backward trajectories and measurements of TGM in Windsor. The results of PSCF were analyzed in conjunction with the Hg emissions inventory of North America (by state/province) to identify regions affecting Windsor. In addition to annual modeling, seasonal PSCF modeling was also conducted. The potential source region was identified between 24–61° N and 51–143° W. Annual PSCF modeling identified major sources southwest of Windsor, stretching from Ohio to Texas. The emissions inventory also supported the findings, as Hg emissions were high in those regions. Results of seasonal PSCF modeling were analyzed to find the combined effects of regional sources, meteorological conditions, and surface reemissions, on intra-annual variability of Hg concentrations. It was found that the summer and winter highs of atmospheric Hg can be attributed to areas where large numbers of coal fired power plants are located in the USA. Weak atmospheric dispersion due to low winds and high reemission from surfaces due to higher temperatures contributed to high concentrations in the summer. In the winter, the atmospheric removal of Hg was slow, but strong winds led to more dispersion, resulting in lower concentrations than the summer. Future studies could use smaller grid sizes and refined emission inventories, for more accurate analysis of source-receptor relationship of atmospheric Hg. Abbreviations of states/provinces: Alabama (AL), Arkansas (AR), British Columbia (BC), Georgia (GA), Iowa (IA), Illinois (IL), Indiana (IN), Kentucky (KY), Louisiana (LA), Manitoba (MB), Michigan (MI), Minnesota (MN), Mississippi (MS), Missouri (MO), Newfoundland and Labrador (NL), Ohio (OH), Ontario (ON), Oregon (OR), Pennsylvania (PA), Tennessee (TN), Texas (TX), West Virginia (WV), Wisconsin (WI).