#### Abstract

Organizational routines are means through which organizations can reutilize best practices and so their replication, i.e., duplicating beneficial routines across context, is a key value-creating strategy. However, it is difficult to map network effects on routine replication. Here, we investigated routine replicating dynamics considering two types of network contexts, namely, (1) connections between different (geographically distributed) units in a decentralized organization and (2) the coupling relation between routines, i.e., a bundle of different routines involved in each unit. By considering routine replication as one kind of template-based activities between different units, we examined interrelations between routines with a NK-based fitness landscape model. Our results show that when there is an appropriate level of absorptive capacities (i.e., when organizations are capable of identifying and acquiring externally generated knowledge), there is an optimal combination of these two types of networks, which is beneficial to routine replicating practices and organization adaptation. Furthermore, we also found that intraorganizational variations, including template-duplicating errors and innovative activities, are instrumental to enhance adaptive changes. Our findings suggest measures to control and manage best practice diffusion across organizations.

#### 1. Introduction

Recent research has shown that organizational routines, i.e., “repetitive, recognizable pattern of interdependent actions” carried out by multiple individual actors [1], are “ubiquitous” in almost every organization [2]. While certain scholars have focused on the so-called “defensive routines”, i.e., actions or policies that “prevent either individuals or segments of the organization from experiencing embarrassment or threat” [3, 4], and explored their negative effect on organizational learning and change [5–7], there is evidence that organizational routines are a key repository of organizational memory, skills and tacit knowledge which enables both organizational stability and change [1, 8–11] and constitutes organizations’ sustained competitive advantages [12–14]. In this sense, by applying routines representing “best practices” to new environments, organizational core capabilities can be enhanced that in turn further sustain competitive advantages and profitability [15]. This is why the replication of organizational routines is considered as an important value-creating strategy and attracted more increasing attention from both researchers and practitioners [16–18].

Firesl and Larty [18] summarized three value-creating mechanisms of replicating routines for which replication is viewed (1) as a strategy for growth; (2) a means to create value by improving organizational efficiency; and (3) an organizational scaffold that may enhance flexibility. These value-creating mechanisms, according to Jenson [19] and D’Adderio [20], are particularly applicable to organizations, such as large or multinational firms, within which multiple similar units perform the same function, or franchise and chain stores, which imply the geographic replication of similar units.

According to Szulanski and Jensen [21], replicating a routine involves the creation of another routine that is ‘similar to the original routine in significant respects.’ However, Howard-Grenville [22] pointed out that organizational routines might perform differently depending on specific organizational environments. Given certain context-dependent characteristics, the replication of organizational routines can be understood as a simple knowledge transfer process. Indeed, any match between these newly created routines and the context they were originally coping with imposes what Winter and Szulanski [15] called a “replication dilemma” problem. This means that besides the exact “cloning” of the original routines, these new routines must be to some extent flexible and adaptive to be suitable to the new, in part similar but also different, organizational context [18, 20].

Therefore, replication of routines is a means through which organizations can reutilize knowledge that is already in use [21, 23]. Some scholars focused on the importance of replication as a “knowledge transfer” across intraorganizational boundaries [23–26] and considered the intraorganizational transferring process between different units, often across different geographic locations, from a network-based perspective [27–29]. Here, there is the important role of templates during the replicating process of routines [21, 30] as a means to stimulate forward and reverse knowledge flows [18].

As suggested by D’Adderio [20] and Friesl and Larty [31], replicating routines is a complex social practice, which involves multiple actors and artifacts. For instance, the capacity of individual actors to understand specific routines is always intertwined with their own behaviors, experiences, and the particular environment within which they live [12]. Therefore, it is likely that individuals have to make context-specific adjustments so as to reproduce “the similar routine” in prospection [11, 32–34]. Unfortunately, as suggested by Friesl and Larty [18, 31], we still know little about the micro processes through which individual actors and actions shape and enact new routines in their new contexts.

Furthermore, research on organizational routines has emphasized the connective nature of routines in that organizations and capabilities are better understood as dependent on some “specific networks of routines,” rather than just on bundles of them [11, 35–37]. This is because routines often do not work alone but are interrelated with each other [38]. Secondly, recent findings indicated that the dynamics of several interrelated routines together, which are often what occurs in practice and to some extent determines organizational outcomes, are definitely different from that of individual routines [39]. If so, routine interdependency can be a fundamental source of both organizational stability and flexibility [40]. However, as recently suggested by Feldman et al. [11] and Kremser and Schreyögg [39], more research is needed to understand ‘implications and consequences of routine interrelatedness.’

Our paper aims to investigate the replicating dynamics of organizational routines within two types of network contexts. First, we hypothesized connections between (geographically) distributed units of an organization [18, 41, 42]. In a second scenario, we hypothesized a routine-related network – i.e., a network of different routines involved in the whole organization [11]. To do so, we first followed the NK and fitness landscape modeling [43] to depict routine networks and built up the fitness landscape of an organization. Then, we translated this into an agent-based model (from now on, ABM) [44–46] to explore complex behaviors of routine replications.

In our model, we focused neither on individual actors nor on artifacts in the enactment of routines in new contexts. Rather, we treated an aggregation of these individuals and artifacts, namely, the unit distributed within an organization, which we called a “unit”, as the basic heterogeneous but interdependent agent. All the units (i.e., agents) randomly interacted with each other, while by using templates specific routines were replicated from unit to unit. Furthermore, we also considered the impact of both internal variability and external dynamics of the organizational environment on routine replications. This included, for example, template-duplicating errors and innovative activities among units and the pace of environmental changes.

