Significance tests for covariate-dependent trends in inhomogeneous spatio-temporal point processes

Modeling and inference for spatial and spatio-temporal point processes is an issue that has been broadly investigated in the last years. Application fields such as forestry, epidemiology and ecology have been the main engine driving such raised interest. The inclusion of spatially varying covariates in the models for the intensity function is becoming of particular interest, but little attention has been paid to testing the significance of such covariates. Testing the significance of covariates is important if one seeks to explain which covariates have an effect in the spatial or spatio-temporal distribution of the point pattern observed. We thus provide practical procedures to build statistical tests of significance for covariates that have an effect on the intensity function of a point pattern. Our approximation focuses on the conditional intensity function, by considering nonparametric kernel-based estimators. We calculate thinning probabilities under the conditions of absence and presence of a covariate and compare them through divergence measures. Based on Monte Carlo experiments, we approximate the statistical properties of our tests under a variety of practical scenarios. An application on testing the significance of a covariate in a spatio-temporal data set on wildfires is also developed.     

A physical interpretation of the deterministic fractal–multifractal method as a realization of a generalized multiplicative cascade

In this study, we attempt to offer a solid physical basis for the deterministic fractal–multifractal (FM) approach in geophysics (Puente, Phys Let A 161:441–447, 1992; J Hydrol 187:65–80, 1996). We show how the geometric construction of derived measures, as Platonic projections of fractal interpolating functions transforming multinomial multifractal measures, naturally defines a non-trivial cascade process that may be interpreted as a particular realization of a random multiplicative cascade. In such a light, we argue that the FM approach is as “physical” as any other phenomenological approach based on Richardson’s eddies splitting, which indeed lead to well-accepted models of the intermittencies of nature, as it happens, for instance, when rainfall is interpreted as a quasi-passive tracer in a turbulent flow. Although neither a fractal interpolating function nor the specific multipliers of a random multiplicative cascade can be measured physically, we show how a fractal transformation “cuts through” plausible scenarios to produce a suitable realization that reflects specific arrangements of energies (masses) as seen in nature. This explains why the FM approach properly captures the spectrum of singularities and other statistical features of given data sets. As the FM approach faithfully encodes data sets with compression ratios typically exceeding 100:1, such a property further enhances its “physical simplicity.” We also provide a connection between the FM approach and advection–diffusion processes.     

Fully probabilistic seismic risk assessment considering local site effects for the portfolio of buildings in Medellín,Colombia

A fully probabilistic seismic risk analysis using a comprehensive approach is conducted for Medellin, the second largest city of Colombia, using a building by building database constructed and complemented from aerial images, considering characteristics such as building use categories, socio-economic levels and replacement values. The seismic hazard used for the analysis corresponds to the most updated study available in the country with the same model that was included in the national building code maps definition. Spectral transfer functions are determined for each of the seismic microzonation zones in order to take into account the dynamic soil response and amplification effects in the risk analysis. Several building types are defined for the city and individual vulnerability functions are assigned to each of them. Risk results are presented in the state of the art metrics such as the loss exceedance curve, probable maximum losses for different return periods, average annual losses and risk maps. The obtained results can be classified by use and socio-economic sectors as well as by structural systems that may help the stakeholders to identify where the risk concentrates.     

Fitting the integrated spectral energy distributions of galaxies

Fitting the spectral energy distributions (SEDs) of galaxies is an almost universally used technique that has matured significantly in the last decade. Model predictions and fitting procedures have improved significantly over this time, attempting to keep up with the vastly increased volume and quality of available data. We review here the field of SED fitting, describing the modelling of ultraviolet to infrared galaxy SEDs, the creation of multiwavelength data sets, and the methods used to fit model SEDs to observed galaxy data sets. We touch upon the achievements and challenges in the major ingredients of SED fitting, with a special emphasis on describing the interplay between the quality of the available data, the quality of the available models, and the best fitting technique to use in order to obtain a realistic measurement as well as realistic uncertainties. We conclude that SED fitting can be used effectively to derive a range of physical properties of galaxies, such as redshift, stellar masses, star formation rates, dust masses, and metallicities, with care taken not to over-interpret the available data. Yet there still exist many issues such as estimating the age of the oldest stars in a galaxy, finer details of dust properties and dust-star geometry, and the influences of poorly understood, luminous stellar types and phases. The challenge for the coming years will be to improve both the models and the observational data sets to resolve these uncertainties. The present review will be made available on an interactive, moderated web page (), where the community can access and change the text. The intention is to expand the text and keep it up to date over the coming years.     

