Prioritizing Institutional and Technological Interventions for Engineering Project Credit-Risk Governance: An Integrated DEMATEL–ISM–MICMAC and Exploratory System Dynamics Approach
Abstract:
Engineering project credit-risk governance requires regulators and project participants to coordinate institutional controls with digital supervision capabilities. Yet decision-makers often lack a structured basis for identifying the factors that should receive priority and for judging how institutional and technological interventions may perform over time. This study investigates the causal structure of engineering project credit risk and examines the policy implications of alternative governance interventions. An online questionnaire collected 86 complete responses covering 33 predefined directional relationships among 13 factors organised under the Technology–Organization–Environment (TOE) framework. Full-precision mean scores were analysed using the Decision-Making Trial and Evaluation Laboratory (DEMATEL), Interpretive Structural Modeling (ISM), and Matrix of Cross-Impact Multiplications Applied to Classification (MICMAC). An exploratory system dynamics (SD) model was then used to compare the baseline, institutional-response, technology, and combined scenarios. Robustness was examined through alternative response coding, threshold sensitivity tests, and 1,000 bootstrap resamples. The results showed that insufficient credit verification by supervision units, environmental and resource compliance risk, and lagging credit-management methods were the three most prominent factors. The ISM analysis placed environmental and resource compliance risk at the root of the four-level hierarchy, while MICMAC classified four factors as independent drivers. All bootstrap samples retained the same three leading factors, and alternative coding preserved the complete prominence ranking (Spearman’s $\rho$ = 1.000). In the exploratory simulation, the technology intervention produced a credit index of 39.09 at time 20, compared with 4.06 under the baseline scenario. The institutional-response intervention showed no clear long-horizon advantage. At time 50, the combined scenario produced a value of 33.73, only slightly higher than the technology-only value of 33.42. These findings indicate that engineering project credit-risk governance should prioritise verifiable supervision, interoperable monitoring, and timely credit-management processes. The integrated framework provides a transparent basis for intervention prioritisation and lifecycle governance, while the simulation results should be interpreted as policy experiments rather than industry forecasts.1. Introduction
Engineering project credit risk arises when project participants fail to meet contractual, regulatory, quality, safety, schedule, cost, or environmental obligations. Institutions establish the formal and informal rules that shape such behaviour [1]. In construction management, institutional theory has been used to examine governance arrangements [2], environmentally responsible adoption [3], responsive multi-party regulation [4], and the interaction between regulation and digital technology [5]. Digital artefacts can carry or initiate institutional change [6], but the digital transformation of construction still depends on the alignment of technological capability, organisational readiness, and external support [7]. For regulators and project organisations, the central decision problem is therefore not whether digital technologies should be adopted in isolation, but how institutional arrangements and smart supervision technologies should be coordinated and prioritised to control credit risk throughout the project lifecycle.
Institutional theory provides a meta-theoretical basis for understanding the regulatory setting in which these decisions are made. North [1] defines institutions as the “rules of the game” in society, covering both formal institutions, such as laws and regulations, and informal institutions, such as industry norms and credit culture. Qiu and Chen [2] systematically reviewed the application of institutional theory in construction project management and distinguished its regulative, normative, and cognitive dimensions. Ayres and Braithwaite [4] proposed a tripartite regulatory framework involving government, enterprises, and public-interest groups, arguing that third-party participation can strengthen regulatory efficacy. However, institutional rules do not automatically produce effective supervision. They must be translated into operational responsibilities, verifiable procedures, information flows, and enforceable responses before they can function as practical regulatory measures [6]. This translation creates a resource-allocation problem because decision-makers must determine which institutional weaknesses require attention and which technological capabilities are most likely to support their implementation.
Recent engineering research demonstrates why institutional and technological measures need to be considered jointly. Studies of blockchain and smart contracts have addressed payment reliability and the reconfiguration of trust [8], [9]; digital twins and connected sensing have been used for real-time condition monitoring [10], [11]; and reviews of artificial intelligence have documented both its potential and its implementation limits in architecture, engineering, and construction [12], [13]. Research on contract claims has also linked delays, cost overruns, and quality failures to performance risk [14], while research on digitally enabled construction management has shown that governance and organisational coordination remain consequential in technologically advanced project environments [15]. These studies identify relevant technological and contractual mechanisms, but they do not establish which institutional, organisational, and technological constraints should be addressed first. Nor do they show whether an intervention that appears influential in a static causal structure will retain its relative value when feedback effects and implementation delays are considered.
The Technology–Organization–Environment (TOE) framework provides a conceptual structure for organising the conditions that influence technology adoption [16]. Related construction research has examined payment automation [17], the causes of schedule and cost overruns [18], dynamic modelling of site risks [19], safety prediction [20], governmental and industry applications of blockchain [21], and the mediating role of relational governance in digital project resilience [22]. Taken together, this literature indicates that smart supervision cannot operate as a stand-alone technical system. Its effectiveness depends on clearly assigned participant responsibilities, interoperable data, organisational capacity, and enforceable institutional rules. In this study, the TOE framework is used to organise the factor set before analysis. It does not determine the causal hierarchy, the relative importance of the factors, or the interventions that should receive priority; these are treated as empirical decision questions.
Smart supervision research examines how emerging technologies alter regulatory processes, information exchange, and the allocation of supervisory responsibilities. Wu et al. [21] reviewed the mechanisms through which blockchain can support efficiency and sustainability in construction. Liu and Wang [22] found that relational governance mediates the relationship between emerging digital technologies and construction-project resilience. Despite this progress, an imbalance remains between technical feasibility and the institutional and organisational work needed to make supervisory information usable in practice. A monitoring platform may collect extensive project data, for example, but those data have limited regulatory value when verification duties are unclear, systems are not interoperable, or identified defaults do not trigger a defined response. The management question is therefore how institutional and technological interventions can be combined within an auditable governance process rather than how a particular technology performs on its own.
Methodologically, integrated decision models have begun to connect structural judgements with dynamic behaviour. Decision-Making Trial and Evaluation Laboratory–System Dynamics (DEMATEL–SD) models have been used to examine project-cost factors [23] and schedule-delay controls [24]; Decision-Making Trial and Evaluation Laboratory–Interpretive Structural Modeling–System Dynamics (DEMATEL–ISM–SD) has been applied to smart-city carbon emissions [25]; and fuzzy DEMATEL–SD has been used to study social sustainability in infrastructure projects [26]. The unresolved methodological issue in engineering project credit-risk governance is not the simple combination of several established techniques. It is the construction of a traceable analytical chain from a clearly bounded set of questionnaire judgements to causal prominence, hierarchical position, driving dependence, and dynamic intervention behaviour. Such a chain is necessary if the results are to support defensible priority setting. It must also disclose which relationships were assessed, how the structural threshold was selected, which factors entered the dynamic model, and which factors were omitted.
Evidence for the factor set was drawn from research on collaborative environmental, social, and governance practices [27], barriers to construction robotics and digital twins [28], [29], governmental supervision technologies [30], construction risk analytics [31], [32], [33], and contractual risk allocation [34]. Additional evidence came from studies of supply-chain collaboration [35], barriers to technology adoption [36], [37], digital tendering [38], studies of China’s green-construction transition and construction-industry digital-transformation policies [39], [40], and data-governance requirements in sustainable smart cities [41]. These sources support the inclusion of institutional, organisational, technological, contractual, environmental, and project-execution factors within a common decision structure. The available study records document the iterative synthesis and consolidation of these factors and the construction of the final 13-factor instrument. They do not, however, contain a defensible count of the initial candidate pool or a formal record of factor deletion and expert consultation. The account of factor development is therefore restricted to procedures that can be verified from the available records rather than reconstructed retrospectively.
Against this background, the study addresses three questions: (1) Which factors have the greatest prominence and causal influence within the traceable direct-relation network? (2) What hierarchy and dependence structure results from a data-adaptive threshold, and how sensitive is that structure to alternative analytical settings? (3) How do selected institutional and technological interventions behave when feedback effects and implementation delays are examined in an exploratory dynamic model? The study makes three contributions. First, it develops a reproducible 33-link direct-relation matrix from full-precision questionnaire means and treats unassessed relationships explicitly as instrument-defined structural non-links. This provides a transparent basis for identifying intervention priorities rather than relying on an unrestricted or retrospectively completed influence matrix. Second, it separates the conceptual organisation supplied by TOE from the causal, hierarchical, and dependence structures obtained through DEMATEL, ISM, and Matrix of Cross-Impact Multiplications Applied to Classification (MICMAC), while testing the stability of the findings through alternative coding, threshold sensitivity analysis, and bootstrap resampling. Third, it connects the structural results to an exploratory system dynamics model limited to factors represented in the Vensim equations. The resulting scenarios are used as policy experiments for comparing institutional, technological, and combined interventions rather than as forecasts of actual industry performance.
2. Conceptual Background and Factor Construction
The smart governance of engineering project credit risk is treated as a management system that connects institutional rules, participant responsibilities, project data, and supervisory action. Institutional arrangements specify who must provide guarantees, record defaults, verify performance, impose consequences, and authorise credit repair. Smart supervision technologies affect the speed, coverage, and consistency with which these obligations can be monitored and coordinated. A digital platform without enforceable responsibilities may simply reproduce existing ambiguities in electronic form, while rules unsupported by interoperable monitoring may remain slow and episodic in practice. The relevant unit of analysis is therefore the combined governance arrangement rather than a stand-alone technology. From a decision-making perspective, the central task is to determine which institutional and technological constraints exert the greatest influence and which interventions should receive priority when supervisory resources are limited.
Engineering project credit risk also develops across the project lifecycle. At the entry stage, decisions concern participant qualifications, performance guarantees, and permissions for accessing and sharing data. During contract performance, project participants assume obligations relating to payment, schedule, cost, quality, safety, and environmental compliance. Construction monitoring determines whether observed signals are verified and recorded as credit events. Default handling then connects the available evidence to sanctions, guarantee claims, corrective measures, appeals, and credit repair. This lifecycle perspective provides the basis for interpreting the causal structure and comparing the institutional and technological interventions after the structural and dynamic results have been obtained.
Table 1 positions the present study as a decision-oriented integration of structural and dynamic analysis, with particular attention to traceability and clearly defined analytical roles. DEMATEL estimates factor prominence and net causal position, ISM arranges the retained relationships into hierarchical levels, and MICMAC classifies the factors according to their driving power and dependence. The SD model then examines the behaviour of selected feedback mechanisms under alternative intervention settings. These outputs address different parts of the decision problem and are used as complementary forms of evidence. They are not combined through an assumed one-to-one transfer of all 13 factors, and the DEMATEL prominence scores are not entered into the SD equations as numerical weights. The integration therefore supports intervention prioritisation without implying that the structural and dynamic models are mathematically identical.
| Study | Method | Domain | Distinguishing Feature |
|---|---|---|---|
| Alsugair et al. [23] | DEMATEL + SD | Project cost | Uses DEMATEL results to inform a dynamic model of project costs. |
| Ajayi and Chinda [24] | DEMATEL + SD | Schedule delay | Connects judgements about delay controls with schedule dynamics. |
| Cheng et al. [25] | DEMATEL + ISM + SD | Smart-city emissions | Links the hierarchy of emission drivers to dynamic emission pathways. |
| Rostamnezhad et al. [26] | Fuzzy DEMATEL + SD | Construction social sustainability | Combines influence judgements with the simulation of feedback relationships. |
| Present study | TOE + DEMATEL + ISM + MICMAC + SD | Engineering project credit-risk governance | Establishes a traceable path from a bounded 33-link instrument to factor prioritisation, hierarchical structuring, dependence classification, robustness assessment, and exploratory intervention comparison. |
The final factor set was developed through an iterative review of the literature and policy documents, followed by the consolidation of overlapping concepts and refinement of the questionnaire. The available study records document the final 13 factors (S1–S13) and their wording but do not contain a recoverable count of the initial candidate factors, a formal deletion log, respondent demographic information, or a documented Delphi process. The study therefore reports only those factor-development procedures that can be verified from the available records. Under the TOE framework, S2 and S10 are assigned to Technology, S4, S5, S6, S7, S8, and S9 to Organization, and S1, S3, S11, S12, and S13 to Environment. S7–S9 are placed under Organization because they originate in project execution and the performance of participating organisations, even though their consequences may extend beyond organisational boundaries. These assignments provide a conceptual structure for the subsequent analysis, while the empirical models determine factor prominence, causal position, hierarchical level, dependence, and relevance to intervention decisions. The final factor system, TOE classification, and operational definitions are summarised in Table 2.
| Code | Factor | TOE | Operational Definition |
| S1 | Lack of credit guarantee mechanism | Environment | Missing project- or market-level performance guarantees, default-risk sharing arrangements, and enforceable credit constraints. |
| S2 | Backward smart supervision technology | Technology | Insufficient platform capability, data connectivity, real-time monitoring, or analytical support for credit supervision. |
| S3 | Cross-subject credit information asymmetry | Environment | Fragmented or unequally accessible credit information among owners, contractors, supervision organisations, and regulators. |
| S4 | Deviant credit behavior of construction units | Organization | Owner or developer conduct that violates contractual commitments, payment duties, or regulatory requirements. |
| S5 | Deviant credit behavior of contractors | Organization | Contractor conduct that violates contract, quality, safety, schedule, payment, or reporting obligations. |
| S6 | Lack of credit verification by supervision units | Organization | Insufficient verification, recording, escalation, or follow-up of participant credit events by supervision organisations. |
| S7 | Quality and safety credit default risk | Organization | Quality and safety non-compliance generated within project execution and attributable to participant performance. |
| S8 | Schedule and cost credit default risk | Organization | Delay, overrun, payment, or cost-control failures arising during contractual delivery. |
| S9 | Environmental and resource compliance credit risk | Organization | Failure to meet environmental, resource-use, or related compliance obligations during project delivery. |
| S10 | Lagging credit management methods | Technology | Manual, fragmented, or delayed credit-management procedures and workflows, distinct from the digital capability itself. |
| S11 | Absence of credit regulatory regulations | Environment | Gaps in external legal rules, regulatory standards, and statutory enforcement applicable to engineering credit. |
| S12 | Lack of integrity culture | Environment | Weak shared norms and expectations supporting honest performance and credible reporting. |
| S13 | Blocked social credit supervision channels | Environment | Limited public or third-party channels for reporting, accessing, and following up engineering credit information. |
Several conceptual boundaries are important when interpreting the factor system. S1 concerns project- or market-level assurance arrangements, including performance guarantees, default-risk sharing, and the financial or contractual cost of non-performance. S11 instead concerns the external legal and regulatory framework that defines applicable standards and enforcement authority. Similarly, S2 refers to the capability of the digital infrastructure, including platforms, connectivity, real-time monitoring, and analytical functions. S10 refers to the procedures and workflows through which credit information is collected, verified, escalated, and used. It is placed under Technology because the factor, as defined in the instrument, concerns the extent to which these workflows remain manual, fragmented, and delayed rather than the broader allocation of organisational responsibilities.
The TOE categories are used to organise the factors conceptually; they do not predetermine the empirical relationships among them. An environmental factor such as S11 may therefore appear as autonomous if it has limited driving power and dependence within the assessed network. Conversely, an organisational risk such as S9 may occupy a deep causal position when its relationships with other factors give it substantial driving power. These empirical positions do not contradict the TOE assignments and should not be altered to produce a more conventional pattern. Instead, they show how the factors function within the questionnaire-defined decision system examined in this study.
3. Materials and Methods
Figure 1 summarises the analytical design. The DEMATEL method was used to quantify the direct and propagated influences among the identified factors [42]. The Maximum Mean De-Entropy (MMDE) algorithm was then applied to determine the threshold for retaining relationships [43]. ISM was used to arrange the retained relationships into hierarchical levels, while the MICMAC was used to classify the factors according to their driving power and dependence. Finally, a scenario-based system dynamics model was developed using stock-flow and feedback principles to compare selected institutional-response and technological interventions [44], [45], [46].
The methods addressed successive parts of the decision problem. DEMATEL identified the factors requiring priority attention, ISM represented their hierarchical relationships, and MICMAC distinguished potential drivers from dependent outcomes. The system dynamics model then examined how selected intervention levers behaved over time under the stated equations and parameter settings. The structural and dynamic analyses were connected through an explicit factor-mapping procedure rather than by transferring all 13 factors directly into the simulation.

