Outlier identification and correction for GRACE aggregated data

Gravity measurements within the Gravity Recovery and Climate Experiment (GRACE) provide a direct measure of monthly changes in mass over the Earth’s land masses. As such changes in mass mainly correspond to water storage changes, these measurements allow to close the continental water balance on large spatial scales and on a monthly time scale within the respective error bounds. When quantifying uncertainties, positive and negative peaks are detected in GRACE aggregated monthly time series (from different data providers) that do not correspond to hydrological or hydro-meteorological signals. These peaks must be interpreted as outliers, which carry the danger of signal degradation. In this paper an algorithm is developed to identify outliers and replace them with hydrologically plausible values. The algorithm is based on a statistical approach in which hydrological and hydro-meteorological signals are used to control the algorithm. The procedure of outlier detection is verified by evaluating catchment based aggregated GRACE monthly signals with ground truth from hydrology and hydro-meteorological signals. The results show improvement in the correlation of GRACE versus hydrometeorological and hydrological signals in most catchments. Also, the noise level is significantly reduced over 255 largest catchments.     

Assessment of agri-spillways as a soil erosion protection measure in Mediterranean sloping vineyards

Suitable vineyard soils enhance soil stability and biodiversity which in turn protects roots against erosion and nutrient losses. There is a lack of information related to inexpensive and suitable methods and tools to protect the soil in Mediterranean sloping vineyards(25° of slope inclination). In the vineyards of the Montes de Málaga(southern Spain), a sustainable land management practice that controls soil erosion is actually achieved by tilling rills in the down-slope direction to canalize water and sediments. Because of their design and use, we call them agri-spillways. In this research, we assessed two agri-spillways(between 10 m and 15 m length, and slopes between 25.8° and 35°) by performing runoff experiments under extreme conditions(a motor driven pump that discharged water flows up to 1.33 l s~(-1) for 12 to 15 minutes: ≈1000 l). The final results showed: i) a great capacity by these rills to canalize large amounts of water and sediments; and, ii) higher water flow speeds(between 0.16 m s-1 and 0.28 m s~(-1)) and sediment concentrationrates(up to 1538.6 g l~(-1)) than typically found in other Mediterranean areas and land uses(such as badlands, rangelands or extensive crops of olives and almonds). The speed of water flow and the sediment concentration were much higher in the shorter and steeper rill. We concluded that agri-spillways, given correct planning and maintenance, can be a potential solution as an inexpensive method to protect the soil in sloping Mediterranean vineyards.     

I Zwicky 1: Decomposition and dynamics of the nearby QSO host

We address current theories of research on morphology, environment and evolution of active galaxies in the form of a detailed case study of the nearby QSO host I Zw 1. This study is based on sub-kpc resolution ¹²CO(1-0) observations with the BIMA mm-interferometer and on near-infrared imaging with ISAAC at the VLT. The ¹²CO(1-0) maps reveal a circumnuclear molecular gas ring with a radius of 0.9 kpc. The imaging data in the J-band are analyzed with respect to a disk-bulge decomposition. Together with a model of the gas rotation curve, the radial profile of the dynamical J-band mass-to-light ratio (M/L) is discussed. The J-band images give new evidence for an interaction between I Zw 1 and its nearby companion. First results from a sequence of N-body simulations for an extensive test of the parameter space of the interaction are presented. This revised version was published online in August 2006 with corrections to the Cover Date.     

On the impact of salinity observations on state estimates in Ems Estuary

