An Mw based earthquake Catalogue for central, northern and northwestern Europe using a hierarchy of magnitude conversions

Data from 25 local catalogues and 30special studies of earthquakes in central,northern and northwestern Europe have beenincorporated into a Databank. The dataprocessing includes discriminating eventtypes, eliminating fake events and dupletsand converting different magnitudes andintensities to M_w if this is not givenby the original source. The magnitudeconversion is a key task of the study andimplies establishment of regressionequations where no local relations exist.The Catalogue contains tectonic events fromthe Databank within the area44°N–72°N,25°W–32°E and the time period1300–1993. The lower magnitude level forthe Catalogue entries is setat M_w == 3.50. The area covered by thedifferent catalogues are associated withpolygons. Within each polygon only datafrom one or a small number of the localcatalogues, supplemented by data fromspecial studies, enter the Catalogue. Ifthere are two or more such catalogues orstudies providing a solution for an event,a priority algorithm selects one entry forthe Catalogue. Then M_w is calculatedfrom one of the magnitude types, or frommacroseismic data, given by the selectedentry according to another priority scheme.The origin time, location, M_w magnitude and reference are specified for eachentry of the Catalogue. So is theepicentral intensity, I₀, if providedby the original source. Following thesecriteria, a total of about 5,000earthquakes constitute the Catalogue.Although originally derived for the purposeof seismic hazard calculation within GSHAP,the Catalogue provides a data base for manytypes of seismicity and seismic hazardstudies.     

The early Caledonian (Finnmarkian) event reassessed in Finnmark: Ar/Ar cleavage age data from NW Varangerhalvøya, N. Norway

The relative significance of early (Finnmarkian) and late (Scandian) Caledonian deformation in N. Norway is uncertain. Early studies suggested pervasive Finnmarkian deformation whilst later results indicated a restricted Finnmarkian domain. The present work suggests it was more widespread than accepted and that inter Finnmarkian–Scandian deformation occurred. ⁴⁰Ar/³⁹Ar dating of 2–6 and 6–11 μm pelitic fractions from the lower to mid-greenschist facies Tanahorn Nappe (five samples; base Middle Allochthon) and the epizone Løkvikfjellet and Barents Sea Groups (three samples; North Varanger Region) in the north Scandinavian Caledonides show slightly discordant spectra. Most spectra from the Tanahorn Nappe preserve possible evidence of an early Caledonian event in the high temperature steps, with recoil/excess Ar effects in the low temperature steps; no pre-Caledonian relict component has been recorded. The results indicate Finnmarkian deformation continued to 460 Ma, with Scandian reactivation at 425–415 Ma. From the North Varanger Region, a strongly crenulated sample yielded plateau ages (444–442 Ma); means of combined young steps from weakly to uncrenulated samples gave 470–450 Ma, suggesting penetrative strike-slip deformation occurred in the late Finnmarkian to inter-Finnmarkian–Scandian period. No Scandian ages were recorded in the North Varanger Region. Reassessment of published data from the Laksefjord Nappe and Gaissa Thrust Belt suggests they were affected by Finnmarkian deformation.     

Geoid determination using adapted reference field, seismic Moho depths and variable density contrast

 The traditional remove-restore technique for geoid computation suffers from two main drawbacks. The first is the assumption of an isostatic hypothesis to compute the compensation masses. The second is the double consideration of the effect of the topographic–isostatic masses within the data window through removing the reference field and the terrain reduction process. To overcome the first disadvantage, the seismic Moho depths, representing, more or less, the actual compensating masses, have been used with variable density anomalies computed by employing the topographic–isostatic mass balance principle. In order to avoid the double consideration of the effect of the topographic–isostatic masses within the data window, the effect of these masses for the used fixed data window, in terms of potential coefficients, has been subtracted from the reference field, yielding an adapted reference field. This adapted reference field has been used for the remove–restore technique. The necessary harmonic analysis of the topographic–isostatic potential using seismic Moho depths with variable density anomalies is given. A wide comparison among geoids computed by the adapted reference field with both the Airy–Heiskanen isostatic model and seismic Moho depths with variable density anomaly and a geoid computed by the traditional remove–restore technique is made. The results show that using seismic Moho depths with variable density anomaly along with the adapted reference field gives the best relative geoid accuracy compared to the GPS/levelling geoid. Received: 3 October 2001 / Accepted: 20 September 2002 Correspondence to: H.A. Abd-Elmotaal     

Intercomparison of Stratospheric Chemistry Models under Polar Vortex Conditions

Several stratospheric chemistry modules from box, 2-D or 3-D models, have been intercompared. The intercomparison was focused on the ozone loss and associated reactive species under the conditions found in the cold, wintertime Arctic and Antarctic vortices. Comparisons of both gas phase and heterogeneous chemistry modules show excellent agreement between the models under constrained conditions for photolysis and the microphysics of polar stratospheric clouds. While the mean integral ozone loss ranges from 4–80% for different 30–50 days long air parcel trajectories, the mean scatter of model results around these values is only about ±1.5%. In a case study, where the models employed their standard photolysis and microphysical schemes, the variation around the mean percentage ozone loss increases to about ±7%. This increased scatter of model results is mainly due to the different treatment of the PSC microphysics and heterogeneous chemistry in the models, whereby the most unrealistic assumptions about PSC processes consequently lead to the least representative ozone chemistry. Furthermore, for this case study the model results for the ozone mixing ratios at different altitudes were compared with a measured ozone profile to investigate the extent to which models reproduce the stratospheric ozone losses. It was found that mainly in the height range of strong ozone depletion all models underestimate the ozone loss by about a factor of two. This finding corroborates earlier studies and implies a general deficiency in our understanding of the stratospheric ozone loss chemistry rather than a specific problem related to a particular model simulation.     

