Climate change and North Sea storm surge extremes: an ensemble study of storm surge extremes expected in a changed climate projected by four different regional climate models

The coastal zones are facing the prospect of changing storm surge statistics due to anthropogenic climate change. In the present study, we examine these prospects for the North Sea based on numerical modelling. The main tool is the barotropic tide-surge model TRIMGEO (Tidal Residual and Intertidal Mudflat Model) to derive storm surge climate and extremes from atmospheric conditions. The analysis is carried out by using an ensemble of four 30-year atmospheric regional simulations under present-day and possible future-enhanced greenhouse gas conditions. The atmospheric regional simulations were prepared within the EU project PRUDENCE (Prediction of Regional scenarios and Uncertainties for Defining EuropeaN Climate change risks and Effects). The research strategy of PRUDENCE is to compare simulations of different regional models driven by the same global control and climate change simulations. These global conditions, representative for 1961–1990 and 2071–2100 were prepared by the Hadley Center based on the IPCC A2 SRES scenario. The results suggest that under future climatic conditions, storm surge extremes may increase along the North Sea coast towards the end of this century. Based on a comparison between the results of the different ensemble members as well as on the variability estimated from a high-resolution storm surge reconstruction of the recent decades it is found that this increase is significantly different from zero at the 95% confidence level for most of the North Sea coast. An exception represents the East coast of the UK which is not affected by this increase of storm surge extremes.     

Sensitivity of transient climate change to tidal mixing: Southern Ocean heat uptake in climate change experiments performed with ECHAM5/MPIOM

We investigate the sensitivity of the transient climate change to a tidal mixing scheme. The scheme parameterizes diapycnal diffusivity depending on the location of energy dissipation over rough topography, whereas the standard configuration uses horizontally constant diffusivity. We perform ensemble climate change experiments with two setups of MPIOM/ECHAM5, one setup with the tidal mixing scheme and the second setup with the standard configuration. Analysis of the responses of the transient climate change to CO₂ increase reveals that the implementation of tidal mixing leads to a significant reduction of the transient surface warming by 9 %. The weaker surface warming in the tidal run is localized particularly over the Weddell Sea, likely caused by a stronger ocean heat uptake in the Southern Ocean. The analysis of the ocean heat budget reveals that the ocean heat uptake in both experiments is caused by changes in convection and advection. In the upper ocean, heat uptake is caused by reduced convection and enhancement of the Deacon Cell, which appears also in isopycnal coordinates. In the deeper ocean, heat uptake is caused by reduction of convective cooling associated with the circulation polewards of 65°S. Tidal mixing leads to stronger heat uptake in the Southern Ocean by causing stronger changes in advection, namely a stronger increase in the Deacon Cell and a stronger reduction in advective cooling by the circulation polewards of 65°S. Counter-intuitively, the relation between tidal mixing and greater heat storage in the deep ocean is an indirect one, through the influence of tidal mixing on the circulation.     

The seasonal cycle in coupled ocean-atmosphere general circulation models

We examine the seasonal cycle of near-surface air temperature simulated by 17 coupled ocean-atmosphere general circulation models participating in the Coupled Model Intercomparison Project (CMIP). Nine of the models use ad hoc “flux adjustment” at the ocean surface to bring model simulations close to observations of the present-day climate. We group flux-adjusted and non-flux-adjusted models separately and examine the behavior of each class. When averaged over all of the flux-adjusted model simulations, near-surface air temperature falls within 2?K of observed values over the oceans. The corresponding average over non-flux-adjusted models shows errors up to ～6?K in extensive ocean areas. Flux adjustments are not directly applied over land, and near-surface land temperature errors are substantial in the average over flux-adjusted models, which systematically underestimates (by ～5?K) temperature in areas of elevated terrain. The corresponding average over non-flux-adjusted models forms a similar error pattern (with somewhat increased amplitude) over land. We use the temperature difference between July and January to measure seasonal cycle amplitude. Zonal means of this quantity from the individual flux-adjusted models form a fairly tight cluster (all within ～30% of the mean) centered on the observed values. The non-flux-adjusted models perform nearly as well at most latitudes. In Southern Ocean mid-latitudes, however, the non-flux-adjusted models overestimate the magnitude of January-minus-July temperature differences by ～5?K due to an overestimate of summer (January) near-surface temperature. This error is common to five of the eight non-flux-adjusted models. Also, over Northern Hemisphere mid-latitude land areas, zonal mean differences between July and January temperatures simulated by the non-flux-adjusted models show a greater spread (positive and negative) about observed values than results from the flux-adjusted models. Elsewhere, differences between the two classes of models are less obvious. At no latitude is the zonal mean difference between averages over the two classes of models greater than the standard deviation over models. The ability of coupled GCMs to simulate a reasonable seasonal cycle is a necessary condition for confidence in their prediction of long-term climatic changes (such as global warming), but it is not a sufficient condition unless the seasonal cycle and long-term changes involve similar climatic processes. To test this possible connection, we compare seasonal cycle amplitude with equilibrium warming under doubled atmospheric carbon dioxide for the models in our data base. A small but positive correlation exists between these two quantities. This result is predicted by a simple conceptual model of the climate system, and it is consistent with other modeling experience, which indicates that the seasonal cycle depends only weakly on climate sensitivity.     