The rest of the paper is structured as follows. In the second section, we provide a brief literature overview on replication of organizational routines. Further, we also summarized the feasibility and effectiveness of adopting an NK-inspired ABM to examine organizational routines. Then, we present a series of concrete scenarios with combinations of input settings and show our simulation results. Finally, the last section summarizes our main findings, discusses limitations, and elaborates on implications.

#### 2. Theoretical Foundations

##### 2.1. Replication of Routines within Organizations

Nelson and Winter [8] and Zollo and Winter [47] considered an organization as “a set of interdependent operational and administrative routines.” These different kinds of routines embody an important part of the organization’s productive knowledge [8] and are considered as one of the key sources for sustained competitive advantages [12, 48, 49]. This suggests that sustaining competitive advantages can also depend on the capability of replicating the most efficacious routines from a unit to another, or even from an organization to another, just as efficacious habits are replicated from an individual to another one [17, 50]. However, in some respects, routine replication can be more difficult than the replication of individual habits. First, after discovering an efficacious practice, managers tend to build organizational barriers against replication, so as to exploit the system’s advantage and avoid imitation by others [51]. Second, knowledge embedded in routines is always to some extent not explicit but tacit and enacted in practice [32, 47]. At the same time, either units or organizations have only limited capacities to discern superior practices from their neighbors [50].

Furthermore, given the “context-dependent” characteristics of routines [22] and the need for adaptive changes [15, 52], for any organization, routine replication has value when it can support appropriate variations that are instrumental to cope with environmental challenges [15, 18, 20].

In this respect, by discussing the connection between intraorganizational replication and interorganizational imitation, Rivkin [51] further pinpointed that replication can be more efficient than imitation in promoting sustained competitive advantage only when the complexity of the task and the environment is moderate. Winter and Szulanski [23] distinguished the term “replication” from similar concepts such as “transfer of knowledge” and “diffusion of innovations” within and/or even between organizations. They suggested that replication involves the creation of “a series of local routines that are quite similar to the original ones in significant respects” [23]. This implies that replication often relates to an organization’s attempts to “reproduce at multiple internal sites the outcome of an existing activity” and so “applies to a broad class of exploitation situations” [8, 23].

Since Nelson and Winter [8], many scholars explored the dynamics of replication, such as routine replication barriers [24, 53], the impact of replicator/replicatee relationship, and the role of institutional contexts [54]. Others looked at the movement of employees or independent experts as a carrier of knowledge and experience transfer from one context to another (e.g., Hodgson [50]). If replication is a “knowledge transferring” process across intraorganizational boundaries [24, 25], this process depends on the fact that the ‘core aspects of a routine’, or even “a whole business format”, can be truly communicated from one unit to another, which then transforms related knowledge into practice [18].

Research suggested that “templates” are key for replicating organizational routines [21, 52]. Here, the term “template,” according to Nelson and Winter [8], referred to “working examples” that contain some critical aspects of the routine. Given that knowledge embedded in routines is always to some extent not explicit, but tacit and enacted in practice [32, 47], templates define the “details of how the work gets done, in what sequence, and how its various components and subroutines are interconnected” [19]. They can help us codifying and explicitly representing the tacit part of knowledge, thereby providing an effective way to overcome stickiness to the original context of knowledge [21]. Meanwhile, by codifying and explicitly representing knowledge inherent in routines, templates can also be helpful to the replicatee units in discerning and adopting the so-called “best practices” [55]. This requires to look at the relationship between the sending and receiving units [24, 56], as the ability and willingness of replicatees to accept knowledge transferred from the sending unit are also prerequisites for the replicating process of routines [57].

##### 2.2. Routine Networks and the NK Fitness Landscape

Organizations can be considered as certain specific aggregations of “bundle of routines” [8, 48, 58, 59], which can form and change through a process of “combinatorics” [12]. However, it is only until recently that scholars began to focus on how routines affect one another in a more complexity-friendly perspective [11, 36, 60]. For instance, Peng et al. [36] emphasized that looking at routine interrelations is key to understand how routines bundled together to generate organizational capabilities – including both incremental and radical innovation. According to Kremser and Schreyögg [39], this interrelatedness can be attributed to the division of labor in that, due to economics of specialization, organizations while dividing work into “multiple, separate routines” are also required to “securing a sufficient integration of efforts” [39]. This suggests the need for looking not only at the division of multiple specialized routines from the complex organizational task, but also at the interrelationship between these routines. This implies that only by specific routine integration organizations can fully accomplish their organizational task.

Here, research suggests a holistic perspective as organizational routines are embedded in systems of routines [14, 36, 61]. This means that they are systemically integrated and reinforce each other so giving rise to a collective effect, which is greater than only the single contribution of a simple set of individual components [62]. In this line, Argote and Ingram [25] considered the “network of routines” as an important form of organizational memories and the repository of organizational knowledge. Salvato [63] argued that ‘the interplay of existing routines’, together with individuals’ everyday actions, should be treated as real “unit of analysis” rather than single factors. More recently, Sele and Grand [60] investigated the role of both human and nonhuman actants in the connection between routines and considered the interaction of routines as a source of generativity for the whole organizational system.