A Case Study Of An On-Ice Air Flow Over The Arctic Marginal Sea-Ice Zone

A case study of warm air advection over the Arctic marginalsea-ice zone is presented, based on aircraft observations with direct flux measurements carriedout in early spring, 1998. A shallow atmospheric boundary layer (ABL) was observed, which wasgradually cooling with distance downwind of the ice edge. This process was mainly connected with astrong stable stratification and downward turbulent heat fluxes of about 10–20 W m^-2, but wasalso due to radiative cooling. Two mesoscale models, one hydrostatic and the other non-hydrostatic,having different turbulence closures, were applied. Despite these fundamental differences betweenthe models, the results of both agreed well with the observed data. Various closure assumptions had amore crucial influence on the results than the differences between the models.Such an assumption was, for example,the parameterization of the surface roughness for momentum (z₀) and heat (z_T). This stronglyaffected the wind and temperature fields not only close to the surface but also within and abovethe temperature inversion layer. The best results were achieved using a formulation for z₀ that took intoaccount the form drag effect of sea-ice ridges together withz_T = 0.1z₀. The stability within theelevated inversion strongly depended on the minimum eddy diffusivity K_min. A simple ad hocparameterization seems applicable, where K_min is calculated as 0.005 timesthe neutral eddy diffusivity. Although the longwave radiative cooling was largest within the ABL, theapplication of a radiation scheme was less important there than above the ABL. This was related to theinteraction of the turbulent and radiative fluxes. To reproduce the strong inversion, it wasnecessary to use vertical and horizontal resolutions higher than those applied in most regional andlarge-scale atmospheric models.     

Transport in a 2-D saturated porous medium: A new method for particle tracking

A new method for solving the transport equation based on the management of a large numbe of particles in a discretized 2-D domain is presented. The method uses numerical variables to represent the number of particles in a given mesh and is more complex than the 1-D problem. The first part of the paper focuses on the specific management of particles in a 2-D problem. The method also would be valid for three dimensions as long as the medium can be modeled similar to a layered system. As the particles are no longer tracked individually, the algorithm is fast and does not depend on the number of particles present. The numerical tests show that the method is nearly numerical dispersion free and permits accurate calculations even for simulations of low-concentration transport. Because each mesh is considered as a closed system between two successive time steps, it is easy to add adsorption phenomenon without any problem of numerical stability. The model is tested under conditions that are extremely demanding for its operating mode and gives a good fit to analytical solutions. The conditions in which it can be used to best advantage are discussed.     

Thorium and uranium isotopes in a manganese nodule from the Peru basin determined by alpha spectrometry and thermal ionization mass spectrometry (TIMS): Are manganese supply and growth related to climate?

Thorium- and uranium isotopes were measured in a diagenetic manganese nodule from the Peru basin applying alpha- and thermal ionization mass spectrometry (TIMS). Alpha-counting of 62 samples was carried out with a depth resolution of 0.4 mm to gain a high-resolution²³⁰Th_excess profile. In addition, 17 samples were measured with TIMS to obtain precise isotope concentrations and isotope ratios. We got values of 0.06–0.59 ppb (²³⁰Th), 0.43–1.40 ppm (²³²Th), 0.09–0.49 ppb (²³⁴U) and 1.66–8.24 ppm (²³⁸U). The uranium activity ratio in the uppermost samples (1–6 mm) and in two further sections in the nodule at 12.5±1.0 mm and 27.3–33.5 mm comes close to the present ocean water value of 1.144±0.004. In two other sections of the nodule, this ratio is significantly higher, probably reflecting incorporation of diagenetic uranium. The upper 25 mm section of the Mn nodule shows a relatively smooth exponential decrease in the²³⁰Th_excess concentration (TIMS). The slope of the best fit yields a growth rate of 110 mm/Ma up to 24.5 mm depth. The section from 25 to 30.3 mm depth shows constant²³⁰Th_excess concentrations probably due to growth rates even faster than those in the top section of the nodule. From 33 to 50 mm depth, the growth rate is approximately 60 mm/Ma. Two layers in the nodule with distinct laminations (11–15 and 28–33 mm depth) probably formed during the transition from isotopic stage 8 to 7 and in stage 5e, respectively. The Mn/Fe ratio shows higher values during interglacials 5 and 7, and lower ones during glacials 4 and 6. A comparison of our data with data from adjacent sediment cores suggests (a) a variable supply of hydrothermal Mn to sediments and Mn nodules of the Peru basin or (b) suboxic conditions at the water sediment interface during periods with lower Mn/Fe ratios.