Next Article in Journal
Stratigraphic and Paleoecological Significance of the Early/Middle Pleistocene Vertebrate Fauna of the Süttő 21 Site
Next Article in Special Issue
Range Contraction and Population Decline of the European Dupont’s Lark Population
Previous Article in Journal
Genetic Conservation and Use of Genetic Resources of 18 Mexican Pine Species
Previous Article in Special Issue
Marked Range Regression and Possible Alteration of Distribution of the Dupont’s Lark Chersophilus duponti in Tunisia: Conservation Consequences of Vanishing Alfa Grass Stipa tenacissima Steppes in North Africa
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Responses of Bird Communities to Habitat Structure along an Aridity Gradient in the Steppes North of the Sahara

by
Juan J. Oñate
1,2,*,
Francisco Suárez
1,2,†,
María Calero-Riestra
3,4,
Jorge H. Justribó
5,
Israel Hervás
1,2,
Eladio L. García de la Morena
1,6,
Álvaro Ramírez
7,
Javier Viñuela
3 and
Jesús T. García
3
1
Terrestrial Ecology Group (TEG-UAM), Department of Ecology, Universidad Autónoma de Madrid, Darwin, 2, 28049 Madrid, Spain
2
Centro de Investigación en Biodiversidad y Cambio Global (CIBC-UAM), Universidad Autónoma de Madrid, Darwin, 2, 28049 Madrid, Spain
3
Instituto de Investigación en Recursos Cinegéticos, IREC (CSIC, UCLM, JCCM), Ronda de Toledo, 12, 13071 Ciudad Real, Spain
4
Instituto Pirenaico de Ecología, Consejo Superior de Investigaciones Científicas (IPE-CSIC), 22700 Zaragoza, Spain
5
Minuartia, 08011 Barcelona, Spain
6
Biodiversity Node S.L. Sector Foresta, Tres Cantos, 28760 Madrid, Spain
7
Departamento de Biodiversidad, Ecología y Evolución, Facultad de Ciencias Biológicas, Universidad Complutense de Madrid, José Antonio Novais, 12, 28040 Madrid, Spain
*
Author to whom correspondence should be addressed.
The author has passed away.
Diversity 2023, 15(6), 737; https://doi.org/10.3390/d15060737
Submission received: 19 April 2023 / Revised: 25 May 2023 / Accepted: 1 June 2023 / Published: 2 June 2023
(This article belongs to the Special Issue Conservation of Farmland Birds)

Abstract

:
We explored the influence of habitat structure on bird density and species richness in the poorly known bird communities in the steppes of Eastern Morocco, along a 200 km long N–S gradient of increasing aridity. The birds were surveyed, and habitat structure was measured in 44 transects regularly distributed along the gradient and during the winter and spring seasons in two consecutive years. After applying a principal component analysis (PCA), five axes were identified, including one related to the latitude–altitude–soil-type gradient and another describing the development of herbaceous vegetation. Generalized linear models were used to explore the relations between bird density and species richness with PCA axes in each season, considering both the entire community and groups of granivorous, insectivorous, and mixed-diet species. More than 90% of the birds were year-round residents, with larks dominating the community in both seasons. We conclude that a distinct multifactorial response can be identified for each functional group of species. In the winter, the community is mainly affected by the structure of the habitat, while aridity (and its assumed relation to primary production) is less influential. In the spring, habitat structure continues to have the greatest explanatory power, but location along the aridity gradient becomes more relevant. These findings reveal the interaction of the negative effects of climatic and anthropogenic changes in the habitat available to these bird communities, with a greater impact expected on birds with diets that include seeds, as well as a general shift of optimal breeding conditions toward more northerly latitudes.

1. Introduction

The structure and composition of bird communities are variable in space and time [1], and the relative importance of the main governing factors continues to be debated among experts. At larger spatial scales, bird communities in arid–semiarid landscapes are thought to be relatively independent of habitat structure, and the role of climatic variables appears to be more influential [2,3,4,5]. The role of vegetation (habitat) structure as a determinant of bird abundance and species richness has been supported by approaches focusing on local to regional geographical scales [6,7,8,9,10], where the relevance of terrain topography has also been invoked [5]. Aridity has been considered the main single structuring factor of steppe bird communities at regional scales, where the uneven distribution of rainfall drives pulses of primary production distributed irregularly both in time and space [11], which may be used by many vagrant and seasonal species [12,13]. In any case, most studies on the factors governing the structure and composition of bird communities in arid–semiarid environments have usually been conducted over a single time period and the variability between seasons (wintering/breeding) or years has rarely been considered (but see [14]).
Both aridity and land cover heterogeneity may have an impact on food resources for birds, ultimately determining the species diversity and richness in a given area (e.g., [15]). This relationship will depend on the use of resources within a community [16], so that the functional relationships manifested in terms of the feeding habits of birds in that particular ecosystem must be taken into account to show how climate and habitat changes influence community dynamics beyond individual species [17].
The steppes of eastern Morocco offer optimal conditions to assess the relative importance of the candidate factors for the explanation of the structure and composition of bird communities in arid and semiarid regions. The “Hauts Plateaux marrocaines” expand over a wide area between latitudes of 34°23′ N and 32°26′ N and longitudes of 1°55′ W and 2°33′ W, comprising a north–south 200 km expanse free of geographical barriers, which could condition the distribution of species [18,19,20]. Along that expanse, rainfall decreases and temperature increases toward the south, determining a pronounced aridity gradient that covary with latitude (see the Methods section). Land cover is varied, with perennial bunch grasses, therophytes, chamaephytes, and nanophanerophytes, more or less degraded after pastoralist activities, cereal crops, and fallow land, and also bare soils with sandy, pebbly, or rocky surfaces.
In this study, we use latitude as a surrogate for aridity variation. However, latitude functions as a surrogate for many environmental factors, including vegetation structure, food availability, and altitude [21,22]. Therefore, we aim (a) to assess the relative importance of vegetation structure and latitude (aridity) as explanatory factors of the structure and composition of bird communities in eastern Morocco, which we expect to be different among bird species depending on their main diet; and (b) to explore the seasonal variation in the influence of the factors considered, which we expect to be different between the winter and the spring.

2. Materials and Methods

2.1. Study Area