In this paper, we tried to incorporate some aspects of these networks of routines by using Kauffman’s [43] NK model to formalize interactions among these multiple organizational routines. Being originally developed to analyze biological organizations’ evolutionary problems [64], Kauffman’s NK modeling and fitness landscape theory was then introduced into organizational studies by Levinthal [65]. It eventually became popular in studies on technological innovation [66], strategic planning [67, 68], and organizational learning [69]. NK modeling now can be considered as one of the crucial toolkits for coping with a series of classic organizational problems of search in multidimensional spaces. According to Ganco and Hoetker [68], NK modeling provides a “relatively simple and replicable methodology” addressing organizational problems that cannot be comprehensively understood empirically.

#### 3. The Model

##### 3.1. A fitness Landscape Approach to Routine Networks

Based on Nelson and Winter’s [8] notion of routines as the “genes” of organizations, Ventresca and Kaghan [70] and Yi et al. [71] suggested to apply the NK model to the study of organizational routines. Through the NK-based fitness landscape, we can build up a description of the linkage between the network of routines and the organizational fit with its environment. For example, let us assume that there are* N* routines involved in a unit, , and for each routine (*i* = 1, 2, …,* N*), each owning a number of different states – or in other words, different kinds of approaches for unit members to accomplish the corresponding organizational subtasks . Thus, according to Kauffman’s [43] theory, this routine’s contribution to the overall unit fitness would depend not only on its own state, but also on the states of some other* K* routines (where and ) with which it is interconnected.

Kauffman [43] described two main methods for assigning relations between dependent parts – namely nearest neighbor influences and randomly assigned influences. He then proved that the distribution of* K* did not dramatically affect the NK model outcomes. In this paper, for simplicity, we adopted the latter strategy and randomly assigned relations between routines [72], as shown in Figure 1.

Thus, we had the following fitness contribution of each routine ,Here, the value of the key parameter* K* can be from 0 to* N* – 1. On the one hand,* K* = 0 represented a situation in which each of the* N* routines involved in the unit was independent. On the other hand, with the increase of the value of* K* (*K* > 0), the change of one routine would influence the performance of other interrelated routines in that any search for improving the configuration of all of them would be more and more difficult [71]. Finally,* K* =* N* – 1 represented another very specific situation in which all routines were interconnected with each other.

The fitness value of the unit with its task environment can be calculated asaccording to which the higher the fitness value* Fitness*, the more adaptive the unit to its living environment.

##### 3.2. Agent Design

In our model, we focused on network characteristics of distributed branches/units within an organization. For the sake of simplicity, we assumed that there was a decentralized organization [73] with* unit*_*num* = 40 units. All these units were interdependent entities with similar characteristics, facing a recurring but common problem. They were randomly located in a 100 × 100 lattice and interconnecting with each other with a given probability* random_p* (0 ≤* random_p* ≤ 1). The units constituting the whole organization faced the same environment defined by a landscape described in formulas (1) and (2). We assumed that every unit (*j* = 1, 2, …,* unit_num*) was considered a heterogeneous agent dealing with a series of organizational subtasks and having a correspondingly specific routine for each task (*i* = 1, 2, …,* N*) [71]. All routines were interrelated with each other, resulting in interrelationships both within and between bundles of routines, as shown in Figure 2.

Furthermore, we assumed that each routine could generate only two different patterns of action, “0” and “1”, respectively. Although a specific routine may have more than two states, a set of routines could also be interpreted “as subroutines that constitute a routine of more than two states” [71].

We assumed that each unit agent had attributes as shown in Table 1.

However, as we mentioned in Section 2.1, unit agents often used templates when replicating routines from neighbors. We assumed that a template could not be the exact copy of the network of routines of the replicator unit agent. Indeed, given that at least a part of knowledge embedded in routines is to some extent tacit and enacted in practice [32, 47], these differences mimicked in difficulties for replicatee units to explicitly encode routines into templates.

More specifically, we assumed that when replicating activities occurred, the replicatee unit agents obtained a template which was duplicated from the network of routines of the replicator unit agents but with some given probability of errors* error_p* (0 ≤* error_p* ≤ 1). This allowed us to mimic errors that may occur during codifying and explicitly representing process of the original set of organizational routines [18, 25]. However, note that in the model we considered the value of* error_p* as constant throughout all simulations. Considering reasons and/or even mechanisms for this kind of errors was out of the scope of our work.

Finally, we considered the* environment* as a regime of possibilities and constraints for unit agents while accomplishing their organizational task [22]. This environment determined both the organizational subtask series and the fitness landscape as shown in formulas (1) and (2). In order to simulate environmental dynamics, in our model, we adopted a Boolean variable* is_env_change?* to represent whether the fitness landscape changed or not. On the one hand, when* is_env_change?* =* TRUE*, the fitness landscape would randomly vary once after every certain number of simulation rounds* env_change_P* so as to mimic environmental changes occurred in practice. Further, for any* env_change_P *≥ 1, we have that the lower the value of this variable was, the more rapid environmental changes happened. On the other hand, when* is_env_change?* =* FALSE*, it meant that there was no change of the organizational environment. Thus, any values of the variable* env_change_P *are meaningless in this situation.

##### 3.3. Rules of the Replicating Activities

According to Hodgson [50], replicating routines from peers is not easy often also due to organizational inertia (e.g., [71, 76]) and bounded rationality [75]. To mimic this, we assumed that unit agents were prone to replicate “superior practices” from neighbors only when they felt competitive pressures, i.e., they had a significantly lower fitness value than the average of neighbors.

