Geoarchaeology of the Cerro La China locality (Buenos Aires,Argentina): Site 2 and site 3

Cerro La China, within the Tandilia Range, has a long archaeological record in which three different human occupations were identified. The geoarchaeological approach allowed us to correlate the archaeosedimentary columns at Site 2 and Site 3 and place the local stratigraphy in a regional context. The archaeologically rich stratigraphic sequence is composed of three episodes of loess deposition. The first bears a 10,600 yr B.P. Paleoindian occupation with Fishtail or Fell's Cave Stemmed points. It is followed by a soil-forming interval which ended about 5000 yr B.P. with an erosive episode. The second event of loess deposition continued through an unknown span of time until a new soil-forming interval started; it includes triangular stemless medium size projectile points corresponding to a different human group. By the sixteenth century another erosive event interrupted the soil-forming interval. After this, the third short-term eolian episode occurred and bears the most recent archaeological occupation marked by small triangular stemless points.     

How much does an extreme rainfall event cost? Material damage and relationships between insurance,rainfall, land cover and urban flooding

This research aims to understand how insurance, rainfall, land cover and urban flooding are related and how these variables influenced the material damage in the Lisbon Metropolitan Area (LMA) during the 2000–2011 period. Correlation coefficients show strong relationships between built-up areas and claims (0.94) and payouts (0.88). Despite no significant relationships being found between rainfall and the amount of material damage per event, three likelihood levels of flooding were determined for hourly rainfall. Unlike the studied period, the number of claims and their spatial distribution during the 2008 extreme rainfall event were strongly dependent on rainfall. Flooding related to the old watercourses assumed greater importance during this extreme event, recovering a more natural/ancient hydrological behaviour. In the LMA, the greatest material damage was the result of high-magnitude/low-probability rainfall events. Lower magnitude events can trigger numerous claims in heavily built-up areas, but they are hardly capable of producing large material damage.     

An analysis of carrier phase differential kinematic GPS positioning using DynaPos

A series of activities have been carried out at the University of New Brunswick in an effort to evaluate advances in long-range marine kinematic differential positioning. These activities involved processing and analysis of GPS carrier phase kinematic data sets. Some of the data was collected by UNB and some was provided by The XYZs of GPS Inc. Data were collected using Trimble 5700 and Ashtech Z-12 receivers. The data sets were processed using the software DynaPos provided by the The XYZs of GPS Inc. The best results obtained in our analysis indicate an agreement of 5 cm RMS for the horizontal component and 12 cm RMS for the vertical component between two ionospheric-delay free solutions, in baselines varying from 40 to 100 km.     

Evaluation of methods of spatial interpolation for monthly rainfall data over the state of Rio de Janeiro,Brazil

Five deterministic methods of spatial interpolation of monthly rainfall were compared over the state of Rio de Janeiro, southeast Brazil. The methods were the inverse distance weight (IDW), nearest neighbor (NRN), triangulation with linear interpolation (TLI), natural neighbor (NN), and spline tension (SPT). A set of 110 weather stations was used to test the methods. The selection of stations had two criteria: time series longer than 20 years and period of data from 1960 to 2009. The methods were evaluated using cross-validation, linear regression between values observed and interpolated, root mean square error (RMSE), coefficient of determination (r ²), coefficient of variation (CV, %), and the Willmott index of agreement (d). The results from different methods are influenced by the meteorological systems and their seasonality, as well as by the interaction with the topography. The methods presented higher precision (r ²) and accuracy (d, RMSE) during the summer and transition to autumn, in comparison with the winter or spring months. The SPT had the highest precision and accuracy in relation to other methods, in addition to having a good representation of the spatial patterns expected for rainfall over the complex terrain of the state and its high spatial variability.     

Multivariate geostatistical application for climate characterization of Minas Gerais State, Brazil

The objective of the present study was to assess for Minas Gerais the cokriging methodology, in order to characterize the spatial variability of Thornthwaite annual moisture index, annual rainfall, and average annual air temperature, based on geographical coordinates, altitude, latitude, and longitude. The climatic element data referred to 39 INMET climatic stations located in the state of Minas Gerais and in nearby areas and the covariables altitude, latitude, and longitude to the SRTM digital elevation model. Spatial dependence of data was observed through spherical cross semivariograms and cross covariance models. Box–Cox and log transformation were applied to the positive variables. In these situations, kriged predictions were back-transformed and returned to the same scale as the original data. Trend was removed using global polynomial interpolation. Universal simple cokriging best characterized the climate variables without tendentiousness and with high accuracy and precision when compared to simple cokriging. Considering the satisfactory implementation of universal simple cokriging for the monitoring of climatic elements, this methodology presents enormous potential for the characterization of climate change impact in Minas Gerais state.     

Estimating Buoyancy Heat Flux Using the Surface Renewal Technique over Four Amazonian Forest Sites in Brazil