The study area was located in the region of the Haut Plateaux in eastern Morocco, along a N–S strip of 200 km length, between the Col of Jerada and the south of Bouarfa (Figure 1).
Detailed meteorological data are not available for the entire area due to lack of meteorological stations, but climate data provided by the European Centre for Medium-Range Weather Forecasts (ECMWF; www.ecmwf.int/, accessed on 25 May 2023) have been used to produce monthly rainfall and average temperature series for four locations in the area, from north to south: Jerada, Aïn Bni Mathar, Tendrara, and Bouarfa (available at www.climate-data.org, accessed on 25 May 2023).
The climate is characteristically continental Mediterranean (altitude between 930 and 1400 m.a.s.l.), with cool winters and hot summers (average minimum and maximum temperatures of the coldest and warmest months, respectively, 1.6 °C and 35.1 °C in Jerada, and 1.1 °C and 36.6 °C in Bouarfa), and rainfall concentrated in autumn and spring (averages of 311 mm/yr in Jerada and 135 mm/yr in Bouarfa). The Martonne aridity index [23] varies in the study area from 14 in the northern fringe of the gradient (semiarid) to 6 in the southern extreme (arid), while the Global Aridity Index [24] ranges from arid values in the north (0.03) to hyper-arid in the south (0.008), showing a clear relationship with latitude.
Clay and silty loam soils tend to predominate toward the north of the area, although they are shallow, poor in organic matter, and have frequent capping crusts, thereby reducing water infiltration and increasing run-off and evaporation. Sandy loam soils are more frequent toward the south, where sand, pebbles, and bare rocky covers are predominant due to increasing erosion toward the transitional Saharan zones [25]. Water resources are limited even in the north, where run-off is harvested for livestock use into dayas (natural ponds) and rdirs (earth dikes) across the oueds (watercourses), and for human use into jboubs (underground cisterns). Underground water is available only in two communes (Bouarfa and Tendrara) from a few shallow aquifers (50 m deep) with low yields [26].
Vegetation varies considerably along the latitudinal gradient. Cereal crops and fallow land are dominant in the north (Jerada), subject to very large annual fluctuations, depending on autumn and early winter rains. Crops are progressively substituted toward the south by alfa grass Stipa tenacissima and white wormwood Artemisia herba-alba steppes around Aïn Bni Mathar; shrublands with Noaea mucronata, Peganum harmala, Atractylis serratuloides, and Anabasis aphylla around Tendrara, and Fredolia aretoides and Haloxylon scoparium in the most meridional limit south of Bouarfa. Occasional plantations of Atriplex nummularia and Pinus halepensis [27] reflect efforts to combat land degradation as a result of the combined effects of climate deterioration and overgrazing [28].

2.2. Sampling Methods

We established 1000 m long transects that were regularly distributed along the latitudinal gradient, separated by about six km from each other as long as flat or undulating relief was predominant and there were no human settlements in the immediate vicinity. The central point of each transect was located with a hand-held GPS (error ± 10 m) in order to be able to return to it in the different surveys. Surveys were conducted during four consecutive periods: first week of January 2005 (Winter1, n = 44 transects), first week of April 2005 (Spring1, n = 43 transects), first week of January 2006 (Winter2, n = 40 transects), and first week of May 2006 (Spring2, n = 40 transects).
Bird abundance and habitat structure were measured in each transect by experienced observers (researchers) trained on the local avifauna and vegetation and in distance estimation and vegetation data collection. Abundance of birds was estimated by the number of auditive and/or visual contacts registered along each transect with a 25 m band at either side of the observer [29]. Since we were primarily interested in studying the total diversity of the bird community, all present bird species were counted, including resident, migrant, and nomadic ones. During spring, surveys took place between one hour after dawn and until noon, avoiding the hottest central hours of the day. Winter surveys were carried out throughout the day, since patterns of detectability do not significantly vary with time in desert habitats [30]. In all cases, surveys were conducted only in good weather conditions, without precipitation, excessive wind, or sand storms. To avoid bias between seasons or sites due to potential differences in bird detection ability between survey participants, site visits were rotated between observers.
Habitat structure was measured in six circles of five-meter radius (sampling points), separated by 200 m along each 1000 m transect. We recorded the following variables at each sampling point: geographical location (UTM coordinates latitude and longitude); elevation (in meters above sea level); estimation of average slope of the terrain (three levels: 0 = flat; 1 = low slope; 2 = moderate slope); total cover of bare ground; total extent of rocks (>10 cm), pebbles (<10 cm), sand, and lime cover; cover of dry cereal crops; and total cover of vegetation. The latter distinguished cover of herbaceous species (including Peganum harmala), shrubs, and alfa grass was measured at three different heights (1, 20, and 40 cm). Vegetation cover was estimated after [31], using the classes <1%; 1–3%; 3–5%; 5–10%; 10–20%; 20–30%; 30–40%; 40–50%; and >50%. Alfa clumps that lacked shoots were also included in vegetation cover and estimated at the same heights. Maximum and modal vegetation heights were also measured.

2.3. Data Analyses

