The MARAS dataset, vegetation and soil characteristics of dryland rangelands across Patagonia

We present the MARAS (Environmental Monitoring of Arid and Semiarid Regions) dataset, which stores vegetation and soil data of 426 rangeland monitoring plots installed throughout Patagonia, a 624.500 km2 area of southern Argentina and Chile. Data for each monitoring plot includes basic climatic and landscape features, photographs, 500 point intercepts for vegetation cover, plant species list and biodiversity indexes, 50-m line-intercept transect for vegetation spatial pattern analysis, land function indexes drawn from 11 measures of soil surface characteristics and laboratory soil analysis (pH, conductivity, organic matter, N and texture). Monitoring plots were installed between 2007 and 2019, and are being reassessed at 5-year intervals (247 have been surveyed twice). The MARAS dataset provides a baseline from which to evaluate the impacts of climate change and changes in land use intensity in Patagonian ecosystems, which collectively constitute one of the world´s largest rangeland areas. This dataset will be of interest to scientists exploring key ecological questions such as biodiversity-ecosystem functioning relationships, plant-soil interactions and climatic controls on ecosystem structure and functioning.


Background & Summary
Drylands, areas where the aridity index -the ratio between precipitation and potential evapotranspiration -is below 0.65, are distributed over 100 of the world's nations, comprise about 45% of the earth's land total area 1 , and are highly prone to land degradation and desertification 2 . In addition to expanding the global area covered by drylands, ongoing climate change and the increases in aridity associated with it 3 are expected to trigger large variations in vegetation cover, biodiversity and key soil properties such as soil carbon content 4 . These ecosystem features largely influence the functioning and capacity of drylands to provide essential ecosystem services, such as biomass production and the maintenance of soil fertility, that sustain the livelihoods of more than 2 billion people living mostly in developing countries 5 .
Programs aiming to describe the status of natural resources and the onset of desertification processes 6 through change in monitoring plots include variables such as plant cover, species richness and indicators of soil fertility 7 , which are sensitive to changes in environmental conditions and human-induced land degradation 2,8,9 . Additional relevant information is provided by approaches such as the Landscape Function Analysis (LFA) 10 , a field methodology using simple and visual indicators strongly linked to fundamental physical, chemical and biological processes within ecosystems. Monitoring these variables and indicators requires the establishment of a network of field-based observations over regional scales encompassing different ecosystem types. Field data is also (a) Location, altitude, stocking density and basic climate variables for 426 monitoring plots. (b) Values for 500 sequential point intercepts, each including up to two vascular plant species or type of soil cover (bare soil, rocks, litter, cryptogams) and synthetic variables: total vegetation cover, absolute cover of each species, Richness and Shannon Wiener biodiversity index. Data for 696 plots (426 initial + 247 second + 23 third assessments) (c) Initial and final points of sequential patches and interpatches along a 50-m line intercept transect including type of patch (Woody, Herbaceous, Standing dead, Fixed litter) or interpatch (Bare soil, Gravel, Rock, Litter) and width and height of each patch and synthetic variables: mean length of patch/interpatch, mean height of patches, number of patches in 10-m. Data is provided for 680 plots (413 initial + 245 second + 22 third assessments that showed recognizable patch/interpatch structure).   www.nature.com/scientificdata www.nature.com/scientificdata/