The questionnaire was administered exclusively online through Questionnaire Star (Wenjuanxing), a web-based survey platform. The archived record contained 86 submissions and 33 directional influence items, although the number of questionnaires or invitations distributed could not be recovered. All 86 records were complete, none were excluded, and all were retained for analysis. No exclusion criteria were applied because the archive contained no verifiable record of mechanical-response screening. Each item used a four-level scale, where 1 represented slight influence, 2 moderate influence, 3 significant influence, and 4 very strong influence. The instrument did not include 0 or a separate “no direct influence” option.
The anonymous dataset did not retain information on occupation, years of experience, organisational affiliation, main project type, engineering credit-supervision experience, or geographic group. Respondent-profile subgroup analyses could therefore not be conducted. The verified questionnaire-administration information and the available sample characteristics are summarised in Table 3.
Item | Verified Value |
|---|---|
Administration mode | Online survey via Questionnaire Star (Wenjuanxing) |
Questionnaires/invitations distributed | Not recoverable from the archived Wenjuanxing record |
Submissions received | 86 |
Submissions excluded | 0 |
Exclusion criteria applied | None; all 86 archived submissions were complete, and no verifiable mechanical-response screening record was available |
Final analytical sample | 86 |
Directional influence items | 33 |
Influence scale | 1–4; no zero/no-influence option |
Missing item responses | 0 |
Respondent profile/geography | Not collected in the anonymous survey dataset |
The absence of profile variables prevents post hoc comparisons of DEMATEL rankings across respondent groups. Bootstrap analysis can assess the stability of the aggregate rankings under resampling, but it cannot replace comparisons between professional, organisational, experiential, or geographic groups. Future questionnaires should therefore collect non-identifying information on respondent roles and experience under an explicit consent and data-minimisation protocol.
Because respondent occupations and professional experience were not retained in the archived dataset, the sample cannot be retrospectively verified as an expert panel. The resulting influence scores are therefore interpreted as aggregate questionnaire judgements rather than as formally validated expert consensus.
The questionnaire completed by the 86 respondents assessed 33 of the 156 possible off-diagonal ordered pairs among the 13 factors. The direct-relation matrix $A$ was constructed from these assessed relationships. Each of the 33 corresponding cells was assigned the full-precision arithmetic mean of the 86 item scores. The remaining 123 off-diagonal cells were set to zero because these relationships had not been evaluated by the questionnaire. The archived study materials did not contain a formal record explaining why these 33 directional relationships were selected from the 156 possible ordered pairs. Accordingly, the present analysis does not retrospectively reconstruct a link-selection procedure that cannot be verified. The resulting causal network should therefore be interpreted as conditional on the coverage of the original questionnaire instrument.
These zero entries were treated as instrument-defined structural non-links for the purpose of matrix construction. They do not represent zero ratings selected by respondents, nor do they demonstrate that the corresponding substantive relationships are absent. This distinction preserves the boundary of the questionnaire-defined network and prevents unassessed relationships from being interpreted as observed judgements. The 33 directional relationships actually evaluated by the questionnaire, together with their full-precision means and response distributions, are reported in Table 4.
| Item | Direction | Mean | Standard Deviation | $\boldsymbol{n(1)}$ | $\boldsymbol{n(2)}$ | $\boldsymbol{n(3)}$ | $\boldsymbol{n(4)}$ |
|---|---|---|---|---|---|---|---|
| 1 | S1 to S9 | 3.453488 | 0.849318 | 6 | 2 | 25 | 53 |
| 2 | S1 to S10 | 3.302326 | 0.736762 | 3 | 5 | 41 | 37 |
| 3 | S2 to S6 | 3.232558 | 0.890134 | 7 | 5 | 35 | 39 |
| 4 | S2 to S7 | 3.395349 | 0.786340 | 4 | 4 | 32 | 46 |
| 5 | S2 to S13 | 3.209302 | 0.704878 | 3 | 5 | 49 | 29 |
| 6 | S3 to S13 | 3.267442 | 0.758085 | 4 | 4 | 43 | 35 |
| 7 | S4 to S3 | 3.360465 | 0.780841 | 4 | 4 | 35 | 43 |
| 8 | S4 to S5 | 3.267442 | 0.726384 | 3 | 5 | 44 | 34 |
| 9 | S4 to S10 | 3.267442 | 0.788512 | 5 | 3 | 42 | 36 |
| 10 | S6 to S1 | 3.255814 | 0.689178 | 2 | 6 | 46 | 32 |
| 11 | S6 to S2 | 3.279070 | 0.697072 | 2 | 6 | 44 | 34 |
| 12 | S6 to S3 | 3.406977 | 0.757362 | 3 | 5 | 32 | 46 |
| 13 | S6 to S4 | 3.244186 | 0.750467 | 4 | 4 | 45 | 33 |
| 14 | S6 to S5 | 3.337209 | 0.776273 | 4 | 4 | 37 | 41 |
| 15 | S6 to S7 | 3.186047 | 0.759437 | 5 | 3 | 49 | 29 |
| 16 | S6 to S10 | 3.348837 | 0.715664 | 2 | 6 | 38 | 40 |
| 17 | S7 to S8 | 3.290698 | 0.700693 | 2 | 6 | 43 | 35 |
| 18 | S7 to S10 | 3.302326 | 0.704102 | 2 | 6 | 42 | 36 |
| 19 | S7 to S12 | 3.255814 | 0.842769 | 7 | 1 | 41 | 37 |
| 20 | S8 to S1 | 3.337209 | 0.806014 | 5 | 3 | 36 | 42 |
| 21 | S8 to S2 | 3.325581 | 0.742680 | 3 | 5 | 39 | 39 |
| 22 | S8 to S5 | 3.325581 | 0.710292 | 2 | 6 | 40 | 38 |
| 23 | S8 to S10 | 3.302326 | 0.768034 | 4 | 4 | 40 | 38 |
| 24 | S9 to S2 | 3.197674 | 0.699912 | 3 | 5 | 50 | 28 |
| 25 | S9 to S4 | 3.290698 | 0.733505 | 3 | 5 | 42 | 36 |
| 26 | S9 to S6 | 3.348837 | 0.715664 | 2 | 6 | 38 | 40 |
| 27 | S9 to S8 | 3.302326 | 0.827040 | 6 | 2 | 38 | 40 |
| 28 | S10 to S1 | 3.348837 | 0.747819 | 3 | 5 | 37 | 41 |
| 29 | S10 to S5 | 3.337209 | 0.679281 | 1 | 7 | 40 | 38 |
| 30 | S11 to S4 | 3.395349 | 0.786340 | 4 | 4 | 32 | 46 |
| 31 | S12 to S9 | 3.325581 | 0.710292 | 2 | 6 | 40 | 38 |
| 32 | S12 to S11 | 3.279070 | 0.791888 | 5 | 3 | 41 | 37 |
| 33 | S13 to S7 | 3.279070 | 0.791888 | 5 | 3 | 41 | 37 |
Because the forced-choice scale did not include a no-influence option, every assessed relationship necessarily received a positive score. The resulting means should therefore not be interpreted as evidence of the frequency with which respondents perceived no direct relationship. Restricting the matrix to the 33 relationships actually assessed prevented the unasked ordered pairs from being assigned artificial positive values. At the same time, it made the resulting network dependent on the coverage of the original instrument. Both consequences are recognised as limitations of the analysis.
Let $a_{ij}$ denote the mean reported influence of factor $i$ on factor $j$ for an evaluated directional item, and let $a_{ij}$ = 0 for an unasked structural non-link. All calculations use float64 values before display rounding.
where $A$ is the traceable direct-relation matrix, $s$ is its maximum row sum (23.0581395349), and $B$ is the normalised direct-relation matrix.
where $T$ is the total-relation matrix and $I$ is the 13 $\times$ 13 identity matrix. The spectral radius of $B$ is 0.3410762589, so the inverse exists for the observed matrix.
where $t_{ij}$ is an element of $T$, $D_i$ is the total influence exerted by factor $i$, and $R_i$ is the total influence received by factor $i$.
where $C_i$ is prominence (centrality) and $N_i$ is net relation. $N_i>$ 0 identifies a net cause and $N_i<$ 0 a net effect.
All substantive analysis was based on the full-precision traceable matrix, and values were rounded only for display. To distinguish the effect of numerical rounding from changes in network topology, the same sparse 33-link matrix was recalculated after the 33 means had been conventionally rounded to one decimal place. The legacy dense matrix was not used as the rounding comparator.
The MMDE procedure is applied to the non-diagonal total-relation values and selects $\lambda$ = 0.168627236856. A directed edge is retained when $t_{ij}$ is at least $\lambda$. The diagonal is set to one for reachability analysis, and transitive closure yields the reachability matrix $M$. For factor $i$, the reachability set contains factors reachable from $i$, while the antecedent set contains factors that can reach $i$. A factor enters the current surface level when the intersection of those sets equals its reachability set; identified factors are removed iteratively until all levels are assigned.
where $m_{ij}$ is an off-diagonal element of the thresholded adjacency matrix and $\lambda$ is the MMDE threshold. Diagonal elements are set to one before transitive closure.
where $P_i$ is MICMAC driving power and $Q_i$ is dependence power in the final reachability matrix. The common classification cutoff is their mean, 2.461538.
Four checks were conducted to assess measurement and model sensitivity. First, the same 33-link topology was recalculated using means rounded to one decimal place rather than full-precision means. The analysis considered both a newly estimated MMDE threshold and the fixed threshold obtained from the full-precision matrix.
Second, the original 1–4 scores were retained as the primary coding, while an alternative coding test mapped scores of 1, 2, 3, and 4 to 0, 1, 2, and 3, respectively. This recoding was used only as a sensitivity test and does not imply that the original questionnaire contained a zero or no-influence response option.
Third, 1,000 bootstrap samples of the 86 questionnaire records were drawn with replacement using seed 12202. DEMATEL prominence rankings and cause-effect classifications were recalculated for every bootstrap sample. This procedure examined whether the principal factor rankings depended on particular combinations of respondents within the available sample.
Fourth, the ISM structure was evaluated at MMDE – 10%, the selected MMDE threshold, MMDE + 10%, and the absolute thresholds of 0.38, 0.40, and 0.42. The absolute thresholds provided a scale-stress test and were deliberately high relative to the values in the sparse total-relation matrix. Together, these checks examined the stability of factor prioritisation, causal classification, and hierarchical structure under alternative numerical treatments.
The Vensim .mdl model defined a stock variable named the industry credit index and a set of feedback relationships involving regulatory pressure, guarantee completeness, smart supervision technology, information symmetry, quality and safety default risk, dishonesty cost, and the rates of credit improvement and decline. The industry credit index was treated as a dimensionless, model-generated representation of engineering-sector credit performance rather than an indicator estimated from observed industry data.
The equations were simulated using Euler integration with a time step of 0.25 over a simulation horizon of 0–50. One unit represented a simulation or policy interval and was not assumed to correspond to a calendar year. Times 5, 10, and 20 were selected as practical comparison points. Times 30–50 were used only to examine the long-horizon behaviour of the modelled trajectories and were not interpreted as forecasts of actual industry conditions.
Four scenarios were examined using the same model structure and equations. The baseline scenario retained the original parameter values. The institutional-response intervention reduced the regulatory response delay from 2.0 to 1.0. The technology intervention increased the technology investment coefficient from 1.0 to 2.0. The combined intervention applied both parameter changes simultaneously. Referring to the second scenario as an institutional-response intervention reflects its specific modelled mechanism and avoids implying that the simulation represents a comprehensive reform of the wider institutional system.
The scenarios were designed to compare the relative dynamic behaviour of selected governance levers under the stated model assumptions. The model was not fitted to historical observations of the engineering-sector credit index and should therefore be treated as an exploratory, scenario-based policy model. Learning-by-doing, diffusion, and feedback modelling supplied the conceptual basis for the dynamic relationships [44], [45], [46], but they did not convert the assumed equations or parameter values into empirically observed relationships. Accordingly, the scenario results were used to examine intervention timing, relative performance, and trajectory stability rather than to predict future industry credit levels.
4. Structural Results
Table 5 reports the full-precision DEMATEL results, including prominence, net relation, cause–effect classification, ranking stability under one-decimal rounding, and bootstrap rank intervals.
Factor | D | R | D + R | D $-$ R | Group | Full-Precision Rank | One-Decimal Rank* | Bootstrap 95% Rank Interval |
|---|---|---|---|---|---|---|---|---|
S6 | 1.420 | 0.434 | 1.854 | 0.986 | Cause | 1 | 1 | 1–1 |
S9 | 1.069 | 0.435 | 1.504 | 0.634 | Cause | 2 | 2 | 2–3 |
S10 | 0.363 | 1.119 | 1.482 | $-$0.756 | Effect | 3 | 3 | 2–3 |
S2 | 0.756 | 0.611 | 1.367 | 0.145 | Cause | 4 | 4 | 4–4 |
S7 | 0.662 | 0.638 | 1.300 | 0.023 | Cause | 5 | 5 | 5–5 |
S8 | 0.811 | 0.439 | 1.250 | 0.371 | Cause | 6 | 6 | 6–6 |
S1 | 0.505 | 0.719 | 1.224 | $-$0.213 | Effect | 7 | 7 | 7–7 |
S4 | 0.506 | 0.580 | 1.086 | $-$0.073 | Effect | 8 | 8 | 8–8 |
S5 | 0.000 | 0.946 | 0.946 | $-$0.946 | Effect | 9 | 9 | 9–9 |
S12 | 0.472 | 0.231 | 0.704 | 0.241 | Cause | 10 | 10 | 10–10 |
S13 | 0.236 | 0.429 | 0.665 | $-$0.192 | Effect | 11 | 11 | 11–11 |
S3 | 0.175 | 0.442 | 0.617 | $-$0.267 | Effect | 12 | 12 | 12–12 |
S11 | 0.222 | 0.175 | 0.397 | 0.047 | Cause | 13 | 13 | 13–13 |
The full-precision prominence order was S6 $>$ S9 $>$ S10 $>$ S2 $>$ S7 $>$ S8 $>$ S1 $>$ S4 $>$ S5 $>$ S12 $>$ S13 $>$ S3 $>$ S11. S6 had the highest prominence value (1.854) and the largest positive net relation (0.986), identifying insufficient credit verification by supervision units as the most prominent net driver within the assessed network. S9 ranked second in prominence and was also classified as a net cause. S10 ranked third but had a strongly negative net relation ($-$0.756), indicating that lagging credit-management methods primarily received propagated influence from other factors.
The combination of high prominence and positive net influence placed S6 and S9 at the head of the initial intervention-priority set. This structural position does not by itself demonstrate that an intervention targeting either factor will produce the largest practical effect. The dynamic implications of selected intervention levers were therefore examined separately in the exploratory system dynamics model. All subsequent structural interpretation was based on the full-precision ordering.
The rounding-only analysis preserved all 13 prominence ranks, with Spearman’s $\rho$ = 1.000, and retained every cause-effect assignment. Re-estimating MMDE from the one-decimal matrix increased the threshold from 0.168627 to 0.183536, reduced the number of retained edges from 13 to 5, and produced a three-level hierarchy in which S7 and S8 moved to the surface level. When the full-precision MMDE threshold was held constant, the one-decimal matrix retained 13 edges and reproduced the original four-level hierarchy. The prominence rankings and net-relation signs were therefore insensitive to rounding. The change in the hierarchy resulted from re-estimation of the threshold rather than from any change in the underlying 33-link topology.
The prominence and net-relation positions of all 13 factors are visualised in Figure 2.

