ALBERT

Description: Perturbations in stratospheric aerosol due to explosive volcanic eruptions are a primary contributor to natural climate variability. Observations of stratospheric aerosol are available for the past decades, and information from ice cores has been used to derive estimates of stratospheric sulfur injections and aerosol optical depth over the Holocene (approximately 10 000 BP to present) and into the last glacial period, extending back to 60 000 BP. Tephra records of past volcanism, compared to ice cores, are less complete but extend much further into the past. To support model studies of the potential impacts of explosive volcanism on climate variability across timescales, we present here an ensemble reconstruction of volcanic stratospheric sulfur injection (VSSI) over the last 140 000 years that is based primarily on terrestrial and marine tephra records. VSSI values are computed as a simple function of eruption magnitude based on VSSI estimates from ice cores and satellite observations for identified eruptions. To correct for the incompleteness of the tephra record, we include stochastically generated synthetic eruptions assuming a constant background eruption frequency from the ice core Holocene record. While the reconstruction often differs from ice core estimates for specific eruptions due to uncertainties in the data used and reconstruction method, it shows good agreement with an ice-core-based VSSI reconstruction in terms of millennial-scale cumulative VSSI variations over the Holocene. The PalVol reconstruction provides a new basis to test the contributions of forced vs. unforced natural variability to the spectrum of climate and the mechanisms leading to abrupt transitions in the palaeoclimate record with low- to high-complexity climate models. The PalVol volcanic forcing reconstruction is available at https://doi.org/10.26050/WDCC/PalVolv1 (Toohey and Schindlbeck-Belo, 2023).

Description: TR17-08, a marine sedimentary core (14.6 m), was collected during 2017 from the Edisto Inlet (Ross Sea, Antarctica), a small fjord near Cape Hallett. The core is characterized by expanded laminated sedimentary sequences making it suitable for studying submillennial processes during the Early Holocene. By studying different well-known foraminifera species (Globocassidulina biora, G. subglobosa, Trifarina angulosa, Nonionella iridea, Epistominella exigua, Stainforthia feylingi, Miliammina arenacea, Paratrochammina bartrami and Portatrochammina antarctica), we were able to identify five different foraminiferal assemblages over the last ∼ 2000 years BP. Comparison with diatom assemblages and other geochemical proxies retrieved from nearby sediment cores in the Edisto Inlet (BAY05-20 and HLF17-1) made it possible to distinguish three different phases characterized by different environmental settings: (1) a seasonal phase (from 2012 to 1486 years BP) characterized by the dominance of calcareous species, indicating a seasonal opening of the inlet by more frequent events of melting of the sea-ice cover during the austral summer and, in general, a higher-productivity, more open and energetic environment; (2) a transitional phase (from 1486 to 696 years BP) during which the fjord experienced less extensive sea-ice melting, enhanced oxygen-poor conditions and carbonate dissolution conditions, indicated by the shifts from calcareous-dominated association to agglutinated-dominated association probably due to a freshwater input from the retreat of three local glaciers at the start of this period; and (3) a cooler phase (from 696 years BP to present) during which the sedimentation rate decreased and few to no foraminiferal specimens were present, indicating ephemeral openings or a more prolonged cover of the sea ice during the austral summer, affecting the nutrient supply and the sedimentation regime.

Description: Published

Description: 95–115

Description: OSA2: Evoluzione climatica: effetti e loro mitigazione

Description: JCR Journal

Description: Abstract. Clouds are assumed to play an important role in the Arctic amplification process. This motivated a detailed investigation of cloud processes, including radiative and turbulent fluxes. Data from the aircraft campaign ACLOUD were analyzed with a focus on the mean and turbulent structure of the cloudy boundary layer over the Fram Strait marginal sea ice zone in late spring and early summer 2017. Vertical profiles of turbulence moments are presented from contrasting atmospheric boundary layers (ABLs) from 4 d. They differ by the magnitude of wind speed, boundary-layer height, stability, the strength of the cloud-top radiative cooling and the number of cloud layers. Turbulence statistics up to third-order moments are presented, which were obtained from horizontal-level flights and from slanted profiles. It is shown that both of these flight patterns complement each other and form a data set that resolves the vertical structure of the ABL turbulence well. The comparison of the 4 d shows that especially during weak wind, even in shallow Arctic ABLs with mixing ratios below 3 g kg-1, cloud-top cooling can serve as a main source of turbulent kinetic energy (TKE).Well-mixed ABLs are generated where TKE is increased and vertical velocity variance shows pronounced maxima in the cloud layer. Negative vertical velocity skewness points then to upside-down convection. Turbulent heat fluxes are directed upward in the cloud layer as a result of cold downdrafts. In two cases with single-layer stratocumulus, turbulent transport of heat flux and of temperature variance are both negative in the cloud layer, suggesting an important role of large eddies. In contrast, in a case with weak cloud-top cooling, these quantities are positive in the ABL due to the heating from the surface. Based on observations and results of a mixed-layer model it is shown that the maxima of turbulent fluxes are, however, smaller than the jump of the net terrestrial radiation flux across the upper part of a cloud due to the (i) shallowness of the mixed layer and (ii) the presence of a downward entrainment heat flux. The mixed-layer model also shows that the buoyancy production of TKE is substantially smaller in stratocumulus over the Arctic sea ice compared to subtropics due to a smaller surface moisture flux and smaller decrease in specific humidity (or even humidity inversions) right above the cloud top. In a case of strong wind, wind shear shapes the ABL turbulent structure, especially over rough sea ice, despite the presence of a strong cloud-top cooling. In the presence of mid-level clouds, cloud-top radiative cooling and thus also TKE in the lowermost cloud layer are strongly reduced, and the ABL turbulent structure becomes governed by stability, i.e., by the surface–air temperature difference and wind speed. A comparison of slightly unstable and weakly stable cases shows a strong reduction of TKE due to increased stability even though the absolute value of wind speed was similar. In summary, the presented study documents vertical profiles of the ABL turbulence with a high resolution in a wide range of conditions. It can serve as a basis for turbulence closure evaluation and process studies in Arctic clouds.

Description: Springtime Arctic mixed-phase convection over open water in the Fram Strait as observed during the recent ACLOUD (Arctic CLoud Observations Using airborne measurements during polar Day) field campaign is simulated at turbulence-resolving resolutions. The first objective is to assess the skill of large-eddy simulation (LES) in reproducing the observed mixed-phase convection. The second goal is to then use the model to investigate how aerosol modulates the way in which turbulent mixing and clouds transform the low-level air mass. The focus lies on the low-level thermal structure and lapse rate, the heating efficiency of turbulent entrainment, and the low-level energy budget. A composite case is constructed based on data collected by two research aircraft on 18 June 2017. Simulations are evaluated against independent datasets, showing that the observed thermodynamic, cloudy, and turbulent states are well reproduced. Sensitivity tests on cloud condensation nuclei (CCN) concentration are then performed, covering a broad range between pristine polar and polluted continental values. We find a significant response in the resolved mixed-phase convection, which is in line with previous LES studies. An increased CCN substantially enhances the depth of convection and liquid cloud amount, accompanied by reduced surface precipitation. Initializing with the in situ CCN data yields the best agreement with the cloud and turbulence observations, a result that prioritizes its measurement during field campaigns for supporting high-resolution modeling efforts. A deeper analysis reveals that CCN significantly increases the efficiency of radiatively driven entrainment in warming the boundary layer. The marked strengthening of the thermal inversion plays a key role in this effect. The low-level heat budget shifts from surface driven to radiatively driven. This response is accompanied by a substantial reduction in the surface energy budget, featuring a weakened flow of solar radiation into the ocean. Results are interpreted in the context of air–sea interactions, air mass transformations, and climate feedbacks at high latitudes.

