Numerical modeling of wave equation by a truncated high-order finite-difference method

Finite-difference methods with high-order accuracy have been utilized to improve the precision of numerical solution for partial differential equations.However, the computation cost generally increases linearly with increased order of accuracy.Upon examination of the finite-difference formulas for the first-order and second-order derivatives, and the staggered finite-difference formulas for the first-order derivative, we examine the variation of finite-difference coefficients with accuracy order and note th...     

The NAFE’06 data set: Towards soil moisture retrieval at intermediate resolution

The National Airborne Field Experiment 2006 (NAFE’06) was conducted during a three week period of November 2006 in the Murrumbidgee River catchment, located in southeastern Australia. One objective of NAFE’06 was to explore the suitability of the area for SMOS (Soil Moisture and Ocean Salinity) calibration/validation and develop downscaling and assimilation techniques for when SMOS does come on line. Airborne L-band brightness temperature was mapped at 1 km resolution 11 times (every 1–3 days) over a 40 by 55 km area in the Yanco region and 3 times over a 40 by 50 km area that includes Kyeamba Creek catchment. Moreover, multi-resolution, multi-angle and multi-spectral airborne data including surface temperature, surface reflectance (green, read and near infrared), lidar data and aerial photos were acquired over selected areas to develop downscaling algorithms and test multi-angle and multi-spectral retrieval approaches. The near-surface soil moisture was measured extensively on the ground in eight sampling areas concurrently with aircraft flights, and the soil moisture profile was continuously monitored at 41 sites. Preliminary analyses indicate that (i) the uncertainty of a single ground measurement was typically less than 5% vol. (ii) the spatial variability of ground measurements at 1 km resolution was up to 10% vol. and (iii) the validation of 1 km resolution L-band data is facilitated by selecting pixels with a spatial soil moisture variability lower than the point-scale uncertainty. The sensitivity of passive microwave and thermal data is also compared at 1 km resolution to illustrate the multi-spectral synergy for soil moisture monitoring at improved accuracy and resolution. The data described in this paper are available at www.nafe.unimelb.edu.au.     

Aeolian fetch distance and secondary airflow effects: the influence of micro‐scale variables on meso‐scale foredune development

Unsuccessful attempts to use process‐scale models to predict long‐term aeolian sediment transport patterns have long been a feature of aeolian research. It has been proposed that one approach to overcome these problems is to identify micro‐scale variables that are important at longer timescales. This paper assesses the contribution of two system variables (secondary airflow patterns and fetch distance) to medium‐term (months to years) dune development. The micro‐scale importance of these variables had been established during previous work at the site (Magilligan Strand, Northern Ireland). Three methods were employed. First, sand drift potentials were calculated using 2 years of regional wind data and a sediment transport model. Second, wind data and large trench traps (2 m length × 1 m width × 1·5 m depth) were used to assess the actual sediment transport patterns over a 2‐month period. Third, a remote‐sensing technique for the identification of fetch distance, a saltation impact sensor (Safire) and wind data were utilized to gauge, qualitatively, sediment transport patterns over a 1‐month period. Secondary airflow effects were found to play a major role in the sediment flux patterns at these timescales, with measured and predicted rates matching closely during the trench trap study. The results suggest that fetch distance is an unimportant variable at this site. Copyright © 2007 John Wiley & Sons, Ltd.     

The tectonic and geomagnetic significance of paleomagnetic observations from volcanic rocks from central Afar, Africa