Next, a* template* was generated through explicitly codifying the knowledge embedded in the network of routines of the replicator unit agent [19, 21]. In this model, we used a vector of “0” or “1” values as follows:

being templates of a copy of the original series of routines, but allowing each of the element values of the vector* tpl*_{k} (0 ≤* k* ≤* N*) to vary with a certain probability of errors* error_p* (0 ≤* error_p* ≤ 1). Namely,Here, the value of the coefficient was either “0” or “1”, which represented the two different states (i.e., patterns of actions) of the* k*th original routine of the replicator unit agent, and .

Furthermore, given that replicatee unit agents were always constrained by a specific level of absorptive capacity [77, 78], we assumed that they could not adopt the whole template but only a certain part. Here, we used a variable* unit_AC* to represent the absorptive capacity of unit agents. Its value numerically equaled the number of routines that replicatee unit agents could change at a time according to the template,* unit_AC* = . Obviously, the higher the value of* unit_AC*, the higher level of the unit agent’s absorptive capacity. In this respect, note that exploring details of agents’ absorptive capacities was beyond the scope of our study. For the sake of simplicity, we assumed that the value of* unit_AC* was assigned at the beginning of the simulation and remained fixed throughout our simulations.

The replicatee unit agents could randomly select a number of routines (which numerically equaled the* unit_AC* value), which owned different states with those from the template – if they existed, of course. Then, the replicatee unit agents changed the relating routines’ states (values) according to the template and generated a new list of routines* R*’.

Finally, when the new collection of routines* R*’ was formed, replicatee unit agents checked if they could perform any innovation to match their work of routines (obviously, here, referring to the list of routines* R*’) and the specific context in which they were living. This was to mimic the “context-dependent” characteristic of organizational routines [22] as well as the interconnection between routines [15, 18, 20, 39]. For example, replicatee units could revise parts of routines while considering both the embeddedness of the system of routines [11] and certain fitness landscapes – namely, the units’ living environment [22, 42]. By doing these innovative activities, replicatee unit agents transformed the temporary list of routines* R*’ into their final network of routines* R* with newly updated status.

In our model, we assumed that replicatee unit agents could innovate and randomly vary the state of one of their routines with a given probability* innovate_p* (0 ≤* innovate_p* ≤ 1) – when this was helpful to improve their fitness values.

Figure 3 provides a description of the routine replicating mechanisms between different unit agents. Figure 4 shows a stylized description of the procedure in which unit agents involved in the model should execute during the whole simulation process.

#### 4. Simulation Experiments, Initialization, and Verification of the Model

##### 4.1. Simulation Experiments

