Laboratory studies about the interaction of ammonia with ice crystals at temperatures between 0 and −20°C

In laboratory experiments the interactions of ammonia with ice crystals were studied within the temperature range between 0 and −20°C. In a first series of experiments dendritic ice crystals were grown from water vapor in presence of ammonia gas in various concentrations between 4 and 400 ppbv. In a second series of experiments pure ice crystals were exposed to a humidified ammonia–air mixture inside a horizontal flow tube. The influence of temperature, ammonia gas concentration (0.6, 1.5, and 10 ppmv), exposure time, and the presence of impurities such as sulfate on the ammonia uptake by the ice surface was investigated by determining the ammonium content in the melt water of the ice crystals by ion chromatography. During the growth of ice crystals significant amounts of ammonia (around 200 μg/l) were taken up even at small gas concentrations. In contrast, even at high gas concentrations the uptake of ammonia by non-growing ice crystals was lower by approximately one order of magnitude. The presence of sulfate on the ice surface affected an enhanced uptake of ammonia by a factor of 5–10. A model is presented which describes the uptake of ammonia by ice considering the chemical processes occurring in the ice surface layer and simultaneous diffusion of ammonia into bulk ice. Even the increased uptake of ammonia by growing ice is rather small compared to the uptake by water droplets; thus, the major process for scavenging of ammonia from the atmosphere via the ice phase might not be the direct uptake by ice crystals but the riming involving super-cooled droplets containing ammonia.     

The Scaling Behaviour of a Turbulent Kinetic Energy Closure Model for Stably Stratified Conditions

We investigate the scaling behaviour of a turbulent kinetic energy (TKE) closure model for stably stratified conditions. The mixing length scale for stable stratification is proportional to the ratio of the square root of the TKE and the local Brunt–Väisälä frequency, which is a commonly applied formulation. We analyze the scaling behaviour of our model in terms of traditional Monin–Obukov Similarity Theory and local scaling. From the model equations, we derive expressions for the stable limit behaviour of the flux–gradient relations and other scaling quantities. It turns out that the scaling behaviour depends on only a few model parameters and that the results obey local scaling theory. The analytical findings are illustrated with model simulations for the second GABLS intercomparison study. We also investigate solutions for the case in which an empirical correction function is used to express the eddy diffusivity for momentum as a function of the Richardson number (i.e. an increasing turbulent Prandtl number with stability). In this case, it seems that for certain parameter combinations the model cannot generate a steady-state solution. At the same time, its scaling behaviour becomes unrealistic. This shows that the inclusion of empirical correction functions may have large and undesired consequences for the model behaviour.     

Astrophysical smooth particle hydrodynamics

The paper presents a detailed review of the smooth particle hydrodynamics (SPH) method with particular focus on its astrophysical applications. We start by introducing the basic ideas and concepts and thereby outline all ingredients that are necessary for a practical implementation of the method in a working SPH code. Much of SPH’s success relies on its excellent conservation properties and therefore the numerical conservation of physical invariants receives much attention throughout this review. The self-consistent derivation of the SPH equations from the Lagrangian of an ideal fluid is the common theme of the remainder of the text. We derive a modern, Newtonian SPH formulation from the Lagrangian of an ideal fluid. It accounts for changes of the local resolution lengths which result in corrective, so-called “grad-h-terms”. We extend this strategy to special relativity for which we derive the corresponding grad-h equation set. The variational approach is further applied to a general-relativistic fluid evolving in a fixed, curved background space-time. Particular care is taken to explicitly derive all relevant equations in a coherent way.     

Submarine landslides at the eastern Sunda margin: observations and tsunami impact assessment

Our analysis of new bathymetric data reveals six submarine landslides at the eastern Sunda margin between central Java and Sumba Island, Indonesia. Their volumes range between 1 km3 in the Java fore-arc basin up to 20 km3 at the trench off Sumba and Sumbawa. We estimate the potential hazard of each event by modeling the corresponding tsunami and its run-up on nearby coasts. Four slides are situated remarkably close to the epicenter of the 1977 tsunamigenic Sumba M _w  = 8.3 earthquake. However, comparison of documented tsunami run-up heights and arrival times with our modeling results neither allows us to confirm nor can we falsify the hypothesis that the earthquake triggered these submarine landslides.     

Predicting riverine dissolved silica fluxes to coastal zones from a hyperactive region and analysis of their first-order controls