Seven factors were classified as net causes: S2, S6, S7, S8, S9, S11, and S12. The remaining six factors—S1, S3, S4, S5, S10, and S13—were classified as net effects. A cause designation refers to the factor’s net position in the assessed influence network; it does not imply moral responsibility, universal causal primacy, or proven intervention effectiveness. S11, for example, had a slightly positive net relation but the lowest prominence of all 13 factors. Its cause status must therefore be interpreted together with its limited connectivity in the threshold and MICMAC results and should not be used alone to support broad regulatory claims.
The effects of threshold selection and numerical rounding on the ISM structure are summarised in Table 6.
Threshold | Value | Edges | Isol. $\boldsymbol{n}$ | Isolated factors | Levels | Surface-to-root partition |
|---|---|---|---|---|---|---|
MMDE $-$10% | 0.151765 | 19 | 3 | S11, S12, S13 | 4 | L1:S3,S5,S10,S11,S12,S13 | L2:S4,S7 | L3:S2 | L4:S1,S6,S8,S9 |
MMDE | 0.168627 | 13 | 3 | S11, S12, S13 | 4 | L1:S1,S3,S4,S5,S10,S11,S12,S13 | L2:S7,S8 | L3:S2,S6 | L4:S9 |
1 d.p., own MMDE | 0.183536 | 5 | 6 | S3, S4, S8, S11, S12, S13 | 3 | L1:S1,S3,S4,S5,S7,S8,S10,S11,S12,S13 | L2:S2,S6 | L3:S9 |
1 d.p., fixed MMDE | 0.168627 | 13 | 3 | S11, S12, S13 | 4 | L1:S1,S3,S4,S5,S10,S11,S12,S13 | L2:S7,S8 | L3:S2,S6 | L4:S9 |
MMDE $+$10% | 0.185490 | 4 | 7 | S1, S3, S4, S8, S11, S12, S13 | 3 | L1:S1,S3,S4,S5,S7,S8,S10,S11,S12,S13 | L2:S2,S6 | L3:S9 |
Absolute 0.38 | 0.380000 | 0 | 13 | S1–S13 | 1 | L1:S1–S13 |
Absolute 0.40 | 0.400000 | 0 | 13 | S1–S13 | 1 | L1:S1–S13 |
Absolute 0.42 | 0.420000 | 0 | 13 | S1–S13 | 1 | L1:S1–S13 |
At the selected MMDE threshold, 13 relationships were retained. The four strongest were S6 to S10 (0.2248), S6 to S5 (0.2064), S2 to S7 (0.1919), and S9 to S2 (0.1867). These relationships connected supervision verification, smart supervision technology, execution risk, and credit-management procedures within the main structural pathway.
At MMDE – 10%, 19 relationships remained, while S11, S12, and S13 were isolated. At MMDE + 10%, only four relationships remained and seven factors were isolated. S11, S12, and S13 remained isolated under every tested threshold. The three absolute thresholds removed all off-diagonal relationships and consequently placed every factor in a single surface level.
Under the formal ISM procedure, an isolated node is assigned to the surface partition because its reachability set contains only itself. A Level 1 position therefore does not necessarily identify an operational outcome, a highly dependent factor, or a factor of little managerial importance. It only describes the position produced by the retained relationships in the questionnaire-defined network. The rapid reduction in retained edges at higher thresholds also shows that the detailed hierarchy was more sensitive to threshold selection than the DEMATEL prominence ranking.
Figure 3 further illustrates how the number of retained edges, isolated factors, and ISM levels changed across the tested thresholds.

