A comparison of shelf observation platforms for assimilation in an eddy-resolving ocean model

An assessment of the likely benefits of assimilating in situ temperature (T) and salinity (S) observations from repeat glider transects and surface velocity observations from high-frequency radar arrays into an eddy-resolving ocean model is presented. The deployment of new shelf observation platforms around Australia is being undertaken through the Australian Integrated Marine Observing System program. In this study, various options for an observing system along the coast of New South Wales, Australia, are assessed for their benefits to an ocean forecast and reanalysis system. The forecast system considered here uses ensemble optimal interpolation (EnOI) for data assimilation. Using error estimates from the EnOI scheme, estimates of the theoretical analysis errors are calculated for different observing systems that include a range of remotely sensed and in situ observations. The results demonstrate that if HF radar observations are assimilated along with the standard components of the global ocean observing system, the analysis errors are likely to reduce by as much as 80% for velocity and 60% for T, S and sea-level in the vicinity of the observations. Owing to the relatively short along-shore decorrelation length-scales for T and S near the shelf, the glider observations are likely to provide the forecast system with a more modest gain.     

Effect of snow cover on pan-Arctic permafrost thermal regimes

This study quantitatively evaluated how insulation by snow depth (SND) affected the soil thermal regime and permafrost degradation in the pan-Arctic area, and more generally defined the characteristics of soil temperature (T_SOIL) and SND from 1901 to 2009. This was achieved through experiments performed with the land surface model CHANGE to assess sensitivity to winter precipitation as well as air temperature. Simulated T_SOIL, active layer thickness (ALT), SND, and snow density were generally comparable with in situ or satellite observations at large scales and over long periods. Northernmost regions had snow that remained relatively stable and in a thicker state during the past four decades, generating greater increases in T_SOIL. Changes in snow cover have led to changes in the thermal state of the underlying soil, which is strongly dependent on both the magnitude and the timing of changes in snowfall. Simulations of the period 2001–2009 revealed significant differences in the extent of near-surface permafrost, reflecting differences in the model’s treatment of meteorology and the soil bottom boundary. Permafrost loss was greater when SND increased in autumn rather than in winter, due to insulation of the soil resulting from early cooling. Simulations revealed that T_SOIL tended to increase over most of the pan-Arctic from 1901 to 2009, and that this increase was significant in northern regions, especially in northeastern Siberia where SND is responsible for 50 % or more of the changes in T_SOIL at a depth of 3.6 m. In the same region, ALT also increased at a rate of approximately 2.3 cm per decade. The most sensitive response of ALT to changes in SND appeared in the southern boundary regions of permafrost, in contrast to permafrost temperatures within the 60°N–80°N region, which were more sensitive to changes in snow cover. Finally, our model suggests that snow cover contributes to the warming of permafrost in northern regions and could play a more important role under conditions of future Arctic warming.     

Transient hysteresis of near-surface permafrost response to external forcing

Estimates of changes in near-surface permafrost (NSP) area S _p relative to change in globally averaged surface air temperature T _g are made by using the global climate model developed at the A.M. Obukhov Institute of Atmospheric Physics RAS (IAP RAS CM). For ensemble of runs forced by scenarios constructed as return-to-preindustrial continuations of the RCP (Representative Concentration Pathways) scenarios family, a possibility of transient hysteresis in dependence of S _p versus T _g is exhibited: in some temperature range which depends on imposed scenario of external forcing, NSP area is larger, at the same global mean surface air temperature, in a warming climate than in a cooling climate. This hysteresis is visible more clearly for scenarios with higher concentration of greenhouse gases in the atmosphere in comparison to those in which this concentration is lower. Hysteresis details are not sensitive to the type of the prescribed continuation path which is used to return the climate to the preindustrial state. The multiple-valued dependence of S _p on T _g arises due to dependence of soil state in the regions of extra-tropical wetlands and near the contemporary NSP boundaries on sign of external climatic forcing. To study the dependence of permafrost hysteresis on amplitude and temporal scale of external forcing, additional model runs are performed. These runs are forced by idealised scenarios of atmospheric CO₂ content varying, depending on run, with periods from 100 to 1,000 year and with different amplitudes. It is shown that the above-mentioned hysteresis is related to the impact of phase transitions of soil water on apparent inertia of the system as well as to the impact of soil state on atmospheric hydrological cycle and radiation transfer in the atmosphere.     

Spherical integral transforms of second-order gravitational tensor components onto third-order gravitational tensor components

New spherical integral formulas between components of the second- and third-order gravitational tensors are formulated in this article. First, we review the nomenclature and basic properties of the second- and third-order gravitational tensors. Initial points of mathematical derivations, i.e., the second- and third-order differential operators defined in the spherical local North-oriented reference frame and the analytical solutions of the gradiometric boundary-value problem, are also summarized. Secondly, we apply the third-order differential operators to the analytical solutions of the gradiometric boundary-value problem which gives 30 new integral formulas transforming (1) vertical-vertical, (2) vertical-horizontal and (3) horizontal-horizontal second-order gravitational tensor components onto their third-order counterparts. Using spherical polar coordinates related sub-integral kernels can efficiently be decomposed into azimuthal and isotropic parts. Both spectral and closed forms of the isotropic kernels are provided and their limits are investigated. Thirdly, numerical experiments are performed to test the consistency of the new integral transforms and to investigate properties of the sub-integral kernels. The new mathematical apparatus is valid for any harmonic potential field and may be exploited, e.g., when gravitational/magnetic second- and third-order tensor components become available in the future. The new integral formulas also extend the well-known Meissl diagram and enrich the theoretical apparatus of geodesy.     