Silicate weathering and resulting transport of dissolved matter influence the global carbon cycle in two ways. First by the uptake of atmospheric/soil CO₂ and second by providing the oceanic ecosystems via the fluvial systems with the nutrient dissolved silica (DSi). Previous work suggests that regions dominated by volcanics are hyperactive or even “hot spots” concerning DSi-mobilization. Here, we present a new approach for predicting DSi-fluxes to coastal zones, emphasizing “first-order” controlling factors (lithology, runoff, relief, land cover and temperature). This approach is applied to the Japanese Archipelago, a region characterized by a high percentage of volcanics (29.1% of surface area). The presented DSi-flux model is based on data of 516 catchments, covering approximately 56.7% of the area of the Japanese Archipelago. The spatial distribution of lithology—one of the most important first order controls—is taken from a new high resolution map of Japan. Results show that the Japanese Archipelago is a hyperactive region with a DSi-yield 6.6 times higher than the world average of 3.3 t SiO₂ km⁻² a⁻¹, but with large regional variations. Approximately 10% of its area exceeds 10 times the world average DSi-yield. Slope constitutes another important controlling factor on DSi-fluxes besides lithology and runoff, and can exceed the influence of runoff on DSi-yields. Even though the monitored area on the Japanese Archipelago stretches from about 31° to 46°N, temperature is not identified as a significant first-order model variable. This may be due to the fact that slope, runoff and lithology are correlated with temperature due to regional settings of the Archipelago, and temperature information is substituted to a certain extent by these factors. Land cover data also do not improve the prediction model. This may partly be attributed to misinterpreted land cover information from satellite images. Implications of results for Earth System and global carbon cycle modeling are discussed.     

In situ Sr/Sr investigation of igneous apatites and carbonates using laser-ablation MC-ICP-MS

In situ Sr isotopic compositions of coexisting apatite and carbonate for carbonatites from the Sarfartoq alkaline complex, Greenland, have been determined by laser-ablation multicollector inductively coupled plasma mass spectrometry. This study is the first to examine the extent of Sr isotopic homogeneity among coexisting igneous minerals containing high Sr (>3000 ppm) and low Rb (?1 ppm) contents within a single ∼50-μm-thick thin-section mount. This technique is capable of producing measured ⁸⁷Sr/⁸⁶Sr values with analytical precision (∼0.005%, 2σ) approaching those obtained by thermal ionization mass spectrometry but in a much shorter interval of time (100 s vs. >1 h, respectively). The combined total analyses (n = 107) of apatite and carbonate yield ⁸⁷Sr/⁸⁶Sr compositions ranging from ∼0.7025 to ∼0.7031. This relatively large variation in Sr isotopic compositions (∼0.0006) is ∼1 order of magnitude larger than the estimated external reproducibility (∼0.00005, 2σ) of the method. The large range in ⁸⁷Sr/⁸⁶Sr values suggests that apatite and carbonate precipitated predominantly under nonequilibrium conditions. The isotopic variations observed within individual hand specimens may therefore reflect larger (regional) scale open-system processes, possibly involving mixing of carbonatitic melts derived from distinct mantle sources or from a common isotopically heterogeneous mantle.     

A Straightforward Random Walk Model for Fast Push‐Pull Tracer Test Evaluation

In this article, we present a straightforward random walk model for fast evaluation of push‐pull tracer tests. By developing an adaptive algorithm, we overcome the problem of manually defining how many particles have to be used to simulate the transport problem. Beside this, we validate the random walk model by evaluating a push‐pull tracer test with drift phase and confirm the results with MT3DMS. The random walk model took less than 1% of computational time of MT3DMS, thus allowing a remarkable faster evaluation of push‐pull tracer tests.     

Comparison of Daily GRACE Gravity Field and Numerical Water Storage Models for De-aliasing of Satellite Gravimetry Observations

Reducing aliasing effects of insufficiently modelled high-frequent, non-tidal mass variations of the atmosphere, the oceans and the hydrosphere in gravity field models derived from the Gravity Recovery and Climate Experiment (GRACE) satellite mission is the topic of this study. The signal content of the daily GRACE gravity field model series (ITG-Kalman) is compared to high-frequency bottom pressure variability and terrestrially stored water variations obtained from recent numerical simulations from an ocean circulation model (OMCT) and two hydrological models (WaterGAP Global Hydrology Model, Land Surface Discharge Model). Our results show that daily estimates of ocean bottom pressure from the most recent OMCT simulations and the daily ITG-Kalman solutions are able to explain up to 40 % of extra-tropical sea-level variability in the Southern Ocean. In contrast to this, the daily ITG-Kalman series and simulated continental total water storage variability largely disagree at periods below 30 days. Therefore, as long as no adequate hydrological model will become available, the daily ITG-Kalman series can be regarded as a good initial proxy for high-frequency mass variations at a global scale. As a second result of this study, based on monthly solutions as well as daily observation residuals, it is shown that applying this GRACE-derived de-aliasing model supports the determination of the time-variable gravity field from GRACE data and the subsequent geophysical interpretation. This leads us to the recommendation that future satellite concepts for determining mass variations in the Earth system should be capable of observing higher frequeny signals with sufficient spatial resolution.