The hydrodynamics of Ems Estuary are dominated by tides and their interaction with buoyancy forcing. Such an environment is challenging for any effort to bring together observations and model results. In this study, we investigate how salinity measurements in the Ems Estuary affect the reconstruction of the salinity field. Similar to the traditional observing system experiments, the impact of specific observational arrays is simulated in the framework of statistical experiments. The experimental algorithm mainly relies on the model covariance matrix. Each experiment results in an estimate of the reconstruction error. The analysed observation configurations involve single and multiple, as well as stationary and non-stationary observing arrays. Generally, the reconstruction of the ocean state improves with increasing the density of observations. It appears that certain locations are more favourable for reconstruction than others. In fact, the regions separating the main dynamical realms resist strongest to the reconstruction effort. Extending the covariance matrix by the temporal cross-covariances between the model grid points enables to evaluate the impact of observations taken from a moving platform. This approach further improves the outcome of the experiments, resulting in reconstruction errors near zero with the exception of the tidal river. The cross-covariance information is able to tackle even the irregular dynamics arising on the border between the different physical regimes.     

On Flare-CME Characteristics from Sun to Earth Combining Remote-Sensing Image Data with <Emphasis Type="Italic">In Situ</Emphasis> Measurements Supported by Modeling

We analyze the well-observed flare and coronal mass ejection (CME) from 1 October 2011 (SOL2011-10-01T09:18) covering the complete chain of effects – from Sun to Earth – to better understand the dynamic evolution of the CME and its embedded magnetic field. We study in detail the solar surface and atmosphere associated with the flare and CME using the Solar Dynamics Observatory (SDO) and ground-based instruments. We also track the CME signature off-limb with combined extreme ultraviolet (EUV) and white-light data from the Solar Terrestrial Relations Observatory (STEREO). By applying the graduated cylindrical shell (GCS) reconstruction method and total mass to stereoscopic STEREO-SOHO (Solar and Heliospheric Observatory) coronagraph data, we track the temporal and spatial evolution of the CME in the interplanetary space and derive its geometry and 3D mass. We combine the GCS and Lundquist model results to derive the axial flux and helicity of the magnetic cloud (MC) from in situ measurements from Wind. This is compared to nonlinear force-free (NLFF) model results, as well as to the reconnected magnetic flux derived from the flare ribbons (flare reconnection flux) and the magnetic flux encompassed by the associated dimming (dimming flux). We find that magnetic reconnection processes were already ongoing before the start of the impulsive flare phase, adding magnetic flux to the flux rope before its final eruption. The dimming flux increases by more than 25% after the end of the flare, indicating that magnetic flux is still added to the flux rope after eruption. Hence, the derived flare reconnection flux is most probably a lower limit for estimating the magnetic flux within the flux rope. We find that the magnetic helicity and axial magnetic flux are lower in the interplanetary space by ～?50% and 75%, respectively, possibly indicating an erosion process. A CME mass increase of 10% is observed over a range of \({\sim}\,4\,\mbox{--}\,20~\mathrm{R}_{\odot }\). The temporal evolution of the CME-associated core-dimming regions supports the scenario that fast outflows might supply additional mass to the rear part of the CME.     

Determination of Mislocation Vectors to Evaluate bias at GSETT-3 primary stations

Using the database provided by the Reviewed Event Bulletins (REBs) for the first 2.5 years of the Group of Scientific Experts Technical Test-3 (GSETT-3) experiment, we compiled mislocation vectors for both arrays and selected three-component stations of the primary network from the published slowness and azimuth information gained through f-k- and polarization analysis. Imposing constraints such as a minimum signal-to-noise ratio (SNR) and number of defining phases, we aim at eliminating location bias as the hypocentral parameters are taken from the REBs. Results from 14 arrays with apertures from about 1 km to more than 20 km are presented as well as from 18 three-component stations, which indicate that the mislocation vectors in many cases can improve location accuracy considerably. If these mislocation vectors are compiled to provide coverage of a sufficient portion of the slowness domain these empirical corrections can easily be applied prior to location processing. In the context of the Comprehensive Nuclear Test-Ban Treaty (CTBT), these mislocation patterns could be essential for providing accurate event location of suspicious low-magnitude events, as these location parameters will be used to pinpoint the area where to conduct an on-site inspection.