A study on the reference frame consistency in recent Earth gravitational models

All gravity field functionals obtained from an Earth gravitational model (EGM) depend on the underlying terrestrial reference frame (TRF), with respect to which the EGM’s spherical harmonic coefficients refer to. In order to maintain a coherent framework for the comparison of current and future EGMs, it is thus important to investigate the consistency of their inherent TRFs, especially when their use is intended for high precision studies. Following the methodology described in an earlier paper by Kleusberg (1980), the similarity transformation parameters between the associated reference frames for several EGMs (including the most recent CHAMP/GRACE models at the time of writing this paper) are estimated in the present study. Specifically, the differences between the spherical harmonic coefficients for various pairs of EGMs are parameterized through a 3D-similarity spatial transformation model that relates their underlying TRFs. From the least-squares adjustment of such a parametric model, the origin, orientation and scale stability between the EGMs’ reference frames can be identified by estimating their corresponding translation, rotation and scale factor parameters. Various aspects of the estimation procedure and its results are highlighted in the paper, including data weighting schemes, the sensitivity of the results with respect to the selected harmonic spectral band, the correlation structure and precision level of the estimated transformation parameters, the effect of the estimated differences of the EGMs’ reference frames on their height anomaly signal, and the overall feasibility of Kleusberg’s formulae for the assessment of TRF inconsistencies among global geopotential models.     

Using gravity and topography-implied anomalies to assess data requirements for precise geoid computation

Many regions around the world require improved gravimetric data bases to support very accurate geoid modeling for the modernization of height systems using GPS. We present a simple yet effective method to assess gravity data requirements, particularly the necessary resolution, for a desired precision in geoid computation. The approach is based on simulating high-resolution gravimetry using a topography-correlated model that is adjusted to be consistent with an existing network of gravity data. Analysis of these adjusted, simulated data through Stokes’s integral indicates where existing gravity data must be supplemented by new surveys in order to achieve an acceptable level of omission error in the geoid undulation. The simulated model can equally be used to analyze commission error, as well as model error and data inconsistencies to a limited extent. The proposed method is applied to South Korea and shows clearly where existing gravity data are too scarce for precise geoid computation.     

A probabilistic methodology to estimate future coastal flood risk due to sea level rise