Integral transformations of gradiometric data onto a GRACE type of observable

Integral transformations of gravitational gradients onto a Gravity Recovery And Climate Experiment (GRACE) type of observable are derived in this article. The gravitational gradients represent components of the gravitational tensor in the local north-oriented frame. The GRACE type of observable corresponds to a difference between two gravitational vectors as projected onto the line of sight between the two GRACE satellites. In total, three integral transformations relating vertical–vertical, vertical–horizontal and horizontal–horizontal gravitational gradients with the GRACE type of observable are provided. Spectral and closed forms of corresponding isotropic kernels are derived for each transformation. Special cases show that the integral transformations are general and relate gravitational gradients to many other quantities of the gravitational field, such as the gravitational vector, and its radial and tangential components. Correctness of the mathematical derivations is validated in a closed-loop simulation using synthetic data.     

The GEOMON network of Czech catchments provides long-term insights into altered forest biogeochemistry: From acid atmospheric deposition to climate change

In 1994, a network of small catchments (GEOMON) was established in the Czech Republic to determine input–output element fluxes in semi-natural forest ecosystems recovering from anthropogenic acidification. The network consists from 16 catchments and the primary observations of elements fluxes were complemented by monitoring of biomass stock, element pools in soil and vegetation, and the main water balance components. Over last three decades, reductions of SO₂, NO_x and NH₃ emissions were followed by sulphur (S) and nitrogen (N) deposition reductions of 75% and 30%, respectively. Steeper declines of strong acid anion concentrations compared to cations (Ca, Mg, Na, K, NH₄) in precipitation resulted in precipitation pH increase from 4.5 to 5.2 in bulk precipitation and from 4.0 to 5.1 in spruce throughfall. Stream chemistry responded to changes in deposition: S leaching declined. However at majority of catchments soils acted as a net source of S to runoff, delaying recovery. Stream pH increased at acidic streams (pH < 6) and aluminium concentration decreased. Stream nitrate (NO₃) concentration declined by 60%, considerably more than N deposition. Stream NO₃ concentration was tightly positively related to stream total dissolved nitrogen to total phosphorus (P) ratio, suggesting the role of P availability on N retention. Trends in dissolved organic carbon fluxes responded to both acidification recovery and to runoff temporal variation. An exceptional drought occurred between 2014 and 2019. Over this recent period, streamflow decreased by ≈ 40% on average compared to 1990s, due to the increases of soil evaporation and vegetation transpiration by ≈ 30% and declines in precipitation by ≈ 15% on average across the elevational gradient. Sharp decreases of stream runoff at catchments <650 m a.s.l. corresponded to areas of recent forest decline caused by bark beetle infestation on drought stressed spruce forests. Understanding of the interactions among legacies of acidification and eutrophication, drought effects on the water cycle and forest disturbance dynamics is requisite for effective management of forested ecosystems under anthropogenic influence.     

Optimal Model for Geoid Determination from Airborne Gravity

Two different approaches for transformation of airborne gravity disturbances, derived from gravity observations at low-elevation flying platforms, into geoidal undulations are formulated, tested and discussed in this contribution. Their mathematical models are based on Green's integral equations. They are in these two approaches defined at two different levels and also applied in a mutually reversed order. While one of these approaches corresponds to the classical method commonly applied in processing of ground gravity data, the other approach represents a new method for processing of gravity data in geoid determination that is unique to airborne gravimetry. Although theoretically equivalent in the continuous sense, both approaches are tested numerically for possible numerical advantages, especially due to the inverse of discretized Fredholm's integral equation of the first kind applied on different data. High-frequency synthetic gravity data burdened by the 2-mGal random noise, that are expected from current airborne gravity systems, are used for numerical testing. The results show that both approaches can deliver for the given data a comparable cm-level accuracy of the geoidal undulations. The new approach has, however, significantly higher computational efficiency. It would be thus recommended for real life geoid computations. Additional errors related to regularization of gravity data and the geoid, and to accuracy of the reference field, that would further deteriorate the quality of estimated geoidal undulations, are not considered in this study.     

On the Integrability of Limited-Area Numerical Weather Prediction Model ALADIN over Extended Time Periods

This note summarizes results of the first integration of regional numerical weather prediction model ALADIN in a climate mode. The ALADIN model, developed in an international cooperation led by Météo France, is operationally used for weather prediction. The grid step of the model is 12 km; the integration domain covers a major part of Europe. A one-month-long run has been performed with this model on observed boundary conditions (represented by assimilations by the global model ARPEGE). It is demonstrated that no excessive error is generated and accumulated in the model during the integration; hence the model is integrable for extended time periods and may serve a basis for a development towards a regional climate model.     

Boundary conditions in the hydromagnetic dynamo problem

au am nu¶rt;, nuau ¶rt;u¶rt;aum ¶rt;ua, a u ¶rt; ma aum u¶rt;uu aam ma¶rt;am . a a au ¶rt; uuauu u nm nmu a u m muna, mm m¶rt; ma a anma n u. uunua muau u m¶rt; nua [4, 5]. n¶rt; am nuam m¶rt; u, u¶rt;u u m¶rt;a a u a nm nu nu¶rt;um au m u m¶rt;.