The basic sample units were 1000 m-length transects. For birds, we used the number of different species and density of each species, calculated as the sum of contacts (within the 50 m band) of each species multiplied by the sampled area (5 ha per transect) and multiplied by 2 to obtain the number of individuals per 10 ha. Data on habitat structure from the six sampling points in each transect were averaged and used as representative for the transect.
A matrix of counts of each species in each season (winter and spring) was created, pooling the counts of the two replicates (surveys) of each season. We pooled counts because there were no significant differences in richness or density between the surveys of each season (ANOVA; period (nested in season): p > 0.18 and p > 0.41 for richness and density, respectively) and because the year effect was not the focus of our study. Based on this matrix, different diversity measures were calculated for each season, considering counts of all the species detected. We calculated richness (S) as the total number of species within the community; Shannon’s diversity index (H′ = −∑pi ln(pi)), where pi is the proportion of individuals belonging to species i; and Simpson’s diversity (D1 = 1 − ∑pi2), Simpson’s dominance (D2 = 1/∑pi2), and Simpson’s evenness index (E = D2/S) (formulas from [32,33,34]. We also calculated species turnover between the two seasons following the framework proposed by [35], where the turnover and nestedness components of beta diversity are decomposed. The spatial turnover component used here is calculated as a Simpson-based dissimilarity index (βSIM) [36]: min(b,c)/a + min(b,c), where a is the number of species common to both seasons, b is the number of species exclusive to the focal season, and c is the number of species exclusive to the other season.
Since the habitat variables measured in the transects were highly correlated to one another, we used a Principal Components Analysis (PCA) to derive a set of uncorrelated, synthetic components. In addition to the variables on soil and vegetation cover, we included other variables in the PCA that can covary with them along the gradient, such as latitude, average slope, and elevation, to obtain a reduced number of (uncorrelated) factors that summarize soil, habitat, topography, and latitude. We then used the synthetic components derived from the PCA as independent variables in subsequent analyses to eliminate multicollinearity while retaining the variation of the environmental variables for the models.
Besides the whole bird community, we also considered functional groups of species (Table 1) that could be influenced differently by environmental features. While we recognize that bird species are rarely exclusive in their use of resources and that their food requirements vary seasonally, we categorized all bird species registered during surveys on the basis of their primary food (diet and feeding habits) as granivores, insectivores, or mixed diet. Allocation of species to functional groups (Table 1) was based on published studies [37] and the authors’ own observations in the study area. For each functional group, bird density (individuals/10 ha) and species richness were derived considering counts of all species in the group.
We used a general linear model (GLM) with a nested design to test whether each of the principal components (PCs) varied between seasons (spring and winter) and between the two surveys within each season (nested effect of survey within season). Generalized linear models (GLZ) were used to test the effect of season (factor) on richness and density of the entire community and of each of the functional groups (dependent variables). We also used GLZ to explore, separately for each season, the relationships between richness and density (dependent variables) and the components derived from the PCA (predictors). Since we have two replicates of each season, we included survey as fixed factor (two levels) in analyses. We used a Poisson error distribution and a log link function (richness) or quasi-Poisson distribution (density) to correct overdispersion [38].
Data were analyzed using Statistica 8.0 [39] and the package lme4 [40] in the statistical environment R version (4.2.1) [41].
Average values and their variations along the text have been quoted as average ± standard error.
Table 1. List of species and their residency status in the study area following [42] (R = resident; S = spring visitor; W = winter visitor; N = nomadic), number of birds (ni) and proportion of the total (pi) recorded in the spring and winter surveys. The functional groups (I = insectivores; G= granivores; M = mixed diet) and the diversity metrics used in this study are also indicated.
Table 1. List of species and their residency status in the study area following [42] (R = resident; S = spring visitor; W = winter visitor; N = nomadic), number of birds (ni) and proportion of the total (pi) recorded in the spring and winter surveys. The functional groups (I = insectivores; G= granivores; M = mixed diet) and the diversity metrics used in this study are also indicated.
SpringWinter
SpeciesStatusDietnipinipi
Cursorius cursorNI30.0000.00
Pterocles orientalisRG60.0100.00
Alaemon alaudipesNI30.00120.01
Chersophilus dupontiRI20.0010.00
Ammomanes cincturaRG360.05680.05
Ammomanes desertiRM10.0000.00
Ramphocoris clotbeyNG140.02570.04
Melanocorhypha calandraRG1080.161810.14
Calandrella brachydactylaSM1060.1600.00
Calandrella rufescensRM1350.201850.14
Eremophila bilophaRG1070.162580.20
Galerida cristataRM350.05570.04
Galerida theklaeRM130.02390.03
Alauda arvensisRM120.02270.02
Anthus pratensisWI00.0030.00
Motacilla albaSI20.0000.00
Oenenthe desertiSI50.0100.00
Oenanthe hispanicaSI30.0000.00
Oenanthe leucuraRI10.0000.00
Oenanthe moestaRI180.0380.01
Oenanthe oenantheSI20.0000.00
Saxicola rubetraWI20.0000.00
Lanius excubitorWI00.0010.00
Lanius senatorSI10.0000.00
Passer domesticusWG00.0050.00
Passer hispaniolensisRM400.0600.00
Bucanetes githagineusNG20.004100.31
Emberiza calandraRG20.0020.00
N° birds 659 1314
N° transects 83 84
Mean density/transect ± s.d. 7.67 ± 11.3 31.28 ± 46.4
Richness (S) 25 16
Mean richness/transect ± s.d. 1.92 ± 1.61 1.78 ± 1.45
Shannon’s diversity (H’) 1.00 0.85
Simpson’s diversity (D1) 0.86 0.81
Simpson’s dominance (D2) 7.53 5.47
Simpson’s evenness (E) 7.49 6.39
Simpson-based dissimilarity (βSIM) 0.19

3. Results

A total of 1973 contacts of 28 different species were registered during all the surveys (Table 1). Some species were very scarcely detected, but we preferred to consider all of them for the subsequent analyses given our primary interest in studying the total diversity of the bird community.
The whole community was dominated by larks in both seasons, accounting for 87% of all contacts in the spring (48% of all species) and 67% in the winter (62.5% of all species). The numerical importance of winter migrant species in the wintering community was low (0.7% of the winter birds; 19% of the species), being somewhat higher than that of spring migrants in the breeding community (18% and 24% of spring birds and species, respectively). Nomadic birds were more abundant in the winter (36%; 19% of all species) than in the spring (3%; 16% of all species). More than half (56%) of the contacted species in the spring were year-round residents (78% of detected birds), while residents accounted for 62.5% of the species in the winter (63% of detected birds). While the mean bird density per transect was clearly larger in the winter, the mean richness remained much more stable between the seasons.
All diversity metrics consistently indicated greater richness, average diversity, and evenness in the spring compared to the winter (Table 1), and the temporal dynamic in the community structure showed a moderate species turnover rate between the two seasons.
The PCA considering soil, habitat, topography, and latitude variables yielded five components with eigenvalues greater than 1, accounting for 73.5% of the total variance (Table 2). The first component (PC1) alone explained 22.2% of the total variance, and the correlation coefficient values show that it has a strong positive correlation with alfa grass cover and vegetation height. PC1 may thus be regarded as a gradient of development of alfa grass, given that this species reaches the maximum height measured among all plants detected in the surveys when fully grown and well preserved. The second component (PC2; 17% of variation) was positively correlated with latitude, silty loam soils, and cereal crops, and negatively correlated with elevation and sand cover. Therefore, this component may be regarded as a geographic (elevation–latitude) gradient in soil texture that determines land suitability for agriculture, from silty loam arable lands at low elevations in the north, to the southern sandy soils in the highest plateau that are less suitable for cultivation.
The third and fourth components (14.2% and 11.6% of variation, respectively) were strongly correlated with lower coverage of shrubs and larger development of herbaceous cover, respectively. The fifth component was mainly positively correlated with the slope, and the cover of rocks and pebbles, in such a way that its highest values correspond to places with stony soils on sloping terrains.
We found a “season” effect for PC4, with significantly higher values in the spring than in the winter, and a “season” (spring > winter) and “survey” effect for PC5 (larger values in the second spring survey as compared to the first one). No significant variations between seasons or between surveys within each season were found for PC1, PC2, or PC3 (Table 3).
The seasonal variation of total species richness was not statistically significant (GLZ: Estimate = 0.044, Wald = 0.62, p = 0.43). Insectivores and mixed-diet groups showed higher richness in the spring than in the winter, although neither reached the significance level (p-values of 0.07 and 0.089, respectively). In contrast, granivores showed the opposite trend, with greater (non-significant, p-value = 0.09) richness in the winter than in the spring. The total bird density per transect was significantly higher in the winter than in the spring (GLZ: Estimate = –0.338, Wald = 7.81, p = 0.0052). The effect of season on bird density was mainly conditioned by the group of granivorous species (GLZ: Estimate = –0.629, Wald = 13.94, p = 0.0001), since the density of the mixed-diet group remained stable between seasons (p = 0.73) and the density of insectivores was almost significantly higher in the spring than in the winter (GLZ: Estimate = 0.265, Wald = 3.51, p = 0.060).
Generalized linear models revealed no significant influence of the principal component 1 (alfa grass cover and vegetation height) on species richness or bird density in any of the two seasons (Table 4 and Table 5). PC5 (stony soils on sloping terrain) only favored granivores density in the spring, and PC3 (low coverage of shrubs) negatively influenced the bird density of the mixed-diet group in the winter. The most consistent influences were identified for PC2 and PC4. In the spring, species richness and bird density, both for all species as a whole and for the mixed-diet group, significantly increased towards northern loam clay arable lands at lower altitude (PC2) and in patches of higher herbaceous cover (PC4). Herbaceous cover also favored species richness and bird density of the granivorous group and bird density of the insectivorous one. In the winter, the extent of herbaceous cover significantly favored species richness and bird density of all species as a whole and of those of the mixed-diet group, as well as bird density of the granivores group. The influence of PC2 was more limited, favoring species richness and bird density of the mixed-diet group and negatively influencing the bird density of the granivorous group. Significant differences were detected between spring surveys for total and mixed-diet species richness (Table 4) and for total, granivorous, and mixed-diet bird densities (Table 5).

4. Discussion

A total of 28 species were counted in the transects. Residents comprised the greatest proportion of detected species (50%) and of individuals counted (68%). Despite the higher detectability due to courtship and breeding behavior [43], the number of individuals counted was higher in the winter than in the spring, which is a reflection of the gregarious nature in the winter of part of the species with the largest bird numbers. Four of the recorded species were nomadic (14%, including Alaemon alaudipes, Bucanetes githagineus, and Cursorius cursor). Spring (24% of the species) and wintering migrants (19%) had a reduced importance in comparison to other Mediterranean environments with greater structural complexity, in particular those forests and shrubs with abundant fruits [44], but it was similar to what happens in the desert areas of Tunisia [45] or in semiarid Spain [14]. Our results corroborate previous knowledge regarding the dominance of larks in the composition of breeding bird communities in the Northern Maghreb, both in species and bird numbers [45,46,47]. The community species turnover was moderate (Simpson-based dissimilarity index, βSIM = 0.19), and there was a high coincidence between seasons in the species with larger quantitative importance, in particular lark species, such as Melanocorhypha calandra, Calandrella rufescens, and Eremophila bilopha. Higher numbers of species found during the spring may reflect the passage or arrival of migrants, although the bird count was higher in the winter, mainly due to the detected flocks and nomadic species.
PCA allowed a comprehensive interpretation of the main factors structuring the bird community, with PC1 capturing the potential influence of alfa grass coverage and development; PC2 integrating the combined role of increasing latitude (less aridity), silty loam soils, and cover of cereal crops in opposition to lower latitude (greater aridity), larger sandy soil cover, and higher elevation; PC3 and PC4 reflecting lower and larger coverage, respectively, of shrub and herbaceous vegetation; and PC5 capturing the dominance of stony soils (rocks and pebbles) on sloping terrains.
The consistency of this interpretation between surveys is large, with no significant variation for the first three PCs, and expected higher values of herbaceous cover development in the spring (PC4). The interpretation of the seasonal variation of PC5 is more doubtful, given that slopes or pebble cover should not vary appreciably between seasons or years. Most probably, the observed variation is an artifact derived from location errors (despite the use of GPS) of the transects between surveys, linked to the inaccurate selection of the exact point and the orientation of the line of progression from that point. However, it could not be discarded that pebble cover could vary between years in relation to seasonal/interannual effects of high, hot, and dry winds, giving rise to significant sandstorms, particularly frequent in the dry summer season in the area [48].
Although the seasonal variation of total richness was not statistically significant, the total density per transect was significantly higher in the winter than in the spring. These larger densities were not observed by [45], who reported larger densities in the spring. However, the communities they studied were dominated by insectivorous species, while in our case, the seasonal differences were mainly due to the higher bird density of granivorous species in the winter. Furthermore, the study in Tunisia focused on a single southern study area at the edge of the desert (e.g., ecologically equivalent to the southern tip of our latitudinal gradient, see Figure 1), while we sampled a considerable latitudinal range in eastern Morocco. In fact, we found a positive effect of PC2 on the overall bird density in the spring (Table 5); that is, bird density increased toward northern areas. Our community seems to be relatively stable in its specific composition (as expected, given the harsh conditions in both periods of the year), but with increasing numbers of wintering bird populations coming from more northerly latitudes.
Habitat (vegetation) structure has a significant influence on the composition and structure of bird communities along the 200 km extension studied in the Moroccan High Plateaux. This general result coincides with other studies at local or regional scales [6,7,8,9,10,43]. Aridity (an inverse correlate of latitude in our study, PC2) was less influential [1,3,4,5].
The extent and development of herbaceous vegetation (PC 4, Table 2) resulted in our study to the main single factor related to higher total species richness and bird density, particularly in the spring, but also in the winter for some functional groups (Table 4 and Table 5). Species richness and bird density of the three functional groups considered were all positively related to this factor in the spring, and also richness and density of the mixed-diet group and density of the granivorous group showed this relation in the winter. In contrast, neither the richness nor density of the insectivorous group showed a relationship to this factor in the winter. Furthermore, bird density of the mixed-diet group appeared to be positively related to shrub cover in the winter.
Two possible explanations may be considered for the role of habitat structure as a factor shaping the studied bird communities. First, although our study expanded along a 200 km N–S strip, it did not reach the desertic areas beyond Bouarfa to the south (<100 mm/yr of highly irregular precipitation) [49]. Therefore, primary productivity pulses in the studied area (with precipitation of 130–319 mm/yr) are more predictable than if we had considered a broader climate gradient, which could favor resident species at the expense of vagrant ones. This fact seems to be reflected in the percentage of nomadic species, which in our case, is less than half that of the Australian deserts (14% in our case; 46% in Australians [13]. In addition, the importance of spring or winter species is also low in our communities, with a quite moderate species turnover between seasons.
A second explanation may be related to the higher spatial variability of habitat structures in the studied area in comparison to the relatively homogeneous environments in which the prevalence of climate-related variables was stated. Although the study area was only 200 km in length, it encompassed a marked diversity of soil types (from silty loam in the north to sandy, pebbly, and rocky soils in the south). Further, human activities linked to pastoralism and cultivation have altered the original alfa grass and white wormwood formations into a variety of degradation facies [25], with heterogeneous composition and coverage. These influences are reflected in a high habitat patchiness at different spatial scales (plot and landscape, mainly), originating a set of environments with contrasting structural characteristics (pers. Obs). This means that the species can use a wide range of structural conditions of the habitat in the study area, producing a clearer differentiation between them than if it were a more homogeneous landscape. Our results are, therefore, consistent with previous findings in suggesting the importance of scale and environmental heterogeneity in the relative importance of resident vs. vagrant species in ecological communities. The larger the scale, the greater the chance of including a suitable habitat in sufficient quantity to support persistent or resident populations (e.g., [50]). In addition, landscapes showing low environmental heterogeneity tend to support communities with low temporal turnover and a higher proportion of resident species—probably because more heterogeneous landscapes are more spatially segmented, effectively reducing the area and resources available per habitat type to support viable resident populations [51,52].
PC2 was the second-most important factor, significantly explaining bird density and species richness in these communities. As mentioned above, this factor reflects the different linked geographical gradients existing in the area, from less arid northern arable lands to arid and sandy lowlands in the south. Total species richness and bird density were related in the spring to the less arid northern croplands, and the mixed-diet group showed higher richness and density in these areas also in the winter, coinciding with lower densities of the granivorous group. The influence of soil type was also reflected in the inverse relation shown by the bird density of the mixed-diet group to the dominance of pebbly and rocky covers and sloping terrains. The importance of soil type in the composition of bird communities or the selection of habitats in semi-desert environments has been frequently highlighted (for North Africa, see [37,42]; see also [53,54]. However, in our case, the typology of soil types also varies in an accentuated way with the same north–south geographic gradient, potentially confounding its effects with those of aridity (latitude) or the abundance of crops. In North American shrub steppes, it has been suggested that the overlap between species of the factors affecting their distribution and density is very high, making it difficult to differentiate their habitat selection [55,56]. Our results suggest a multiple and differential response between species to the different parameters that determine habitat structure, including the type of soil, as expected, considering the existence of interspecific differences in their biology.
Our results offer a valuable first insight on the primary factors influencing the poorly studied bird communities in the High Plateaux of eastern Morocco, especially considering the accelerated process of degradation suffered by their natural environments. The combined effect of climatic change and land use intensification is driving widespread desertification in the arid and semi-arid environments of North Africa [57]. In the High Plateaux of eastern Morocco, continued overgrazing and the substitution of grassland by itinerant cropland have reduced the alfa grass and white wormwood formations from 2 million hectares in the 1970s to less than 420,000 hectares in less than 50 years [28]. Consequently, moderately to severely degraded and very severely degraded vegetation classes have become dominant (1,800,000 ha and 650,000 ha, respectively) [28]. In the inventories of [58], one of the most abundant species was the shrub warbler Scotocerca inquieta, which was not contacted in our transects and was only observed once throughout the study. We could not establish an historical comparison with the work of these authors since they only sampled alfa grass stands. However, these formations have suffered severe degradation, particularly around Aïn Bni Mathar [59], resulting in the disappearance of one of the bird species more strongly associated with this specific vegetation, the Dupont’s lark [60], also conditioning corvids’ and raptors’ abundance [61]. However, other passerines must also have been affected. In the face of these accelerated changes, detailed monitoring of bird communities in these areas is urgently needed, especially considering that quantitative records are almost non-existent.
Four essential conclusions are derived from our study. First, the explanatory variables considered affect the species richness and density of individuals in the studied bird communities in a differential way, making it possible to establish multifactorial behavior for each functional group of species that includes geographical, soil, and vegetation structure factors. Second, the total density of the bird community is strongly affected by the development of herbaceous vegetation, independently of the prevailing aridity gradient N–S, probably reflecting the erratic rainfall in these semi-desert areas, as well as the vagrant behavior of many bird species in these communities. Third, species richness along the studied N–S 2000 km gradient was affected in the winter by factors related to habitat structure and not by variables related to the aridity gradient. Fourthly, the rapid changes in habitat structure that seem to be occurring in the area, associated with anthropic activities and climate degradation, suggest that they could have important effects on these singular bird communities. Our data reveal that ongoing climate change may exert negative impacts on the community, especially on birds with mixed or granivorous diets and during the breeding season. Further northward shifts in environmental favorability are predicted for the study species due to the interaction between climate change and human-induced changes in the habitat. However, to better understand the combined effects of habitat and climate changes as well as determine gains and losses in favorability across the range, more long-term studies are necessary.

Author Contributions

Conceptualization, F.S., J.T.G. and J.J.O.; methodology F.S., J.T.G. and J.J.O.; data collection, all authors; formal analysis, J.T.G.; writing—original draft preparation, J.J.O.; writing—review and editing, all authors. All authors have read and agreed to the published version of the manuscript.

Funding

The authors did not receive any specific financial support for the research and used their own funds and resources (e.g., vehicles) in some expeditions, as well as other resources available to the IREC and the UAM.

Institutional Review Board Statement

Ethical review and approval were waived for this study, since birds were not manipulated and only visually contacted.

Data Availability Statement

The data presented in this study are available upon request from the corresponding author. The data are not publicly available due to privacy restrictions as they are being used in further studies for upcoming publications.

Acknowledgments

We will forever be grateful to the late Quico Suárez, who started and led this research. We wish to thank the Moroccan government (Le Secretaire General du Departement des Eaux et Forets) for fieldwork permits. Meg A. Cefalo kindly provided us with warm accommodation during our stays in Aïn Beni Mathar. We thank V. Garza for his help during fieldtrips. The Instituto de Investigación en Recursos Cinegéticos, IREC (CSIC-UCLM-JCCM), provided vehicles for the expeditions. This paper contributes to project REMEDINAL TE-CM(P2018/EMT4338).

Conflicts of Interest

The authors have no relevant financial or non-financial interest to disclose.

References

  1. Wiens, J.A. The Ecology of Bird Communities; Cambridge University Press: Cambridge, UK, 1989. [Google Scholar]
  2. Wiens, J.A. Habitat heterogeneity and avian community structure in North American Grasslands. Am. Midl. Nat. 1974, 91, 195–213. [Google Scholar] [CrossRef]
  3. Rotenberry, J.T. Components of avian diversity along a multifactorial climatic gradient. Ecology 1978, 59, 693–699. [Google Scholar] [CrossRef]
  4. Naranjo, L.G.; Raitt, R.J. Breeding bird distribution in Chihuahuan desert habitats. Southwest. Nat. 1993, 38, 43–51. [Google Scholar] [CrossRef]
  5. Kaboli, M.; Guillaumet, A.; Prodon, R. Avifaunal gradients in two arid zones of central Iran in relation to vegetation, climate and topography. J. Biogeogr. 2006, 33, 133–144. [Google Scholar] [CrossRef]
  6. Cody, M.L. Mulga bird communities. I. Species composition and predictability across Australia. Aust. J. Ecol. 1974, 19, 206–219. [Google Scholar] [CrossRef]
  7. Tomoff, C.S. Avian species diversity in desert scrub. Ecology 1974, 55, 396–403. [Google Scholar] [CrossRef]
  8. Rotenberry, J.T. The role of habitat in avian community composition: Physiognomy or floristics? Oecologia 1985, 67, 213–217. [Google Scholar] [CrossRef]
  9. Recher, H.F.; Davis, W.E., Jr. Foraging ecology of a mulga bird community. Wildl. Res. 1997, 24, 27–43. [Google Scholar] [CrossRef]
  10. Pavel, C.R.; Nano, C.E.M. Bird assemblages of arid Australia: Vegetation patterns have a greater effect than disturbance and resource pulses. J. Arid. Environ. 2009, 73, 634–642. [Google Scholar] [CrossRef]
  11. Le Houérou, H.N.; Bingham, R.L.; Skerbek, W. Relationship between the variability of primary production and the variability of annual precipitation in world arid lands. J. Arid. Environ. 1988, 15, 1–18. [Google Scholar] [CrossRef]
  12. Dean, W.R.J.; Hockey, P.A.R. An ecological perspective of lark (Alaudidae) distribution and diversity in the southwest-arid zone of Africa. Ostrich 1989, 60, 27–34. [Google Scholar] [CrossRef]
  13. Dean, W.R.J. Nomadic Desert Birds; Springer: Berlin/Heidelberg, Germany, 2004. [Google Scholar]
  14. Tellería, J.L.; Santos, T.; Suárez, F. Bird communities of the Iberian shrubsteppes: Seasonality and structure along a climatic gradient. Holartic Ecol. 1988, 11, 171–177. [Google Scholar]
  15. Lorenzón, R.E.; Beltzer, A.H.; Olguin, P.F.; Ronchi-Virgolini, A.L. Habitat heterogeneity drives bird species richness, nestedness and habitat selection by individual species in fluvial wetlands of the Paraná River, Argentina. Austral Ecol. 2016, 41, 829–841. [Google Scholar] [CrossRef]
  16. Ricklefs, R.E. Evolutionary diversification, coevolution between populations and their antagonists, and the filling of niche space. Proc. Natl. Acad. Sci. USA 2010, 4, 1265–1272. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. González-Salazar, C.; Martínez-Meyer, E.; López-Santiago, G. A hierarchical classification of trophic guilds for North American birds and mammals. Rev. Mex. Biodivers. 2014, 85, 931–941. [Google Scholar] [CrossRef]
  18. García, J.T.; Suárez, F.; Garza, V.; Calero-Riestra, M.; Hernández, J.; Pérez-Tris, J. Genetic and phenotypic variation among geographically isolated populations of the globally threatened Dupont’s lark Chersophilus duponti. Mol. Phylogenet. Evol. 2008, 46, 237–251. [Google Scholar] [CrossRef] [PubMed]
  19. Guillaumet, A.; Crochet, P.A.; Godelle, B. Phenotipic variation in Galerida larks in Morrocco: The role of history and natural selection. Mol. Ecol. 2004, 14, 3809–3821. [Google Scholar] [CrossRef]
  20. Guillaumet, A.; Pons, J.M.; Godelle, B.; Crochet, P.A. History of the Crested Lark in the Mediterranean region as reveales by mtDNA sequences and morphology. Mol. Phylogenet. Evol. 2006, 39, 645–656. [Google Scholar] [CrossRef]
  21. Meiri, S.; Dayan, T. On the validity of Bergmann’s rule. J. Biogeogr. 2003, 30, 331–351. [Google Scholar] [CrossRef]
  22. Rahbek, C. The role of spatial scale and the perception of large-scale species-richness patterns. Ecol. Lett. 2005, 8, 224–239. [Google Scholar] [CrossRef]
  23. De Martonne, E. Une nouvelle fonction climatologique: L’ indice d’aridité. Météorologie 1926, 2, 449–458. [Google Scholar]
  24. Zomer, R.J.; Xu, J.; Trabucco, A. Version 3 of the Global Aridity Index and Potential Evapotranspiration Database. Sci. Data 2002, 9, 409. [Google Scholar] [CrossRef]
  25. Bechchari, A.; El Aich, A.; Mahyou, H.; Baghdad, B.; Bendaou, M. Analyse de l’évolution du système pastoral du Maroc oriental. Rev. Élev. Méd. Vét. Pays Trop. 2015, 67, 151–162. [Google Scholar] [CrossRef] [Green Version]
  26. El Harradji, A. Aménagement, érosion et désertification sur les Hauts-Plateaux du Maroc oriental. Méditerranée 1997, 86, 15–23. [Google Scholar] [CrossRef]
  27. Ben El Mostafa, S.; Haloui, B.; Berrichi, A. Contribution à l’étude de la végétation steppique du Maroc oriental: Transect Jerrada—Figuig. Acta Bot. Malacit. 2001, 6, 295–301. [Google Scholar] [CrossRef]
  28. INRA-ONUDI. Etude sur la Situation de Référence au Niveau des Hauts Plateaux du Maroc Oriental. Rapport Final. Projet de Lutte Participative Contre la Désertification et de Réduction de la Pauvreté dans les Écosystèmes Arides et Semi Arides des Hauts Plateaux du Maroc Oriental; Centre Régional de La Recherche Agronomique d’Oujda: Oujda, Morocco, 2012. [Google Scholar]
  29. Bibby, C.J.; Burgess, N.D.; Hill, D. A. Bird Census Techniques; Academic Press: London, UK, 1992. [Google Scholar]
  30. Weathers, W.W.; Mayhew, W.W. Time of day and desert bird censuses. West. Birds 1981, 12, 157–172. [Google Scholar]
  31. Prodon, R.; Lebreton, J.D. Breeding avifauna of a Mediterranean succession: The holm oak and cork oak series in eastern Pyrenees. 1. Analysis and modelling of the structure gradient. Oikos 1981, 37, 21–38. [Google Scholar] [CrossRef]
  32. Morris, E.K.; Caruso, T.; Buscot, F.; Fischer, M.; Hancock, C.; Maier, T.S.; Meiners, T.; Müller, C.; Obermaier, E.; Prati, D.; et al. Choosing and using diversity indices: Insights for ecological applications from the German Biodiversity Exploratories. Ecol. Evol. 2014, 4, 3514–3524. [Google Scholar] [CrossRef] [Green Version]
  33. Shannon, C. A mathematical theory of communication. Bell Syst. Technol. J. 1948, 27, 379–423. [Google Scholar] [CrossRef] [Green Version]
  34. Simpson, E.H. Measurement of diversity. Nature 1949, 163, 688. [Google Scholar] [CrossRef]
  35. Baselga, A. Partitioning the turnover and nestedness components of beta diversity. Global Ecol. Biogeogr. Let. 2010, 19, 134–143. [Google Scholar] [CrossRef]
  36. Lennon, J.J.; Koleff, P.; Greenwood, J.J.D.; Gaston, K.J. The geographical structure of British bird distributions: Diversity, spatial turnover and scale. J. Anim. Ecol. 2001, 70, 966–979. [Google Scholar] [CrossRef] [Green Version]
  37. De Juana, E.; Suárez, F.; Ryan, P. Family Alaudidae (Larks). In Handbook of the Birds of the World, Cotingas to Pipits and Wagtails; Del Hoyo, J., Elliott, A., Christie, D., Eds.; Lynx Edicions: Barcelona, Spain, 2005; Volume 9, pp. 496–601. [Google Scholar]
  38. Crawley, M.J. The R Book; John Wiley & Sons: Chichester, UK, 2012. [Google Scholar] [CrossRef]
  39. StatSoft, Inc. STATISTICA (Data Analysis Software System), Version 8.0. 2007. Available online: www.statsoft.com (accessed on 25 May 2023).
  40. Bates, D.; Mächler, M.; Bolker, B.; Walker, S. Fitting linear mixed-effects models Usinglme4. J. Stat. Softw. 2015, 67, 1–48. [Google Scholar] [CrossRef]
  41. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2022; Available online: https://www.R-project.org/ (accessed on 25 May 2023).
  42. Thévenot, M.; Vernon, R.; Bergier, P. The Birds of Morocco; BOU Checklist nº 20; British Ornithologist’ Union-British Ornithologists’ Club: Tring, UK, 2003. [Google Scholar]
  43. Heezik, Y.V.; Seddon, P.J. Effects of season and habitat on bird abundance and diversity in a steppe desert, northern Saudi Arabia. J. Arid. Environ. 1999, 43, 310–317. [Google Scholar] [CrossRef]
  44. Santos, T.; Tellería, J.L. Patrones generales de la distribución invernal de passeriformes en la Península Ibérica. Ardeola 1985, 32, 17–30. [Google Scholar]
  45. Heatwole, H.; Muir, R. Population densities, biomass and trophic relations of birds in the pre Saharan steppe of Tunisia. J. Arid. Environ. 1982, 5, 145–167. [Google Scholar] [CrossRef]
  46. Blondel, J. Donnés écologiques sur l’avifaune des Monts des Ksours (Sahara septentrional). Terre Vie 1962, 16, 209–251. [Google Scholar]
  47. Suárez, F. Introducción al estudio de las comunidades de aves reproductoras de los espartales norteafricanos. Bol. Estac. Cent. Ecol. 1985, 28, 29–34. [Google Scholar]
  48. Mahyou, H.; Tychon, B.; Balaghi, R.; Louhaichi, M.; Mimouni, J. A knowledge-based approach for mapping land degradation in the arid rangelands of North Africa. Land Degrad. Dev. 2016, 27, 1574–1585. [Google Scholar] [CrossRef]
  49. Le Houérou, H.N. Classification écoclimatique des zones arides (s.l.) de l’Afrique du Nord. Ecol. Mediterr. 1989, 15, 95–146. [Google Scholar] [CrossRef]
  50. Jenkins, M.F.; White, E.P.; Hurlbert, A.H. The proportion of core species in a community varies with spatial scale and environmental heterogeneity. PeerJ 2018, 6, e6019. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Stegen, J.C.; Freestone, A.L.; Crist, T.O.; Anderson, M.J.; Chase, J.M.; Comita, L.S. Vellend M. Stochastic and deterministic drivers of spatial and temporal turnover in breeding bird communities. Glob. Ecol. Biogeogr. 2013, 22, 202–212. [Google Scholar] [CrossRef]
  52. Gaston, K.J.; Davies, R.G.; Orme, C.; David, L.; Olson, V.A.; Thomas, G.H.; Ding, T.S.; Rasmussen, P.C.; Lennon, J.J.; Bennett, P.M.; et al. Spatial turnover in the global avifauna. Proc. R. Soc. B 2007, 274, 1567–1574. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Kantrud, H.A.; Kologiski, R.L. Avian associations of the northern Great plains grasslands. J. Biogeogr. 1983, 10, 331–350. [Google Scholar] [CrossRef]
  54. Van der Haegen, W.M.; Dobler, F.C.; Price, D.J. Shrubsteppe bird reponse to habitat and landscape variables in Eastern Washington, USA. Conserv. Biol. 2000, 14, 1145–1160. [Google Scholar] [CrossRef] [Green Version]
  55. Rotenberry, J.T.; Wiens, J.A. Habitat structure, patchiness, and avian communities in North American steppe vegetation: A multivariate analysis. Ecology 1980, 61, 1228–1250. [Google Scholar] [CrossRef]
  56. Wiens, J.A.; Rotenberry, J.T. Habitat associations and community structure of birds in shrubsteppe environments. Ecol. Monogr. 1981, 51, 21–41. [Google Scholar] [CrossRef]
  57. Le Houérou, H.N. Biogeography of the arid steppeland north of the Sahara. J. Arid. Environ. 2001, 48, 103–128. [Google Scholar] [CrossRef]
  58. Suárez, F.; Fernández, A.; de Lope, M.J. Note sur les effects de l’aridité sur la structure et la composition des communautes de Passeriformes des hauts-plateaux a alpha (Stipa tennacissima) au Maroc. Bull. Inst. Sci. Rabat 1986, 10, 185–192. [Google Scholar]
  59. Bechchari, A.; El Aich, A.; Mahyou, H.; Baghdad, B.; Bendaou, M. Study of the degradation of steppic rangelands in Béni Mathar and Maâtarka communes (northeastern of Morocco). J. Mater. Environ. Sci. 2014, 5, 2572–2583. [Google Scholar]
  60. García, J.T.; Suárez, F.; Garza, V.; Justribó, J.H.; Oñate, J.J.; Hervás, I.; Calero, M.; García de la Morena, E.L. Assessing the Distribution, Habitat, and Population Size of the Threatened Dupont’s Lark Chersophilus Duponti in Morocco: Lessons for Conservation. Oryx 2008, 42, 592–599. [Google Scholar] [CrossRef] [Green Version]
  61. Brosset, A. L’évolution récente de l’avifaune du nord-est marocain: Perte et gains depuis 35 ans. Rev. Écol. 1990, 45, 237–245. [Google Scholar] [CrossRef]
Figure 1. (A) Geographic location of the study area in Morocco (shaded area). (B) Details of the explored area, showing the Morocco–Algeria border (dashed red line), and the main villages (red circles) and routes (yellow lines). (Na: Nador, Ta: Taourirt, Ou: Oujda, Je: Jerada, Aï: Aïn Bni Mathar, Te: Tendrara, Bo: Bouarfa, Fi: Figuig, Tal: Talsint. (C) Elevation profile (derived from Google earth) of the latitudinal gradient along which the bird surveys were carried out. IDM: De Martone aridity index (see text).
Figure 1. (A) Geographic location of the study area in Morocco (shaded area). (B) Details of the explored area, showing the Morocco–Algeria border (dashed red line), and the main villages (red circles) and routes (yellow lines). (Na: Nador, Ta: Taourirt, Ou: Oujda, Je: Jerada, Aï: Aïn Bni Mathar, Te: Tendrara, Bo: Bouarfa, Fi: Figuig, Tal: Talsint. (C) Elevation profile (derived from Google earth) of the latitudinal gradient along which the bird surveys were carried out. IDM: De Martone aridity index (see text).
Diversity 15 00737 g001
Table 2. Results of the principal component analysis (PCA), showing the score values obtained for each variable for each factor (PCs). Vegetation cover was measured at three different heights (1, 20, and 40 cm). Values with the highest positive (>0.65) and negative (<−0.65) significant correlations with each PC are highlighted in bold.
Table 2. Results of the principal component analysis (PCA), showing the score values obtained for each variable for each factor (PCs). Vegetation cover was measured at three different heights (1, 20, and 40 cm). Values with the highest positive (>0.65) and negative (<−0.65) significant correlations with each PC are highlighted in bold.
VariablePC1PC2PC3PC4PC5
Latitude0.130.78−0.01−0.11−0.26
Elevation0.05−0.660.11−0.020.24
Slope0.120.000.150.150.67
Rocks cover−0.10−0.060.000.040.79
Pebble cover−0.07−0.300.10−0.210.65
Sand cover−0.05−0.770.08−0.14−0.30
Silty loam cover−0.080.80−0.030.00−0.20
Herb cover1−0.170.190.000.710.18
Herb cover200.010.040.030.94−0.03
Herb cover400.05−0.05−0.010.90−0.07
Shrub cover0−0.040.10−0.92−0.06−0.04
Shrub cover200.000.08−0.94−0.05−0.10
Shrub cover400.01−0.03−0.740.05−0.08
Alfa cover00.93−0.010.14−0.120.09
Alfa cover 200.950.030.12−0.120.07
Alfa cover 400.920.160.11−0.100.04
Vegetation heightmax0.80−0.21−0.240.14−0.11
Vegetation heightmode0.87−0.08−0.150.08−0.17
Crops cover−0.190.660.120.380.17
Eigenvalue4.223.262.712.211.57
Total variance (%)22.2017.1714.2511.648.26
Cumulative variance22.2039.3753.6265.2573.52
Table 3. Nested ANOVA table, testing for the effect of season (winter vs. spring) and survey (nested within each season) on each component resulting from the PCA (see Table 2).
Table 3. Nested ANOVA table, testing for the effect of season (winter vs. spring) and survey (nested within each season) on each component resulting from the PCA (see Table 2).
ComponentEffectsdfMSFp-Value
PC1Season10.00030.0000.98
Survey (Season)20.68130.6740.51
Standard error1631.0100
PC2Season11.50781.5170.21
Survey (Season)21.28781.2950.27
Standard error1630.9939
PC3Season10.14740.1460.70
Survey (Season)20.74120.7340.48
Standard error1631.0086
PC4Season15.03005.1760.02
Survey (Season)21.36151.4010.25
Standard error1630.9716
PC5Season15.76747.3230.007
Survey (Season)216.14720.50<0.001
Standard error1630.7874
Table 4. Results of the generalized linear model (GLZ) for the effect of the five components derived from the PCA (Table 2) and the survey factor (two levels) on species richness of the different functional groups and total species richness in spring and winter. The reference levels for Survey factor are “spring 2” and “winter 2”.
Table 4. Results of the generalized linear model (GLZ) for the effect of the five components derived from the PCA (Table 2) and the survey factor (two levels) on species richness of the different functional groups and total species richness in spring and winter. The reference levels for Survey factor are “spring 2” and “winter 2”.
SpringWinter
GroupEstimateSEZ-ValueEstimateSEZ-Value
Granivores
Intercept−0.2070.180−1.15 0.1460.1530.95
Survey−0.3850.323−1.19 −0.2840.240−1.18
PC1−0.2340.2201.06 −0.2890.171−1.68^
PC2−0.0380.131−0.29 −0.1750.120−1.46
PC30.0010.1770.00 −0.0450.093−0.48
PC40.1830.0722.53*0.3420.1871.83^
PC50.0750.1400.54 0.1940.1031.89^
Mixed diet
Intercept−0.1170.177−0.66 −1.0360.282−3.66***
Survey−0.7030.337−2.08*0.3020.3140.96
PC1−0.1270.177−0.71 0.1680.0961.74^
PC20.5100.1333.82***0.6930.1883.68***
PC3−0.2370.154−1.54 −0.1670.088−1.89^
PC40.2540.0713.59***0.4970.2042.44*
PC50.0490.1550.32 0.2130.1511.40
Insectivores
Intercept−1.0220.270−3.78***−1.6340.388−4.21***
Survey0.0090.4510.02 −0.0370.463−0.08
PC1−0.2510.318−0.78 −0.0330.211−0.15
PC20.1090.1780.61 0.1360.2280.59
PC3−0.3250.207−1.56 −0.0080.185−0.04
PC40.1730.0941.83^−0.4760.655−0.72
PC5−0.3160.244−1.29 −0.3200.359−0.89
Total richness
Intercept0.7750.1116.97***0.6760.1155.86***
Survey−0.4440.205−2.16*−0.1430.171−0.83
PC1−0.1960.127−1.53 −0.0210.076−0.27
PC20.2120.0802.63**0.1180.0861.37
PC3−0.1580.101−1.57 −0.0750.059−1.26
PC40.2090.0444.76***0.3130.1332.35*
PC5−0.0120.095−0.13 0.1540.0841.84^
^ p < 0.10, * p < 0.05, ** p < 0.01, *** p < 0.001.
Table 5. Results of the generalized linear model (GLZ) presenting the effect of the five components derived from the PCA (Table 2) and the survey factor (two levels) on bird density of the different functional groups and total bird density in spring and winter. The reference levels for survey factor are “spring 2” and “winter 2”.
Table 5. Results of the generalized linear model (GLZ) presenting the effect of the five components derived from the PCA (Table 2) and the survey factor (two levels) on bird density of the different functional groups and total bird density in spring and winter. The reference levels for survey factor are “spring 2” and “winter 2”.
SpringWinter
GroupEstimateSEZ-ValueEstimateSEZ-Value
Granivores
Intercept1.9960.2079.64***3.0940.28510.85***
Survey−0.8230.379−2.17*0.0820.3980.20
PC1−0.3420.278−1.22 −0.0790.223−0.35
PC20.2850.1541.84^−0.4480.210−2.13*
PC30.0270.2240.12 0.3380.2801.20
PC40.3090.0783.95***0.6250.2712.30*
PC50.3520.1432.46*0.2170.1591.36
Mixed diet
Intercept2.1140.2657.96***1.4390.3953.64***
Survey−1.4790.536−2.76**0.1400.4180.33
PC1−0.3630.340−1.06 0.2620.1361.91^
PC20.8370.1964.26***0.7520.2602.89**
PC3−0.3960.206−1.91^−0.2440.113−2.15*
PC40.3640.1003.62***0.7600.2802.71**
PC5−0.1010.269−0.37 −0.3450.361−0.95
Insectivores
Intercept−0.2670.283−0.94 −0.5850.368−1.58
Survey0.3140.4370.72 −0.1500.475−0.31
PC1−0.3510.353−0.99 −0.0490.218−0.22
PC20.0940.1770.53 0.1140.2340.49
PC3−0.3740.211−1.77^−0.1050.157−0.66
PC40.1830.0892.05*−0.3970.618−0.64
PC5−0.4270.248−1.72^−0.1020.319−0.32
Total density
Intercept2.9230.15918.33***3.5230.22215.83***
Survey−1.1010.318−3.46***−0.0290.333−0.09
PC1−0.3880.229−1.69^0.0420.1480.28
PC20.5450.1224.43***−0.1380.168−0.82
PC3−0.2340.151−1.55 0.0350.1480.24
PC40.3240.0645.06***0.5820.2282.55*
PC50.1420.1381.02 0.1460.1560.93
^ p < 0.10, * p < 0.05, ** p < 0.01, *** p < 0.001.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Oñate, J.J.; Suárez, F.; Calero-Riestra, M.; Justribó, J.H.; Hervás, I.; de la Morena, E.L.G.; Ramírez, Á.; Viñuela, J.; García, J.T. Responses of Bird Communities to Habitat Structure along an Aridity Gradient in the Steppes North of the Sahara. Diversity 2023, 15, 737. https://doi.org/10.3390/d15060737

AMA Style

Oñate JJ, Suárez F, Calero-Riestra M, Justribó JH, Hervás I, de la Morena ELG, Ramírez Á, Viñuela J, García JT. Responses of Bird Communities to Habitat Structure along an Aridity Gradient in the Steppes North of the Sahara. Diversity. 2023; 15(6):737. https://doi.org/10.3390/d15060737

Chicago/Turabian Style

Oñate, Juan J., Francisco Suárez, María Calero-Riestra, Jorge H. Justribó, Israel Hervás, Eladio L. García de la Morena, Álvaro Ramírez, Javier Viñuela, and Jesús T. García. 2023. "Responses of Bird Communities to Habitat Structure along an Aridity Gradient in the Steppes North of the Sahara" Diversity 15, no. 6: 737. https://doi.org/10.3390/d15060737

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop