Abstract
Road construction is a catalyst of land-use change. In this analysis, drawing on the growth machine framework, the investigators assess whether there is any reciprocal feedback in this process, specifically asking if there is a bidirectional association between road construction and land development. To answer that question, the investigators use recently released satellite imagery from the National Land Cover Database’s Impervious Descriptor dataset. The analysis covers the years 2010 to 2024, during a period of sustained growth for the 3,203 census tracts in the four major metropolitan areas of the “Texas Triangle.” The investigators compare results from two types of longitudinal models using first differences: cross-lagged structural equation models and spatial regression models. Results from these models reveal that the impact of road construction on land development is proportionately greater than the reciprocal impact of land development on road construction, providing an example of an asymmetric reciprocal socioecological association. From the growth machine perspective, the investigators argue that road construction and land development do not function as equal, twin engines of growth but rather as primary and subordinate catalysts of urban expansion, implying an infrastructure-led strategy of the growth coalition elite.
In this study, we set out to test whether there is an asymmetric reciprocal association between road construction and land development across the 3,203 census tracts in the four major metropolitan areas of Texas, covering the years 2010 to 2024, during a period of sustained urban growth in the state. Over this 15-year period, the population of the metropolitan areas of Austin, Dallas–Fort Worth, Houston, and San Antonio grew by more than 30 percent, a growth rate greater than four times the national increase of approximately 8 percent during the same time frame. These four metropolitan areas are part of an emerging macro-metropolitan system, sometimes referred to as the “Texas Triangle” (Fulton et al. 2021). This region represents a unique experience in the recent history of urban growth, with projections indicating that it will soon be home to four of the nation’s 10 most populous cities, with a collective economic output that would rank as the fifteenth largest global economy (Fulton et al. 2021). Meanwhile, Texas’s metropolitan growth has been facilitated by relatively permissive land-use regulations and state-level policies that enable aggressive land development (e.g., Brannan 2010; Martinez 2025). Yet this high-magnitude growth has resulted in various urban infrastructure challenges; for instance, the pace of road construction has not been able to keep up with the demand that arises from aggressive land development, resulting in more traffic congestion and increased accidents (Barfield 2022; Saginor et al. 2011). According to the Texas Department of Transportation (2024a), between 1975 and 2015, the state’s number of registered vehicles grew by 172 percent, but its highway capacity (in lane miles) grew by only 19 percent. And, this infrastructure challenge continues into the present, particularly for the “Texas Triangle”; as of 2024, according to a standardized ranking of 2,100 congested roadways across the state, measured in terms of delay per mile, “94 of the top 100 . . . are in Austin, Dallas-Fort Worth, Houston and San Antonio” (TxDOT 2024b).
Thus, in Texas’s major metropolitan areas, because of all this growth, there has been a mismatch between infrastructure capacity and the demand for infrastructure. In this study, using recently released land cover data from satellite imagery (USGS 2024), we draw from growth machine theory to assess this mismatch, asking whether there is an asymmetric association between road construction and land development. As discussed later, growth machine theory provides a comprehensive framework to study not only whether this association is asymmetric but also whether it is reciprocal. Thus, we ask whether the association between these two socioecological dynamics of urban growth is bidirectional and asymmetric, such that the impact of road construction on land development is proportionately unequal to the reciprocal impact of land development on road construction. These questions draw from theoretical and empirical research in the growth machine literature. Indeed, for empirical studies of the growth machine, Logan, Whaley, and Crowder (1997) emphasized the importance of “studying effects over time, and estimating causal reciprocal relationships” (p. 624). We take up this call, following in the path of more recent empirical research on the topic using satellite imagery on land cover (e.g., Clement et al. 2023; see also Aljoufie et al. 2011; Tong and Qiu 2020). For example, Clement et al. (2023) examined the reciprocal association between core features of the growth machine, including land development. Yet although Clement et al. (2023) claimed that the reciprocal association was asymmetric, they did not report the results of any statistical techniques to test for asymmetry, for example, using a likelihood ratio test to compare unconstrained versus equality-constrained structural equation models (e.g., Ambrona and López-Pérez 2014; Kenny 1979).
For this study, drawing from the reciprocal dynamics of the growth machine theory, we conduct formal, statistical tests of the mismatch between road construction and land development, the results of which would signal whether infrastructure-led development, as a core dynamic of the growth machine model, is unbalanced or asymmetric. In this context, we assess whether there is a bidirectional association between road construction and land development and whether this association is asymmetric. By evaluating the potential for an asymmetric reciprocal association between these two key socioecological dimensions of urban growth, and by identifying which of these forces might exert a dominant influence on the other, we can help clarify the degree of symmetry or asymmetry in the operation of the growth machine. Moreover, from the growth machine perspective, given the significant roles played by private real estate interests as well as the state and local governments in Texas, testing whether this association is symmetric or asymmetric can also help clarify the functional relationship between public and private actors of the growth coalition, in terms of their support for road construction and land development. Specifically, this analysis can provide evidence about whether public subsidies and private investment function equally as twin engines of growth, or whether one serves as the primary catalyst that triggers urban expansion.
Theoretical Framework
In sociology, growth machine theory offers a multipurpose sociological lens for addressing the socioecological dimensions of community change, from population growth and changes in unemployment to the physical reconfiguration of the landscape (Logan and Molotch 2007; Molotch 1976). With the notion that growth is a “syndrome of associated events” (Molotch 1976:310), the growth machine draws from the human ecological tradition (Molotch 1967), adopting the view that the socioecological dimensions of community change are connected through interdependent, reciprocal processes (see also Elliott and Frickel 2015). To this human ecological idea, the growth machine adds a political economy critique, claiming that “any given piece of real estate has both a use value and an exchange value” (Logan and Molotch 2007:1). Yet it is the pursuit of exchange value that serves as the primary motive of urban growth; indeed, as Logan and Molotch argued, urban growth functions as a mechanism for increasing real estate value, an objective that unites a broad coalition of stakeholders. The exchange value potential of real estate generates a widely shared interest in land-use intensification (e.g., Gerber and Phillips 2003); all the same, despite the broad support it may receive, the economic gains of land development remain concentrated among the wealthy and powerful of the growth coalition elite (e.g., Dye and McMillen 2007; Guerrieri, Hartley, and Hurst 2013). In this light, land development, and the way it economically benefits real estate developers and wealthy residents, represents a core dynamic in the growth machine perspective, making it a fundamental, interdependent part of urban growth.
Despite the theory’s emphasis on the physical intensification of land, foundational empirical work by Molotch (1976) and subsequent quantitative research often relied on more accessible demographic and economic measures, such as population growth or unemployment (e.g., Fodor 2012). Although scholars produced insightful growth machine studies on land use, this research tended to focus on the demographic, economic, and political correlates of land development, rather than relying on direct measures of land cover (e.g., Logan and Zhou 1989; Warner and Molotch 2000). To be sure, qualitative sociologists have provided detailed historical accounts of landscape transformation from the growth machine perspective (Rudel 2009). Nevertheless, of relevance to the present analysis, the expanding availability of high-resolution satellite imagery has enabled quantitative sociologists to develop concrete measures of land cover, which has contributed to a flourishing body of research for testing growth machine propositions on landscape transformation (Clement and Podowski 2013; Ergas, Clement, and McGee 2016; Landis 2017; Pravin 2019; Smiley and Hakkenberg 2020; Steele and Luloff 2003; Wongso and Göçmen 2025). Meanwhile, scholars have long pointed out that road construction is a primary catalyst of land development; the construction of roads penetrates undeveloped areas, enabling the flow of human activity and facilitating the transformation of this land (e.g., Bernier and McGlinchy 2023; Greenberg, Wishart, and Danielsen 2024; Rudel 2009). Indeed, Logan and Molotch (2007) argued that “transportation doesn’t just serve growth, it creates it” (p. 74). As such, viewed through the lens of the growth machine, road construction can be framed as the catalytic mechanism for the “increasing intensification of the land use of the area” (Molotch 1976:309), signaling its role in enhancing exchange value by fostering land development. Yet in the growth machine framework, land development serves a dual purpose: it maximizes exchange value for private real estate interests while simultaneously providing the state with the fiscal capacity to fund further road construction, which effectively socializes the costs of the infrastructure required for continued urban expansion (e.g., Harvey 1982). Indeed, through various channels, such as local taxes, developer impacts fees, and tax increment financing (TIF), more land development generates more public revenue, which is dedicated specifically for more road construction (e.g., Briffault 2010; Jeong and Feiock 2006; Kirkpatrick 2016; Saginor et al. 2011; Sroka 2024). In this way, although road construction is a primary catalyst of land-use change, the subsequent land development then reciprocally funds new road construction, thereby reigniting the catalyst of urban growth.
Yet as mentioned in the introduction, the continued challenges of congestion and traffic accidents in Texas’s metropolitan areas suggest an imbalance or asymmetry in the association between road construction and land development, an observation that informs our research questions, which we evaluate in the following analysis.
To develop a testable idea of an unbalanced bidirectional relationship, we use the framework of asymmetry as established in the quantitative social science literature (e.g., Allison 2019; Lieberson 1985; York and Light 2017). In examining the link between road construction and land development, we specifically test for reciprocal asymmetry, which we treat as a distinct concept from directional asymmetry. Indeed, the notion of reciprocal asymmetry (the impact of X on Y ≠ the impact of Y on X) is related to the notion of directional asymmetry (the impact of ↑X on Y ≠ the impact of ↓X on Y), but they ultimately offer different insights. On one hand, directional asymmetry not only can be helpful in identifying differences in the magnitude of impact, but it can also reveal how the mechanism through which ↑X affects Y is different than ↓X affecting Y. On the other hand, the concept of reciprocal asymmetry is primarily helpful in identifying differences in the magnitude of impact rather than revealing new information about the mechanisms of influence (e.g., Ambrona and López-Pérez 2014).
To demonstrate this idea for our study, we will use the X variable for road construction and the Y variable for land development. So, with reciprocal asymmetry, because road construction (X) and land development (Y) are different variables, the mechanism through which X affects Y is plainly different from the mechanism through which Y affects X. Thus, a test of reciprocal asymmetry will not reveal any new information about the mechanisms through which X and Y influence each other; but, with the appropriate statistical techniques, an assessment of reciprocal asymmetry reveals differences in the proportional influences of these two variables on each other. In this context, reciprocal asymmetry points to an imbalance in the bidirectional system of change, indicating that one variable acts as a primary stimulus of change while the other variable is still a significant yet subordinate force of feedback.
On that note, if the interdependent factors within the growth machine model show reciprocal asymmetry, it means their proportional influences on each other are unequal, reflecting an overall imbalance in the operating forces of urban growth. Specifically, with respect to the association between road construction and land development, reciprocal asymmetry would point to which variable is the primary instrument and which variable is the subordinate feedback force in catalyzing local, urban growth. Such a finding can contribute to a more nuanced discussion about the theoretical and practical implications of urban growth as well as alternatives to the current form of development in the United States, especially given its car-centric history. For instance, if there is an asymmetric reciprocal relationship between road construction and land development, then we argue that road construction serves as a primary instrument of land-use change, but land development still operates as a significant, yet secondary, force of feedback. This reciprocal asymmetry would point to a systemic pattern where road construction stimulates infrastructure-led land development, which then in turn subsequently feeds back into, albeit at a lower magnitude, more road construction, maintaining an unbalanced system of car-centric urban growth.
In other words, reciprocal asymmetry, consistent with the growth machine perspective, would help clarify the state’s role in the growth coalition. Through state support, road construction acts as a lead but interdependent socioecological force that contributes to subsequent private development, which then reciprocally contributes to road construction as a form of urban infrastructure that is designed around the private automobile (e.g., Angelo and Wachsmuth 2020; see also Harvey 1982). By providing the infrastructure that precedes private investment, the state does more than facilitate urban growth; it effectively helps socialize the initial costs and risks of urban expansion. If there is a reciprocal association between road construction and land development, then this would create a cycle where the public subsidizes the infrastructure which then ensures the intensification of private land development. In growth machine terms, if road construction is found to be a primary catalyst of expansion, conceptually, it is helping maximize exchange value and, theoretically, road construction emerges as a structural cornerstone of the growth machine perspective. Moreover, if urban growth is calibrated for infrastructure-led development, then a finding of asymmetry could also be seen as making a contribution to the vast literature describing how public transit has been marginalized in urban America (e.g., Lee and Sener 2017; Mattioli et al. 2020). With evidence of asymmetry, the reciprocal forces of urban growth could be described as structurally designed to favor projects that maximize exchange value through prioritizing private automobile access, instead of public transit, even if road construction projects, in terms of traffic congestion and accidents, come at a higher social and environmental cost than public transit (e.g., Warner and Molotch 2000).
To that end, we test for reciprocal asymmetry between road construction and land development across 3,203 census tracts in Texas’s major metropolitan areas, evaluating whether these dynamics function in a balanced feedback loop or in an unbalanced system where public infrastructure serves as the primary catalyst for private exchange value. In an unbalanced system, such a relationship would, by extension, imply that alternative development models are marginalized in favor of a car-centric landscape.
Data and Analytic Techniques
For our analysis, we focus on the 3,203 census tracts in the four major metropolitan areas of Texas: Austin, Dallas–Fort Worth, Houston, and San Antonio. As mentioned, the study covers the years 2010 to 2024, during a period of sustained growth in the state of Texas. To account for the multiyear construction time frame required to complete infrastructure projects (Glaeser and Poterba 2021), we use four waves of data from 2010, 2015, 2019, and 2024; furthermore, these four- and five-year intervals generally align with the windows for major infrastructure planning in Texas (e.g., NCTCOG 2022). Variable names, descriptions, and sources are provided in Table 1.
Variable Names, Descriptions, and Sources (n = 3,203 Census Tracts).
The data for the road construction and land development variables come from the National Land Cover Database’s Impervious Descriptor satellite imagery (USGS 2024). With a tract-level shapefile, based on 2010 census boundaries, we use the Tabulate Area tool in ArcGIS to quantify, for each census tract, and for each year 2010, 2015, 2019, and 2024, the area (in square meters) covered by roads and the area (in square meters) covered by all other nonroad, human-built impervious surfaces, including residential structures, commercial centers, industrial facilities, and the like. For the analysis, both variables are turned into percentages according to the total area (in square meters) of the census tract; the road construction variable equals the percentage of the total area covered by all road types, which includes everything from unpaved rural roads to residential streets to interstate freeways. 1 The land development variable equals the percentage of the total area covered by all nonroad, human-built impervious surfaces. The use of percentages for the analysis standardizes the slope estimates, which indicate the percentage point change in the dependent variable for every 1 percentage point change in the predictor variable, holding the rest of the equation constant. 2
Drawing from previous social science research on land-use change using satellite imagery (e.g., Clement et al. 2023; Smiley and Hakkenberg 2020), we also include three variables to control for local demographics and socioeconomic status: population density, percentage white, and median household income. Population density is the total resident population of the census tract divided by the tract’s total land area (in 1,000 m2); percentage white is the percentage of the tract’s population that is non-Hispanic white; and median household income is the middle income of all households in the census tract, adjusted for inflation, in 2010 constant dollars. Before computing first differences, we calculate the natural logarithm for population density and median household income after adding a constant of 1. 3 These variables are used as controls to assess whether their impact may confound the reciprocal association between road construction and land development; on that note, for the control variables, we use only data from 2010, 2015, and 2019 as the lagged predictors of road construction and land development.
To test whether there is an asymmetric reciprocal association between road construction and land development, we use two types of modeling techniques: structural equation models and spatial regression models. For both, we estimate models using two waves of first-differenced values from the four time periods: 2010, 2015, 2019, and 2024; because the variables are percentages, these first differences are interpreted as percentage point changes. Using lagged first differences from four times periods of data yields three waves of first differences: Δ2010–2015, Δ2015–2019, Δ2019–2024, from which we can discern whether slope estimates change over time. All structural equation and spatial regression models are estimated using the robust variance-covariance estimator to account for potential heteroskedasticity (the “robust” vce option in Stata) (White 1980).
For the structural equation models, we use a cross-lagged structure with autoregressive parameters. See Figure 1, which does not include the lagged predictors for the control variables. The cross-lagged associations are between road construction and land development: ΔRoadt−1 → ΔLand t and ΔLandt−1 → ΔRoad t . The series of cross-lagged models we estimate involve comparing results from an unconstrained model, where the cross-lagged estimates are free to vary, and an equality-constrained model, where the cross-lagged estimates are set to be equal. Because the only parameters that are different between the two nested models (i.e., unconstrained and constrained models) are the cross-lagged slopes, then comparing their model fits with a likelihood ratio test would tell us whether we should treat the cross-lagged effects as equal or unequal; if the likelihood ratio test is not significant, then the model with equal slopes exhibits better fit, which provides evidence that the reciprocal association is symmetric; if the likelihood ratio test is significant, then the model with unequal slopes exhibits better fit, which provides evidence that the reciprocal association is asymmetric (e.g., Ambrona and López-Pérez 2014; Kenny 1979).

Cross-lagged first-differenced structural equation model.
Covering the years 2010, 2015, 2019, and 2024, using two waves of first-differenced values, we estimate two models: an unconstrained versus a constrained model. Following are the generic equations for these models.
Unconstrained structural equation model:
with:
Equality-constrained structural equation model:
with:
In these equations, ΔYit and ΔXit are the first-differenced values of Y and X, and ΔYi,t−1 and ΔXi,t−1 are the lagged first-differenced values of Y and X. β1 and β3 are the autoregressive slope estimates. In the unconstrained model, β2 and β4 are the cross-lagged slope estimates for ΔYi,t−1 and ΔXi,t−1, which are allowed to vary; in the equality-constrained model, the cross-lagged slope estimates β2equal are made equal. The terms u* it and v* it capture the unexplained variation in ΔYit and ΔXit. The term σ uv is the covariance between the unexplained residuals in ΔYit and ΔXit.
For the spatial regression analysis, our analysis begins with a temporally lagged first-differenced model structure, on the basis of the approach recommended by Allison (2009). Allison (2009) demonstrated that with three time periods of data, two separate ordinary least squares models of lagged first differences can be run, effectively creating a reciprocal association fixed-effects model to assess bidirectional relationships. We applied this strategy using first-differenced data from three time periods covering the years 2010, 2015, 2019, and 2024, which yielded two models where road construction is the dependent variable (and land development is the independent variable) and two models where land development is the dependent variable (and road construction is the independent variable), giving us a total of four cross-lagged models. We then added a spatial lag term ρ and spatial error term λ to these reciprocal association fixed-effects models. With both ρ and λ, the model is a spatial autoregressive model with spatial autoregressive disturbances or a spatial autoregressive combined model (see Anselin and Florax 1995). The generic equations for these models are as follows:
Spatial regression models:
In these equations, β1 and β2 represent the (temporally) lagged slope estimates;
Results
Table 2 reports univariate statistics for 2010, 2015, 2019, and 2024 for the road construction and land development variables, and for 2010, 2015, and 2019 for the control variables. We use the first differences for the analysis. Looking at the first differences, between 2010, 2015, 2019, and 2024, there were consistent increases for both road construction and land development (all first differences were positive; p < 0.001, paired t tests). In terms of the control variables, the results show inconsistent variation over time. On one hand, population density, on average, increased between 2010, 2015, and 2019 (as log differences, these numbers are interpreted, roughly, as 4.1 percent and 2.7 percent increases in population density). On the other hand, the typical tract became more white and less affluent between 2010 and 2015; but between 2015 and 2019, we observe the opposite: the typical tract became less white and more affluent.
Univariate Statistics (n = 3,203 Census Tracts).
Note. The values for population density and median household income were log transformed before calculating the first differences.
p < .001 (one-tailed paired t tests).
To assess whether road construction and land development are reciprocally associated, we begin with the first-differenced cross-lagged structural equation modeling results in Tables 3 and 4. The unconstrained results are also displayed in Figure 1, which does not display the control variables.
Unconstrained and Constrained Cross-Lagged Structural Equation Models (n = 3,2023 Census Tracts).
p < .10, *p < .05, **p < .01, and ***p < .001 (two-tailed tests).
Cross-Lagged Structural Equation Model Fit Indices and Asymmetry Test (n = 3,2023 Census Tracts).
Note. The model fit indices are derived from structural equation models that do not include population density as an additional control. Although adding population density as a control does not change the substantive findings of the analysis, including the results of the asymmetry tests, it does seriously degrade model fit in terms of the TLI, which decreases from 0.947 to 0.838 after including population density. AIC = Akaike information criterion; BIC = Bayesian information criterion; CD = coefficient of determination; CFI = confirmatory fit index; RMSEA = root mean square error of approximation; SRMR = standardized root mean square residual; TLI = Tucker-Lewis index.
p < .001 (two-tailed tests for slope equality and residual spatial autocorrelation).
Using the 2015–2019 interval for the endogenous outcomes, the significant slope estimates from both the unconstrained and equality-constrained models show that there was a positive feedback loop between road construction and land development. Because the variables are percentages, the slope estimates are interpreted as a percentage point change; for instance, in Table 3, looking at model 1, across the 3,203 census tracts, for every 1 percentage point increase in the rate of change of road construction between 2010 and 2015, the rate of change of land development between 2015 and 2019 increased by 0.159 percentage points (p < .01), all else being equal (including the prior rate of change of land development, which is based on the significant slope estimate for the land development parameter in Table 4). Likewise, for every 1 percentage point increase in the rate of change of land development between 2010 and 2015, the rate of change of road construction between 2015 and 2019 increased by 0.026 percentage points (p < .05), all else being equal (including the prior rate of change of road construction).
However, using the 2019–2024 interval for the endogenous outcomes, there was not a significant bidirectional association between road construction and land development. In fact, only the unidirectional impact of road construction on land development was significant (b = 0.368, p < .05).
For the unconstrained model in Table 3, the slope estimates for road construction were larger than the slope estimates for land development. Thus, in Table 4, to assess whether these differences were significant and whether there is evidence of asymmetry, we ran two types of tests: a likelihood ratio test to check for differences in model fit, and two z tests to check if the differences between the slope estimates were significant for both first-differenced waves. The results for both the likelihood ratio test and the two z tests were statistically significant (p < .001). Taken together with the slope estimates, these results highlight three points for our study. First, that allowing the cross-lagged estimates to vary freely significantly improved model fit, providing evidence of asymmetry. Second, the slope estimates for road construction were significantly greater than the slope estimates for land development; moreover, this difference in slope estimates, on average, increased between the 2010–2015 and 2015–2019 first-differenced waves; this is evidence that the asymmetry became more pronounced over time. Third, according to the structural equation modeling results, the relationship between road construction and land development was asymmetric yet still reciprocal between 2010 and 2019; however, between 2015 and 2024, there was no bidirectional association between these two variables. This suggests a shift in the growth dynamic, where the early interdependence between these variables evolved into a more linear process in the most recent time interval.
Looking at Table 4, although the traditional structural equation model fit indices (i.e., χ2, root mean square error of approximation, Akaike information criterion, Bayesian information criterion, Tucker-Lewis index, confirmatory fit index, standardized root mean square residual, and coefficient of determination) all show that the unconstrained model exhibits better fit than the constrained model, 4 we also test whether the residuals from these models are spatially autocorrelated. If they are, this spatial dependence might be biasing slope estimates, suggesting the need for a spatial regression approach. At the bottom of Table 4, we report these results; we find that the residuals from the primary endogenous outcome variables (i.e., road construction and land development) do exhibit highly significant, positive spatial autocorrelation. Figure 2 displays a map of the residuals from the unconstrained model when land development Δ2015–2019 is the endogenous outcome. We do this only for Harris County, given that it would become difficult to see tract boundaries if we presented residuals for all census tracts in all the counties across all four metropolitan areas. In this map, consistent with the estimate for the spatial autocorrelation of the residuals (ρ = 0.334, p < .001), we observe a discernible spatial pattern in the clustering of the residuals. Although some tracts outside the urban core did exhibit lower than predicted land development, the spatial clustering of higher residuals outside the urban core reflects extra development the model could not account for, which is likely the above-average pace of land development characteristic of suburban sprawl. On that note, to check whether this spatial dependency might be biasing the slope estimates, we next ran a series of cross-lagged spatial regression models.