Table 7 summarises the ISM level, driving power, dependence power, and MICMAC classification of each factor.
| Code | Factor | TOE | ISM | Driving | Depend. | MICMAC Class |
| S1 | Lack of credit guarantee mechanism | Environment | Level 1 | 1 | 4 | Dependent |
| S2 | Backward smart supervision technology | Technology | Level 3 | 3 | 2 | Independent |
| S3 | Cross-subject credit information asymmetry | Environment | Level 1 | 1 | 3 | Dependent |
| S4 | Deviant credit behavior of construction units | Organization | Level 1 | 1 | 2 | Autonomous |
| S5 | Deviant credit behavior of contractors | Organization | Level 1 | 1 | 4 | Dependent |
| S6 | Lack of credit verification by supervision units | Organization | Level 3 | 6 | 2 | Independent |
| S7 | Quality and safety credit default risk | Organization | Level 2 | 2 | 4 | Dependent |
| S8 | Schedule and cost credit default risk | Organization | Level 2 | 4 | 1 | Independent |
| S9 | Environmental and resource compliance credit risk | Organization | Level 4 | 9 | 1 | Independent |
| S10 | Lagging credit management methods | Technology | Level 1 | 1 | 6 | Dependent |
| S11 | Absence of credit regulatory regulations | Environment | Level 1 | 1 | 1 | Autonomous |
| S12 | Lack of integrity culture | Environment | Level 1 | 1 | 1 | Autonomous |
| S13 | Blocked social credit supervision channels | Environment | Level 1 | 1 | 1 | Autonomous |
The ISM analysis produced four levels extending from surface to root. Level 1 contained S1, S3, S4, S5, S10, S11, S12, and S13. Level 2 contained S7 and S8, while Level 3 contained S2 and S6. S9 was the sole factor at Level 4.
Within the assessed network, the hierarchy indicated a pathway in which environmental and resource compliance credit risk was connected to smart supervision technology and verification capacity. These factors were subsequently linked to project-execution risks and several observable governance deficiencies. The position of S9 at the root of the hierarchy applied specifically to the questionnaire-defined 33-link network. It should not be interpreted as evidence that S9 is the universal root cause of engineering project credit risk in all institutional or project settings.
Figure 4 presents the transitive reduction of the MMDE-thresholded directed acyclic graph rather than all 13 retained relationships. The direct S6-to-S10 edge was omitted because S6 and S10 were also connected through the S6-to-S7-to-S10 pathway. Removing this redundant edge left 12 irreducible relationships. Irreducible edges extending across non-adjacent levels were retained.