The surface renewal (SR) method was applied for the first time to measurements of air temperature over four Amazonian forest sites and different seasons in order to obtain estimates of buoyancy heat flux. The required calibration of this method against eddy covariance resulted in a value for a specific parameter that is close to the range reported in other studies, contributing to the generalization of the SR method to different kinds of canopies. The comparison with fluxes obtained using the eddy covariance technique revealed a good match between the two methods for different sites, heights and seasons. Sites with high levels of non-stationarity in the signals of temperature and wind speed presented higher scatter in the regression with fluxes from eddy covariance. For a particular site with previously reported influence of low-frequency motions, the regression was only satisfactory, i.e., slope parameter close to unity and small offset, when oscillations with periods longer than $\approx $ ≈ 13 min were filtered out. The SR method has a great potential due to the simplicity of the instrumentation required. However, care should be taken when measuring under the influence of mesoscale motions, which can lead to high levels of non-stationarity, compromising the fundamental concepts of the SR theory.     

On the Nature of the Transition Between Roll and Cellular Organization in the Convective Boundary Layer

Both observational and numerical studies of the convective boundary layer (CBL) have demonstrated that when surface heat fluxes are small and mean wind shear is strong, convective updrafts tend to organize into horizontal rolls aligned within 10–20\(^\circ \) of the geostrophic wind direction. However, under large surface heat fluxes and weak to negligible shear, convection tends to organize into open cells, similar to turbulent Rayleigh-Bénard convection. Using a suite of 14 large-eddy simulations (LES) spanning a range of \(-z_i/L\) between zero (neutral) and 1041 (highly convective), where \(z_i\) is the CBL depth and L is the Obukhov length, the transition between roll- and cellular-type convection is investigated systematically for the first time using LES. Mean vertical profiles including velocity variances and turbulent transport efficiencies, as well the “roll factor,” which characterizes the rotational symmetry of the vertical velocity field, indicate the transition occurs gradually over a range of \(-z_i/L\); however, the most significant changes in vertical profiles and CBL organization occur from near-neutral conditions up to about \(-z_i/L \approx \) 15–20. Turbulent transport efficiencies and quadrant analysis are used to characterize the turbulent transport of momentum and heat with increasing \(-z_i/L\). It is found that turbulence transports heat efficiently from weakly to highly convective conditions; however, turbulent momentum transport becomes increasingly inefficient as \(-z_i/L\) increases.     

A Hybrid Spectral/Finite-Volume Algorithm for Large-Eddy Simulation of Scalars in the Atmospheric Boundary Layer

Pseudospectral methods are frequently used in the horizontal directions in large-eddy simulation of atmospheric flows. However, the same approach often creates unphysical oscillations for scalar fields if there are horizontal heterogeneities in the sources and/or sinks, as is usual in air pollution problems. A hybrid approach is developed to combine the use of pseudospectral representation of the velocity field and bounded finite-volumes for the scalar concentration. An interpolation scheme that yields a divergence-free interpolated velocity field is derived and implemented, and its importance is illustrated by two sample applications.     

Assessment of troposphere mapping functions using three-dimensional ray-tracing

The troposphere delay is an important source of error for precise GNSS positioning due to its high correlation with the station height parameter. It has been demonstrated that errors in mapping functions can cause sub-annual biases as well as affect the repeatability of GNSS solutions, which is a particular concern for geophysical studies. Three-dimensional ray-tracing through numerical weather models (NWM) is an excellent approach for capturing the directional and daily variation of the tropospheric delay. Due to computational complexity, its use for positioning purposes is limited, but it is an excellent tool for evaluating current state-of-the-art mapping functions used for geodetic positioning. Many mapping functions have been recommended in the past such as the Niell Mapping Function (NMF), Vienna Mapping Function 1 (VMF1), and the Global Mapping Function (GMF), which have been adopted by most IGS analysis centers. A new Global Pressure Temperature model (GPT2) has also been developed, which has been shown to improve upon the original atmospheric model used for the GMF. Although the mapping functions mentioned above use the same functional formulation, they vary in terms of their atmospheric source and calibration approach. A homogeneous data set of three-dimensional ray-traced delays is used to evaluate all components of the mapping functions, including their underlying functional formulation, calibration, and compression method. Additionally, an alternative representation of the VMF1 is generated using the same atmospheric source as the truth data set to evaluate the differences in ray-tracing methods and their effect on the end mapping function. The results of this investigation continue to support the use of the VMF1 as the mapping function of choice when geodetic parameters are of interest. Further support for the GPT2 and GMF as reliable back-ups when the VMF1 is not available was found due to their high consistency with the NWM-derived mapping function. Additionally, a small latitude-dependent bias in station height was found in the current mapping functions. This bias was identified to be due to the assumption of a constant radius of the earth and was largest at the poles and at the equator. Finally, an alternative version of the VMF1 is introduced, namely the UNB-VMF1 which provides users with an independent NWM-derived mapping function to support geodetic positioning.