Preprocessing of gravity gradients at the GOCE high-level processing facility

One of the products derived from the gravity field and steady-state ocean circulation explorer (GOCE) observations are the gravity gradients. These gravity gradients are provided in the gradiometer reference frame (GRF) and are calibrated in-flight using satellite shaking and star sensor data. To use these gravity gradients for application in Earth scienes and gravity field analysis, additional preprocessing needs to be done, including corrections for temporal gravity field signals to isolate the static gravity field part, screening for outliers, calibration by comparison with existing external gravity field information and error assessment. The temporal gravity gradient corrections consist of tidal and nontidal corrections. These are all generally below the gravity gradient error level, which is predicted to show a 1/f behaviour for low frequencies. In the outlier detection, the 1/f error is compensated for by subtracting a local median from the data, while the data error is assessed using the median absolute deviation. The local median acts as a high-pass filter and it is robust as is the median absolute deviation. Three different methods have been implemented for the calibration of the gravity gradients. All three methods use a high-pass filter to compensate for the 1/f gravity gradient error. The baseline method uses state-of-the-art global gravity field models and the most accurate results are obtained if star sensor misalignments are estimated along with the calibration parameters. A second calibration method uses GOCE GPS data to estimate a low-degree gravity field model as well as gravity gradient scale factors. Both methods allow to estimate gravity gradient scale factors down to the 10⁻³ level. The third calibration method uses high accurate terrestrial gravity data in selected regions to validate the gravity gradient scale factors, focussing on the measurement band. Gravity gradient scale factors may be estimated down to the 10⁻² level with this method.     

Preindustrial, historical, and fertilization simulations using a global ocean carbon model with new parameterizations of iron limitation, calcification, and N2 fixation

The Canadian Model of Ocean Carbon (CMOC) has been developed as part of a global coupled climate carbon model. In a stand-alone integration to preindustrial equilibrium, the model ecosystem and global ocean carbon cycle are in general agreement with estimates based on observations. CMOC reproduces global mean estimates and spatial distributions of various indicators of the strength of the biological pump; the spatial distribution of the air-sea exchange of CO₂ is consistent with present-day estimates. Agreement with the observed distribution of alkalinity is good, consistent with recent estimates of the mean rain ratio that are lower than historic estimates, and with calcification occurring primarily in the lower latitudes. With anthropogenic emissions and climate forcing from a 1850-2000 climate model simulation, anthropogenic CO₂ accumulates at a similar rate and with a similar spatial distribution as estimated from observations. A hypothetical scenario for complete elimination of iron limitation generates maximal rates of uptake of atmospheric CO₂ of less than 1 PgC y⁻¹, or about 11% of 2004 industrial emissions. Even a ‘perfect’ future of sustained fertilization would have a minor impact on atmospheric CO₂ growth. In the long term, the onset of fertilization causes the ocean to take up an additional 77 PgC after several thousand years, compared with about 84 PgC thought to have occurred during the transition into the last glacial maximum due to iron fertilization associated with increased dust deposition.     

Development and validation of a finite element morphological model for shallow water basins

Recognising the importance of understanding sediment dynamics to evaluate the status of a coastal lagoon environment, this work has been focused on the investigation of the hydrodynamic and sediment transport processes occurring in such basins. In order to describe the lagoon system, a modelling approach combining hydrodynamics, waves and sediment dynamics has been developed. The framework of the numerical model consists of a finite element hydrodynamic model, a third generation finite element spectral wave model and a sediment transport and morphodynamic model for both cohesive and non-cohesive sediments. The model adopts the finite element technique for spatial integration, which has the advantage to describe more accurately complicated bathymetry and irregular boundaries for shallow water areas. The developed model has been applied to test cases and to a very shallow tidal lagoon, the Venice Lagoon, Italy. Numerical results show good agreement with water level, waves and turbidity measurements collected in several monitoring stations inside the Lagoon of Venice. Such a model represents an indispensable tool in analysing coastal problems and assessing morphological impacts of human interference.     

A high-resolution geophysical investigation of sediment distribution controlled by catchment size and tides in a multi-basin turbid outwash fjord: Simpson Bay,Prince William Sound,Alaska