Description: The presence of reactive bromine in polar regions is a widespread phenomenon that plays an important role in the photochemistry of the Arctic and Antarctic lower troposphere, including the destruction of ozone, the disturbance of radical cycles, and the oxidation of gaseous elemental mercury. The chemical mechanisms leading to the heterogeneous release of gaseous bromine compounds from saline surfaces are in principle well understood. There are, however, substantial uncertainties about the contribution of different potential sources to the release of reactive bromine, such as sea ice, brine, aerosols, and the snow surface, as well as about the seasonal and diurnal variation and the vertical distribution of reactive bromine. Here we use continuous long-term measurements of the vertical distribution of bromine monoxide (BrO) and aerosols at the two Antarctic sites Neumayer (NM) and Arrival Heights (AH), covering the periods of 2003–2021 and 2012–2021, respectively, to investigate how chemical and physical parameters affect the abundance of BrO. We find the strongest correlation between BrO and aerosol extinction (R=0.56 for NM and R=0.28 for AH during spring), suggesting that the heterogeneous release of Br2 from saline airborne particles (blowing snow and aerosols) is a dominant source for reactive bromine. Positive correlations between BrO and contact time of air masses, both with sea ice and the Antarctic ice sheet, suggest that reactive bromine is not only emitted by the sea ice surface but by the snowpack on the ice shelf and in the coastal regions of Antarctica. In addition, the open ocean appears to represent a source for reactive bromine during late summer and autumn when the sea ice extent is at its minimum. A source–receptor analysis based on back trajectories and sea ice maps shows that main source regions for BrO at NM are the Weddell Sea and the Filchner–Ronne Ice Shelf, as well as coastal polynyas where sea ice is newly formed. A strong morning peak in BrO frequently occurring during summer and that is particularly strong during autumn suggests a night-time build-up of Br2 by heterogeneous reaction of ozone on the saline snowpack in the vicinity of the measurement sites. We furthermore show that BrO can be sustained for at least 3 d while travelling across the Antarctic continent in the absence of any saline surfaces that could serve as a source for reactive bromine.

Description: Stable water isotopologues of snow, firn and ice cores provide valuable information on past climate variations. Yet single profiles are generally not suitable for robust climate reconstructions. Stratigraphic noise, introduced by the irregular deposition, wind-driven erosion and redistribution of snow, impacts the utility of high-resolution isotope records, especially in low-Accumulation areas. However, it is currently unknown how stratigraphic noise differs across the East Antarctic Plateau and how it is affected by local environmental conditions. Here, we assess the amount and structure of stratigraphic noise at seven sites along a 120 km transect on the plateau of Dronning Maud Land, East Antarctica. Replicated oxygen isotope records of 1 m length were used to estimate signal-To-noise ratios as a measure of stratigraphic noise at sites characterised by different accumulation rates (43-64 mm w.e. a-1), snow surface roughnesses and slope inclinations. While we found a high level of stratigraphic noise at all sites, there was also considerable variation between sites. At sastrugi-dominated sites, greater stratigraphic noise coincided with stronger surface roughnesses, steeper slopes and lower accumulation rates, probably related to increased wind speeds. These results provide a first step to modelling stratigraphic noise and might guide site selection and sampling strategies for future expeditions to improve high-resolution climate reconstructions from low-Accumulation regions.

Description: The Global Ocean Data Analysis Project (GLODAP) is a synthesis effort providing regular compilations of surface-to-bottom ocean biogeochemical bottle data, with an emphasis on seawater inorganic carbon chemistry and related variables determined through chemical analysis of seawater samples. GLODAPv2.2022 is an update of the previous version, GLODAPv2.2021 (Lauvset et al., 2021). The major changes are as follows: data from 96 new cruises were added, data coverage was extended until 2021, and for the first time we performed secondary quality control on all sulphur hexafluoride (SF6) data. In addition, a number of changes were made to data included in GLODAPv2.2021. These changes affect specifically the SF6 data, which are now subjected to secondary quality control, and carbon data measured onboard the RV Knorr in the Indian Ocean in 1994–1995 which are now adjusted using CRM measurements made at the time. GLODAPv2.2022 includes measurements from almost 1.4 million water samples from the global oceans collected on 1085 cruises. The data for the now 13 GLODAP core variables (salinity, oxygen, nitrate, silicate, phosphate, dissolved inorganic carbon, total alkalinity, pH, CFC-11, CFC-12, CFC-113, CCl4, and SF6) have undergone extensive quality control with a focus on systematic evaluation of bias. The data are available in two formats: (i) as submitted by the data originator but converted to World Ocean Circulation Experiment (WOCE) exchange format and (ii) as a merged data product with adjustments applied to minimize bias. For the present annual update, adjustments for the 96 new cruises were derived by comparing those data with the data from the 989 quality controlled cruises in the GLODAPv2.2021 data product using crossover analysis. SF6 data from all cruises were evaluated by comparison with CFC-12 data measured on the same cruises. For nutrients and ocean carbon dioxide (CO2) chemistry comparisons to estimates based on empirical algorithms provided additional context for adjustment decisions. The adjustments that we applied are intended to remove potential biases from errors related to measurement, calibration, and data handling practices without removing known or likely time trends or variations in the variables evaluated. The compiled and adjusted data product is believed to be consistent to better than 0.005 in salinity, 1 % in oxygen, 2 % in nitrate, 2 % in silicate, 2 % in phosphate, 4 μmol kg-1 in dissolved inorganic carbon, 4 μmol kg-1 in total alkalinity, 0.01–0.02 in pH (depending on region), and 5 % in the halogenated transient tracers. The other variables included in the compilation, such as isotopic tracers and discrete CO2 fugacity (fCO2), were not subjected to bias comparison or adjustments.