Harris County residuals from the unconstrained structural equation model in Table 3, with land development (2015–2019) as the endogenous outcome.
Table 5 reports the results from these spatial regression models. We first note that the estimates for ρ and λ are significant in all models (at least p < .05), and the residuals from these models no longer exhibit significant spatial autocorrelation. Again, only for Harris County, Figure 3 displays the residuals from model 1 of Table 5, in which land development Δ2015–2019 is the dependent variable. Now, in this spatial regression model, we do not observe a discernible spatial clustering in the residuals; instead, we see more randomly distributed residuals. Thus, although suburban sprawl is occurring, the inclusion of the spatial parameters successfully minimized any potential bias in the slope estimates (of the predictor variables) that stems from such spatial clustering processes. Indeed, looking at the slope estimates in Table 5, we do observe changes in the slope estimates. First, in the spatial regression models, the cross-lagged slope estimates for both road construction and land development are now all positive and significant, which refines the results from the structural equation models and now suggests a significant reciprocal association for all three first-differenced waves: Δ2010–2015, Δ2015–2019, and Δ2019–2024. Second, in Table 5, the results from the Z-tests of the equality of slope estimates indicate that the estimates for road construction were significantly greater than the estimates for land development; yet this asymmetric association did not intensify over time like it did in the structural equation modeling analysis. Ultimately, these findings suggest that while the intensity of the imbalance persisted, the bidirectional link between road construction and land development is more enduring than the structural equation modeling analysis initially indicated.
Cross-Lagged Spatial Regression Results (n = 3,203 Census Tracts).
Note. The slopes for the cross-lagged and control variables are the direct effects from the spatial regression model.
p < .10, *p < .05, **p < .01, and ***p < .001 (two-tailed tests).

Harris County residuals from spatial regression model 1 in Table 5, with land development (2015–2019) as the dependent variable.
Before we discuss these results and make some concluding comments, we briefly report the results for the control variables. Across the structural equation and spatial regression models, the slope estimates for the control variables were not consistent. For the structural equation modeling analysis, we will focus on the unconstrained model because it exhibited significantly better fit than the constrained model. In the unconstrained structural equation model, an increase in the percentage of the white population was significantly associated with lower rates of land development in the earlier first-differenced wave, but then, in the later first-differenced wave, percentage white was positively associated with road construction, along with median household income. Meanwhile, population density was negatively associated with road construction in the later first-differenced wave in the unconstrained structural equation model. However, we caution making any conclusive generalizations about the control variables given that, in the spatial regression models, the slope estimates for the control variables were no longer significant.
Discussion and Conclusion
Drawing from growth machine theory, the foregoing analysis of 3,203 census tracts across the four major metropolitan areas of the “Texas Triangle” (Fulton et al. 2021) reveals that road construction and land development are characterized by a persistent asymmetric association. Although initial structural equation modeling results suggested a shift in the later waves toward a unidirectional relationship, the spatial regression models suggest that there is a persistent bidirectional link that is significant across the entire study period. Although only the spatial regression model shows that the reciprocal relationship persists across the entire study period, the association between these two socioecological processes is consistently asymmetric in both the structural equation modeling and spatial regression analyses, as the slope estimates for road construction significantly exceed those for land development for all first-differenced waves. Thus, the spatial regression results identify a feature of reciprocal asymmetry in the process of urban growth, which helps refine growth machine theory: while both road construction and land development are socioecological dynamics that function as twin engines of urban growth, these twins are not equal; rather, infrastructure provision acts as the primary catalyst of change. Ultimately, in describing these key socioecological dynamics of urban growth in the “Texas Triangle,” we argue that road construction dictated the pace of land development but in the presence of a subordinate yet significant feedback from land development.
To elaborate on these findings and connect them to the actual on-the-ground experience of urban growth in Texas, we display two sets of images (Figures 4 and 5) based on an analysis of the residuals from the spatial regression models 1 and 2. Again, focusing on Harris County as an example, we highlight two census tracts in that county, where the rates of land development and road construction were far higher than predicted on the basis of the equations from models 1 and 2. In Figure 4, we display imagery for census tract 48201343200, which is in southeast Harris County, containing much of the Fairmont Park neighborhood in La Porte. Although development in Fairmont Park began in the late 1950s, and the neighborhood was relatively built out by the 1980s, the area, like much of greater Houston, has continued to experience development. For example, this tract experienced a rate of land development between Δ2015–2019 that was 13.652 percentage points higher than predicted on the basis of the equation used to estimate model 1. Then, in Figure 5, we display imagery for census tract 48201241103, which is in north central Harris County, containing a portion of the Timber Lane area of Spring. Spurring the development of the Breckenridge Forest neighborhood, this tract experienced a rate of road construction between Δ2015–2019 that was 10.653 percentage points higher than predicted based on the equation used to estimate model 2. In Figure 5, we can see the sequence of infrastructure-led development, where prior road construction penetrates undeveloped land, opening it up for later development, which would be the Breckenridge Forest neighborhood.

