Multi-species chemical data assimilation with the Danish Eulerian hemispheric model: system description and verification

Satellite retrievals of atmospheric composition provide a wealth of data on a global scale. These complement results from atmospheric chemistry-transport models (CTMs), and can be combined using data assimilation. We present two assimilation schemes coupled to the Danish Eulerian Hemispheric Model (DEHM), a three-dimensional, off-line CTM with full photochemistry: a variant on the ensemble Kalman filter and the three-dimensional variational scheme. The aim of this paper is to describe the two schemes and present an initial assessment of their impacts on model skill. Retrievals of multiple atmospheric trace gases are assimilated, namely: NO₂ tropospheric column densities, CH₄ total column densities, and partial column concentrations of O₃, CO and CH₄; these data are retrieved from four satellite sensors. Data for each species are assimilated independently of one another, and other species are only adjusted indirectly via the model’s chemistry and dynamics. Assimilation results are compared with measurements from surface monitoring stations and other satellite retrievals, and preliminary validation results are presented.Reference simulations (without assimilation) grossly underestimate surface CO concentrations, and both assimilation schemes eliminate this large and systematic model bias. The assimilation improves the spatial correlation of modelled CO with surface observations, and improves the spatial correlation between forecasts and retrievals for CO, NO₂ and O₃. Results for CH₄ show a loss of skill due to a mismatch in model bias between two assimilated CH₄ data-sets. Finally, we discuss differences in methodology and results between this paper and a recent study on multi-species chemical data assimilation. Joint optimisation of initial conditions and emission rates offers a promising direction for improving modelled boundary-layer concentrations.     

Tracer experiments with turbulently dispersed air ions

Unipolar air ions released into the wind constitute a tracer which can be measured with high resolution. An ion source produces a cloud with homogeneous charge density, insensitive to source strength, dependent on time since formation only. It is well suited for tracing concentration changes due to turbulence, less suited for cloud size tracing. A tight array of 8 sensors has been used to examine turbulently dispersed ions. High-resolution records are presented and discussed. The highest concentrations measured could be undiluted source material. The frequency distribution within a plume did not differ from that in a multitude of puffs. The distribution seems to be log-normal with a geometric standard deviation of about 2.45. The time resolution used corresponds to volume resolutions of 40, 225, and 650 cm³. Sample size had no apparent effect.     

Oxygen and nutrient dynamics within mats of the filamentous macroalga Chaetomorpha linum

Concentration profiles of O₂, NH₄ ⁺, NO₃ ⁻, and PO₄ ³⁻ were measured at high spatial resolution in a 12-cm thick benthic mat of the filamentous macroalga Chaetomorpha linum. Oxygen and nutrient concentration profiles varied depending on algal activity and water turbulence. High surface irradiance stimulated O₂ production in the surface layers and introduced O₂ to deeper parts of the mat while the bottom layers of the mat and the underlying sediment were anoxic. Nutrient concentrations were highest in the bottom layers of the mat directly above the sediment nutrient source and decreased towards the surface layers due to algal assimilation and enhanced mixing with the overlying water column. Increased turbulence during windy periods resulted in more homogeneous oxygen and nutrient concentration profiles and shifted the oxic-anoxic interface downward. Denitrification within the mat, as measured by the isotope pairing technique on addition of ¹⁵NO₃ ⁻, was found to take place directly below the oxic-anoxic interface. Denitrification activity was always due to coupled nitrification-denitrification, whereby nitrifiers in the mat utilize NH₄ ⁺ diffusing from below and O₂ diffusing from above. The denitrification rate in the mat ranged from 22 μmol m⁻² h⁻¹ to 28 μmol m⁻² h⁻¹, approximately equivalent to that measured in the surrounding nonvegetated sediment. Although sediment denitrification is suppressed when the sediment surface is covered by a dense macroalgal mat, the denitrification zone may migrate up into the mat. In eutrophic estuaries with a large area of macroalgal cover, the physical structure and growth stage of algal mats may thus play an important role in the regulation of nitrogen removal by denitrification.     

An inflationary cosmological model in Dunn's scalar-tensor theory of gravitation

The gravitational field equations in Dunn's scalar-tensor theory of gravitation are generalized by including a cosmological constant. The resulting equations are solved for a Robertson-Walker line-element with flat three-space. The solution represents a cosmological model that develops into an inflationary era.     

Wide field monitoring of the X-ray sky using Rotation Modulation Collimators