Confirmation of beach accretion by grain-size trend analysis: Camposoto beach, Cádiz, SW Spain

An application of the grain size trend analysis (GSTA) is used in an exploratory approach to characterize sediment transport on Camposoto beach (Cádiz, SW Spain). In May 2009 the mesotidal beach showed a well-developed swash bar on the upper foreshore, which was associated with fair-weather conditions prevailing just before and during the field survey. The results were tested by means of an autocorrelation statistical test (index I of Moran). Two sedimentological trends were recognized, i.e. development towards finer, better sorted and more negatively skewed sediment (FB–), and towards finer, better sorted and less negatively or more positively skewed sediment (FB+). Both vector fields were compared with results obtained from more classical approaches (sand tracers, microtopography and current measurements). This revealed that both trends can be considered as realistic, the FB+ trend being identified for the first time in a beach environment. The data demonstrate that, on the well-developed swash bar, sediment transported onshore becomes both finer and better sorted towards the coast. On the lower foreshore, which exhibits a steeper slope produced by breaking waves, the higher-energy processes winnow out finer particles and thereby produce negatively skewed grain-size distributions. The upper foreshore, which has a flatter and smoother slope, is controlled by lower-energy swash-backwash and overwash processes. As a result, the skewness of the grain-size distributions evolves towards less negative or more positive values. The skewness parameter appears to be distributed as a function of the beach slope and, thus, reflects variations in hydrodynamic energy. This has novel implications for coastal management.

Spectral expressions for modelling the gravitational field of the Earth’s crust density structure

We derive expressions for computing the gravitational field (potential and its radial derivative) generated by an arbitrary homogeneous or laterally varying density contrast layer with a variable depth and thickness based on methods for a spherical harmonic analysis and synthesis of gravity field. The newly derived expressions are utilised in the gravimetric forward modelling of major known density structures within the Earth’s crust (excluding the ocean density contrast) beneath the geoid surface. The gravitational field quantities due to the sediments and crust components density contrasts, shown in numerical examples, are computed using the 2 × 2 arc-deg discrete data from the global crustal model CRUST2.0. These density contrasts are defined relative to the adopted value of the reference crustal density of 2670 kgm⁻³. All computations are realised globally on a 1 × 1 arc-deg geographical grid at the Earth’s surface. The maxima of the gravitational signal due to the sediments density contrast are mainly along continental shelf regions with the largest sedimentary deposits. The corresponding maxima due to the consolidated crust components density contrast are over areas of the largest continental crustal thickness with variable geological structure.     

Topographic predictors of susceptibility to alluvial fan flooding,Southern Apennines

The flooding susceptibility of alluvial fans in the Southern Apennines has long been neglected. To partly address this oversight, we focus on the region of Campania which contains highly urbanized piedmont areas particularly vulnerable to flooding. Our findings are based on stratigraphic analysis of the fans and morphometric analysis of the basin‐fan systems. Using geomorphological analysis we recognized active alluvial fans while stratigraphic analysis together with statistical analysis of the morphometric variables was used to classify the fans in terms of the transport process involved. The results indicate that in the geological context examined, the best discrimination between debris flow (Df) and water flood (Wf) processes is achieved by means of two related variables, one for the basin (feeder channel inclination, Cg) and one for the fan (fan length, Fl). The probability that an unclassified fan belongs to group Wf is computed by applying a logistic function in which a P value exceeding 0.5 indicates that a basin/fan system belongs to group Wf. This important result led to the classification of the entire basin/fan system data. As regards process intensity, debris flow‐dominated fans are susceptible to the occurrence of flows with high viscosity and hence subject to more severe events than water flood‐dominated fans. Bearing this in mind, the data gathered in this study allow us to detect where alluvial fan flooding might occur and give information on the different degrees of susceptibility at a regional scale. Regrettably, urban development in recent decades has failed to take the presence of such alluvial fans into account due to the long recurrence time (50–100 years) between floods. This paper outlines the distribution of such susceptibility scenarios throughout the region, thereby constituting an initial step to implementing alluvial fan flooding control and mitigation. Copyright © 2012 John Wiley & Sons, Ltd.     