Comparison of land development over time in census tract 48201343200.

Comparison of road construction over time in census tract 48201241103.
Before we offer some concluding thoughts on the theoretical and practical implications of the study, we highlight some limitations of the foregoing analysis, with an intention to stimulate future scholarship.
First, although the National Land Cover Database’s Impervious Descriptor satellite imagery distinguishes between road area and nonroad area, it does not distinguish between the different types of roads (e.g., unpaved rural roads, residential streets, highways) and the different types of land development (e.g., residential, commercial, industrial). As such, this study is not able to test whether the road construction-land development reciprocal association varies depending on the type of road or the type of land development. To be fair, although deriving both variables from a single satellite imagery dataset ensures consistency in how land cover categories are defined, future research might consider whether merging different data products is feasible.
Second, for the purpose of this study, which was focused on whether there was an asymmetric reciprocal association between two key socioecological processes of urban growth, the control variables were not treated as endogenous outcomes but rather simply as predictors of road construction and land development. For instance, as we discussed, adding population density to the SEM models, simply as an exogenous control variable, resulted in a deteriorated model fit, in terms of a lower Tucker-Lewis index; however, there is scholarship showing that both land development and road construction can reciprocally influence changes in population size (Baum-Snow 2007; Clement et al. 2023). Future research, in more elaborate models, may assess whether these key socioecological processes of urban growth feed back into the demographic and socioeconomic variables we examined here.
Third, although the spatial regression models successfully accounted for residual autocorrelation, the structural equation modeling results exhibited spatial clustering of high residual values in suburban tracts outside the urban core, which is a pattern reflecting the rapid pace of development characteristic of suburban sprawl. Consequently, given the highly significant, positive spatial lag estimates (ρ) for land development in the spatial regression models, all these results taken together suggest that this sprawling development is driven significantly by a spatial contagion process, where the intensification of land use in one suburban tract significantly increases the probability of development in adjacent tracts. Thus, although our spatial regression models successfully accounted for this contagion, future research might directly address this spatial process, perhaps asking how networks of growth coalition elites coordinate land development strategies across suburban tracts.
Fourth, our findings about census tracts in the major metropolitan areas of Texas are situated within a unique period of urban growth, between 2010 and 2024. Although this time frame provides a clear window into the operations of a high-magnitude form of urban growth, it may also reflect specific land use policies and economic conditions that are not universal. Future scholarship can explore whether and how the finding of reciprocal asymmetry may manifest in different political contexts or in areas of contraction or stagnation.
To conclude, we elaborate on the theoretical and practical implications of the results of our analysis. First, this analysis adds to a burgeoning body of research using satellite imagery on land cover to test growth machine propositions (Clement and Podowski 2013; Ergas et al. 2016; Landis 2017; Pravin 2019; Smiley and Hakkenberg 2020; Steele and Luloff 2003; Wongso and Göçmen 2025). Furthermore, although Clement et al. (2023) claimed that the reciprocal association they observed was asymmetric, they did not use any formal statistical tests for asymmetry, such as a likelihood ratio test comparing unconstrained versus equality-constrained structural equation models; our study redresses this gap (see also Ambrona and López-Pérez 2014; Kenny 1979). Thus, our study, using these techniques, showed that road construction created arterial corridors through which human activity moved and then expanded its impervious footprint (see Figure 5), and reciprocally, land development fed back into more road construction, although its proportional influence was smaller. These results are consistent with but add a reciprocal-asymmetry nuance to how the growth machine framework describes infrastructure-led development (e.g., Bernier and McGlinchy 2023; Greenberg et al. 2024; Rudel 2009).
As previously mentioned, in the state of Texas, through various direct and indirect channels, land development can spur further road construction. With land development comes more property taxes; while state and federal taxes fund major highway construction, property taxes are used for local road construction. Municipalities in Texas also have specific policy mechanisms they can use to fund road construction through land development; as two examples: (1) developers pay impact fees on new projects, creating dedicated road construction funds, and (2) tax increment reinvestment zones create revenue to fund public improvements, including major road construction (e.g., Briffault 2010; Jeong and Feiock 2006; Kirkpatrick 2016; Sroka 2024). In the first instance, although developer impact fees require private interests to contribute to the cost of infrastructure for new projects, these fees generally cover only a fraction of the total expenditure, leaving much of the financial burden for road expansion to be subsidized by the public through broader tax revenues. In the second instance, TIF allows a portion of the future property tax revenue generated by a specific development to be captured and reinvested directly back into that area’s infrastructure. This mechanism essentially uses anticipated public funds to pay for the immediate road construction necessary to make the private development viable, thereby allowing developers to realize the benefits of increased land values while the public assumes the long-term financial risk.
On that note, this study further reinforces the supportive role of the government in the growth coalition. Not only does the government fund road construction, which, as a primary catalyst, opens up land to new development, but it also can then generate revenue from that newly developed land, reciprocally, to fund more road construction, although this feedback has a proportionally smaller impact. Practically, as mentioned, developers’ impact fees do go to road construction, but they represent a smaller source of funding for road construction than public tax dollars. Yet consistent with the political-economy of the growth machine, developers and growth coalition boosters also openly advocate for TIF. For instance, Near Southside, Inc. (NSI), is a nonprofit organization that works with the real estate industry to promote the further development of neighborhoods immediately south of downtown Fort Worth. During the time period of this study, NSI was contracted with the city to administer a TIF district, working with various stakeholders to support “public infrastructure constructed in conjunction with dozens of private redevelopment projects, providing an incentive that defrays development costs for projects that advance district goals and, without the TIF’s assistance, wouldn’t be financially feasible” (NSI 2026). Thus, developers and growth coalition boosters, like NSI, advocate for state support, like the creation of TIF districts, to bridge the financing gaps for the public infrastructure necessary for their projects. In this way, with state assistance, the real estate industry can have risk transferred to the public, yet they can still capture the profit from the subsequent increase in land value that comes from the development (e.g., Briffault 2010; Kirkpatrick 2016). Whether through large-scale greenfield projects or the development of smaller parcels, this state-led infrastructure investment would lower the barrier to land intensification, pointing to a growth coalition strategy to incentivize development and maximize exchange value. And, our study clarifies that road construction, as car-centric public infrastructure, plays a significant role in this process.
Drawing from the growth machine framework, our finding of reciprocal asymmetry adds detail to how local land development strategies are unequally intertwined together with car-centric infrastructure projects, validating arguments that the ascendance of the private automobile has happened at the expense of public transit (e.g., Angelo and Wachsmuth 2020; Lee and Sener 2017; Mattioli et al. 2020). In other words, the government is socializing the risk and cost of road construction reciprocally in service of private land development and the private automobile. This dynamic reveals the asymmetric reciprocity that links the socialization of public cost to private, car-centric development. Although this gives us a socioecological example of an asymmetric reciprocal association, conceptually, it demonstrates that road construction, as infrastructure-led development, is a functional yet unbalanced component in the interdependent framework of the growth machine.
Footnotes
Funding
The author received no financial support for the research, authorship, and/or publication of this article.
Declaration of Conflicting Interests
The author declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
1
The National Land Cover Database Impervious Descriptor satellite imagery does not distinguish among the various road types; the raster file simply signals whether the cell pixel is covered by a road. We discuss this as a potential limitation in the conclusion of the study.
2
Among 3,203 census tracts, 2 have no roads and 3 have no (nonroad) developed land, representing 0.06 percent and 0.09 percent of the sample size, respectively.
3
Nine census tracts had no resident population between 2010 and 2024, representing about 0.28 percent of the sample size. These include wildlife refuges, industrial areas, and airports.
4
With the exception of the χ2 statistic, which is likely the result of the large sample size. As noted in
, the model fit indices are derived from structural equation models that do not include population density as an additional control. Although adding population density as a control does not change the substantive findings of the analysis, including the results of the asymmetry tests, it does seriously degrade model fit in terms of the Tucker-Lewis index, which decreases from 0.947 to 0.838 after including population density.