The MICMAC analysis classified S2, S6, S8, and S9 as independent factors; S1, S3, S5, S7, and S10 as dependent factors; and S4, S11, S12, and S13 as autonomous factors. No linkage factor appeared at the mean cutoff. In MICMAC terminology, “independent” refers to factors with relatively high driving power and low dependence; it does not mean statistical independence.
The combined ISM and MICMAC results identified S9 as both a root-level factor and an independent driver, while S6 occupied a deep level and also had substantial driving power. S2 and S8 were likewise classified as independent, although their prominence values were lower than those of S6 and S9. These differences show why factor priority should not be determined from a single indicator. Prominence describes overall involvement in the network, net relation distinguishes causes from effects, ISM represents hierarchical position, and MICMAC describes driving power and dependence.
S11 illustrates the distinction between conceptual classification and empirical connectivity. TOE classified S11 as an environmental condition because it concerns the external regulatory framework, whereas MICMAC classified it as autonomous because it had low driving power and low dependence in the assessed network. The two classifications therefore answer different questions and are not contradictory.
The resulting driving–dependence positions and MICMAC classifications are shown in Figure 5.

The alternative 0–3 coding produced the same prominence order as the primary 1–4 coding, with Spearman’s $\rho$ = 1.000. The factor ranking was therefore unchanged by a uniform one-point shift in the response coding. This test could not, however, reproduce a genuine no-influence response because such an option had not been included in the original questionnaire.
Across the 1,000 bootstrap resamples, S6 ranked first in every draw. S9 and S10 occupied ranks 2 and 3 in every draw, and each of the three factors appeared in the top three with a frequency of 1.000. Cause-effect membership remained unchanged for all factors except S7, which was classified as a net cause in 99.6% of the bootstrap draws.
These findings indicated that the principal prominence ordering was stable under resampling of the observed responses. The narrow rank intervals should nevertheless be interpreted in light of the fixed 33-link topology. The bootstrap procedure represented sampling variation in the collected item scores, but it did not capture uncertainty about the 123 directional relationships that were not assessed by the questionnaire. The robustness results therefore support the stability of factor prioritisation within the observed network rather than the completeness of the network itself.
5. Exploratory System Dynamics Simulation
Table 8 reports how each factor from the DEMATEL–ISM–MICMAC analysis was represented in the system dynamics model. S1–S7 and S11 were included either directly, through an explicitly identified composite variable, or as an exogenous model input. S8, S9, S10, S12, and S13 were not simulated because the source model contained no corresponding equations. The mapping did not transfer DEMATEL prominence values into the simulation as coefficients, nor did it assume that the static hierarchy had been dynamically validated. Instead, the system dynamics model examined a narrower set of feedback relationships associated with regulatory response and technology investment.
| Code | SD Variable | Type | Representation | Sign | Equation/Parameter |
| S1 | Credit guarantee completeness; dishonesty cost | Auxiliary/composite | Direct/composite | $+$ | Guarantee smooth; dishonesty-cost equation |
| S2 | Smart-supervision technology level | Auxiliary | Direct | $+$ | Technology smooth; investment coefficient |
| S3 | Information symmetry degree | Auxiliary | Direct | $+$ | Technology-adoption lookup |
| S4 | Developer misconduct index | Exogenous/step | Direct | $-$ | Credit-decline step at time 5 |
| S5 | Contractor misconduct index | Exogenous | Direct | $-$ | Quality/safety risk equation |
| S6 | Supervision efficacy decay coefficient | Parameter | Composite proxy | $-$ | Quality/safety risk equation |
| S7 | Quality and safety default risk | Auxiliary | Direct | $-$ | Risk and decline equations |
| S8 | None | Not included | Not simulated | n/a | No source-model equation |
| S9 | None | Not included | Not simulated | n/a | No source-model equation |
| S10 | None | Not included | Not simulated | n/a | No source-model equation |
| S11 | Regulation completeness | Exogenous | Direct | $+$ | Dishonesty-cost conditional |
| S12 | None | Not included | Not simulated | n/a | No source-model equation |
| S13 | None | Not included | Not simulated | n/a | No source-model equation |
This distinction is important for interpreting the intervention results. S9 occupied the root position in the ISM hierarchy, and S10 ranked third in DEMATEL prominence, but neither factor was represented in the dynamic model. S6 was represented only indirectly through a composite proxy for the decay of supervision efficacy. Accordingly, the simulation did not test whether interventions directed specifically at S9 or S10 would outperform the modelled technology-investment and regulatory-response interventions. It provided a separate exploratory comparison of the intervention levers available in the source model.
Figure 6 shows the stock-flow structure and the principal feedback relationships used in the exploratory system dynamics model.