Deformation throughout Afar over the past 2 myr has been characterized by widespread and intense crustal fragmentation that results from inhomogeneous extension across the region. In eastern Afar, this situation has evolved to localized extension associated with the westward propagation of the Gulf of Aden/Gulf of Tadjurah seafloor spreading system into the Asal–Ghoubbet Rift. During the gradual process of rift propagation and localization, crustal blocks in eastern Afar sustained clockwise rotations of 11°. To better understand the processes of rift propagation and localization and how they affect the rest of Afar, we have collected and analyzed over 400 oriented paleomagnetic samples from 67 lava flows from central and southern Afar. Unlike eastern Afar, the mean paleomagnetic direction from central Afar indicates that vertical-axis rotations are statistically insignificant (3.6°±4.4°), though small clockwise rotations (<8°) are permitted. Thus, propagation and localization in central Afar have not had the same influence in causing crustal block rotations or, perhaps more likely, have not reached the same stage of evolution as seen in eastern Afar. In addition, several of the lava flows record intriguing geomagnetic field behavior associated with polarity transitions, excursions, or large secular variation events. Interestingly, the transitional or anomalous virtual geomagnetic poles (VGPs) tend to cluster in two nearly antipodal regions, one in the northern Pacific Ocean and the other in the southwest Indian Ocean. One lava flow has recorded both of the antipodal transitional components, with the two components residing in magnetic minerals with unblocking temperatures above and below 500°C, respectively. Reheating and partial remagnetization by the overlying flow cannot explain either of the transitional directions because both differ significantly from that of the reversely magnetized overlying flow. The high-temperature component gives a VGP in the northern Pacific, whereas the lower-temperature component gives a nearly antipodal VGP south of Cape Town, South Africa. Hence, the configuration of the geomagnetic field appears to have jumped nearly instantaneously from a northern-hemisphere transitional state to a southern-hemisphere one during this normal-to-reverse polarity transition.     

Simple automatic time-stepping for improved simulation of groundwater hydrographs

An automatic time-stepping algorithm is presented, based on the intensity of driving groundwater recharge, that improves the simulation of groundwater level fluctuations in regional models while maintaining model run-times. The algorithm is implemented in the ZOOMQ3D finite difference groundwater flow code and controls the discretization of time using two user-defined criteria: a maximum time-step length and a maximum recharge per time-step. Daily recharge is accumulated in time until either of these criteria is violated when the model then calculates a solution. The efficiency and accuracy of the algorithm is tested using an idealized groundwater model and an existing regional groundwater model of a UK aquifer. The approach is illustrated using simulations of high groundwater levels and associated groundwater flood events in a responsive, high-diffusivity aquifer. Simulations using the automatic time-stepping technique are presented that reduce the maximum absolute error in groundwater level by 45% and the run-time by 51% compared to models using conventional, a priori defined stress-periods and time-steps.     

The Australia Telescope 20-GHz (AT20G) Survey: the Bright Source Sample

The Australia Telescope 20-GHz (AT20G) Survey is a blind survey of the whole southern sky at 20 GHz (with follow-up observations at 4.8 and 8.6 GHz) carried out with the Australia Telescope Compact Array from 2004 to 2007.
The Bright Source Sample (BSS) is a complete flux-limited subsample of the AT20G Survey catalogue comprising 320 extragalactic     ) radio sources south of  δ=−15°  with      Jy. Of these, 218 have near simultaneous observations at 8 and 5 GHz.
In this paper we present an analysis of radio spectral properties in total intensity and polarization, size, optical identifications and redshift distribution of the BSS sources. The analysis of the spectral behaviour shows spectral curvature in most sources with spectral steepening that increases at higher frequencies (the median spectral index α, assuming   S ∝ν^α  , decreases from  α^8.6_4.8= 0.11  between 4.8 and 8.6 GHz to  α²⁰_8.6=−0.16  between 8.6 and 20 GHz), even if the sample is dominated by flat spectra sources (85 per cent of the sample has  α²⁰_8.6 > −0.5)  . The almost simultaneous spectra in total intensity and polarization allowed us a comparison of the polarized and total intensity spectra: polarized fraction slightly increases with frequency, but the shapes of the spectra have little correlation. Optical identifications provided an estimation of redshift for 186 sources with a median value of 1.20 and 0.13, respectively, for QSO and galaxies.     

Contrasting geomorphological storm response from two adjacent shorefaces