Methods
Field work was performed by six INTA workgroups based on the cities of Rio Gallegos, Esquel, Trelew, Bariloche, Viedma and Santa Rosa (Argentina), and by one INIA team based on Punta Arenas (Chile). Field crews always included one or more trained botanists, and all of them received joint training and cross-calibration in project workshops that have taken place yearly.
Monitoring plots were chosen to represent the main landscape units of biozones 20 , which were the regional monitoring units, under their usual management: grazed by sheep, cattle or goats, or otherwise free from domestic animals but often grazed by Lama guanicoe (guanacos, a wild camelid that is the only large herbivore present in Patagonia 22 ). Mean stocking density present at each monitor was established by interviews to the land managers or to local extension officers. It is expressed in Ewe Equivalents (EE), that in south Patagonia is a 49-kg ewe that raises a 20 kg lamb, and consumes about 500 kg of herbage Dry Matter 23 . Patagonian EE is equal to New Zealand's Ewe Equivalent 24 , 1.51 UGO, the "Unidad ganadera ovina patagonica", a dry sheep equivalent used in Northern Patagonia 25 , 1.54 Dry sheep units, a similar Australian equivalent 26 . Other domestic livestock, when present, were converted as 1 cow = 6.4 EE, 1 horse = 3.63 EE and 1 goat = 1 EE 27 . Stocking density is also given in the usual cattle equivalent AU Animal Units Year −1 that refers to a 454-kg cow that consumes 3200 kg of herbage dry matter 28 where 1 AU = 0.12 EE.
Permanent monitoring plots were placed in uniform areas with a dominant vegetation type and a representative management (typically ewe paddocks in sheep extensive stations), at least 500 m away from water sources and roads. Wetlands or other azonal vegetation types were not sampled. Monitoring plots were geo-referenced using two GPS points and following a standardized method with a single protocol. A manual describing in depth the procedure used to select sites and locate the monitoring plots is available in Spanish 19 and English 21 (downloadable from Figshare 29 ). Briefly, each monitoring plot consisted in three 50-m transects (Fig. 3) positioned   Table 2. Number of MARAS monitoring plots installed in Patagonia by February 10 th 2020 (Initial assessment) in each Biozone 20 (regional monitoring unit), number of monitoring plots that were re-evaluated once (Second Assessment), or re-evaluated a second time (Third Assessment) and total number of assessments performed. (2020) 7:327 | https://doi.org/10.1038/s41597-020-00658-0 www.nature.com/scientificdata www.nature.com/scientificdata/ between permanent poles in which vegetation and soil variables were recorded and reassessed every five years. Vegetation was sampled using two 50-m line transects with 250 interception points each, and patch structure was sampled with 50-m line intercept transects for interpatches (areas that lose resources, with a minimum length of 5 cm) and patches (resource sink areas, with a minimum length of 10 cm). Soil stability, infiltration and nutrient cycling were assessed using 11 soil superficial condition indicators recorded in ten bare soil patches (located in interpatch areas) following a modified version of the LFA 30 methodology adjusted to the cover and litter values found in Patagonia (Table 3). Two 0-10 cm depth composite soil samples were obtained, one from patch and one from interpatch areas, which were analyzed in laboratory for Organic Carbon, N, texture, pH and conductivity. Initially soil samples were discarded, but those obtained in subsequent reassessments have been archived at INTA and INIA node laboratories.
The MARAS dataset was developed by the Development and Technical Assistance for Third Parties Program (PRODAT) from the National University of Austral Patagonia (UNPA) as part of the ARGENINTA grant (2015-2016), and was initially hosted in the INTA Bariloche servers in 2015. In 2019, within the Ley Ovina 2017-19 grant, the dataset was upgraded to provide a public access viewer https://maras.inta.gob.ar/portal/app/ that generates reports of change in monitoring plots, land units, biozones, provinces and departments and hosting was migrated to the central INTA servers in Buenos Aires. All field data is collected in standardized paper forms (downloadable from Figshare 29 ) that are deposited in each node. Information that remained in separate spreadsheet datasets in the initial GEF PNUD 2008-2014 stage was incorporated via data-entry contracts curated by a responsible technician in each node starting in 2015. Available information up to February 10 th 2020 was completely transferred to the dataset and is now deposited in Figshare 29 .  www.nature.com/scientificdata www.nature.com/scientificdata/ Characteristics of the field survey. MARAS sampling manuals 19,21 describe the layout of a MARAS monitoring plot, which is similar to that of the WARMS system developed in West Australia 12 . A photographic point ( Fig. 3) was fixed and a 72-m central line was laid following the main resource flux direction (wind or water flow). Three poles were fixed at 8.5 m along this line, separated perpendicularly from each other by 2.5 m, and two additional sets of three poles were set at 22 m and 72 m from the photographic point, separated from each other by 6.5 m. Photos were obtained from the photographic pole at 2-m height and from six fixed positions that show a higher detail on the photographic plot and the three lines used to evaluate vegetation and soil. All of them were uploaded to the dataset with a 2000 × 3000-pixel resolution or higher. Positions of pole 1 and 9 were registered with a GPS (Datum WGS84). Three graduated 50-m tapes were placed, two for vegetation surveys and one for patch structure sampling and LFA sampling. All the poles remain in situ in order to facilitate the reassessment of the monitoring plots.
Four basic field methods to estimate vegetation and soil properties were used at each monitoring plot: (1) Point intercept lines, with 500 sequential point intercepts 31 at 20-cm intervals along two 50 m graduated lines ( Fig. 4). At each point, a needle was set and perennial vascular plants were identified. In case of foliage superposition, the two higher plants were recorded. Plants that could not be recognized in the field were collected and voucher specimens were classified in the laboratory. All plants were identified to the level of species or pseudo species (identified to genus). Once the point intercept lines were finished, a thorough visual search was made along the central point intercept line including the photographic plot and additional plant species were listed. Starting in 2015 this was done following the point and flexible area methodology 32 . The number of individuals and approximate dimensions (length × width) of new species within 1 m at each side of the central transect were recorded. Plants that were farther away were recorded as number, dimensions and modal distance from the central transect. These data are not yet available in the dataset. A total of 636 vascular plants 33 have been identified across all MARAS monitoring plots surveyed so far. Values for cover of each species, Rocks, Bare soil, Litter, Ephemerals, Standing dead, Cryptogams were estimated as number of strikes on each category in the point intercept line divided by the total number of intercept points (500), and are expressed in %. Total vegetation cover cannot be estimated by sum of cover of each species because double vegetation strikes are accepted in the point intercept. It was estimated instead as the complement of non-vegetated point intercepts: 100-(Bare soil % + Rock % + Litter % + Standing dead% + Cryptogam %). This data is available for all the monitoring plots.
(2) In drylands, some portions of the landscape are enriched by trapping resources that are transferred from bare ground areas devoid of perennial vegetation ("interpatches") to sink areas or "patches" (typically discrete perennial vegetation and/or litter/wood accumulation patches 34 ). Interpatches are therefore nutrient-poor areas where resources (soil, propagules, nutrients) are lost, while patches are nutrient-rich areas that retain them, a pattern that takes place at different scales. In the MARAS protocol they are evaluated at site-scale (patterns generated by accumulation by shrubs or grasses) using line intercept transects that record patch structure along the soil 50-m line. The criteria used in order to define patches and interpatches is detailed in the MARAS manual 21 . The length of successive patches (minimum length 10 cm) and interpatches (minimum length 5 cm) was noted in each transect. All patches were classified as herbaceous, woody or standing dead. The height and width of each patch, measured perpendicularly at the middle point of each patch with a ruler, were also recorded. Interpatches were classified as bare soil, litter, or desert pavement (Fig. 5). This data is available for 681 plots (413 initial, 246 second and 22 third assessments) that showed patch/interpatch structure.
(3) Soil surface condition was assessed using the LFA 30 methodology along the line-intercept transect (Fig. 6). Plots were outlined over the first ten interpatches that exceeded 40 cm in length, a modification of the original www.nature.com/scientificdata www.nature.com/scientificdata/ methodology 10,30 where plots are systematically placed at intervals along a line. At each sampling point, which was 20-cm wide and had variable length, eleven indicators of soil surface conditions were visually estimated: (1) Aerial cover for rain interception, (2) Basal cover of patches, (3) Litter cover, origin and degree of incorporation, (4) Cryptogram cover, (5) Erosion features: hummocks, desert pavements, microrills or rills, (6) Deposited materials, (7) Microtopography, (8) Soil crust type and degree to which it is disturbed, (9) Surface crust resistance, (10) Slake test: time that soil aggregates retain integrity in water and (11) Texture. We adapted the range of indicators in each class of the original LFA methodology 30 to the soil and vegetation characteristics of Patagonia (Table 3). The sum of ratings obtained with these indicators for a particular monitoring plot was divided by the sum of the maximum values of the indicators (that represents 100%) and expressed as percentage. Three LFA indexes were estimated based on combinations of these indicators: Stability, Infiltration and Nutrient Cycling, which were calculated accordingly to the MARAS manual 21 . Data is provided for the 677 plots (412 initial, 243 second and 22 third assessments) that showed interpatches.
(4) Two composite 0-10 cm depth soil samples (consisting on five subsamples in patches and five in interpatches) were collected in patches (areas that accumulate resources, usually vegetated) and interpatches (areas that lose resources, usually bare soil). These samples were analyzed in the laboratory for pH 1:2.5, conductivity  www.nature.com/scientificdata www.nature.com/scientificdata/ (dS/m), Organic matter following the Walkley -Black 35 method (%), organic carbon, approximately 58% of organic matter, was estimated as organic matter/1,724 36 , N content by modified Kjeldahl procedure (%) and texture was obtained using a using Bouyoucos hydrometer. This data is available for 397 monitoring plots (295 initial + 87 second + 15 third assessments).