The industry credit index, $CI$, was specified as a dimensionless endogenous stock with an initial value of 50. This value represented a reference state within the model rather than an observed or officially defined industry credit score. The stock was not bounded because the source equation integrated improvement minus decline without imposing a lower limit. A negative value therefore represented a modelled credit deficit relative to the reference state; it did not represent a negative administrative credit score. This formulation exposed the effect of adverse feedback within the model but restricted direct interpretation of the numerical values in real-world terms.
The core equations and the interpretation of the main system dynamics variables are provided in Table 9.
| Variable | Source-Model Equation | Meaning |
|---|---|---|
| Credit index $CI$ | $\mathrm{INTEG}(\text{improvement} - \text{decline},50)$ | Dimensionless stock; negative values denote deficit relative to reference. |
| Improvement | $0.5 \times (\text{guarantee})^{0.8} \times (1 + \text{technology gain}) \times \text{coordination}$ | Positive inflow to $CI$. |
| Decline | $0.3 \times (\text{quality risk})^{0.8} + 0.2 \times \max(0,100 - \text{dishonesty cost})^{0.8} + \mathrm{STEP}(0.1 \times \text{developer misconduct},5)$ | Negative outflow from $CI$. |
| Regulatory pressure | $\mathrm{SMOOTHI}(0.6 \times \text{credit gap} + 0.4 \times \text{crisis frequency}, \text{response delay},50)$ | Lagged response to credit deficit and crisis frequency. |
| Quality risk | $50 \times (1 - \text{information symmetry}/100)^{0.8} \times (1 + 0.3 \times \text{contractor misconduct}/100) \times \text{supervision decay}$ | Composite quality/safety default risk. |
| Dishonesty cost | $\min(100, \text{guarantee} \times 1.2 \times \text{regulation multiplier})$ | Constraint on default-related decline; multiplier is 1.5 only when regulation completeness $> 50$. |
| Technology level | $\mathrm{SMOOTH}(\text{regulatory pressure} \times \text{technology investment} \times 0.6, \text{technology delay})$ | Lagged smart-supervision capability. |
| Information symmetry | $\mathrm{TECH\_ADOPT\_LOOKUP}(\text{technology level} \times 0.7/50)$ | Nonlinear conversion of technology into information symmetry. |
Table 10 reports the baseline values, uncertainty ranges, and provenance of the principal model parameters.
Parameter | Baseline | Range/Treatment | Basis | Calibration Target |
Improvement rate coefficient | 0.5 | $+/-$20% MC | Archived source-model value (retained) | No historical fit |
Institution-technology coordination | 1.0 | Fixed | Archived source-model value (retained) | No historical fit |
Decline quality weight | 0.3 | $+/-$20% MC | Archived source-model value (retained) | No historical fit |
Decline dishonesty weight | 0.2 | Fixed | Archived source-model value (retained) | No historical fit |
Regulatory response delay | 2.0 | 1.6–2.4 scan; $+/-$20% MC | Archived source-model value (retained) | No historical fit |
Technology response delay | 1.5 | Fixed | Archived source-model value (retained) | No historical fit |
Regulation credit/crisis weights | 0.6/0.4 | Fixed | Archived structural weights (retained) | No historical fit |
Information symmetry coefficient | 0.7 | $+/-$20% MC | Archived source-model value (retained) | No historical fit |
Baseline risk level | 50 | Fixed | Archived reference value (retained) | No historical fit |
Dishonesty cost coefficient | 1.2 | Fixed | Archived source-model value (retained) | No historical fit |
Technology investment coefficient | 1.0 | 0.8–1.2 scan; $+/-$20% MC | Archived source-model value (retained) | No historical fit |
Crisis quality/default weights | 0.4/0.3 | Fixed | Archived structural weights (retained) | No historical fit |
Regulation completeness | 50 | Fixed | Archived reference value (retained) | No historical fit |
Developer misconduct index | 50 | Fixed; step at time 5 | Archived reference value (retained) | No historical fit |
Contractor misconduct index | 50 | Fixed | Archived reference value (retained) | No historical fit |
Supervision efficacy decay | 1.0 | Fixed | Archived reference value (retained) | No historical fit |
Nonlinearity exponent | 0.8 | Fixed | Archived source-model value (retained) | No historical fit |
Euler time step | 0.25 | Fixed | Archived Vensim setting (retained) | Equation reproduction |
The equations were simulated using Euler integration with a time step of 0.25 over a horizon of 0–50. One unit represented a simulation or policy interval and was not assumed to correspond to a calendar year. Times 5, 10, and 20 were used as practical comparison points, while time 50 was used only as a long-horizon diagnostic. None of these points should be interpreted as a forecast date.
Under the baseline settings, $CI$ started at 50 and decreased to 42.02 at time 5 and 4.98 at time 10. It reached its lowest reported value of $-$5.27 at time 14 before recovering to 4.06 at time 20. The trajectory subsequently approached a low oscillatory regime, reaching 6.38 at time 30 and 5.64 at time 50.
Figure 7 presents the corresponding trajectories of regulatory pressure, credit guarantee completeness, technology level, and $CI$. The temporary negative value resulted from the unbounded stock equation and represented a deficit relative to the model reference state. It should not be interpreted as an observed engineering-industry credit score.

The technology-investment scan showed a monotonic pattern over the range of 0.8–1.2. At time 20, CI increased from $-$7.31 at a coefficient of 0.8 to $-$1.54 at 0.9, 4.06 at 1.0, 9.05 at 1.1, and 13.54 at 1.2. Within this local range, the technology-investment coefficient affected both the timing of recovery and the value reached at the later simulation checkpoint. The corresponding sensitivity trajectories are shown in Figure 8.

The regulatory-response-delay scan produced the opposite time-20 pattern. The corresponding $CI$ values were 5.39, 4.91, 4.06, 2.84, and 1.29 for delays of 1.6, 1.8, 2.0, 2.2, and 2.4, respectively. A shorter response delay produced a higher time-20 value within this local scan. However, the transient trajectory was nonlinear, and the separate response-delay intervention did not produce a higher time-50 value than the baseline scenario. Figure 9 shows the corresponding sensitivity of $CI$ to changes in regulatory response delay.

A 1,000-draw Monte Carlo analysis was conducted using independent triangular distributions for five parameters, with the lower bound, mode, and upper bound set at $-$20%, the baseline value, and +20%, respectively. At time 18, the mean $CI$ was $-$0.30 and the median was 0.50. The empirical 95% interval ranged from $-$14.38 to 11.93, and 51.7% of the simulated values were positive. Because this interval crossed zero, the timing of the baseline recovery should not be treated as precise. The resulting distribution of $CI$ at simulation time 18 is shown in Figure 10.