Improved magnetic anomalies of the Antarctic lithosphere from satellite and near-surface data

The Antarctic magnetic anomaly map compiled marine and airborne surveys collected south of 60°S through 1999 and used Magsat data to help fill in the regional gaps between the surveys. Ørsted and CHAMP satellite magnetic observations with greatly improved measurement accuracies and temporal and spatial coverage of the Antarctic, have now supplanted the Magsat data. We combined the new satellite observations with the near-surface survey data for an improved magnetic anomaly map of the Antarctic lithosphere. Specifically, we separated the crustal from the core and external field components in the satellite data using crustal thickness variations estimated from the terrain and the satellite-derived free-air gravity observations. Regional gaps in the near-surface surveys were then filled with predictions from crustal magnetization models that jointly satisfied the near-surface and satellite crustal anomalies. Comparisons in some of the regional gaps that also considered newly acquired aeromagnetic data demonstrated the enhanced anomaly estimation capabilities of the predictions over those from conventional minimum curvature and spherical harmonic geomagnetic field models. We also noted that the growing number of regional and world magnetic survey compilations involve coverage gaps where these procedures can contribute effective near-surface crustal anomaly estimates.     

Structure of Palaeogene sediments in east Ellesmere Island: Constraints on Eurekan tectonic evolution and implications for the Nares Strait problem

The “Nares Strait problem” represents a debate about the existence and magnitude of left-lateral movements along the proposed Wegener Fault within this seaway. Study of Palaeogene Eurekan tectonics at its shorelines could shed light on the kinematics of this fault. Palaeogene (Late Paleocene to Early Eocene) sediments are exposed at the northeastern coast of Ellesmere Island in the Judge Daly Promontory. They are preserved as elongate SW–NE striking fault-bounded basins cutting folded Early Paleozoic strata. The structures of the Palaeogene exposures are characterized by broad open synclines cut and displaced by steeply dipping strike-slip faults. Their fold axes strike NE–SW at an acute angle to the border faults indicating left-lateral transpression. Weak deformation in the interior of the outliers contrasts with intense shearing and fracturing adjacent to border faults. The degree of deformation of the Palaeogene strata varies markedly between the northwestern and southeastern border faults with the first being more intense. Structural geometry, orientation of subordinate folds and faults, the kinematics of faults, and fault-slip data suggest a multiple stage structural evolution during the Palaeogene Eurekan deformation: (1) The fault pattern on Judge Daly Promontory is result of left-lateral strike-slip faulting starting in Mid to Late Paleocene times. The Palaeogene Judge Daly basin formed in transtensional segments by pull-apart mechanism. Transpression during progressive strike-slip shearing gave rise to open folding of the Palaeogene deposits. (2) The faults were reactivated during SE-directed thrust tectonics in Mid Eocene times (chron 21). A strike-slip component during thrusting on the reactivated faults depends on the steepness of the fault segments and on their obliquity to the regional stress axes.Strike-slip displacement was partitioned to a number of sub-parallel faults on-shore and off-shore. Hence, large-scale lateral movements in the sum of 80–100 km or more could have been accommodated by a set of faults, each with displacements in the order of 10–30 km. The Wegener Fault as discrete plate boundary in Nares Strait is replaced by a bundle of faults located mainly onshore on the Judge Daly Promontory.     

