Performance of Jacobi preconditioning in Krylov subspace solution of finite element equations

This paper examines the performance of the Jacobi preconditioner when used with two Krylov subspace iterative methods. The number of iterations needed for convergence was shown to be different for drained, undrained and consolidation problems, even for similar condition number. The differences were due to differences in the eigenvalue distribution, which cannot be completely described by the condition number alone. For drained problems involving large stiffness ratios between different material zones, ill‐conditioning is caused by these large stiffness ratios. Since Jacobi preconditioning operates on degrees‐of‐freedom, it effectively homogenizes the different spatial sub‐domains. The undrained problem, modelled as a nearly incompressible problem, is much more resistant to Jacobi preconditioning, because its ill‐conditioning arises from the large stiffness ratios between volumetric and distortional deformational modes, many of which involve the similar spatial domains or sub‐domains. The consolidation problem has two sets of degrees‐of‐freedom, namely displacement and pore pressure. Some of the eigenvalues are displacement dominated whereas others are excess pore pressure dominated. Jacobi preconditioning compresses the displacement‐dominated eigenvalues in a similar manner as the drained problem, but pore‐pressure‐dominated eigenvalues are often over‐scaled. Convergence can be accelerated if this over‐scaling is recognized and corrected for. Copyright © 2002 John Wiley & Sons, Ltd.     

Multi-wavelength observations of Asteroid 2100 Ra-Shalom

We observed near-Earth asteroid (NEA) 2100 Ra-Shalom over a six-year period, obtaining rotationally resolved spectra in the visible, near-infrared, thermal-infrared, and radar wavelengths. We find that Ra-Shalom has an effective diameter of Deff=2.3±0.2 km, rotation period P=19.793±0.001 h, visual albedo pv=0.13±0.03, radar albedo , and polarization ratio μc=0.25±0.04. We used our radar observations to generate a three-dimensional shape model which shows several structural features of interest. Based on our thermal observations, Ra-Shalom has a high thermal inertia of ∼10³ J m⁻² s^−0.5 K⁻¹, consistent with a coarse or rocky surface and the inferences of others [Harris, A.W., Davies, J.K., Green, S.F., 1998. Icarus 135, 441-450; Delbo, M., Harris, A.W., Binzel, R.P., Pravec, P., Davies, J.K., 2003. Icarus 166, 116-130]. Our spectral data indicate that Ra-Shalom is a K-class asteroid and we find excellent agreement between our spectra and laboratory spectra of the CV3 meteorite Grosnaja. Our spectra show rotation-dependent variations consistent with global variations in grain size. Our radar observations show rotation-dependent variations in radar albedo consistent with global variations in the thickness of a relatively thin regolith.     

Tropical forcing of Australian extreme low minimum temperatures in September 2019

We explore the causes and predictability of extreme low minimum temperatures (T_min) that occurred across northern and eastern Australia in September 2019. Historically, reduced T_min is related to the occurrence of a positive Indian Ocean Dipole (IOD) and central Pacific El Niño. Positive IOD events tend to locate an anomalous anticyclone over the Great Australian Bight, therefore inducing cold advection across eastern Australia. Positive IOD and central Pacific El Niño also reduce cloud cover over northern and eastern Australia, thus enhancing radiative cooling at night-time. During September 2019, the IOD and central Pacific El Niño were strongly positive, and so the observed T_min anomalies are well reconstructed based on their historical relationships with the IOD and central Pacific El Niño. This implies that September 2019 T_min anomalies should have been predictable at least 1–2 months in advance. However, even at zero lead time the Bureau of Metereorolgy ACCESS-S1 seasonal prediction model failed to predict the anomalous anticyclone in the Bight and the cold anomalies in the east. Analysis of hindcasts for 1990–2012 indicates that the model's teleconnections from the IOD are systematically weaker than the observed, which likely stems from mean state biases in sea surface temperature and rainfall in the tropical Indian and western Pacific Oceans. Together with this weak IOD teleconnection, forecasts for earlier-than-observed onset of the negative Southern Annular Mode following the strong polar stratospheric warming that occurred in late August 2019 may have contributed to the T_min forecast bust over Australia for September 2019.

     

Geochemistry of soils of King George Island, South Shetland Islands, West Antarctica: Implications for pedogenesis in cold polar regions

Fine fractions of soils on the Barton Peninsula, King George Island, West Antarctica have been forming during the last 6000 yr since the last deglaciation. Texturally, they are mostly composed of mineral and rock fragments with some volcanic ashes, which are also indicated by geochemical compositions representing for the nonclay silicate minerals and low values of chemical index of alteration. No significant changes are observed in major- and trace element abundances. Such geochemical characteristics suggest that chemical weathering of bedrocks on the Barton Peninsula seems insignificant and that the soils are composed of physically weathered mineral and rock fragments which are mixed with eolian additions of volcanic ashes and Patagonian dusts. Chondrite-normalized rare earth element (REE) distribution patterns of the Barton Peninsula soils are slightly different from those of bedrocks, indicating that the REE abundances and characteristics were influenced by eolian additions. Mixing calculations, which mass-balance the REEs, suggest that volcanic ashes blown from Deception Island were the major eolian contributor, followed by atmospheric dusts sourced from Patagonia, South America. Even in the warmer and humid climatic conditions in the maritime Antarctic region, the chemical weathering of bedrocks appears to be insignificant, probably due to the relatively short duration of weathering since the last deglaciation.     

An agent-based evacuation model for the 2011 Brisbane City-scale riverine flood

In this study, an agent-based model is proposed in order to provide new insights into the policy analysis and strategy assessment of city-scale evacuation management. The proposed model is suitable for assessment of the influence of different departure times and communications among peer evacuees on the number of residents at risk who arrive at official shelters. A case study is applied to build a simulation model for the coastal city of Brisbane in Australia. The Brisbane River catchment experiences regular flooding almost every year; the second severest flood since the twentieth century occurred in 2011. During that event, over 15,000 properties were inundated and around 3600 households evacuated in metropolitan Brisbane alone. Making use of high-performance computing clusters, the evacuation simulation was coupled with results from a validated hydrodynamic model to test a variety of escaping scenarios based on the 2011 flood situation. This case study demonstrates the proposed model’s capacity to represent the dynamic evacuation process and also shows that the model is able to help develop flood emergency plans and evaluate response measures through exploring key elements in a range of scenarios.