At time 18, the improvement-rate coefficient had the strongest rank correlation with $CI$ ($\rho$ = 0.702), followed by the technology-investment coefficient ($\rho$ = 0.560). The remaining correlations were $-$0.224 for the decline-quality weight, $-$0.172 for regulatory response delay, and 0.137 for the information-symmetry coefficient. These values describe model sensitivity under the stated uncertainty design and are not empirical estimates of policy or causal effects. Figure 11 presents the Monte Carlo trajectory envelope together with the rank-based parameter sensitivity results.

Four scenarios were compared using the same equations and initial conditions. The baseline scenario retained all original parameter values. The institutional-response intervention reduced regulatory response delay from 2.0 to 1.0. The technology intervention increased the technology-investment coefficient from 1.0 to 2.0. The combined intervention applied both parameter changes.
The technology-investment value of 2.0 represented a discrete policy experiment and lay outside the local sensitivity range of 0.8–1.2 examined in Section 5.4. The local scan assessed behaviour around the baseline parameter, whereas the scenario analysis examined the model response to a larger hypothetical change. The results of the technology scenario should therefore be interpreted as stress-test evidence under the model equations rather than as an empirically calibrated estimate of a feasible investment programme.
The credit-index values and differences from baseline at the selected comparison points are summarised in Table 11.
Scenario | $\boldsymbol{t}$ = 5 | Delta | $\boldsymbol{t}$ = 10 | Delta | $\boldsymbol{t}$ = 20 | Delta | $\boldsymbol{t}$ = 50* | Delta |
|---|---|---|---|---|---|---|---|---|
Baseline | 42.02 | 0.00 | 4.98 | 0.00 | 4.06 | 0.00 | 5.64 | 0.00 |
Institutional | 37.69 | $-$4.32 | 7.73 | $+$2.75 | 5.43 | $+$1.37 | 5.52 | $-$0.13 |
Technology | 67.61 | $+$25.59 | 25.85 | $+$20.86 | 39.09 | $+$35.03 | 33.42 | $+$27.77 |
Combined | 56.54 | $+$14.53 | 28.78 | $+$23.79 | 34.22 | $+$30.16 | 33.73 | $+$28.09 |
The technology intervention produced the highest $CI$ values at times 5 and 20. Under the model equations, doubling the investment coefficient directly increased the technology stock, information symmetry, and the technology-related contribution to credit improvement. The combined intervention produced the highest value at time 10 but remained below the technology-only scenario at time 20.
At time 50, the combined scenario produced a value of 33.73, compared with 33.42 under the technology-only scenario. The difference of 0.31 was small and should not be interpreted as evidence that the combined intervention was materially more stable or effective. It only shows that the combined scenario produced a slightly higher value at this long-horizon diagnostic point.
The institutional-response intervention produced a value 4.32 below the baseline at time 5, 2.75 above the baseline at time 10, and 1.37 above the baseline at time 20. At time 50, it was 0.13 below the baseline. Reducing the response delay changed the timing of regulatory pressure and guarantee adjustment but did not introduce a new long-run target or a separate enforcement-capacity stock. Under the stated parameterisation, the institutional-response intervention moderated the downturn and improved selected intermediate outcomes but did not produce a clear long-run advantage.
These comparisons indicate that the relative performance of the modelled interventions depended on both the intervention mechanism and the selected time point. They do not establish that technology investment is universally preferable to institutional reform, particularly because the institutional-response scenario represented only a change in response delay and the model did not explicitly include S8, S9, S10, S12, or S13. The scenario results therefore support conditional intervention comparison within the stated model rather than a general ranking of all available engineering project credit-risk policies. The trajectories of the four scenarios are shown in Figure 12.