Tectonic significance of Triassic mafic rocks in the June Complex,Sanandaj–Sirjan zone,Iran

This study concentrates on the petrological and geochemical investigation of mafic rocks embedded within the voluminous Triassic June Complex of the central Sanandaj–Sirjan zone (Iran), which are crucial to reconstruct the geodynamics of the Neotethyan passive margin. The Triassic mafic rocks are alkaline to sub-alkaline basalts, containing 43.36–49.09 wt% SiO₂, 5.19–20.61 wt% MgO and 0.66–4.59 wt% total alkalis. Based on MgO concentrations, the mafic rocks fall into two groups: cumulates (Mg# = 51.61–58.94) and isotropic basaltic liquids (Mg# = 24.54–42.66). In all samples, the chondrite-normalized REE patterns show enrichment of light REEs with variable (La/Yb)_N ratios ranging from 2.48 to 9.00, which confirm their amalgamated OIB-like and E-MORB-like signatures. Enrichment in large-ion lithophile elements and depletion in high field strength elements (HFSE) relative to the primitive mantle further support this interpretation. No samples point to crustal contamination, all having undergone fractionation of olivine + clinopyroxene + plagioclase. Nevertheless, elemental data suggest that the substantial variations in (La/Sm)_PM and Zr/Nb ratios can be explained by variable degrees of partial melting rather than fractional crystallization from a common parental magma. The high (Nb/Yb)_PM ratio in the alkaline mafic rocks points to the mixing of magmas from enriched and depleted mantle sources. Abundant OIB alkaline basalts and rare E-MORB appear to be linked to the drifting stage on the northern passive margin of the Neotethys Ocean.     

Geochronological and isotopic records of crustal storage and assimilation in the Wolverine Creek–Conant Creek system,Heise eruptive centre,Snake River Plain

Understanding the processes of differentiation of the Yellowstone–Snake River Plain (YSRP) rhyolites is typically impeded by the apparent lack of erupted intermediate compositions as well as the complex nature of their shallow interaction with the surrounding crust responsible for their typically low O isotopic ratios. A pair of normal-δ¹⁸O rhyolitic eruptions from the Heise eruptive centre in eastern Idaho, the Wolverine Creek Tuff and the Conant Creek Tuff, represent unique magmatic products of the Yellowstone hotspot preserving abundant vestiges of the intermediate differentiation steps leading to rhyolite generation. We address both shallow and deep processes of magma generation and storage in the two units by combining high-precision ID–TIMS U–Pb zircon geochronology, trace element, O and Hf isotopic studies of zircon, and Sr isotopic analyses of individual high-Mg# pyroxenes inherited from lower- to mid-crustal differentiation stages. The zircon geochronology confirms the derivation of both tuffs from the same rhyolitic magma reservoir erupted at 5.5941 ± 0.0097 Ma, preceded by at least 92 ± 14 ky of continuous or intermittent zircon saturation approximating the length of pre-eruptive magma accumulation in the upper crust. Some low-Mg# pyroxenes enclosing zircons predate the eruption by at least 45 ± 27 ky, illustrating the co-crystallisation of major and accessory phases in the near-liquidus rhyolitic melts of the YSRP over a significant period of time. Coeval zircon crystals are isotopically heterogeneous (two populations at εHf ~?5 and ?13), requiring the assembly of isotopically distinct melt pockets directly prior to, or during, the eruption. The primitive Mg# 60–90 pyroxenes are out of isotopic equilibrium with the host rhyolitic melt (⁸⁷Sr/⁸⁶Sr_i = 0.70889), covering a range of ⁸⁷Sr/⁸⁶Sr_i = 0.70705–0.70883 corresponding to ratios typical of the most radiogenic YSRP basalts to the least radiogenic YSRP rhyolites. Together with the low εHf in zircon, the Sr isotopic ratios illustrate limited assimilation dominated by radiogenic Archean crustal source materials incorporated into variably evolved YSRP melts as they progress towards rhyolitic compositions by assimilation–fractional crystallisation.