In this paper we present a methodology to estimate the probability of future coastal flooding given uncertainty over possible sea level rise. We take as an example the range of sea level rise magnitudes for 2100 contained in the IPCC Third Assessment Report [Church, J.A., Gregory, J.M., Huybrechts, P., Kuhn, M., Lambeck, K., Nhuan, M.T., Qin, D., Woodworth, P.L., Anisimov, O.A., Bryan, F.O., Cazenave, A., Dixon, K.W., Fitzharris, B.B., Flato, G.M., Ganopolski, A., Gornitz, V., Lowe, J.A., Noda, A., Oberhuber, J.M., O'Farrell, S.P., Ohmura, A., Oppenheimer, M., Peltier, W.R., Raper, S.C.B., Ritz, C., Russell, G.L., Schlosser, E., Shum, C.K., Stocker, T.F., Stouffer, R.J., van de Wal, R.S.W., Voss, R., Wiebe, E.C., Wild, M., Wingham, D.J. and Zwally, H.J., 2001. Changes in sea level. In Houghton, J.T. et al. (eds), Climate Change 2001: The Scientific Basis. Contribution of Working Group I to the Third Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge University Press, Cambridge, United Kingdom, 881pp.] and infer a plausible probability distribution for this range. We then use a Monte Carlo procedure to sample from this distribution and use the resulting values as an additional boundary forcing for a two-dimensional model of coastal inundation used to simulate a 1 in 200 year extreme water level event. This yields an ensemble of simulations for an event of this magnitude occurring in 2100, where each member represents a different possible scenario of sea level rise by this time. We then develop a methodology to approximate the probability of flooding in each model grid cell over the ensemble and by combining these hazards maps with maps of land use values (consequence) we are able to estimate spatial contributions to flood risk that can aid planning and investment decisions. The method is then applied to a 32 km section of the UK coast in Somerset, South-West England and used to estimate the monetary losses and risk due a 1 in 200 year recurrence interval event under: (a) current conditions; (b) with the IPCC's most plausible value for sea level rise by 2100 (0.48 m) and (c) using the above methodology to fully account for uncertainty over possible sea level rise. The analysis shows that undertaking a risk assessment using the most plausible sea level rise value may significantly underestimate monetary losses as it fails to account for the impact of low probability, high consequence events. The developed method provides an objective basis for decisions regarding future defence spending and can be easily extended to consider other sources of uncertainty such as changing event frequency–magnitude distribution, changing storm surge conditions or model structural uncertainty, either singly or in combination as joint probabilities.     

Features of coastal upwelling regions that determine net air-sea CO<Subscript>2</Subscript> flux

The influence of the coastal ocean on global net annual air-sea CO₂ fluxes remains uncertain. However, it is well known that air-sea pCO₂ disequilibria can be large (ocean pCO₂ ranging from ∼400 μatm above atmospheric saturation to ∼250 μatm below) in eastern boundary currents, and it has been hypothesized that these regions may be an appreciable net carbon sink. In addition it has been shown that the high productivity in these regions (responsible for the exceptionally low surface pCO₂) can cause nutrients and inorganic carbon to become more concentrated in the lower layer of the water column over the shelf relative to adjacent open ocean waters of the same density. This paper explores the potential role of the winter season in determining the net annual CO₂ flux in temperate zone eastern boundary currents, using the results from a box model. The model is parameterized and forced to represent the northernmost part of the upwelling region on the North American Pacific coast. Model results are compared to the few summer data that exist in that region. The model is also used to determine the effect that upwelling and downwelling strength have on the net annual CO₂ flux. Results show that downwelling may play an important role in limiting the amount of CO₂ outgassing that occurs during winter. Finally data from three distinct regions on the Pacific coast are compared to highlight the importance of upwelling and downwelling strength in determining carbon fluxes in eastern boundary currents and to suggest that other features, such as shelf width, are likely to be important.     

Seasonal changes in community composition and trophic structure of fish populations of five salt marshes along the Essex coastline, United Kingdom

European intertidal salt marshes are important nursery sites for juvenile fish and crustaceans. Due to the increasing threat of habitat loss, the seasonal changes of salt marsh fish communities need to be understood in order to appreciate the ecological and economic importance of the saltmarsh habitat. This study was the first in Great Britain to investigate the seasonal changes of salt marsh fish communities and the variation in community structure between closely located marsh habitats. Between February 2007 and March 2008, five marshes on three estuaries of the Essex coastline were sampled using flume nets to block off intertidal creeks at high tide. Fourteen fish species were caught. The community overall was dominated by three species that made up 91.6% of the total catch: the common goby Pomatoschistus microps (46.2% of the total catch), juvenile herring Clupea harengus (24.3%), and juvenile and larval sea bass Dicentrarchus labrax (21.2%). Cluster analysis demonstrated clear seasonal patterns, with some community structures unique to specific marshes or estuaries. The marsh fish community shifts from a highly diverse community during spring, to a community dominated by D. labrax and P. microps in autumn, and low diversity during winter months. Gravimetric stomach content analysis of fish community identified three main trophic guilds; macroinvertivores, planktivores and omnivores. The macroinvertivore feeding guild contained D. labrax and P. microps, the two most frequently occurring species. This investigation demonstrates the importance of British salt marshes as nursery habitats for commercial fish species.     

Stable isotope ratios in Cape gannets around the southern coasts of Africa reveal penetration of biogeographic patterns in oceanic signatures

The southern coasts of Africa are influenced by two major oceanic currents, leading to biogeographic patterns in inshore and offshore species assemblages, and in the stable isotope signatures of suspended particulate matter and filter-feeding mussels. We used the stable isotope ratios of carbon (¹³C/¹²C) and nitrogen (¹⁵N/¹⁴N) from the blood and feathers of adult and chick Cape gannets (Morus capensis) to investigate whether the geographic differences observed at the lower levels in the marine communities are deep penetrating effects that reach top predators. Additionally, we evaluated whether trophic segregation occurs between adult and reared chick gannets, and whether a shift to wintering habitat occurs in adults. The study was conducted during the 2006 breeding season on Bird Island in the Agulhas system, and on Malgas and Ichaboe Islands, in the south and north Benguela respectively. Our results showed significant differences in the isotope ratios of members of different colonies, but no intra-colony differences between tissues or age groups. These results indicate that there is neither age-related nor temporal segregation in the diet of members of the same colony. Feather isotopic values suggest that adults remain all year round in the same habitats, and do not undertake long migration after reproduction. Since all gannets tend to target similar prey, we attributed among-colony differences in isotope signatures mostly to the oceanic conditions experienced by the main prey of birds rather than substantial differences in diet composition. Overall, isotopic signatures segregate the two current systems, with depleted carbon values in the Agulhas and enriched nitrogen values in the upwelled waters of the Benguela. Within the Benguela birds from Ichaboe in the north had higher δ¹⁵N values than those from Malgas in the south, which we attributed to differences in the functioning of the upwelling cells in the vicinity of the two colonies. Finally, slight variation in the proportion of main prey and discards from fisheries may contribute to the variation in the stable-isotope signatures between colonies in the Benguela.     

Shapes of chondrules determined from the petrofabric of the CR2 chondrite NWA 801

The approximately spherical shapes of chondrules has long been attributed to surface tension acting on ~1 mm melt droplets that formed and cooled in the microgravity field of the solar nebula. However, chondrule shapes commonly depart significantly from spherical. In this study, 109 chondrules in a sample of CR2 chondrite NWA 801 were imaged by X-ray computed tomography and best-fitted to ellipsoids. The analysis confirms that many chondrules are indeed not spherical, and also that the chondrules’ collective shape fabric records a definite 13% compaction in the host meteorite. Dehydration of phyllosilicates within chondrules may account for that strain. However, retro-deforming all chondrules shows that a large majority were already far from spherical prior to accretion. Possible models for these initial shapes include prior deformation of individual chondrules in earlier hosts, and, as suggested by previous authors, rotation of chondrules as they were solidifying, and/or “streaming” of molten chondrules by their differential velocities with their gaseous hosts after melting. More in situ 3-D work such as this study on a variety of unequilibrated chondrites, combined with detailed structural petrography, should help further constrain these models and refine our understanding of chondrite formation.     

Experimental impact cratering: A summary of the major results of the MEMIN research unit

This paper reviews major findings of the Multidisciplinary Experimental and Modeling Impact Crater Research Network (MEMIN). MEMIN is a consortium, funded from 2009 till 2017 by the German Research Foundation, and is aimed at investigating impact cratering processes by experimental and modeling approaches. The vision of this network has been to comprehensively quantify impact processes by conducting a strictly controlled experimental campaign at the laboratory scale, together with a multidisciplinary analytical approach. Central to MEMIN has been the use of powerful two-stage light-gas accelerators capable of producing impact craters in the decimeter size range in solid rocks that allowed detailed spatial analyses of petrophysical, structural, and geochemical changes in target rocks and ejecta. In addition, explosive setups, membrane-driven diamond anvil cells, as well as laser irradiation and split Hopkinson pressure bar technologies have been used to study the response of minerals and rocks to shock and dynamic loading as well as high-temperature conditions. We used Seeberger sandstone, Taunus quartzite, Carrara marble, and Weibern tuff as major target rock types. In concert with the experiments we conducted mesoscale numerical simulations of shock wave propagation in heterogeneous rocks resolving the complex response of grains and pores to compressive, shear, and tensile loading and macroscale modeling of crater formation and fracturing. Major results comprise (1) projectile–target interaction, (2) various aspects of shock metamorphism with special focus on low shock pressures and effects of target porosity and water saturation, (3) crater morphologies and cratering efficiencies in various nonporous and porous lithologies, (4) in situ target damage, (5) ejecta dynamics, and (6) geophysical survey of experimental craters.