Dataset
The main MARAS database was built using open source software tools PHP 5.3 (www.php.net) and PostgreSQL (www.postgresql.org). The source code has a GNU General Public License (GPL), and is hosted in the main INTA servers https://maras.inta.gob.ar/app/. It is accessible with passwords, but a public viewer is available at https:// maras.inta.gob.ar/. This free browser delivers maps, photographs, basic monitoring plot information and change reports at monitoring plot and regional scales. It is currently available only in Spanish, but automatic translation to English is acceptable in most web browsers. It includes a GIS map-based walkthrough with image layers of Open Street Map (https://www.openstreetmap.org). The MARAS dataset is frequently updated through six nodes in Argentina and one in Chile, and the database program code also receives upgrades. A complete backup of the information, code and images of the dataset updated to February 2020 has been deposited in Figshare 29 . This dataset will be updated annually in Figshare to reflect data additions and updates to the code.

Data Records
The dataset presented in Figshare 29 is a copy of the dataset as it was on 10 th February 2020. By this date, a total number of 426 monitoring plots were setup in Patagonia. A single, initial assessment is presented for 156 monitoring plots, 247 of them also show data of a second assessment after a mean period of 5.9 years and 23 of them have a third, additional assessment after a mean period of 8.6 years since their set up (Table 1). A total number of 696 monitoring plot readings are in this way available in the MARAS Figshare files. Monitoring plots are distributed in six provinces of Argentina and two provinces (or regions) of Chile (Table 1), and political division of land also includes 52 departments. Biophysical classification includes 12 Biozones 20 ( Table 2) and 35 Landscape units 37 . The dataset stores additional information on the 381 farms including owners or traditional landholders, and managers. Stocking data is provided in the.csv and.xls files.
"MARAS monitoring plots csv format" is a.csv file, with a similar file available in.xls format that present the location and main vegetation, soil, land function and floristic variables for 426 MARAS monitoring plots, including 247 second assessments and 23 third assessments. Variables of this file are described in Online only Table 1. Each one of the 687 rows represents a monitoring plot reading with location, name, date of assessment, mean temperature and annual precipitation, aridity index and stocking rate data. Vegetation variables are cover, diversity, bare soil, ephemerals, standing dead and cryptogams. Vegetation structure includes patch and interpatch length, basal cover, patch height and width. Land Function Analysis includes indexes of Stability, Infiltration and Nutrient recycling. Soil laboratory results for two samples (patch and interpatch) include pH 1:2.5, conductivity (dS/m), Organic matter (%), organic carbon, Nitrogen (%), clay, silt and sand content (%). The rest of the columns are absolute cover values for all the vascular perennial species detected. Nomenclature follows "Flora del Cono Sur" 33 . Note that the sum of all absolute species cover can be >100% given that the species may be superimposed in two strata (up to two species strikes are recorded per point).
"MARAS SQL Database" is a.sql file that contains the dataset information, which includes the sequence of strikes per species in the line transects, the sequence of patches and interpatches and the individual indicator values for the LFA plots, as well as other information such as sheep or cattle numbers. This information is accessible through the MARAS program that can be restored and installed in a server via the maras-files.tgz file or browsed using other SQL software. It requires PostgreSQL.
"MARAS program and photographs" is a 12.3 Gb compressed.tgz file that upon restoration generates a complete version of the system written in PHP 5.3 (www.php.net) and PostgreSQL (www.postgresql.org) program code. In order to access the database a password must be set up in the file/exports/production/maras-www/bd/ Bd.Main.Configuraciones.Class.php. Additionally, the .tgz file will generate a photograph file folder tree with the structure: /exports/production/maras-www/imagenes/fotos_monitores/id.monitoring_plot/id.observation. The "id.monitoring plot" folders are code numbers of each MARAS monitoring plot, as listed in Online only Table 1. The "id.observation" folder is the number of assessment of the monitoring plot. Each "id.observation" folder has seven files: /diag.der; /diag.izq; /foto-grupo; /poste_central; /transecta_central; /transecta_lateral; / transecta_suelo. These files store the photographs of monitoring plots and observation teams in jpg format, with the original resolution and taken from different positions as detailed in the MARAS manual 21 . For most ecological applications it is not necessary to interrogate the SQL database or restore the compressed files, as the location of monitoring plots including climatic data, stocking density and main vegetation, patch structure, land function, soil analysis and vascular plant species cover have been summarized together in the .csv and .xls files "MARAS monitoring plots". More intensive geostatistical uses may require information on the sequence of point intercepts or patches, and this can be recovered from the SQL database directly or by means of the restored database program.
"Manual de monitores MARAS" is a Spanish (original) version of the installation manual 19 in .pdf format. "MARAS manual English version" is an English version of the MARAS installation manual 21 in .pdf format. "MARAS field worksheets_June 2020" and "Planillas de campo MARAS version 2020" are .pdf files with the worksheets used during field data recording. "Shapefile_Scidata" is a .zip file with Administrative units, biozones and land units used in shape format.
www.nature.com/scientificdata www.nature.com/scientificdata/ technical Validation Errors associated to MaRaS estimations. The MARAS system uses a single methodology and sampling effort (500 transect points along two 50-m lines, 50-m line-intercept transects and 10 LFA plots) to estimate cover, patch-interpatch structure and LFA indicators over a wide range of vegetation types, from shrublands to grasslands or semi deserts. Whether the estimations are representative of the real mean value of the monitoring plot or not depends heavily on the grain of the vegetation heterogeneity in relation to the length and resolution of sampling lines. We assessed errors associated with these estimations using the minimum number of samples equation 39 : where n = number of samples Z α/2 = False-change Type I error rate σ = Standard deviation E = Error in absolute terms.
Where σ (intra plot standard deviation) was estimated in 5 monitors randomly selected for each biozone. The standard deviation for point intercepts was obtained dividing 50-m lines in 10 subsamples with 50 points each (n = 10). In the case of interpatch length (n = 25 to 50 interpatch length measures) and LFA plots (n = 10 measures), standard deviation was estimated directly. False-change Type I error rate Z α/2 was set at 1.96 (0.05 probability). Given that at the plot scale the sampling effort is fixed (500 points, 25-50 patch-interpatch pairs and 10 LFA plots), this equation was modified to estimate the error.
Errors associated to monitoring plot estimations using Equation 2 are shown in Table 4. Lines with 500 intercept points estimate plot total vegetation cover within an error of 4.5%. The line intercept transects consisting of at least 50 patch-interpatch pairs provided interpatch length estimations within a 26 cm error, and the 10 LFA Stability index observations provided estimation errors within 4 units. However, coarse-grained vegetation types such as Austral Monte Shrublands have a higher error (5.8%, 71 cm and 4.2 units, respectively).
This error analysis could not be applied to diversity estimations at the monitoring plot scale, as it is well known that the number of species in a point intercept line does not stabilize with increasing sampling effort, as new, rare species keep appearing with added point intercepts 40 . The number of species detected with the MARAS 500-point line protocol is therefore a sub estimation of total plant biodiversity. To evaluate the precision of the protocol´s estimation of richness as a proxy of monitoring plot biodiversity we analyzed a subset of 160 monitoring plots, where species richness obtained in the MARAS monitoring plots (MR: species detected in the 500 points) was correlated with the total number of species identified by a dedicated and thorough visual search of the whole MARAS monitoring plot including the photographic plot and the area between the 50 m transects (R: species detected by thorough inspection). The linear regression of total species count with the number of species detected by point intercept (y = 1,03x + 5,77) has an R² = 0,8154 (P < 0.01) 17 . In this way, although the species count of MARAS lines underestimated total richness (and missed approximately 6 species in each monitoring plot), the high R² value and a slope close to 1 indicates that it is an effective estimator of the α-diversity of each plot. Oriental Monte Shrubland n/a n/a n/a Península Valdez Region n/a n/a n/a  Table 4. Errors associated to site means estimations using the prescribed sampling effort of MARAS (500 intercept points, 50 patch-interpatch pairs in line-intercept transects and 10 LFA plots) estimated from five monitoring plots in each of the main Biozones of Patagonia using Equation 1 (n/a= errors not assessed). (2020) 7:327 | https://doi.org/10.1038/s41597-020-00658-0 www.nature.com/scientificdata www.nature.com/scientificdata/ Minimum sample size for biozones. Inter-plot differences arise within the biozones surveyed due to climatic or soil heterogeneity, and to differences in grazing management. The minimum sampling effort is the number of monitors necessary to achieve a single estimation within an acceptable error 16 . It was estimated using Equation 1 by analyzing the standard deviation of four main variables between monitoring plots at each biozone. Sample size estimation at regional scale (Table 5) indicates that the 426 monitoring plots installed by December 2019 were enough to estimate the vegetation cover, species richness and LFA mean within the 10% error target (equivalent to ±5% cover, ±2 species of richness, and ±5 LFA units of the general mean). Not enough monitors were in place to estimate Interpatch length within 10% error (±12 cm), but it nevertheless provides an estimation of this variable with a 15% error (±19 cm). This analysis holds for most biozones, where the existing number of monitoring plots are enough to satisfy the 10% and 15% error targets, except for some heterogeneous biozones such as Subandean grasslands and Golfo San Jorge shrublands, where additional monitoring plots should be established to satisfy them. The underlying problem in some Biozones such as the Subandean grasslands may be that they encompass too much environmental heterogeneity and should be subdivided further to adequately detect changes in the future. The limitations of the analysis of vegetation structure in high cover areas through line intercept transects becomes evident in the Humid Magellan Steppe, where interpatches are few and mostly associated to walking trails of domestic stock. In this conditions Interpatch length (and its complementary variable Patch length) varies substantially and cannot be evaluated with an acceptable error. We have kept the line intercept transects mandatory for these habitats in the prospect that a recognizable patch/interpatch structure may arise with degradation processes in the future.