Shorefaces play a critical role in cross‐shore sediment transport between the beach and inner shelf, particularly during storm conditions. A comparison and examination of storm‐driven sedimentary changes on two adjacent shorefaces in Northern Ireland, located only 5 km apart, revealed significantly different geomorphological responses. The steeper shoreface at West Strand responded with extensive sediment deposition across almost the entire shoreface, in contrast with the more dissipative and quasi‐linear shoreface at Portstewart, which mostly showed nearshore bar changes. Results from the two sites, which have similar wave/wind characteristics and seabed sediments, suggest that: (i) cross‐shore morphology, (ii) immediately previous (antecedent) shoreface morphodynamic behaviour and (iii) the presence, or lack of, offshore sand appear to be the primary controls on storm‐driven sedimentary changes attributed to the high‐energy event. Copyright © 2015 John Wiley & Sons, Ltd.     

Major element chemistry of ocean island basalts — Conditions of mantle melting and heterogeneity of mantle source

We estimate average compositions of near-primary, ‘reference’ ocean island basalts (OIBs) for 120 volcanic centers from 31 major island groups and constrain the depth of lithosphere–asthenosphere boundary (LAB) at the time of volcanism and the possible depth of melt–mantle equilibration based on recently calibrated melt silica activity barometer. The LAB depth versus fractionation corrected OIB compositions (lava compositions, X, corrected to Mg# 73, X_OIB^#73, i.e., magmas in equilibrium with Fo₉₀, if olivine is present in the mantle source) show an increased major element compositional variability with increasing LAB depths. OIBs erupted on lithospheres < 40 km thick approach the compositions (e.g. SiO₂^#73, TiO₂^#73, [CaO/Al₂O₃]^#73) of primitive ridge basalts and are influenced strongly by depth and extent of shallow melting. However, X_OIB^#73 on thicker lithospheres cannot be explained by melt–mantle equilibration as shallow as LAB. Melt generation from a somewhat deeper (up to 50 km deeper than the LAB) peridotite source can explain the OIB major element chemistry on lithospheres ≤ 70 km. However, deeper melting of volatile-free, fertile peridotite is not sufficient to explain the end member primary OIBs on ≥ 70 km thick lithospheres. Comparison between X_OIB^#73 and experimental partial melts of fertile peridotite indicates that at least two additional melt components need to be derived from OIB source regions. The first component, similar to that identified in HIMU lavas, is characterized by low SiO₂^#73, Al₂O₃^#73, [Na₂O/TiO₂]^#73, and high FeO?^#73, CaO^#73, [CaO/Al₂O₃]^#73. The second component, similar to that found in Hawaiian Koolau lavas, is characterized by high SiO₂^#73, moderately high FeO?^#73, and low CaO^#73 and Al₂O₃^#73. These two components are not evenly sampled by all the islands, suggesting a heterogeneous distribution of mantle components that generate them. We suggest that carbonated eclogite and volatile-free, silica-excess eclogite are the two most likely candidates, which in conjunction with fertile mantle peridotite, give rise to the two primitive OIB end members.     

Using the ETAS Model for Catalog Declustering and Seismic Background Assessment

The concept of background seismicity is strictly related to the identification of spontaneous and triggered earthquakes. The definition of foreshocks, main shocks and aftershocks is currently based on procedures depending on parameters whose values are notoriously assumed by subjective criteria. We propose a method for recognizing the background and the induced seismicity statistically. Rather than using a binary distinction of the events in these two categories, we prefer to assign to each of them a probability of being independent or triggered. This probability comes from an algorithm based on the ETAS model. A certain degree of subjectivity is still present in this procedure, but it is limited by the possibility of adjusting the free parameters of the algorithm by rigorous statistical criteria such as maximum likelihood. We applied the method to the seismicity of southern California and analyzed the sensitivity of the results to the free parameters in the algorithm. Finally, we show how our statistical declustering algorithm may be used for mapping the background seismicity, or the moment rate in a seismic area.