Description: Arctic river deltas and deltaic near-shore zones represent important land–ocean transition zones influencing sediment dynamics and nutrient fluxes from permafrost-affected terrestrial ecosystems into the coastal Arctic Ocean. To accurately model fluvial carbon and freshwater export from rapidly changing river catchments as well as assess impacts of future change on the Arctic shelf and coastal ecosystems, we need to understand the sea floor characteristics and topographic variety of the coastal zones. To date, digital bathymetrical data from the poorly accessible, shallow, and large areas of the eastern Siberian Arctic shelves are sparse. We have digitized bathymetrical information for nearly 75 000 locations from large-scale (1:25 000–1:500 000) current and historical nautical maps of the Lena Delta and the Kolyma Gulf region in northeastern Siberia. We present the first detailed and seamless digital models of coastal zone bathymetry for both delta and gulf regions in 50 and 200 m spatial resolution. We validated the resulting bathymetry layers using a combination of our own water depth measurements and a collection of available depth measurements, which showed a strong correlation (r〉0.9). Our bathymetrical models will serve as an input for a high-resolution coupled hydrodynamic–ecosystem model to better quantify fluvial and coastal carbon fluxes to the Arctic Ocean, but they may be useful for a range of other studies related to Arctic delta and near-shore dynamics such as modeling of submarine permafrost, near-shore sea ice, or shelf sediment transport. The new digital high-resolution bathymetry products are available on the PANGAEA data set repository for the Lena Delta (https://doi.org/10.1594/PANGAEA.934045; Fuchs et al., 2021a) and Kolyma Gulf region (https://doi.org/10.1594/PANGAEA.934049; Fuchs et al., 2021b), respectively. Likewise, the depth validation data are available on PANGAEA as well (https://doi.org/10.1594/PANGAEA.933187; Fuchs et al., 2021c).

Description: Data from the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition allowed us to investigate the temporal dynamics of snowfall, snow accumulation and erosion in great detail for almost the whole accumulation season (November 2019 to May 2020). We computed cumulative snow water equivalent (SWE) over the sea ice based on snow depth and density retrievals from a SnowMicroPen and approximately weekly measured snow depths along fixed transect paths. We used the derived SWE from the snow cover to compare with precipitation sensors installed during MOSAiC. The data were also compared with ERA5 reanalysis snowfall rates for the drift track. We found an accumulated snow mass of 38 m SWE between the end of October 2019 and end of April 2020. The initial SWE over first-year ice relative to second-year ice increased from 50 % to 90 % by end of the investigation period. Further, we found that the Vaisala Present Weather Detector 22, an optical precipitation sensor, and installed on a railing on the top deck of research vessel Polarstern, was least affected by blowing snow and showed good agreements with SWE retrievals along the transect. On the contrary, the OTT Pluvio2 pluviometer and the OTT Parsivel2 laser disdrometer were largely affected by wind and blowing snow, leading to too high measured precipitation rates. These are largely reduced when eliminating drifting snow periods in the comparison. ERA5 reveals good timing of the snowfall events and good agreement with ground measurements with an overestimation tendency. Retrieved snowfall from the ship-based Ka-band ARM zenith radar shows good agreements with SWE of the snow cover and differences comparable to those of ERA5. Based on the results, we suggest the Ka-band radar-derived snowfall as an upper limit and the present weather detector on RV Polarstern as a lower limit of a cumulative snowfall range. Based on these findings, we suggest a cumulative snowfall of 72 to 107 m and a precipitation mass loss of the snow cover due to erosion and sublimation as between 47 % and 68 %, for the time period between 31 October 2019 and 26 April 2020. Extending this period beyond available snow cover measurements, we suggest a cumulative snowfall of 98-114 m.

Description: Authorship conflicts are a common occurrence in academic publishing, and they can have serious implications for the careers and well-being of the involved researchers as well as the collective success of research organizations. In addition to not inviting relevant contributors to co-author a paper, the order of authors as well as honorary, gift, and ghost authors are all widely recognized problems related to authorship. Unfair authorship practices disproportionately affect those lower in the power hierarchies – early career researchers, women, researchers from the Global South, and other minoritized groups. Here we propose an approach to preparing author lists based on clear, transparent, and timely communication. This approach aims to minimize the potential for late-stage authorship conflicts during manuscript preparation by facilitating timely and transparent decisions on potential co-authors and their responsibilities. Furthermore, our approach can help avoid imbalances between contributions and credits in published papers by recording planned and executed responsibilities. We present authorship guidelines which also include a novel authorship form along with the documentation of the formulation process for a multidisciplinary and interdisciplinary center with more than 250 researchers. Other research groups, departments, and centers can use or build on this template to design their own authorship guidelines as a practical way to promote fair authorship practices.

Description: The Arctic is exposed to even faster temperature changes than most other areas on Earth. Constantly increasing temperature will lead to thawing permafrost and changes in the methane (CH4) emissions from wetlands. One of the places exposed to those changes is the Abisko–Stordalen Mire in northern Sweden, where climate and vegetation studies have been conducted since the 1970s. In our study, we analyzed field-scale methane emissions measured by the eddy covariance method at Abisko–Stordalen Mire for 3 years (2014–2016). The site is a subarctic mire mosaic of palsas, thawing palsas, fully thawed fens, and open water bodies. A bimodal wind pattern prevalent at the site provides an ideal opportunity to measure mire patches with different permafrost status with one flux measurement system. The flux footprint for westerly winds was dominated by elevated palsa plateaus, while the footprint was almost equally distributed between palsas and thawing bog-like areas for easterly winds. As these patches are exposed to the same climatic and weather conditions, we analyzed the differences in the responses of their methane emission for environmental parameters. The methane fluxes followed a similar annual cycle over the 3 study years, with a gentle rise during spring and a decrease during autumn, without emission bursts at either end of the ice-free season. The peak emission during the ice-free season differed significantly for the two mire areas with different permafrost status: the palsa mire emitted 19 mg-C m−2 d−1 and the thawing wet sector 40 mg-C m−2 d−1. Factors controlling the methane emission were analyzed using generalized linear models. The main driver for methane fluxes was peat temperature for both wind sectors. Soil water content above the water table emerged as an explanatory variable for the 3 years for western sectors and the year 2016 in the eastern sector. The water table level showed a significant correlation with methane emission for the year 2016 as well. Gross primary production, however, did not show a significant correlation with methane emissions. Annual methane emissions were estimated based on four different gap-filing methods. The different methods generally resulted in very similar annual emissions. The mean annual emission based on all models was 3.1 ± 0.3 g-C m−2 a−1 for the western sector and 5.5 ± 0.5 g-C m−2 a−1 for the eastern sector. The average annual emissions, derived from these data and a footprint climatology, were 2.7 ± 0.5 and 8.2 ± 1.5 g-C m−2 a−1 for the palsa and thawing surfaces, respectively. Winter fluxes were relatively high, contributing 27 %–45 % to the annual emissions.

Description: This work presents the integration of a gas-phase and particulate atmospheric emission inventory (AEI) for Argentina in high spatial resolution (0.025∘×0.025∘; approx. 2.5 km×2.5 km) considering monthly variability from 1995 to 2020. The new inventory, called GEAA-AEIv3.0M, includes the following activities: energy production, fugitive emissions from oil and gas production, industrial fuel consumption and production, transport (road, maritime, and air), agriculture, livestock production, manufacturing, residential, commercial, and biomass and agricultural waste burning. The following species, grouped by atmospheric reactivity, are considered: (i) greenhouse gases (GHGs) – CO2, CH4, and N2O; (ii) ozone precursors – CO, NOx (NO+NO2), and non-methane volatile organic compounds (NMVOCs); (iii) acidifying gases – NH3 and SO2; and (iv) particulate matter (PM) – PM10, PM2.5, total suspended particles (TSPs), and black carbon (BC). The main objective of the GEAA-AEIv3.0M high-resolution emission inventory is to provide temporally resolved emission maps to support air quality and climate modeling oriented to evaluate pollutant mitigation strategies by local governments. This is of major concern, especially in countries where air quality monitoring networks are scarce, and the development of regional and seasonal emissions inventories would result in remarkable improvements in the time and space chemical prediction achieved by air quality models. Despite distinguishing among different sectoral and activity databases as well as introducing a novel spatial distribution approach based on census radii, our high-resolution GEAA-AEIv3.0M shows equivalent national-wide total emissions compared to the Third National Communication of Argentina (TNCA), which compiles annual GHG emissions from 1990 through 2014 (agreement within ±7.5 %). However, the GEAA-AEIv3.0M includes acidifying gases and PM species not considered in TNCA. Temporal comparisons were also performed against two international databases: Community Emissions Data System (CEDS) and EDGAR HTAPv5.0 for several pollutants; for EDGAR it also includes a spatial comparison. The agreement was acceptable within less than 30 % for most of the pollutants and activities, although a 〉90 % discrepancy was obtained for methane from fuel production and fugitive emissions and 〉120 % for biomass burning. Finally, the updated seasonal series clearly showed the pollution reduction due to the COVID-19 lockdown during the first quarter of year 2020 with respect to same months in previous years. Through an open-access data repository, we present the GEAA-AEIv3.0M inventory as the largest and more detailed spatial resolution dataset for the Argentine Republic, which includes monthly gridded emissions for 12 species and 15 stors between 1995 and 2020. The datasets are available at https://doi.org/10.17632/d6xrhpmzdp.2 (Puliafito et al., 2021), under a CC-BY 4 license.

Description: Icequakes are the result of processes occurring within the ice mass or between the ice and its environment. Studying icequakes provides a unique view on ice dynamics, specifically on the basal conditions. Changes in conditions due to environmental or climate changes are reflected in icequakes. Counting and characterizing icequakes is thus essential to monitor them. Most of the icequakes recorded by the seismic station at the Belgian Princess Elisabeth Antarctica Station (PE) have small amplitudes corresponding to maximal displacements of a few nanometres. Their detection threshold is highly variable because of the rapid and strong changes in the local seismic noise level. Therefore, we evaluated the influence of katabatic winds on the noise measured by the well-protected PE surface seismometer. Our purpose is to identify whether the lack of icequake detection during some periods could be associated with variations in the processes generating them or simply with a stronger seismic noise linked to stronger wind conditions. We observed that the wind mainly influences seismic noise at frequencies greater than 1 Hz. The seismic noise power exhibits a bilinear correlation with the wind velocity, with two different slopes at a wind velocity lower and greater than 6 m s−1 and with, for example at a period of 0.26 s, a respective variation of 0.4 dB (m −1 s) and 1.4 dB (m −1 s). These results allowed a synthetic frequency and wind-speed-dependent noise model to be presented that explains the behaviour of the wind-induced seismic noise at PE, which shows that seismic noise amplitude increases exponentially with increasing wind speed. This model enables us to study the influence of the wind on the original seismic dataset, which improves the observation of cryoseismic activity near the PE station.

Description: Improvements in our capability to reconstruct ancient surface-ocean conditions based on organic-walled dinoflagellate cyst (dinocyst) assemblages from the Southern Ocean provide an opportunity to better establish past position, strength and oceanography of the subtropical front (STF). Here, we aim to reconstruct the late Eocene to early Miocene (37–20 Ma) depositional and palaeoceanographic history of the STF in the context of the evolving Tasmanian Gateway as well as the potential influence of Antarctic circumpolar flow and intense waxing and waning of ice. We approach this by combining information from seismic lines (revisiting existing data and generating new marine palynological data from Ocean Drilling Program (ODP) Hole 1168A) in the western Tasmanian continental slope. We apply improved taxonomic insights and palaeoecological models to reconstruct the sea surface palaeoenvironmental evolution. Late Eocene–early Oligocene (37–30.5 Ma) assemblages show a progressive transition from dominant terrestrial palynomorphs and inner-neritic dinocyst taxa as well as cysts produced by heterotrophic dinoflagellates to predominantly outer-neritic/oceanic autotrophic taxa. This transition reflects the progressive deepening of the western Tasmanian continental margin, an interpretation supported by our new seismic investigations. The dominance of autotrophic species like Spiniferites spp. and Operculodinium spp. reflects relatively oligotrophic conditions, like those of regions north of the modern-day STF. The increased abundance in the earliest Miocene of Nematosphaeropsis labyrinthus, typical for modern subantarctic zone (frontal) conditions, indicates a cooling and/or closer proximity of the STF to the site . The absence of major shifts in dinocyst assemblages contrasts with other records in the region and suggests that small changes in surface oceanographic conditions occurred during the Oligocene. Despite the relatively southerly (63–55∘ S) location of Site 1168, the rather stable oceanographic conditions reflect the continued influence of the proto-Leeuwin Current along the southern Australian coast as Australia continued to drift northward. The relatively “warm” dinocyst assemblages at ODP Site 1168, compared with the cold assemblages at Antarctic Integrated Ocean Drilling Program (IODP) Site U1356, testify to the establishment of a pronounced latitudinal temperature gradient in the Oligocene Southern Ocean.

Description: The importance of Antarctic sea ice and Southern Ocean warming has come into the focus of polar research during the last couple of decades. Especially around West Antarctica, where warm water masses approach the continent and where sea ice has declined, the distribution and evolution of sea ice play a critical role in the stability of nearby ice shelves. Organic geochemical analyses of marine seafloor surface sediments from the Antarctic continental margin allow an evaluation of the applicability of biomarker-based sea-ice and ocean temperature reconstructions in these climate-sensitive areas. We analysed highly branched isoprenoids (HBIs), such as the sea-ice proxy IPSO25 and phytoplankton-derived HBI-trienes, as well as phytosterols and isoprenoidal glycerol dialkyl glycerol tetraethers (GDGTs), which are established tools for the assessment of primary productivity and ocean temperatures respectively. The combination of IPSO25 with a phytoplankton marker (i.e. the PIPSO25 index) permits semi-quantitative sea-ice reconstructions and avoids misleading over- or underestimations of sea-ice cover. Comparisons of the PIPSO25-based sea-ice distribution patterns and TEX86L- and RI-OH′-derived ocean temperatures with (1) sea-ice concentrations obtained from satellite observations and (2) instrument measurements of sea surface and subsurface temperatures corroborate the general capability of these proxies to determine oceanic key variables properly. This is further supported by model data. We also highlight specific aspects and limitations that need to be taken into account for the interpretation of such biomarker data and discuss the potential of IPSO25 as an indicator for the former occurrence of platelet ice and/or the export of ice-shelf water.

Description: Aeromagnetic exploration is an important method of geophysical exploration. We study the compensation method of the towed bird system and establish the towed bird interference model. Due to the geomagnetic gradient changing greatly, the geomagnetic gradient is considered in the towed bird interference model. In this paper, we model the geomagnetic field gradient and analyze the influence of the towed bird system on the aeromagnetic compensation results. Finally, we apply the ridge regression method to solve the problem. We verify the feasibility of this compensation method through actual flight tests and further improve the data quality of the towed bird interference.

Description: Two feature-based verification methods, thus far only used for the diagnostic evaluation of atmospheric models, have been applied to compare ∼7 km resolution pre-operational analyses of chlorophyll-a (Chl-a) concentrations to a 1 km gridded satellite-derived Chl-a concentration product. The aim of this study was to assess the value of applying such methods to ocean models. Chl-a bloom objects were identified in both data sets for the 2019 bloom season (1 March to 31 July). These bloom objects were analysed as discrete (2-D) spatial features, but also as space–time (3-D) features, providing the means of defining the onset, duration and demise of distinct bloom episodes and the season as a whole. The new feature-based verification methods help reveal that the model analyses are not able to represent small coastal bloom objects, given the coarser definition of the coastline, also wrongly producing more bloom objects in deeper Atlantic waters. Model analyses' concentrations are somewhat higher overall. The bias manifests itself in the size of the model analysis bloom objects, which tend to be larger than the satellite-derived bloom objects. The onset of the bloom season is delayed by 26 d in the model analyses, but the season also persists for another month beyond the diagnosed end. The season was diagnosed to be 119 d long in the model analyses, compared to 117 d from the satellite product. Geographically, the model analyses and satellite-derived bloom objects do not necessarily exist in a specific location at the same time and only overlap occasionally.

Description: Non-Gaussian forecast error is a challenge for ensemble-based data assimilation (DA), particularly for more nonlinear convective dynamics. In this study, we investigate the degree of the non-Gaussianity of forecast error distributions at 1 km resolution using a 1000-member ensemble Kalman filter, and how it is affected by the DA update frequency and observation number. Regional numerical weather prediction experiments are performed with the SCALE (Scalable Computing for Advanced Library and Environment) model and the LETKF (local ensemble transform Kalman filter) assimilating phased array radar observations every 30 s. The results show that non-Gaussianity develops rapidly within convective clouds and is sensitive to the DA frequency and the number of assimilated observations. The non-Gaussianity is reduced by up to 40 % when the assimilation window is shortened from 5 min to 30 s, particularly for vertical velocity and radar reflectivity.

Description: The mass of the Greenland ice sheet is declining as mass gain from snow accumulation is exceeded by mass loss from surface meltwater runoff, marine-terminating glacier calving and submarine melting, and basal melting. Here we use the input–output (IO) method to estimate mass change from 1840 through next week. Surface mass balance (SMB) gains and losses come from a semi-empirical SMB model from 1840 through 1985 and three regional climate models (RCMs; HIRHAM/HARMONIE, Modèle Atmosphérique Régional – MAR, and RACMO – Regional Atmospheric Climate MOdel) from 1986 through next week. Additional non-SMB losses come from a marine-terminating glacier ice discharge product and a basal mass balance model. From these products we provide an annual estimate of Greenland ice sheet mass balance from 1840 through 1985 and a daily estimate at sector and region scale from 1986 through next week. This product updates daily and is the first IO product to include the basal mass balance which is a source of an additional ∼24 Gt yr−1 of mass loss. Our results demonstrate an accelerating ice-sheet-scale mass loss and general agreement (coefficient of determination, r2, ranges from 0.62 to 0.94) among six other products, including gravitational, volume, and other IO mass balance estimates. Results from this study are available at https://doi.org/10.22008/FK2/OHI23Z (Mankoff et al., 2021).

Description: China has experienced dramatic changes in emissions since 2010, which accelerated following the implementation of the Clean Air Action program in 2013. These changes have resulted in significant air quality improvements that are reflected in observations from both surface networks and satellite observations. Air pollutants, such as PM2.5, surface ozone, and their precursors, have long enough lifetimes in the troposphere to be easily transported downwind. Emission changes in China will thus not only change the domestic air quality but will also affect the air quality in other regions. In this study, we use a global chemistry transport model (CAM-chem) to simulate the influence of Chinese emission changes from 2010 to 2017 on both domestic and foreign air quality. We then quantify the changes in air-pollution-associated (including both PM2.5 and O3) premature mortality burdens at regional and global scales. Within our simulation period, the population-weighted annual PM2.5 concentration in China peaks in 2011 (94.1 µg m−3) and decreases to 69.8 µg m−3 by 2017. These estimated national PM2.5 concentration changes in China are comparable with previous studies using fine-resolution regional models, though our model tends to overestimate PM2.5 from 2013 to 2017 when evaluated with surface observations. Relative to 2010, emission changes in China increased the global PM2.5-associated premature mortality burdens through 2013, among which a majority of the changes (∼ 93 %) occurred in China. The sharp emission decreases after 2013 generated significant benefits for human health. By 2017, emission changes in China reduced premature deaths associated with PM2.5 by 108 800 (92 800–124 800) deaths per year globally, relative to 2010, among which 92 % were realized in China. In contrast, the population-weighted, annually averaged maximum daily 8 h ozone concentration peaked in 2014 and did not reach 2010 levels by 2017. As such, O3 generated nearly 8500 (6500–9900) more premature deaths per year in 2017 compared to 2010. Downwind regions, such as South Korea, Japan, and the United States, generally experienced O3 improvements following 2013 due to the decreased export of ozone and its precursors. Overall, we conclude that the sharp emission reductions in China over the past decade have generated substantial benefits for air quality that have reduced premature deaths associated with air pollution at a global scale.

Description: In this paper, we present a new version of the chemistry–climate model SOCOL-AERv2 supplemented by an iodine chemistry module. We perform three 20-year ensemble experiments to assess the validity of the modeled iodine and to quantify the effects of iodine on ozone. The iodine distributions obtained with SOCOL-AERv2-I agree well with AMAX-DOAS observations and with CAM-chem model simulations. For the present-day atmosphere, the model suggests that the iodine-induced chemistry leads to a 3 %–4 % reduction in the ozone column, which is greatest at high latitudes. The model indicates the strongest influence of iodine in the lower stratosphere with 30 ppbv less ozone at low latitudes and up to 100 ppbv less at high latitudes. In the troposphere, the account of the iodine chemistry reduces the tropospheric ozone concentration by 5 %–10 % depending on geographical location. In the lower troposphere, 75 % of the modeled ozone reduction originates from inorganic sources of iodine, 25 % from organic sources of iodine. At 50 hPa, the results show that the impacts of iodine from both sources are comparable. Finally, we determine the sensitivity of ozone to iodine by applying a 2-fold increase in iodine emissions, as it might be representative for iodine by the end of this century. This reduces the ozone column globally by an additional 1.5 %–2.5 %. Our results demonstrate the sensitivity of atmospheric ozone to iodine chemistry for present and future conditions, but uncertainties remain high due to the paucity of observational data of iodine species.

Description: Frequently occurring mega-droughts under current global climate change have attracted broad social attention. A paleoclimatic perspective is needed to increase our understanding of the causes and effects of droughts. South-western (SW) China has been threatened by severe seasonal droughts. Our current knowledge of millennial-scale dry and wet phases in this region is primarily based on the variability of the Indian summer monsoon. However, water availability over land does not always follow patterns of monsoonal precipitation but also depends on water loss from evaporation and transpiration. Here, we reconstructed precipitation intensity, lake hydrological balance and the soil water stress index (SWSI) for the last 27 000 years. Grain size, geochemical and pollen records from Yilong Lake reveal the long-term relationships and inconsistencies of dry–wet patterns in meteorological, hydrological and soil systems in the central Yunnan region, SW China. Our results show that the long-term trends among precipitation, hydrological balance and soil moisture varied through time. The hydrological balance and soil moisture were primarily controlled by temperature-induced evaporation change during periods of low precipitation such as the Last Glacial Maximum and Younger Dryas. During periods of high precipitation (the early to late Holocene), intensified evaporation from the lake surface offset the effects of increased precipitation on the hydrological balance. However, abundant rainfall and the dense vegetation canopy circumvented a soil moisture deficit that might have resulted from rising temperature. In conclusion, the hydrological balance in the central Yunnan region was more sensitive to temperature change while soil moisture could be further regulated by vegetation changes over millennial timescales. Therefore, under future climate warming, the surface water shortage in the central Yunnan region may become even more serious. Our study suggests that reforestation efforts may provide some relief to soil moisture deficits in this region.

Description: This paper presents a new technique to derive thermospheric temperature from space-based disk observations of far ultraviolet airglow. The technique, guided by findings from principal component analysis of synthetic daytime Lyman–Birge–Hopfield (LBH) disk emissions, uses a ratio of the emissions in two spectral channels that together span the LBH (2,0) band to determine the change in band shape with respect to a change in the rotational temperature of N2. The two-channel-ratio approach limits representativeness and measurement error by only requiring measurement of the relative magnitudes between two spectral channels and not radiometrically calibrated intensities, simplifying the forward model from a full radiative transfer model to a vibrational–rotational band model. It is shown that the derived temperature should be interpreted as a column-integrated property as opposed to a temperature at a specified altitude without utilization of a priori information of the thermospheric temperature profile. The two-channel-ratio approach is demonstrated using NASA GOLD Level 1C disk emission data for the period of 2–8 November 2018 during which a moderate geomagnetic storm has occurred. Due to the lack of independent thermospheric temperature observations, the efficacy of the approach is validated through comparisons of the column-integrated temperature derived from GOLD Level 1C data with the GOLD Level 2 temperature product as well as temperatures from first principle and empirical models. The storm-time thermospheric response manifested in the column-integrated temperature is also shown to corroborate well with hemispherically integrated Joule heating rates, ESA SWARM mass density at 460 km, and GOLD Level 2 column O/N2 ratio.

Description: The observing system design of multidisciplinary field measurements involves a variety of considerations on logistics, safety, and science objectives. Typically, this is done based on investigator intuition and designs of prior field measurements. However, there is potential for considerable increases in efficiency, safety, and scientific success by integrating numerical simulations in the design process. Here, we present a novel numerical simulation–environmental response function (NS–ERF) approach to observing system simulation experiments that aids surface–atmosphere synthesis at the interface of mesoscale and microscale meteorology. In a case study we demonstrate application of the NS–ERF approach to optimize the Chequamegon Heterogeneous Ecosystem Energy-balance Study Enabled by a High-density Extensive Array of Detectors 2019 (CHEESEHEAD19). During CHEESEHEAD19 pre-field simulation experiments, we considered the placement of 20 eddy covariance flux towers, operations for 72 h of low-altitude flux aircraft measurements, and integration of various remote sensing data products. A 2 h high-resolution large eddy simulation created a cloud-free virtual atmosphere for surface and meteorological conditions characteristic of the field campaign domain and period. To explore two specific design hypotheses we super-sampled this virtual atmosphere as observed by 13 different yet simultaneous observing system designs consisting of virtual ground, airborne, and satellite observations. We then analyzed these virtual observations through ERFs to yield an optimal aircraft flight strategy for augmenting a stratified random flux tower network in combination with satellite retrievals. We demonstrate how the novel NS–ERF approach doubled CHEESEHEAD19's potential to explore energy balance closure and spatial patterning science objectives while substantially simplifying logistics. Owing to its modular extensibility, NS–ERF lends itself to optimizing observing system designs also for natural climate solutions, emission inventory validation, urban air quality, industry leak detection, and multi-species applications, among other use cases.

Description: The source apportionment of aerosol iron (Fe), including natural and combustion Fe, is an important issue because aerosol Fe can enhance oceanic primary production in the surface ocean. Based on our previous finding that combustion Fe emitted by evaporation processes has Fe isotope ratios (δ56Fe) that are approximately 4 ‰ lower than those of natural Fe, this study aimed to distinguish aerosol Fe sources over the northwestern Pacific using two size-fractionated marine aerosols. The δ56Fe values of fine and coarse particles from the eastern or northern Pacific were found to be similar to each other, ranging from 0.0 ‰ to 0.4 ‰. Most of them were close to the crustal average, suggesting the dominance of natural Fe. On the other hand, particles from the direction of East Asia demonstrated lower δ56Fe values in fine particles (−0.5 ‰ to −2.2 ‰) than in coarse particles (on average −0.02 ± 0.12 ‰). The correlations between the δ56Fe values and the enrichment factors of lead and vanadium suggested that the low δ56Fe values obtained were due to the presence of combustion Fe. The δ56Fe values of the soluble component of fine particles in this region were lower than the total, indicating the preferential dissolution of combustion Fe. In addition, we found a negative correlation between the δ56Fe value and the fractional Fe solubility in air masses from the direction of East Asia. These results suggest that the presence of combustion Fe is an important factor in controlling the fractional Fe solubility in air masses from the direction of East Asia, whereas other factors are more important in the other areas. By assuming typical δ56Fe values for combustion and natural Fe, the contribution of combustion Fe to the total (acid-digested) Fe in aerosols was estimated to reach up to 50 % of fine and 21 % of bulk (coarse + fine) particles in air masses from the direction of East Asia, whereas its contribution was small in the other areas. The contribution of combustion Fe to the soluble Fe component estimated for one sample was approximately twice as large as the total, indicating the importance of combustion Fe as a soluble Fe source despite lower emissions than the natural. These isotope-based estimates were compared with those estimated using an atmospheric chemical transport model (IMPACT), in which the fractions of combustion Fe in fine particles, especially in air masses from the direction of East Asia, were consistent with each other. In contrast, the model estimated a relatively large contribution from combustion Fe in coarse particles, probably because of the different characteristics of combustion Fe that are included in the model calculation and the isotope-based estimation. This highlights the importance of observational data on δ56Fe for size-fractionated aerosols to scale the combustion Fe emission by the model. The average deposition fluxes of soluble Fe to the surface ocean were 1.4 and 2.9 nmol m−2 d−1 from combustion and natural aerosols, respectively, in air masses from the direction of East Asia, which suggests that combustion Fe could be an important Fe source to the surface seawater among other Fe sources. Distinguishing Fe sources using the δ56Fe values of marine aerosols and seawater is anticipated to lead to a more quantitative understanding of the Fe cycle in the atmosphere and surface ocean.

Description: This article summarises the results of an analysis of solar radio bursts (SRBs) detected by the Compound Astronomical Low-cost Low-frequency Instrument for Spectroscopy and Transportable Observatory (CALLISTO) spectrometer hosted by the University of Rwanda. The data analysed were detected during the first year (2014–2015) of the instrument operation. Using quick plots provided by the e-CALLISTO website, a total of 201 intense and well-separated solar radio bursts detected by the CALLISTO station located in Rwanda, are found consisting of 4 type II, 175 type III and 22 type IV radio bursts. It is found that all analysed type II and ∼ 37 % of type III bursts are associated with impulsive solar flares, while the minority (∼ 13 %) of type IV radio bursts are associated with solar flares. Furthermore, all type II radio bursts are associated with coronal mass ejections (CMEs), ∼ 44 % of type III bursts are associated with CMEs, and the majority (∼ 82 %) of type IV bursts were accompanied by CMEs. With aid of the atmospheric imaging assembly (AIA) images on board the Solar Dynamics Observatory (SDO), the location of open magnetic field lines of non-flare-associated type III radio bursts are shown. The same images are used to show the magnetic loops in the solar corona for type IV radio bursts observed in the absence of solar flares and/or CMEs. Findings from this study indicate that analysis of SRBs that are observed from the ground can provide a significant contribution to the early diagnosis of solar transients phenomena, such as solar flares and CMEs, which are major drivers of potential space weather hazards.

Description: Recent observational and modeling studies show that variations of stratospheric ozone and the resulting interaction between ozone and the stratospheric circulation play an important role in surface weather and climate. However, in many cases, computationally expensive coupled chemistry models have been used to study these effects. Here, we demonstrate how a much simpler idealized general circulation model (GCM) can be used for studying the impact of interactive stratospheric ozone on the circulation. The model, named Simplified Chemistry-Dynamical Model (SCDM V1.0), is constructed from a preexisting idealized GCM, into which a simplified linear ozone scheme and a parameterization for the shortwave radiative effects of ozone are implemented. The distribution and variability of stratospheric ozone simulated by the new model are in good agreement with the MERRA2 reanalysis, even for extreme circulation events such as Arctic stratospheric sudden warmings. The model thus represents a promising new tool for the study of ozone–circulation interaction in the stratosphere and its associated effects on tropospheric weather and climate.

Description: Forest fires modify soil organic carbon and suppress soil respiration for many decades after the initial disturbance. The associated changes in soil autotrophic and heterotrophic respiration from the time of the forest fire, however, are less well characterized. The FireResp model predicts soil autotrophic and heterotrophic respiration parameterized with a novel dataset across a fire chronosequence in the Yukon and Northwest Territories of Canada. The dataset consisted of soil incubation experiments and field measurements of soil respiration and soil carbon stocks. The FireResp model contains submodels that consider a Q10 (exponential) model of respiration compared to models of heterotrophic respiration using Michaelis–Menten kinetics parameterized with soil microbial carbon. For model evaluation we applied the Akaike information criterion and compared predicted patterns in components of soil respiration across the chronosequence. Parameters estimated with data from the 5 cm soil depth had better model–data comparisons than parameters estimated with data from the 10 cm soil depth. The model–data fit was improved by including parameters estimated from soil incubation experiments. Models that incorporated microbial carbon with Michaelis–Menten kinetics reproduced patterns in autotrophic and heterotrophic soil respiration components across the chronosequence. Autotrophic respiration was associated with aboveground tree biomass at more recently burned sites, but this association was less robust at older sites in the chronosequence. Our results provide support for more structured soil respiration models than standard Q10 exponential models.

Description: The photo-oxidation of myrcene, a monoterpene species emitted by plants, was investigated at atmospheric conditions in the outdoor simulation chamber SAPHIR (Simulation of Atmospheric PHotochemistry In a Large Reaction Chamber). The chemical structure of myrcene consists of one moiety that is a conjugated π system (similar to isoprene) and another moiety that is a triple-substituted olefinic unit (similar to 2-methyl-2-butene). Hydrogen shift reactions of organic peroxy radicals (RO2) formed in the reaction of isoprene with atmospheric OH radicals are known to be of importance for the regeneration of OH. Structure–activity relationships (SARs) suggest that similar hydrogen shift reactions like in isoprene may apply to the isoprenyl part of RO2 radicals formed during the OH oxidation of myrcene. In addition, SAR predicts further isomerization reactions that would be competitive with bimolecular RO2 reactions for chemical conditions that are typical for forested environments with low concentrations of nitric oxide. Assuming that OH peroxy radicals can rapidly interconvert by addition and elimination of O2 like in isoprene, bulk isomerization rate constants of 0.21 and 0.097 s−1 (T=298 K) for the three isomers resulting from the 3′-OH and 1-OH addition, respectively, can be derived from SAR. Measurements of radicals and trace gases in the experiments allowed us to calculate radical production and destruction rates, which are expected to be balanced. The largest discrepancies between production and destruction rates were found for RO2. Additional loss of organic peroxy radicals due to isomerization reactions could explain the observed discrepancies. The uncertainty of the total radical (ROx=OH+HO2+RO2) production rates was high due to the uncertainty in the yield of radicals from myrcene ozonolysis. However, results indicate that radical production can only be balanced if the reaction rate constant of the reaction between hydroperoxy (HO2) and RO2 radicals derived from myrcene is lower (0.9 to 1.6×10-11 cm3 s−1) than predicted by SAR. Another explanation of the discrepancies would be that a significant fraction of products (yield: 0.3 to 0.6) from these reactions include OH and HO2 radicals instead of radical-terminating organic peroxides. Experiments also allowed us to determine the yields of organic oxidation products acetone (yield: 0.45±0.08) and formaldehyde (yield: 0.35±0.08). Acetone and formaldehyde are produced from different oxidation pathways, so that yields of these compounds reflect the branching ratios of the initial OH addition to myrcene. Yields determined in the experiments are consistent with branching ratios expected from SAR. The yield of organic nitrate was determined from the gas-phase budget analysis of reactive oxidized nitrogen in the chamber, giving a value of 0.13±0.03. In addition, the reaction rate constant for myrcene + OH was determined from the measured myrcene concentration, yielding a value of (2.3±0.3)×10-10 cm3 s−1.

Description: A new luminescence erosion meter has huge potential for inferring erosion rates on sub-millennial scales for both steady and transient states of erosion, which is not currently possible with any existing techniques capable of measuring erosion. This study applies new rock luminescence techniques to a well-constrained scenario provided by the Beinn Alligin rock avalanche, NW Scotland. Boulders in this deposit are lithologically consistent and have known cosmogenic nuclide ages and independently derived Holocene erosion rates. We find that luminescence-derived exposure ages for the Beinn Alligin rock avalanche were an order of magnitude younger than existing cosmogenic nuclide exposure ages, suggestive of high erosion rates (as supported by field evidence of quartz grain protrusions on the rock surfaces). Erosion rates determined by luminescence were consistent with independently derived rates measured from boulder edge roundness. Inversion modelling indicates a transient state of erosion reflecting the stochastic nature of erosional processes over the last ∼4.5 kyr in the wet, temperate climate of NW Scotland. Erosion was likely modulated by known fluctuations in moisture availability and to a lesser extent temperature, which controlled the extent of chemical weathering of these highly lithified rocks prior to erosion. The use of a multi-elevated temperature, post-infra-red, infra-red stimulated luminescence (MET-pIRIR) protocol (50, 150 and 225 ∘C) was advantageous as it identified samples with complexities that would not have been observed using only the standard infra-red stimulated luminescence (IRSL) signal measured at 50 ∘C, such as that introduced by within-sample variability (e.g. surficial coatings). This study demonstrates that the luminescence erosion meter can infer accurate erosion rates on sub-millennial scales and identify transient states of erosion (i.e. stochastic processes) in agreement with independently derived erosion rates for the same deposit.

Description: The earthquake early warning systems (EEWSs) in China have achieved great progress, with warning alerts being successfully delivered to the public in some regions. We examined the performance of the EEWS in China's Sichuan Province during the 2019 Changning earthquake. Although its technical effectiveness was tested with the first alert released 10 s after the earthquake, we found that a big gap existed between the EEWS's message and the public's response. We highlight the importance of EEWS alert effectiveness and public participation for long-term resiliency, such as delivering useful alert messages through appropriate communication channels and training people to understand and properly respond.

Description: The global distribution of cropping intensity (CI) is essential to our understanding of agricultural land use management on Earth. Optical remote sensing has revolutionized our ability to map CI over large areas in a repeated and cost-efficient manner. Previous studies have mainly focused on investigating the spatiotemporal patterns of CI ranging from regions to the entire globe with the use of coarse-resolution data, which are inadequate for characterizing farming practices within heterogeneous landscapes. To fill this knowledge gap, in this study, we utilized multiple satellite data to develop a global, spatially continuous CI map dataset at 30 m resolution (GCI30). Accuracy assessments indicated that GCI30 exhibited high agreement with visually interpreted validation samples and in situ observations from the PhenoCam network. We carried out both statistical and spatial comparisons of GCI30 with six existing global CI estimates. Based on GCI30, we estimated that the global average annual CI during 2016–2018 was 1.05, which is close to the mean (1.09) and median (1.07) CI values of the existing six global CI estimates, although the spatial resolution and temporal coverage vary significantly among products. A spatial comparison with two satellite-based land surface phenology products further suggested that GCI30 was not only capable of capturing the overall pattern of global CI but also provided many spatial details. GCI30 indicated that single cropping was the primary agricultural system on Earth, accounting for 81.57 % (12.28×106 km2) of the world's cropland extent. Multiple-cropping systems, on the other hand, were commonly observed in South America and Asia. We found large variations across countries and agroecological zones, reflecting the joint control of natural and anthropogenic drivers on regulating cropping practices. As the first global-coverage, fine-resolution CI product, GCI30 is expected to fill the data gap for promoting sustainable agriculture by depicting worldwide diversity of agricultural land use intensity. The GCI30 dataset is available on Harvard Dataverse: https://doi.org/10.7910/DVN/86M4PO (Zhang et al., 2020).

Description: McMurdo Sound sea ice can generally be partitioned into two regimes: (1) a stable fast-ice cover, forming south of approximately 77.6∘ S around March–April and then breaking out the following January–February, and (2) a more dynamic region north of 77.6∘ S that the McMurdo Sound and Ross Sea polynyas regularly impact. In 2019, a stable fast-ice cover formed unusually late due to repeated break-out events. We analyse the 2019 sea-ice conditions and relate them to a modified storm index (MSI), a proxy for southerly wind events. We find there is a strong correlation between the timing of break-out events and several unusually large MSI events.

Description: The Arabian Sea (AS) hosts one of the most intense oxygen minimum zones (OMZs) in the world. Observations suggest a decline in O2 in the northern AS over the recent decades accompanied by an intensification of the suboxic conditions there. Over the same period, the local sea surface temperature has risen significantly, particularly over the Arabian Gulf (also known as Persian Gulf, hereafter the Gulf), while summer monsoon winds may have intensified. Here, we simulate the evolution of dissolved oxygen in the AS from 1982 through 2010 and explore its controlling factors, with a focus on changing atmospheric conditions. To this end, we use a set of eddy-resolving hindcast simulations forced with winds and heat and freshwater fluxes from an atmospheric reanalysis. We find a significant deoxygenation in the northern AS, with O2 inventories north of 20∘ N dropping by over 6 % per decade between 100 and 1000 m. These changes cause an expansion of the OMZ volume north of 20∘ N at a rate of 0.6 % per decade as well as an increase in the volume of suboxia and the rate of denitrification by 14 and 15 % per decade, respectively. We also show that strong interannual and decadal variability modulate dissolved oxygen in the northern AS, with most of the O2 decline taking place in the 1980s and 1990s. Using a set of sensitivity simulations we demonstrate that deoxygenation in the northern AS is essentially caused by reduced ventilation induced by the recent fast warming of the sea surface, including in the Gulf, with a contribution from concomitant summer monsoon wind intensification. This is because, on the one hand, surface warming enhances vertical stratification and increases Gulf water buoyancy, thus inhibiting vertical mixing and ventilation of the thermocline. On the other hand, summer monsoon wind intensification causes a rise in the thermocline depth in the northern AS that lowers O2 levels in the upper ocean. Our findings confirm that the AS OMZ is strongly sensitive to upper-ocean warming and concurrent changes in the Indian monsoon winds. Finally, our results also demonstrate that changes in the local climatic forcing play a key role in regional dissolved oxygen changes and hence need to be properly represented in global models to reduce uncertainties in future projections of deoxygenation.

Description: The thermokarst lakes of permafrost regions play a major role in the global carbon cycle. These lakes are sources of methane to the atmosphere although the methane flux is restricted by an ice cover for most of the year. How methane concentrations and fluxes in these waters are affected by the presence of an ice cover is poorly understood. To relate water body morphology, ice formation and methane to each other, we studied the ice of three different water bodies in locations typical of the transition of permafrost from land to ocean in a continuous permafrost coastal region in Siberia. In total, 11 ice cores were analyzed as records of the freezing process and methane composition during the winter season. The three water bodies differed in terms of connectivity to the sea, which affected fall freezing. The first was a bay underlain by submarine permafrost (Tiksi Bay, BY), the second a shallow thermokarst lagoon cut off from the sea in winter (Polar Fox Lagoon, LG) and the third a land-locked freshwater thermokarst lake (Goltsovoye Lake, LK). Ice on all water bodies was mostly methane-supersaturated with respect to atmospheric equilibrium concentration, except for three cores from the isolated lake. In the isolated thermokarst lake, ebullition from actively thawing basin slopes resulted in the localized integration of methane into winter ice. Stable δ13C-CH4 isotope signatures indicated that methane in the lagoon ice was oxidized to concentrations close to or below the calculated atmospheric equilibrium concentration. Increasing salinity during winter freezing led to a micro-environment on the lower ice surface where methane oxidation occurred and the lagoon ice functioned as a methane sink. In contrast, the ice of the coastal marine environment was slightly supersaturated with methane, consistent with the brackish water below. Our interdisciplinary process study shows how water body morphology affects ice formation which mitigates methane fluxes to the atmosphere.

Description: Methane emissions from boreal and arctic wetlands, lakes, and rivers are expected to increase in response to warming and associated permafrost thaw. However, the lack of appropriate land cover datasets for scaling field-measured methane emissions to circumpolar scales has contributed to a large uncertainty for our understanding of present-day and future methane emissions. Here we present the Boreal–Arctic Wetland and Lake Dataset (BAWLD), a land cover dataset based on an expert assessment, extrapolated using random forest modelling from available spatial datasets of climate, topography, soils, permafrost conditions, vegetation, wetlands, and surface water extents and dynamics. In BAWLD, we estimate the fractional coverage of five wetland, seven lake, and three river classes within 0.5 × 0.5∘ grid cells that cover the northern boreal and tundra biomes (17 % of the global land surface). Land cover classes were defined using criteria that ensured distinct methane emissions among classes, as indicated by a co-developed comprehensive dataset of methane flux observations. In BAWLD, wetlands occupied 3.2 × 106 km2 (14 % of domain) with a 95 % confidence interval between 2.8 and 3.8 × 106 km2. Bog, fen, and permafrost bog were the most abundant wetland classes, covering ∼ 28 % each of the total wetland area, while the highest-methane-emitting marsh and tundra wetland classes occupied 5 % and 12 %, respectively. Lakes, defined to include all lentic open-water ecosystems regardless of size, covered 1.4 × 106 km2 (6 % of domain). Low-methane-emitting large lakes (〉10 km2) and glacial lakes jointly represented 78 % of the total lake area, while high-emitting peatland and yedoma lakes covered 18 % and 4 %, respectively. Small (〈0.1 km2) glacial, peatland, and yedoma lakes combined covered 17 % of the total lake area but contributed disproportionally to the overall spatial uncertainty in lake area with a 95 % confidence interval between 0.15 and 0.38 × 106 km2. Rivers and streams were estimated to cover 0.12  × 106 km2 (0.5 % of domain), of which 8 % was associated with high-methane-emitting headwaters that drain organic-rich landscapes. Distinct combinations of spatially co-occurring wetland and lake classes were identified across the BAWLD domain, allowing for the mapping of “wetscapes” that have characteristic methane emission magnitudes and sensitivities to climate change at regional scales. With BAWLD, we provide a dataset which avoids double-accounting of wetland, lake, and river extents and which includes confidence intervals for each land cover class. As such, BAWLD will be suitable for many hydrological and biogeochemical modelling and upscaling efforts for the northern boreal and arctic region, in particular those aimed at improving assessments of current and future methane emissions. Data are freely available at https://doi.org/10.18739/A2C824F9X (Olefeldt et al., 2021).

Description: Lakes in permafrost regions are dynamiclandscape components and play an important role for climatechange feedbacks. Lake processes such as mineralizationand flocculation of dissolved organic carbon (DOC), oneof the main carbon fractions in lakes, contribute to thegreenhouse effect and are part of the global carbon cycle.These processes are in the focus of climate research, butstudies so far are limited to specific study regions. Inour synthesis, we analyzed 2167 water samples from 1833lakes across the Arctic in permafrost regions of Alaska,Canada, Greenland, and Siberia to provide first pan-Arcticinsights for linkages between DOC concentrations andthe environment. Using published data and unpublisheddatasets from the author team, we report regional DOCdifferences linked to latitude, permafrost zones, ecoregions,geology, near-surface soil organic carbon contents, andground ice classification of each lake region. The lakeDOC concentrations in our dataset range from 0 to1130 mg L−1(10.8 mg L−1median DOC concentration).Regarding the permafrost regions of our synthesis, wefound median lake DOC concentrations of 12.4 mg L−1(Siberia), 12.3 mg L−1(Alaska), 10.3 mg L−1(Greenland),and 4.5 mg L−1(Canada). Our synthesis shows a significantrelationship between lake DOC concentration and lakeecoregion. We found higher lake DOC concentrationsat boreal permafrost sites compared to tundra sites. Wefound significantly higher DOC concentrations in lakesin regions with ice-rich syngenetic permafrost deposits(yedoma) compared to non-yedoma lakes and a weak butsignificant relationship between soil organic carbon contentand lake DOC concentration as well as between ground icecontent and lake DOC. Our pan-Arctic dataset shows that theDOC concentration of a lake depends on its environmentalproperties, especially on permafrost extent and ecoregion, aswell as vegetation, which is the most important driver of lakeDOC in this study. This new dataset will be fundamental toquantify a pan-Arctic lake DOC pool for estimations of theimpact of lake DOC on the global carbon cycle and climatechange.

Description: 〈jats:p〉Abstract. The risk of carbon emissions from permafrost ground is linked to ground temperature and thus in particular to thermal insulation by vegetation and organic soil layers in summer and snow cover in winter. This ground insulation is strongly influenced by the presence of large herbivorous animals browsing for food. In this study, we examine the potential impact of large herbivore presence on the ground carbon storage in thermokarst landscapes of northeastern Siberia. Our aim is to understand how intensive animal grazing may affect permafrost thaw and hence organic matter decomposition, leading to different ground carbon storage, which is significant in the active layer. Therefore, we analysed sites with differing large herbivore grazing intensity in the Pleistocene Park near Chersky and measured maximum thaw depth, total organic carbon content and decomposition state by δ13C isotope analysis. In addition, we determined sediment grain size composition as well as ice and water content. We found the thaw depth to be shallower and carbon storage to be higher in intensively grazed areas compared to extensively and non-grazed sites in the same thermokarst basin. The intensive grazing presumably leads to a more stable thermal ground regime and thus to increased carbon storage in the thermokarst deposits and active layer. However, the high carbon content found within the upper 20 cm on intensively grazed sites could also indicate higher carbon input rather than reduced decomposition, which requires further studies. We connect our findings to more animal trampling in winter, which causes snow disturbance and cooler winter ground temperatures during the average annual 225 days below freezing. This winter cooling overcompensates ground warming due to the lower insulation associated with shorter heavily grazed vegetation during the average annual 140 thaw days. We conclude that intensive grazing influences the carbon storage capacities of permafrost areas and hence might be an actively manageable instrument to reduce net carbon emission from these sites. 〈/jats:p〉

Description: 〈jats:p〉Abstract. The combined effect of hot and dry extremes can have disastrous consequences for the society, economy, and the environment. While a significant number of studies have been conducted regarding the variability of the individual hot or dry extremes in Romania, the evaluation of the combined effect of these extremes (e.g. compound effect) is still lacking for this region. Thus, in this study we have assessed the spatio-temporal variability and trends of hot and dry summers in the eastern part of Europe, focusing on Romania, between 1950 and 2020 and we have analyzed the relationship between the frequency of hot summers and the prevailing large-scale atmospheric circulation. The length, spatial extent and frequency of HWs in Romania has increased significantly over the last 70 years, while for the drought conditions no significant changes have been observed. The rate of increase in the frequency and spatial extent of HWs has accelerated significantly after the 1990’s, while the smallest number of HWs was observed between 1970 and 1985. The hottest years, in terms of heatwave duration and frequency, were 2007, 2012, 2015, and 2019. One of the key drivers of hot summers, over our analyzed region, is the prevailing large-scale circulation, featuring an anticyclonic circulation over the central and eastern parts of Europe and enhanced atmospheric blocking activity associated with positive temperature anomalies underneath. We conclude that our study can help improve our understanding of the spatio-temporal variability of hot and dry summers, especially at the regional scale, as well as their driving mechanisms which might lead to a better predictability of these extreme events. 〈/jats:p〉

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: In this study, the first fully continuous monitoring of water vapour isotopic composition at Neumayer Station III, Antarctica, during the 2-year period from February 2017 to January 2019 is presented. Seasonal and synoptic-scale variations in both stable water isotopes H182O and HDO are reported, and their links to variations in key meteorological variables are analysed. In addition, the diurnal cycle of isotope variations during the summer months (December and January 2017/18 and 2018/19) has been examined. Changes in local temperature and specific humidity are the main drivers for the variability in δ18O and δD in vapour at Neumayer Station III, on both seasonal and shorter timescales. In contrast to the measured δ18O and δD variations, no seasonal cycle in the Deuterium excess signal (d) in vapour is detected. However, a rather high uncertainty in measured d values especially in austral winter limits the confidence of this finding. Overall, the d signal shows a stronger inverse correlation with specific humidity than with temperature, and this inverse correlation between d and specific humidity is stronger for the cloudy-sky conditions than for clear-sky conditions during summertime. Back-trajectory simulations performed with the FLEXPART model show that seasonal and synoptic variations in δ18O and δD in vapour coincide with changes in the main sources of water vapour transported to Neumayer Station III. In general, moisture transport pathways from the east lead to higher temperatures and more enriched δ18O values in vapour, while weather situations with southerly winds lead to lower temperatures and more depleted δ18O values. However, on several occasions, δ18O variations linked to wind direction changes were observed, which were not accompanied by a corresponding temperature change. Comparing isotopic compositions of water vapour at Neumayer Station III and snow samples taken in the vicinity of the station reveals almost identical slopes, both for the δ18O–δD relation and for the temperature–δ18O relation.

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: 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.