Surficial sediment distribution within Simpson Bay is a function of antecedent bedrock and recently deposited glacial geology, as well as active physical processes both within Simpson Bay and Prince William Sound (PWS). Simpson Bay is a turbid, outwash fjord located in northeastern PWS, Alaska. Freshwater from heavy precipitation, and the melting of high alpine glaciers enter the bay through bay head rivers and small shoreline creeks. The catchment has a high watershed/basin surface area ratio (∼8:1), and easily erodible bedrock that contribute to high sediment loads. The system can be divided into three discrete basins, each with specific morphologic and circulatory characters. Side scan sonar, swath bathymetry, and seismic profiles reveal that bathymetric highs are areas of outcropping glacial surfaces. High backscatter coupled with surface grab samples reveal these surfaces to be composed of coarse sediment and bedrock outcrops. Bathymetric lows are areas of low backscatter, and grab samples reveal these areas to be ponded deposits of organic-rich estuarine muds. The data provide evidence of terminal morainal bank systems, and glacial grounding line deposits at the mouth of the bay and rocky outcrops were identified as subsurface extensions of aerial rocky promontories. Radioisotope analyses of short cores reveal that the bay has an average accumulation rate of approx. 0.5 cm year⁻¹, but that this varies in function of the watershed/basin surface area ratios of the different basins. The interaction of tidal currents and sediment source drives sediment distribution in Simpson Bay. Hydrographic data reveal high spatial variability in surface and bottom currents throughout the bay. Subsurface currents are tide dominated, but generally weak (5–20 cm s⁻¹), while faster currents are found along shorelines, outcrops, and bathymetric highs. Bathymetric data reveal steep slopes with little to no modern sediment throughout the bay, suggesting lack of deposition due to tidal currents.     

MedFlux: Investigations of particle flux in the Twilight Zone

The MedFlux project was devised to determine and model relationships between organic matter and mineral ballasts of sinking particulate matter in the ocean. Specifically we investigated the ballast ratio hypothesis, tested various commonly used sampling and modeling techniques, and developed new technologies that would allow better characterization of particle biogeochemistry. Here we describe the rationale for the project, the biogeochemical provenance of the DYFAMED site, the international support structure, and highlights from the papers published here. Additional MedFlux papers can be accessed at the MedFlux web site (http://msrc.sunysb.edu/MedFlux/).     

The First Solar Seeing Profile Measurement with Two Apertures and Multiple Guide Regions

The sky brightness is a critical parameter for estimating the coronal observation conditions for a solar observatory. As part of a site-survey project in Western China, we measured the sky brightness continuously at the Lijiang Observatory in Yunnan province in 2011. A sky brightness monitor (SBM) was adopted to measure the sky brightness in a region extending from 4.5 to 7.0 apparent solar radii based on the experience of the Daniel K. Inouye Solar Telescope (DKIST) site survey. Every month, the data were collected manually for at least one week. We collected statistics of the sky brightness at four bandpasses located at 450, 530, 890, and 940 nm. The results indicate that aerosol scattering is of great importance for the diurnal variation of the sky brightness. For most of the year, the sky brightness remains under 20 millionths per airmass before local Noon. On average, the sky brightness is less than 20 millionths, which accounts for 40.41% of the total observing time on a clear day. The best observation time is from 9:00 to 13:00 (Beijing time). The Lijiang Observatory is therefore suitable for coronagraphs investigating the structures and dynamics of the corona.     

Fourier Power Spectra of Solar Noise Storms

We analyzed three noise storms recorded on 200?–?400 MHz Trieste Callisto radio spectra on 2 July 2012, 8 July 2012, and 16 July 2012 by the Fourier method. We divided intervals of the noise storms into five-minute intervals, and in these intervals we computed the mean Fourier spectra as a function of the wave numbers in the frequency and height-scale spaces. We found that these Fourier spectra, where the spectrum from the quiet-activity interval was subtracted, are power-law spectra. The mean power-law index of these spectra in the range \(\ln(k_{z}) = [1.8, 2.9]\) (where \(k_{z}\) is the wave number in the height-scale space) is \(-1.7\pm0.14\), \(-1.6\pm0.14\), and \(-1.5 \pm0.12\) for the 2 July 2012, the 8 July 2012, and the 16 July 2012 noise storms, respectively. It appears that as the number of Type-I bursts in the studied interval increases, the power-law index becomes closer to \(-5/3\); this is known as the Kolmogorov spectral index. The power-law index of the noise storms is very similar to that of the narrowband dm-spikes found in our previous studies. Furthermore, we found a break in the power spectra at \(\ln(k_{z}) \approx2.9\), and the mean power-law index values above this break are \(-2.9\pm0.46\), \(-3.1\pm0.65\), and \(-3.4\pm0.98\), respectively.     

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.