6. Discussion
The structural and dynamic analyses addressed related aspects of engineering project credit-risk governance, but they did not operate with identical variables or provide interchangeable evidence. In the structural analysis, S6 was the most prominent factor and had the largest positive net relation. This result indicates that supervision organisations occupy a central position in the assessed network because project information becomes usable for governance only after relevant events have been verified, recorded, and escalated. S9 occupied the root level of the ISM hierarchy and had the highest MICMAC driving power. Within the questionnaire-defined network, environmental and resource compliance risk therefore functioned as an upstream condition connecting project data, supervisory capacity, and organisational obligations. S10 was prominent but dependent, suggesting that improvements in credit-management procedures may rely on prior changes in verification, technological capability, and the project-execution controls that supply information to those procedures.
The dynamic analysis examined a narrower set of mechanisms. Under the selected scenario settings, the technology intervention produced higher credit-index values than the institutional-response intervention at times 5 and 20. This pattern arose because technology investment entered two positive pathways in the model: it increased information symmetry, which reduced quality risk, and it increased the technology-related contribution to credit improvement. The institutional-response intervention changed the regulatory response delay but left the long-run target structure, regulation-completeness reference, and enforcement-capacity structure unchanged. A shorter delay could therefore alter the timing of regulatory pressure and guarantee adjustment without necessarily producing a higher long-horizon value.
These scenario differences should not be interpreted as evidence that technology investment is generally more effective than institutional reform. The two interventions changed different parameters by different magnitudes, and the technology scenario used a coefficient of 2.0, which lay outside the local sensitivity range of 0.8–1.2. In addition, S9 and S10 were not explicitly included in the system dynamics model, while S6 was represented only through a composite proxy. The simulation therefore compared the available levers under the stated equations rather than validating the complete structural ranking.
At time 50, the combined scenario produced a credit-index value of 33.73, compared with 33.42 under the technology-only scenario. This difference of 0.31 was too small to establish that the combined intervention was materially more stable or more effective. It only indicates that the combined setting produced a slightly higher value at the selected long-horizon diagnostic point. Demonstrating greater stability would require separate evidence on convergence, oscillation amplitude, variance, or another defined stability measure.
At project entry, regulators and project owners should establish qualification requirements, performance guarantees, shared identifiers, and data-access permissions before contract award. These actions relate primarily to S1, S3, S11, and S13. The required management output is not simply a connected database. It is an auditable entry decision that identifies the responsible organisation, the applicable guarantee instrument, the authority to access and verify data, the conditions that trigger further action, and the available appeal route.
During contract performance, project owners and contractors should translate the risks represented by S4, S5, and S8 into measurable obligations concerning payment, schedule, cost, disclosure, and change control. Credit consequences should be linked to verified events and proportionate remedies. This distinction is necessary because delays, cost variations, and contractual changes do not automatically constitute misconduct. A documented assessment process is required to distinguish an ordinary project variation from material non-performance and to connect a verified default with guarantee activation, corrective action, or subsequent credit repair.
During construction monitoring, supervision organisations and regulators should coordinate the factors represented by S2, S6, S7, S9, and S10. A smart supervision platform should retain source evidence and provide traceable records, but the supervision organisation remains responsible for verification, exception classification, and escalation. Building information modelling, connected sensors, distributed ledgers, and analytical systems may support these tasks where their functions match the required evidence chain. However, the present model did not estimate the separate effect of any individual technology and does not support a ranking of specific technological solutions.
Environmental and resource compliance requires particular attention because S9 occupied the root of the ISM hierarchy and had the highest MICMAC driving power. In operational terms, environmental and resource obligations should be converted into measurable project requirements, assigned to responsible participants, monitored through verifiable records, and connected to defined responses. This recommendation follows from the structural analysis rather than the system dynamics model because S9 was not explicitly simulated.
During default handling, regulators, project owners, contractors, supervision organisations, and guarantee providers should connect the verified event record to notification, correction, sanction, guarantee claims, appeals, and credit repair. S1 supplies the project- or market-level assurance mechanism, S11 represents the external regulatory setting, and S3 and S13 affect access to credit information. The autonomous MICMAC positions of S11 and S13 in this dataset call for careful interpretation. Their inclusion in governance arrangements is supported by their legal and operational functions rather than by strong connectivity in the assessed network.
The empirical hierarchy did not reproduce the TOE categories, nor was it expected to do so. TOE describes the nature of each condition, whereas DEMATEL, ISM, and MICMAC describe how the factor was positioned within the assessed relationship network. S7–S9 remained organisational factors because they arose from project execution and participant performance. S11 remained an environmental factor because it concerned external laws, regulations, and enforcement authority, even though it had low prominence, was isolated at the selected MMDE threshold, and was classified as autonomous by MICMAC.
These results demonstrate why conceptual classification and empirical network position should be interpreted separately. A factor may be institutionally important while having limited connectivity in a particular questionnaire-defined network. Conversely, an organisational risk may occupy a root position when the assessed relationships give it substantial driving power. Neither result requires the factor to be reassigned to another TOE category.
The absolute thresholds of 0.38–0.42 retained no off-diagonal edges in the sparse matrix. This was a consequence of scale and sparsity rather than evidence that all 13 factors were substantively independent. The local MMDE sensitivity analysis was more informative for evaluating the selected hierarchy. It showed that the central pathway connecting S9 with S2 and S6 persisted around the selected threshold, while the placement of peripheral factors was more sensitive to small threshold changes.
S11, S12, and S13 remained isolated across the tested thresholds. Recommendations concerning these factors should therefore be framed as legal, cultural, and informational conditions that permit the governance system to operate, rather than as empirically dominant causal levers. Their isolated positions also reflect the limited coverage of the original 33-link instrument and should not be treated as proof that they have no meaningful relationships outside the assessed network.
Priority 1—Verification responsibility and lifecycle coordination. The first priority is to assign clear responsibility for entry checks, contract-event recording, construction verification, exception classification, escalation, sanction, and credit repair. Performance dashboards should distinguish the participant that creates a risk from the organisation responsible for verifying and responding to it. This distinction is particularly important for S6, which had the highest prominence and the largest positive net relation. Verification records should identify the underlying evidence, the responsible organisation, the date of verification, the classification decision, and the action taken. These requirements turn project information into an auditable basis for operational decisions.
Priority 2—Data interoperability and smart supervision. Project owners, contractors, supervision organisations, and regulators should establish shared identifiers, minimum data fields, verification timestamps, and accountable interfaces across their information systems. Technology selection should follow the requirements of the evidence and decision chain rather than determine them. The dynamic model supported investment in a composite smart supervision capability under its stated equations, but it did not separately quantify the effects of blockchain, building information modelling, the Internet of Things, distributed sensing, or big-data analytics. S9 should also be incorporated into this information architecture through measurable environmental and resource-compliance records, although its intervention effect was not examined in the dynamic model.
Priority 3—Credit guarantees, default control, and procedural repair. Project-entry arrangements should establish proportionate performance-guarantee requirements and define the circumstances under which a guarantee may be activated. Default decisions should use explicit evidence thresholds and connect verified non-performance with notice, correction, sanction, appeal, and credit-repair procedures. This priority addresses S1 and S10 by linking assurance mechanisms to workable credit-management routines. It also prevents guarantee requirements from becoming detached from actual project events or being applied automatically to ordinary schedule changes and contractual variations.
These priorities combine structural ranking with lifecycle feasibility. They should not be interpreted as a universal sequence for every engineering project. Regulators and project organisations may need to adjust their order according to project type, institutional setting, available data infrastructure, and the severity of the risks being managed. Their purpose is to translate the analytical results into a transparent set of operational and strategic choices.
The evidence is subject to six main limitations. First, the sample consisted of 86 anonymous online submissions. The archived dataset did not contain information on occupation, professional experience, organisational affiliation, project type, engineering credit-supervision experience, or geographic location. Respondent-group and regional comparisons could therefore not be conducted. The factor-development process drew partly on China-specific engineering and regulatory policy documents. This institutional context may limit the direct transfer of the resulting factor structure to regulatory environments that were not examined. The absence of professional-profile information also prevents retrospective verification of the respondents’ domain expertise.
Second, the questionnaire assessed only 33 directional relationships and used a forced-choice scale from 1 to 4 without a zero or no-influence option. Treating the unassessed relationships as structural non-links preserved the traceability of the matrix, but it made the results conditional on the coverage of the original instrument. The zero entries cannot be interpreted as respondent judgements that no substantive relationship existed.
Third, the bootstrap analysis measured score-sampling variation within the collected responses. It did not capture differences among occupational or regional groups, and it could not assess uncertainty concerning relationships that were not included in the questionnaire. The narrow bootstrap rank intervals therefore support stability within the observed 33-link topology rather than the completeness of the network.
Fourth, the ISM levels and MICMAC classifications depended on threshold selection. The local sensitivity tests retained the main deep pathway near the selected MMDE threshold but showed increasing network fragmentation at MMDE + 10%. The absolute thresholds of 0.38–0.42 were too high relative to the sparse total-relation matrix to provide locally informative cutoffs.
Fifth, the system dynamics model represented S1–S7 and S11 either directly, through composite variables, or as an exogenous model input, whereas S8, S9, S10, S12, and S13 were omitted because the source model contained no corresponding equations. S6 was represented only by a composite proxy. The dynamic analysis therefore complemented the 13-factor structural model rather than validating it. In particular, the simulation could not compare direct interventions targeting S9 or S10, even though these factors occupied prominent positions in the structural results.
Sixth, several parameter values were inherited assumptions from the archived source model, and the system dynamics model had no historical calibration target. The scenario differences, Monte Carlo intervals, and trajectories over times 30–50 describe behaviour under the specified equations and parameter settings; they are not predictions of future engineering-industry credit outcomes. The use of different parameter changes in the institutional-response and technology scenarios also means that their numerical results should not be interpreted as a controlled comparison of equally intensive interventions.
Future research should collect non-identifying information on respondent roles, professional experience, project types, and regions; assess a more complete directional network using a genuine no-influence option; and examine whether the factor structure varies across institutional and project settings. An expanded system dynamics model should include S8, S9, S10, S12, and S13 explicitly and should be calibrated against longitudinal indicators where suitable data are available. Further work should also define measurable intervention costs and constraints so that governance alternatives can be compared in terms of effectiveness, timing, feasibility, and resource requirements.
7. Conclusions
The full-precision analysis identified insufficient credit verification by supervision units, environmental and resource compliance credit risk, and lagging credit-management methods as the three most prominent factors in the questionnaire-defined network. The MMDE-based ISM structure placed S9 at the root level and S2 and S6 at the next level, while MICMAC classified S2, S6, S8, and S9 as independent factors with relatively high driving power and low dependence.
The exploratory system dynamics model examined only the institutional-response and technology-related \sloppy mechanisms represented in its equations. Under the selected parameter settings, the technology intervention produced the highest credit-index value at time 20. Reducing the regulatory response delay moderated the downturn and improved selected intermediate outcomes but did not produce a clear long-horizon advantage. At time 50, the combined scenario produced a value of 33.73, only 0.31 above the technology-only value of 33.42. This small difference is a long-horizon model diagnostic and does not establish greater stability or predict performance over a real-world 50-year period.
The findings indicate that engineering project credit-risk governance should be organised as an auditable lifecycle rather than as a stand-alone technology programme. Entry arrangements should establish guarantees, participant identities, and data permissions; contractual obligations should be translated into verifiable project events; supervision organisations should be assigned explicit responsibility for verification and escalation; and confirmed defaults should be connected to proportionate correction, sanction, appeal, and credit-repair procedures. Smart supervision technology supports this process by improving information connectivity and monitoring capacity, but its value depends on the institutional responsibilities and decision rules within which it operates.
The study provides a traceable basis for identifying factor priorities and comparing selected intervention pathways, while its conclusions remain bounded by the 86-response sample, the questionnaire-defined 33-link network, and the exploratory, partially mapped system dynamics model. The results support decision-making within these stated boundaries and should not be interpreted as a universal causal structure or as a forecast of engineering-industry credit performance.
Conceptualization, Z.H.Z.; methodology, S.B.Z. and Z.H.Z.; software, S.B.Z.; validation, Y.Q.H.; formal analysis, S.B.Z.; investigation, Y.Q.H. and W.S.Z.; data curation, Y.Q.H. and W.S.Z.; writing—original draft preparation, S.B.Z.; writing—review and editing, Z.H.Z., Y.Q.H., and W.S.Z.; visualization, S.B.Z.; supervision, Z.H.Z.; funding acquisition, Z.H.Z. All authors have read and agreed to the published version of the manuscript.
The data used to support the research findings are available from the corresponding author upon request. Individual questionnaire responses are not publicly shared because they contain respondent-level information and are subject to privacy considerations.
The authors declare no conflicts of interest.