We constructed our model using NetLogo (https://ccl.northwestern.edu/netlogo/). All the main input settings were shown in Table 2. By applying the unit testing method [79] to each program module and to the whole simulation model, we ensured that our model was comprehensively verified.

The aim of this paper was to investigate the replicating dynamics of routines under the constraints of both the intraorganizational network of unit agents and the routine networks. Furthermore, we considered four other factors, namely, the absorptive capacity of unit agents, the probability of errors occurring during the template-duplicating processes and that of unit agents’ innovative activities, and the dynamic characteristic of the organizational environment.

To do so, we first considered one of the simplest cases (“Scenario I”). By assuming that there was no environmental change (in this model, we let the Boolean variable* is_env_change?* =* FALSE*, referring to no changes occurred of the fitness landscape) and that there are neither template copying errors nor innovative activities (i.e.,* error_p* =* innovate_p* = 0), we purely investigated the influence of both of the two different networks – namely, one is the intraorganizational network of unit agents (determined by the model parameter* random_p*) and the other is the network of routines (represented by the model parameter* K_routine*) – on routines’ replication dynamics.

Then, in “Scenario II”, we added variability. On the one hand, we introduced template-duplicating errors (represented by a series of values of the parameter* error_p*). On the other hand, we assumed that success probability of innovative activities by different replicatee units (represented by variations of the parameter values* innovate_p*) could result in variations of the routine system.

Next, in “Scenario III”, we introduced uncertainty of the organizational environment, by varying model parameter values* is_env_change?* =* TRUE*, and* env_change_P*, so varying the changing rate of the fitness landscape function. This allowed us to examine routine replication within the dynamic organizational environment.

Furthermore, in “Scenario IV”, we combined Scenarios II and III to analyze routine replication in a more complex environment, in which internal sources of variability, i.e., template-duplicating errors and innovative activities, were coupled with external sources, i.e., uncertainty of the organizational environment. It is worth noting that these two sources of variations are typical of complex, dynamic environments in the current business landscape, as suggested by Goldstein, et al. [80].

Finally, we developed four aggregate measures to estimate the impact of scenario manipulations. The first one, called* steadyTime*, measured the number of simulation rounds required by the organization system to arrive at a steady state (e.g., every unit agent found the best bundle of routines that matched its own landscape, and thus there was no replication activity left), if such a specific steady state eventually existed, of course. This is crucial to understand the performance of routine replications as it indicates the replicating speed of routines between different units within the organization. The second one is the average fitness values of all the unit agents. This represented the organizational efficiency that every unit could gain by routine replication, so in principle benefiting the whole organization.

We then also measured the number of replicating and innovative activities occurred among all the unit agents (here, denoted as* replicate_num* and* innovate_num*, respectively). When copying routines from their neighbors, unit agents bore some operating cost, i.e., the expenses or resources required to replicating or innovative activities. In general, we hypothesized that the more these activities were performed by unit agents, the higher the operating cost was. In this sense, the number of replicating and innovative activities occurred together could reflect the level of operating cost required by the organization to perform its internal replicating activities.

##### 4.2. Initialization and Verification of the Model

For verifying the correctness of the model, we first considered a very basic scenario in which there were neither variations from template-duplicating errors and innovative activities nor changes in the organizational environment. We let the parameter value* unit*_*num* = 20, 40, 60, and 80 and keep other variables in the default settings as shown in Table 2. This led to four corresponding but definitely different parameter configurations. In order to determine the ideal number of simulation run for each parameter configuration, we applied the statistical power analysis suggested by Secchi and Seri [81, 82]. Given the significance level *α* = 0.01, the power 1 - *β* = 0.95, and the effect size with a small-level* ES* = 0.1 according to discussions in [81], we found that 569 runs were sufficient. Then, using the ‘smallest effect size of interest’ (SESOI) approach, we set simulation runs to 600 times per parameter configuration, which ensured an estimated expected power of 0.96.

On the other hand, after some pilot runs to study convergence and sensitivity, we decided to let each simulation running for 3000 rounds, each of which representing one routine performance. This allowed us to find a convergence of the routine system to a given steady state, at the same time enabling us to measure one of the most important outputs, i.e., the variable* steadyTime*.

We then collected output data and calculated the mean values of each variable in each simulation run. Figure 5 shows a summary of simulation results. Figure 5(a) indicates that with the increasing of the number of units involved in the organization, there was a continuous variation of the organization’s states, i.e., the values of the three main outputs,* steadyTime*,* replicate_num,* and* avg_Fitness*. It is worth noting here that we did not consider the variable* innovate_num* as that* innovate_p* = 0 and no innovative activities occurred. First, when connections between units were kept at a very lower level (*random_P* = 0.1) and there were only a few units involved (i.e.,* unit_num* = 20), replicating activities occurred randomly but with a very small mean value (see Figure 5(c)). The organization took a short time to reach a steady state (see Figure 5(b)), and the* avg_Fitness* values followed a very standard normal distribution with mean value approximately equaling 0.5 (see Figure 5(d)). This corroborated our idea of initializing routine networks of all the units randomly.

**(a)**

**(b)**

**(c)**

**(d)**

Furthermore, the increasing number of units obviously increased their interactive opportunities and thus led to the emergence of certain specific patterns of replicating activities. More specifically, the distribution of the variable* replicate_num* obtained a higher mean value but with a smaller range of variance (see Figure 5(c)). Then, the* steadyTime* values increased, and the distribution of the variable* avg_Fitness* followed a partial peak distribution with a bigger mean value (see Figures 5(b) and 5(d)). Again, all these results confirmed the validity of our logics in building and parameterizing the model.

#### 5. Simulation Results

Table 2 shows the set of conditions after a series of analyses performed on a testing scenario (as mentioned above) and additional “trial-and-error” experiments. We decided to report our simulation experiments in four subsections corresponding to the four scenarios included in Section 4. To avoid either under- or overpowering our simulation results [81], we let the significance level *α* = 0.01, the power 1 - *β* = 0.95, while assuming the effect size with a micro level* ES* = 0.05 due to a pilot study and previous testing experiments. For each of our four scenarios, we followed a full factorial design and calculated the ideal number of simulation runs (denoted as* i_NofRun*) needed for each scenario – as shown in Table 3. Furthermore, following the SESOI approach [81, 82], we set certain running times, see* r_NofRun*, for each parameter configuration, resulting in an estimated expected power (*EEP*) of 0.96. Moreover, we kept each simulation run at 3000 rounds to ensure that the organization system reached a steady state – if this steady state definitely existed.

*(1) Routine Replications within a Static Environment. *We assumed* unit_num* = 40 and* N_routine* = 8 based on the pilot experimental results and assigned* random_p* = 0.1, 0.3, 0.5, 0.7 and 0.9, and* K_routine* = 1, 3, 5 and 7, to represent different density levels of both the network of units and that of routines, respectively. We let* unit_AC* = 1, 2 and 3 to represent different levels of the units’ absorptive capacity. For each configuration of these three parameters and with default settings of the other parameters as shown in Table 2 (here, we named this as Scenario I), we collected all the 400 simulation running results and calculated mean values of the three outputs, i.e.,* steadyTime*,* replicate_num,* and* avg_Fitness* but excluded the variable* innovate_num* because we removed any possibility of innovative activities due to the parameter* innovate_p* = 0 in this scenario.

Figure 6 shows simulation results of the other three outputs – namely,* steadyTime*,* replicate_num,* and* avg_Fitenss*. We found that the network density of both units and routines had an impact on routine replication processes, so leading to different performance of the organization. However, their effects were also determined by another important variable, namely, the units’ absorptive capacity* unit_AC*. First, when the parameter* unit_AC* = 1, units were limited within a very narrow space when searching target neighbors for replication. With this constraint, even increasing the unit network density (*random_p* from 0.1 to 0.9) could increase the probability for units to link with each other, and all the three outputs* steadyTime*,* replicate_num,* and* avg_Fitness* could increment only minimally. On the other hand, increasing the routine network density (*K_routine* from 1 to 7) had no significant influence on the outputs* steadyTime*. However, it improved both* replicate_num* and* avg_Fitness* with a limited margin – except at a very medium value* K_routine = 3* when the unit network density* random_p* = 0.7, which resulted in a decrease of the* avg_Fitness* value (see the left part in Figures 6(a), 6(b), and 6(c)).

**(a)**

**(b)**

**(c)**

Next, when the parameter* unit_AC* = 2 and 3, we found that the organization system spent a longer time to arrive at a steady state (i.e., a higher* steadyTime* value, see the medium and right parts in Figure 6(a)). More replicating activities occurred (i.e., a high* replicate_num* value, see the medium and right parts in Figure 6(b)), and so a higher level of fitness was reached (i.e., a high* avg_Fitness* value, see the medium and right parts in Figure 6(c)). In this situation, the lower the value of the variable* random_p *was, the greater the effects of its increments to the organizational system were.

On the other hand, although variation of* K_routine* values had a weak effect when the variable* random_p* took low values (e.g.,* random_p* = 0.1 and 0.3 when* unit_AC* = 2; or* random_p* = 0.1 when* unit_AC* = 3), achieving an optimal fitness when taking a medium value (here, for example,* K_routine* = 3 or 5) and the variable* random_p* with higher values (i.e.,* random_p* = 0.5, 0.7 and 0.9) was extremely beneficial for the organization. Obviously, in this case, reaching a steady state required more time.

However, if we increased* K_routine* to a significantly higher value (e.g., in our model,* K_routine* = 7), the system took less time to reach a steady state, but the increasing values of* replicate_num* and the decreasing values of* avg_Fitness* indicate that the organization undertook more replication activities. This permitted to overcome resistance against change due to the strong ties between different routines, yet with no obvious benefits for organization adaptiveness.

To sum up, if unit agents had only very limited absorptive capacity, variations of both the unit and routine networks’ density do not improve organizational performance. If we consider the cost for unit agents to manage their relationship with partners, the contribution of the replicating process of organizational routines to the organization performance depends on an optimal configuration of the two types of networks – i.e., the network of units in the organization represented by the variable* random_p* and that of routines represented by the variable* K_routine*.

*(2) Routine Replications with Template-Duplicating Errors and Units’ Innovative Activities*. In this section, we considered two important factors of variation. One is the probability of errors that may occur when we duplicate routines from the template. The other one is the probability of innovation for unit agents while they replicate routines from their partners. We assumed the values of the parameter* error_p* = 0.1, 0.3, 0.5, 0.7, and 0.9 and another one* innovate_p* = 0.2, 0.4 and 0.6, representing different levels of the two kinds of events’ occurrence probabilities, respectively. According to simulation results of Scenario I (as shown in Figure 6), we selected one of the most representative but not extreme cases by letting* random_p* = 0.5 and* unit_AC* = 2 and kept the other default settings as Scenario I. Then, we separately investigated the impact of variation from both template-duplicating errors and innovative activities of unit agents.

Figure 7 shows that these induced variations created opportunities for replicatee units with significant benefit to the fitness of the whole organization. However, this came to the cost of more replicating and innovative activities required and so more time is required by the organizational system to reach a steady state. First, increasing network density of routines (i.e.,* K_routine*) meant a tighter coupling relation between these routines. Therefore, the organizations spent more time and were more involved in replicating and innovative activities (see Figures 7(a), 7(b), and 7(c)). However, we found that when the variable* error_p* equaled low values (e.g.,* error_p* = 0.1 and 0.3), a medium value of* K_routine* (i.e.,* K_routine* = 3) still existed which was beneficial to the organization’s adaptiveness (see Figure 7(d)).

**(a)**

**(b)**

**(c)**

**(d)**

Second, in case of variation from template-duplicating errors, on the one hand, the higher was the value of the variable* error_p* (e.g.,* error_p* = 0.9), the more randomness occurred during the replicating activities, which could result in higher heterogeneity of routine networks between unit agents. On the other hand, very low values of* error_p* (e.g.,* error_p* = 0.1) caused less fluctuations to the whole routine system. Therefore, these two cases were both beneficial for the organization to economize on the time required to reach a steady state (as shown in Figure 7(a)). Further, increasing* innovate_p* values mainly led to more innovative activities among unit agents (see Figure 7(c)). This had a significant influence both on* steadyTime* and* replicate_num* only when* error_p* equaled low values (see Figures 7(a) and 7(b)). It also influenced the* avg_Fitness* when* error_p* did not equal medium values – i.e.,* error_p* = 0.1 or 0.9 (see Figure 7(d)).

Furthermore, variations from both template-duplicating errors and unit innovations led to an increase of replicating and innovative activities at the same time reducing average fitness values. However, when* K_routine* values increased from 1 to 7, especially when* K_routine* = 7, the influence of these variations decreased. In these cases, for different* error_p* values (from 0.1 to 0.9), the three outputs* steadyTime*,* innovate_num,* and* avg_Fitness* tended to converge (see Figures 7(a), 7(c), and 7(d)).

*(3) Routine Replications in Dynamic Environments. *In this section, we investigated routine replicating processes in a dynamic environment by letting the variables* is_env_change?* =* TRUE* and* env_change_P* = 200 and 1000, respectively. This was to mimic different changing rates of the organizational environment. For the rest, we kept the other input settings as in Scenario I. In this case, given that environmental changes were pervasive, we expected that the organization system as a whole could not arrive at any steady state. Given that the output value* steadyTime* had no meaning here and that by assuming the variable* innovate_p* = 0 the output values* innovate_num* always equaled 0, we reported here only two simulation outputs including* replicate_num* and* avg_Fitness*.

Figure 8 shows that when the fitness landscape changed, we found that the absorptive capacity level of unit agents has a strong effect on routine replicating activities so determining the average fitness value of all the units involved.

**(a)**

**(b)**

First, when absorptive capacity was minimal (e.g.,* unit_AC* = 1), unit agents were very constrained to search for appropriate routine candidates from their neighbors. Independently from the degree of connections between routines (i.e., the* K_routine* values) or that between unit agents (i.e., the* random_p* values), and the changing rate of the organizational environment (i.e., the different* env_change_P* values), the average number of unit agents involved in replicating activities (i.e.,* replicate_num*) and the average fitness of all the unit agents (i.e.,* avg_Fitness*) did not have a significant effect. In addition, we also found that environmental changes influenced both these two outputs even more with the connection density between units increased (e.g., increasing the* random_p* values from 0.1 to 0.9) (see the left part of Figures 8(a) and 8(b)).

Secondly, higher levels of absorptive capacity (e.g.,* unit_AC* = 3) implied that every unit could easily find its ideal routine targets for replication regardless of the pace of environmental changes. To some extent, this explained why variations of* env_change_P* values had less impact on the* replicate_num* values (see the right part of Figure 8(a)). Further, although both* replicate_num* and* avg_Fitness* were very sensitive to variations of the connection degree between units in case of low values of* random_p* (e.g.,* random_p* = 0.1), the variation of* random_p* had no effect when unit agents were highly connected with each other (i.e.,* random_p* = 0.5, 0.7 and 0.9).

Furthermore, when the variable* unit_AC* equaled a medium value (e.g.,* unit_AC* = 2), we found a correlation between the increasing rapidity of environmental changes (e.g., decreasing* env_change_P* from 1000 to 200) and less replicating activities (see the middle part of Figure 8(a)). This was especially true when the network density of unit agents and routines was high (e.g.,* K_routine* = 5 and 7 and* random_p* = 0.7 and 0.9).

Finally, Figure 8(b) shows that when* unit_AC* = 2 and 3, the organization’s fitness (i.e., the* avg_Fitness* values) decreased when organizational environment varied more rapidly (e.g., decreasing* env_change_P* from 1000 to 200) (see the middle and right parts of Figure 8(b)). In addition, we found a positive relationship between* K_routine* and* avg_Fitness*, except when the unit network density was low (e.g.,* random_p* = 0.1). In this cases, there was an optimal* avg_Fitness* when* K_routine* had a medium value (i.e.,* K_routine* = 3), similarly to what occurred in a static environment.

*(4) Routine Replications with Internally and Externally Coupled Uncertainties. *In this section, we combined Scenarios II and III to examine the influence of coupled internal variations and external uncertainties. To do this, we set up* random_p* = 0.5;* unit_AC* = 2;* error_p* = 0.1, 0.3, 0.5, 0.7, and 0.9;* innovate_p* = 0.2, 0.4 and 0.6;* is_env_change?* =* TRUE*; and* env_change_P* = 200 and 1000, while the other input parameters were the same default settings of scenarios I, II, and III.

Figure 9 shows that within a dynamic environment, more replicating and innovative activities were required to cope with more coupling of routines (i.e., routines’ connection density parameter* K_routine*). Moreover, we found that although more template-duplicating errors (i.e., the variable* error_p*) had less impact on unit agents’ innovative activities, this caused more randomness to the routines system, and thus required more routine replicating activities. On the other hand, while increasing* innovate_p* values could only promote innovative activities of unit agents, this determined less influence of the number of unit agents involving in routine replicating activities (see Figures 9(a) and 9(b)). This was especially true when errors in template-duplication were low (e.g.,* error_p* = 0.1 and 0.3).

**(a)**

**(b)**

**(c)**

Secondly, when the organizational environment changed rapidly, density variations in unit agents and routine networks did not have any significant influence on the level of innovative activities, except when less template-duplicating errors (e.g.,* error_p* = 0.1) were matched with a slowly changing environment (e.g.,* env_change_P* = 1000). In this situation, any increase of routine coupling (e.g.,* K_routine* values varied from 1 to 7) led to more innovative activities. However, differences of both replicating and innovative activities caused by variation of either* K_routine* or* error_p* values decreased with the increasing rapidity of environmental changes.

Figure 9(c) shows that environmental changes had a negative effect on organization adaptiveness. Furthermore, in case internal and external variations were not particularly prominent (i.e.,* error_p* = 0.1 and 0.3, and* env_change_P* = 1000), medium values of* K_routine* (e.g.,* K_routine* = 3) were still beneficial to the optimal fitness of the organization. Otherwise, we found a positive relation between* K_routine* and* avg_Fitness* values when* K_routine* = 1, but this relation was not significant when* K_routine* value was high (e.g.,* K_routine* = 5).

Furthermore, in a slowly changing environment (e.g.,* env_change_P *= 1000), with the increase of template-duplicating errors (i.e., increasing* error_p* from 0.1 to 0.9), the average fitness values decreased. However, this negative relation decreased when the speed of environmental changes was higher (see the bottom part of Figure 9(c)). On the other hand, no matter whether the organizational environment changed slowly or rapidly, higher values of* innovate_p* only resulted in more innovative activities (i.e., the increase of* innovate_num* values) as a means to cope with the randomness caused by both internal variations and external dynamics, and this did not have significant effects on the output* avg_Fitness* (see the medium and right parts of Figure 9(c)).

Finally, by comparing the bottom and upper parts of Figures 9(a) and 9(c), we found when the organizational environment was moderate (e.g.,* env_change_P* = 1000), the routine system could reach a higher level of fitness (i.e., the* avg_Fitness* values) but with less replicating activities in contrast to that within a dramatically changing environment (e.g.,* env_change_P* = 200). These findings are consistent with Rivkin’s arguments in [51].

#### 6. Discussions and Conclusions

In this paper, we considered replication of organizational routines as some specific “template-based” activities between different units within an organization [19, 21, 30] and developed an agent-based model to simulate the routine replicating dynamics from the “bottom-up” approach [83]. We considered two distinguished types of network contexts, the first one referred to connections between (geographically) distributed units [24] and the other one was the coupling relation between different routines [11, 36, 60, 63]. By depicting the former as a random network model [74] and the latter as a NK-based fitness landscape [43], we constructed an agent-based model that allowed us to zoom in the routine replicating process in detail. At the same time, the model permitted to consider the impact of structural factors, such as the two different types of network contexts, the internal variations from both template-duplicating errors and innovative activities of unit agents involved, and the uncertainty of environmental changes (more detailed panel regression results can be found in the Appendix).

Our findings showed that in a relatively static environment, a functional replication of routines within a decentralized organization depends on absorptive capacity of unit agents [18, 24]. Only by an appropriate level of absorptive capacity, each unit can have opportunities to replicate the “best practices” from their neighbors. In this case, an optimal combination of the two different types of networks – namely, both the connections between different units and the coupling relations between routines – can be achieved which is beneficial to routine replicating practices that help the organization to adapt to its living environment.

Then, we also found that certain degree of internal variations either from template-duplicating errors or innovative activities of units involved could create opportunities for unit agents to develop “performative drifts” of organizational routines [32, 84]. While this is beneficial to the adaptiveness of the whole organization, it also comes at operating expenses. However, it is worth noting that a higher level of variations could be harmful to routine replication dynamics [12, 85].

Next, we found that although environmental changes require more routine replicating activities, these processes can be detrimental to the whole organization’s adaptiveness [59, 86]. When the environment is dynamic and uncertain, both routine replicating processes and the performance of the organization are sensitive to smaller connecting density values either of the network of unit agents or that of routines.

This implies that increasing connections between different units are not always beneficial to the organization. Managing unit relations always requires balancing costs and benefits [42, 87]. On the other hand, our findings show that the coupling relation of routines is key to cope with the changing pace of the organizational environment [11].

Finally, we found that a certain level of variation from either template-duplicating errors or innovative activities can help organizations to cope with environmental changes better [18]. Here, developing platforms to exploit a certain degree of redundancy and uncertainty from the bottom-up can increase the capacity of an organization to leverage on emergent processes and so cope with turbulent environments better (e.g., [88]). Obviously, this depends on the extent of internal variations, as these could increase operating costs without benefiting organizational performance.

This said our study has important limitations. First, our agent-based model of routine replication dynamics is based on a series of pioneer work (e.g., [15, 17–21]) that inherited Nelson and Winter’s (1982) seminal book [8]. This implies that we focused on routines that are beneficial for organizations’ survival and development; i.e., by reconfiguring these routines, creating new ones, and using them consistently across organizational sites, organizations can adapt to a changing environment and obtain a high-level of performance [47–49]. However, “defensive routines” as mentioned in [5–7] are beyond the scope of our discussion. Second, we only limited our concentration on the network characteristic of distributed units within a decentralized organization – i.e., the organization within which activities (particularly those regarding decision-making) are distributed or delegated away from a central, authoritative leadership – rather than considering factors such as leadership and culture [73, 80]. These factors are, of course, very crucial ones for organization studies, and we would discuss their effects on routine dynamics in detail in another forthcoming paper. Third, the template-based replication of routines between different units is often linked with some kinds of artifacts – i.e., “physical manifestations” of the performance of organizational routines [89]. These artifacts play an important role in routine dynamics [20, 90, 91]. In this sense, template-duplicating activities can incorporate very complicated interactions between human actors and their material contexts, which are channelled by relating artifacts. Investigating the role of different kinds of artifacts in routine replicating dynamics could complement and enrich our findings. Last but not least, our agent-based model is only a* typification* and so should be verified with context-specific details [92, 93].

#### Appendix

#### Panel Regression Results

In order to evaluate if the fixed was preferred to the random coefficient panel regressions, we performed a series of Hausman tests using the R package “plm” [94], as shown in Table 4. These tests had significant results with a* p-value* < 0.01, which indicated that the fixed effects were better fit to the data.

#### Data Availability

No empirical data were involved to support this study. We only got input settings of the simulation based on a series of simple specific scenarios testing exercises and repeated “trial-and-error” experiments. For each scenario mentioned in the paper, we repeated the simulation running for hundreds of times according to the power analysis results as shown in Table 3 and discussed the simulation results in detail. All the simulation outputs data used in this paper are available from the corresponding author upon request. The NetLogo code is available at openABM Archive https://www.comses.net/codebase-release/459505f3-45be-4316-bb61-91233a04e332/.

#### Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

#### Acknowledgments

This research has been supported by Shandong Provincial Natural Science Foundation (under Grant Nos. ZR2016GB06 and ZR2016GQ07), Shandong Technology and Business University’s (SDTBU) Doctoral Foundation (under Grant No. BS201606), and National Natural Science Foundation of China (under Grant No. 71501113).