Wide field monitoring is of particular interest in X-ray astronomy due to the strong time-variability of most X-ray sources. Not only does the time-profiles of the persistent sources contain characteristic signatures of the underlying physical systems, but, additionally, some of the most intriguing sources have long periods of quiesense in which they are almost undetectable as X-ray sources, interspersed with relatively brief periods of intense outbursts, where we have unique opportunities of studying dynamical effects, in, for instance, the evolution of accretion discs. Another question for which wide field monitors may provide key information, is the origin and nature of the cosmic gamma ray bursts.Rotation Modulation Collimators (RMC's) were originally introduced in X-ray astronomy to provide accurate source localizations over extended fields. This role has since been taken over by the grazing incidence telescope systems. The potential of the RMC's as wide field monitors have recently been demonstrated by the WATCH instruments on GRANAT and EURECA. It now appears likely, that for use on large, 3-axis stabilized spacecraft, a pinhole camera system may provide better sensitivity than an RMC-system of corresponding physical dimensions. But due to its simplicity, low data rate, and ability to work on spin stabilized (micro)satellites, the RMC wide field monitor may still have a role to play in the X-ray astronomy of the future.     

A New method for identifying possible causal relationships between CO2, total solar irradiance and global temperature change

The study compares two formulas for calculating the daily evapotranspiration ET₀ for a reference crop. The first formula was proposed by Allen et al. (AL), while the second one was proposed by Katerji and Perrier with the addition of the carbon dioxide (CO₂) effect on evapotranspiration (KP). The study analyses the impact of the calculation by the two formulas on the irrigation requirement (IR). Both formulas are based on the Penman-Monteith equation but adopt different approaches for parameterising the canopy resistance r _c . In the AL formula, r _c is assumed constant and not sensitive to climate change, whereas in the KP formula, r _c is first parameterised as a function of climatic variables, then ET₀ is corrected for the air CO₂ concentration. The two formulas were compared in two periods. The first period involves data from two sites in the Mediterranean region within a measured climate change period (1981–2006) when all the input climatic variables were measured. The second period (2070–2100) involves data from a future climate change period at one site when the input climatic variables were forecasted for two future climate scenarios (A2 and B2). The annual cumulated values of ET₀ calculated by the AL formula are systematically lower than those determined by the KP formula. The differences between the ET₀ estimation with the AL and KP formulas have a strong impact on the determination of the IR for the reference crop. In fact, for the two periods, the annual values of IR when ET₀ is calculated by the AL formula are systematically lower than those calculated by the KP formula. For the actual measured climate change period, this reduction varied from 26 to 28 %, while for the future climate change period, it varied based on the scenario from 16 % (A2) to 20 % (B2).     

Variability in inorganic and organic nitrogen uptake associated with riverine nutrient input in the Gulf of Riga, Baltic Sea

Concentrations and rates of uptake of dissolved organic nitrogen (DON, free amino acids, and urea) and inorganic nitrogen (DIN, nitrate, and ammonium) were measured along two transects in the Gulf of Riga, a sub-basin of the Baltic Sea, during May and July 1996. Concentrations of total dissolved nitrogen (TDN) were 23±3 μg-at N 1⁻¹ in the northern region (mouth) and 41±5 μg-at N 1⁻¹ in the southern region (head) of the Gulf. Rates of nitrogen uptake, determined with¹⁵N-labeled substrates, reflected differences in TDN concentration between the regions. In May, uptake of DIN+DON measured 0.17 and 0.43 μg-at N 1⁻¹ h⁻¹ in the northern and southern parts of the Gulf, respectively. In July, DIN+DON uptake measured 0.38 and 0.68 μg-at N 1⁻¹ h⁻¹ in the north and south, respectively. Most of the variability in total nitrogen flux between the northern and southern regions was due to heterogeneity of DON utilization. Uptake of urea and dissolved free amino acid were up to 6 and 3 times greater in the south compared to the north. As evidenced by size-fractionation, plankton size structure appeared to play a role in the uptake of DON. The community in the southern part was largely composed of cells <5 μm, while up to 67% of the community in the northern part was composed of cells >5 μm. Our results indicate that DON was a major source of nitrogen to phytoplankton, particularly in the southern part of the Gulf.     

New Isotopic data from davidites and sulfides in the bidjovagge gold-copper deposit,Finnmark, Northern Norway

Summary The Bidjovagge gold-copper deposit (69°17'N, 22°29'E) occurs in the north-south trending Proterozoic Kautokeino greenstone belt. The greenstone belt consists of shallow marine sediments and mafic to ultramafic volcanites, which are intruded by diabase sills and granitoids. The sequence can be correlated with the 1930 to 2200 Ma Kiruna greenstone belt. Svecokarelian intrusions took place in two phases between 1880–1890 Ma and between 1860-1870 Ma.The deposit occurs in the stratigraphically lower part of the belt along a north-south trending anticlinal structure on the border between greenschist and amphibolite grade metamorphism within the Baltic-Bothnia megashear zone. The deposit is related to shear zones and hosted by albitic felsite, a metasomatic rock formed along the contact with the diabase sills.The only ore minerals of economic significance are native gold and chalcopyrite. Other common ore minerals are pyrite and pyrrhotite. Marcasite, magnetite, ilmenite, hematite, tellurides, rutile, sphalerite, galena, davidite and pentlandite occur in accessory amounts. A total of 10 samples of davidite, 14 sulfides and 4 samples of albitic felsite were analyzed for their isotopic compositions. Our best estimate for the age of the davidites (1885 ± 18 Ma) is based on the upper intersection of a U/Pb concordia diagram and is identified as the time of davidite (and gold) mineralization. This age is confirmed by a Sm/Nd date of 1886 ± 88 Ma on the same material. At 1339 ± 8 Ma, a metamorphic or metasomatic event disturbed the U/Pb system in some of the davidite samples and completely reset the U/Pb ratios in the albitic felsites. The sulfide minerals, containing trace U and Pb, were strongly disturbed during galena formation at about 500–550 Ma but apparently were initially formed at the same time as, or shortly after, the davidite. Based on an estimated time of formation of galena during the development of the early Paleozoic peneplain in the area, the Pb/Pb data are interpreted as a two-stage isochron with an initial age of 1876 ± 15 Ma.The relation between the mineralization and shear zones, the correspondence between the age of the mineralization and the Svecokarelian orogeny and the initial ratios of the radiogenic isotopes of the ore forming solutions, all support a metamorphic/metasomatic origin for the gold and copper mineralization.

---

Neue Isotopendaten von Davidit und Sulfiden aus der Bidjovagge Gold-Lagerstätte, Finnmark, N.-Norwegen

Zusammenfassung Die Bidjovagge Gold-Kupfer-Lagerstätte (69°17'N, 22°29'E) liegt in dem N-S streichenden proterozoischen Kautokeino Grüngesteinsgürtel. Er besteht aus marinen Flachwassersedimenten, sowie mafischen und ultramafischen Vulkaniten, die von Diabas-Sills und von Granitoiden intrudiert werden. Diese Gesteinsabfolge kann mit dem 1930–2200 Ma alten Kiruna Grüngesteinsgürtel verglichen werden. Sveko-karelische Intrusionen erfolgten zweiphasig, zwischen 1880–1890 Ma und zwischen 1860 und 1870 Ma.Die Lagerstätte liegt in den stratigraphisch tieferen Anteilen diese Grüngesteinsgürtel entlang einer N-S streichenden Antiklinalstruktur im Übergangsbereich zwischen der Grünschiefer und Amphibolitfazies innerhalb der baltisch-bothnischen Mega-Scherzone.Sie ist an Scherzonen und an Albit-Felsite, einem metasomatisch am Kontakt mit den Diabas-Sills gebildeten Gestein, gebunden.Gediegen Gold und Kupferkies sind die wirtschaftlich wichtigsten Erzminerale. Pyrit und Magnetkies sind häufig. Markasit, Magnetit, Ilmenit, Hämatit, Telluride, Rutil, Zinkblende, Bleiglanz, Davidit und Pentlandit kommen akzessorisch vor. 10 Proben von Davidit, 14 Sulfide und 4 Proben von den Albit-Felsiten wurden auf ihre Isotopenzusammensetzung hin untersucht.Der obere Schnittpukt im U/Pb Concordiadiagramm ergab die beste Altersabschätzung für die Davidite (1885 ± 18 Ma). Dieses Alter wird als Zeitpunkt der Davidit (und Gold)-Mineralisation betrachtet. Ein Sm/Nd-Alter von 1886 ± 88 Ma bestätigt dieses. Ein metamorphes oder metasomatisches Ereignis, datiert mit 1339 ± 8 Ma störte das U/Pb System in einigen Daviditproben und führte zu einer vollständigen Reequilibrierung der U/Pb Verhältnisse in den Albit-Felsiten. Die Isotopenverhältnisse der Sulfide, die Spuren von U und Pb enthalten, wurden während der Bildung von Bleiglanz um etwa 500–550 Ma ebenfalls massiv gestört; die Sulfide wurden aber offensichtlich gleichzeitig, oder aber kurz nach Davidit, gebildet. Basierend auf einer Abschätzung des Bleiglanz-Bildungsalters im Zuge der Entwicklung einer altpaläozoischen Peneplain in diesem Gebiet, werden die Pb/Pb Daten als Zweistufen-Isochrone mit einem Ursprungsalter von 1876 ± 15 Ma interpretiert.Die Bindung der Mineralisation an Scherzonen, die Übereinstimmung des Mineralisationsalters mit dem der Sveco-karelischen Orogenese und die Isotopeninitiale der erzbildenden Lösungen belegen einen metamorph/metasomatischen Ursprung der Gold-und Kupfervererzungen.

---

---

With 9 Figures     

Radioelement (U,Th, Rn) concentrations in Norwegian bedrock groundwaters

Twenty-eight samples of groundwater from bedrock boreholes in three distinct Norwegian geological provinces have been taken and analyzed for content of Rn, U, and Th, together with a wide variety of minor and major species. Median values of 290 Bq/1, 7.6 g/1, and 0.02 g/1 were obtained for Rn, U, and Th respectively, while maximum values were 8500 Bq/1, 170 g/1, and 2.2 g/1. Commonly suggested drinking water limits range from 8 to 1000 Bq/1 for radon and 14 to 160 g/1 for uranium. Radioelement content was closely related to lithology, the lowest concentrations being derived from the largely Caledonian rocks of the Trøndelag area, and the highest from the Precambrian Iddefjord Granite of southeast Norway (11 boreholes) where median values of 2500 Bq/1, 15 g/1 and 0.38 g/1, respectively, were obtained. The Iddefjord Granite is not believed to be unique in Norway in yielding high dissolved radionuclide contents in groundwaters, and several other granitic aquifers warrant further investigation in this respect.