Foundation conditions of buildings in undermined areas: An example of evaluation

Evaluation of a set of individual geological and petrophysical parameters is necessary in many geoscientific studies or those with constructional purposes. However, correlation and mutual evaluation of these parameters are very complicated. Our methodology, which is part of the expert evaluation, is based on compilation of four significant effects: local subsurface geology, groundwater table (and accumulation of surface water), deformation of the surface due to mining (subsidence), and landslides. Different values of specified significances were assigned to all four selected effects. The class of foundation conditions was selected according to the summary of values of specified significances for a given point. This value describes how the selected parameters influence vibration effect on the surface and resistance of buildings to this seismic loading. The pilot map of foundation conditions, which was elaborated as part of the evaluation of Stonava area in 2008, is presented.     

Coral communities and reefs from Guerrero,Southern Mexican Pacific

Corals in the Eastern Pacific extend south from the Gulf of California to Ecuador and oceanic Chile, and west from Colombia to Clipperton Atoll. Nevertheless, large stretches of the Mexican Pacific remain fundamentally unstudied. Therefore, to assess the current conditions of coral communities, a coastal fringe ～300 km long (17°40′ N, 101°39′ W to 16°46′ N, 99°49′ W) was surveyed within the Southern Mexican Pacific, between 2005 and 2009. Fifteen stony coral species were identified at 13 coral communities and six Pocillopora‐dominated fringing reefs, with Pocillopora verrucosa and Pocillopora damicornis the primary contributing taxa. Reef development was identified in embayments or behind rocks or islands that offered shelter from northern and northwestern winds. Observations of Pocillopora effusus, Pocillopora inflata, Porites lobata, Pavona clavus, and Pavona varians expanded the species known geographic ranges by several degrees of latitude, suggesting reef building fauna comprised a mixture of widespread and relatively rare Eastern Pacific corals. Results indicated greater live coral cover in the Ixtapa‐Zihuatanejo area (15–73%) than in the Acapulco localities, which had high algal dominance; the reefs in the latter region exhibited high erosion. Regional differences are likely the result of long‐standing anthropogenic pressures around Acapulco since 1950, when it became an important tourist destination. This paper is the first detailed report of ecologically stressed corals and coral reefs from the state of Guerrero on the Mexican Southern Pacific coast.     

Flood risk assessment in a coastal lagoon under present and future scenarios: Ria de Aveiro case study

Iranian strong motion records as well as detailed conditions of their instrument sites and the characteristics of their causative seismic sources are compiled and processed. The dataset consists of 2286 three-component records from 461 Iranian earthquakes with at least two high-quality records having moment magnitude from 3.9 to 7.3. These records are about 20% of the Iranian database and are suitable for seismic hazard analysis and engineering applications. Perhaps for the first time in the literature, the distance to the surface projection of the fault is reported for a great number of records corresponding to earthquakes with M > 6.0. The raw accelerations are processed using the wavelet de-noising method. Having corrected and filtered these raw data, the pseudospectral accelerations are calculated for each of the three components of time series, separately. In addition to the ground motion parameters, a large and comprehensive list of metadata characterizing the recording conditions of each record is also developed. Moreover, careful revision of the characteristics of the earthquakes such as location, magnitude, style of faulting and fault rupture plane geometry, if available, is carried out using the best available information in a scientifically sound manner. Finally, we also focus on special ground motion records including records with peak ground acceleration (PGA) >300 cm/s² and distances less than 30 km. These are “exceptional” records in the Iranian dataset and include less than 2% of the selected dataset.