Repeated measures. Changes in vegetation and soil induced by climate, management or natural events
can be tracked with descriptive statistics and maps, but the system should be able to detect the statistical significance of these changes at regional and Biozone scales. The power to detect change depends on the number of monitoring plots deployed by region or Biozone, but also on the observational errors, which are inevitable and variable according to the techniques applied. In the MARAS sampling design, plots and Biozones are not sampled randomly each time, as sites are re-localized and measuring tapes set between poles permanently fixed in the monitoring plots. Even in this case, errors arise due to random factors such as misalignment of reading tapes and because observers are likely to change, and each have their own skills, criteria and biases. Monitoring plots may be treated in this experimental design as individuals in paired T-tests or repeated-measures ANOVA, and the number of monitors necessary to detect change with a given error in consecutive estimations may be estimated from the standard deviation of the differences between samples as follows 39 : where: sdiff = Standard deviation of the differences between paired samples Z α = Z-coefficient for the false-change (Type I) error rate Z β = Z-coefficient for the missed-change (Type II) error rate. MDC = Minimum detectable change size in absolute terms.
Observational errors in repeated estimations using the MARAS protocol were estimated using a set of monitoring plots that were assessed by different teams of technicians yearly from 2012 to 2015. This is not the usual 5-year assessment interval and was done in order to fulfill particular monitoring requirements of the Vanguardia Mining site 48° 20′ 30.7′′ S; 68° 8′ 50.6′′ W. These eight monitoring plots were set on desertic perennial dwarf Oriental Monte Shrubland 21 n/a n/a n/a n/a Peninsula Valdez 2 n/a n/a n/a n/a  www.nature.com/scientificdata www.nature.com/scientificdata/ shrublands of the Central Plateaus, and located in a 10-km radius with similar topographic positions and no grazing. In these conditions climatic and management differences between monitoring plots were minimal and vegetation showed small changes and, most importantly, similar interannual trends. We assumed that in these conditions most of the variability in change estimation between plots (Sdiff) was attributable to observational errors of three different evaluation teams that performed the assessments. The variability of the differences was used therefore to estimate MDC at monitoring unit scale in Equation 3. In this analysis Z α and Z β were set at 1.96 and 1.64 respectively (5% error rate).
Paired sample interannual differences (Diff) between monitoring plots in the Vanguardia site between 2012 and 2015 (Table 6) were small: +/−0.2% of vegetation cover, +/−0.4 species richness, +/−7 cm of interpatch length and +/−2.6 units of LFA stability index. Standard deviation of these difference values (Sdiff) were 3-10 times higher in all cases, indicating that annual variation in monitoring plots was lower than the observational errors due to imprecision and bias of the observers, especially in the interpatch and stability LFA index measures, which showed yearly variations of 20 cm and 9 units respectively. Minimum detectable change analysis showed that a set of 10 monitoring plots would be able to detect small changes in cover (2.2%) and species richness (1.8 species), but increased observational variability in readings of Interpatch size and Stability LFA index reduce the power of change detection to relatively high levels of variation of 23 cm and 10% LFA units in consecutive readings (Table 6).
Possible use of this data. The MARAS data allow for evaluations of ecosystem change with a precision that has not been previously possible at the regional scale in South America. The MARAS monitoring plots are thus useful for a wide variety of floristic, ecological and biogeographic studies, and have been used to describe biodiversity patterns in Patagonia 41 , to quantify the relative importance of biotic and abiotic factors as drivers of regional variations in plant productivity 42 and soil organic carbon 4 , to validate information from satellite data on the ground 43 , to assess how biodiversity modulates ecosystem responses to drought 44 and to explore how aridity and overgrazing affect the structure and functioning of drylands 45 . From the land manager perspective, the MARAS system has been used to compare the effects of different grazing systems 46,47 and in certification schemes such as sustainable management of grasslands 48 , Responsible Wool 49 or Organic Production. MARAS is also being used to monitor vegetation restoration programs set up by the gold mining industry and to assess the effects of a large-scale hydroelectric project on the vegetation of the Santa Cruz river basin http://represaspatagonia.com. ar/index.php/en/home. The data provided by MARAS are particularly helpful to interpret changes in "slow" variables 6 that have lengthy turnover times and are related to ecosystem attributes linked to climate change and variations in land use intensity, two key components of on-going global environmental change. They also provide relevant data to fulfill monitoring requirements of UN convention to Combat Desertification (UNCCD) and related conventions, the United Nations Framework Convention on Climate Change (UNFCCC) and the Convention on Biological Diversity (CBD). Therefore, the MARAS data is of interest to scientists, administrators and land managers alike, as well as to those aiming to setup ecosystem monitoring programs in drylands worldwide.

Future directions.
Most of the future effort of the MARAS project will be invested in completing the re-sampling of established plots in Patagonia and keeping the dataset and the international repository updated. There will be a decreasing emphasis on new plots over time, and the methodology will be promoted as a reliable way to certify sustainable management within extensive grazing production protocols (a way to densify the network of monitoring plots and to cover some of the costs of maintaining the MARAS system). Work will be directed towards securing the long-time funding of the program and expanding the spatial coverage of key Argentine and Chilean arid and semi-arid regions, including Puna and Monte, and integrating it with complimentary monitoring systems developed for the xeric woodlands of Chaco and Espinal 50 . This will require new monitoring nodes, probably based on agreements with other institutions from Argentina and Chile.  Table 6. Yearly paired-sample mean differences and standard deviation of the differences (Sdiff) for vegetation cover, species richness, interpatch length and LFA stability index for five monitoring plots of the Vanguardia site in Santa Cruz. MDC (Minimum detectable change) using a sample of n=10 or 5 monitors to estimate change using Equation 2.