Loading [MathJax]/jax/output/SVG/jax.js
Research article Special Issues

Urinary DNA lesions as a biomarker for assessing male reproductive health

  • Received: 26 March 2015 Accepted: 02 June 2015 Published: 15 June 2015
  • The study aimed to examine whether urinary 8-hydroxy-2'-deoxyguanosine (8-OHdG) could serve as a biomarker for assessing sperm DNA integrity. Urine and semen samples were collected from 124 coke-oven workers, who had chronically been exposed to polycyclic aromatic hydrocarbons (PAHs), from a steel plant in Kaohsiung, Taiwan. The coke-oven workers were divided into two subgroups: topside-oven workers and side-oven workers. Sperm DNA integrity was assessed in terms of DNA fragmentation and 8-oxo-7,8-dihydro-2'-deoxyguanosine (8-oxodGuo). Urine samples were used to detect 1-hydroxypyrene (1-OHP) and urinary 8-OHdG, which served to assess exposure to PAHs and the whole body oxidative stress, respectively. Urinary 1-OHP concentrations were significantly higher in the topside-oven workers than the side-oven workers (p = 0.02). Sperm 8-oxodGuo concentrations were correlated with percentages of sperm fragmentation (p = 0.044), but urinary 8-OHdG levels failed to correlate with sperm 8-oxodGuo levels and with percentages of fragmentation. In conclusion, exposure to PAHs was linked to increased oxidative stress on sperm. However, urinary 8-OHdG may not be a suitable biomarker for examining sperm DNA damage associated with oxidative stress.

    Citation: Hueiwang Anna Jeng, Chih-Hung Pan, Chao Mu-Rong. Urinary DNA lesions as a biomarker for assessing male reproductive health[J]. AIMS Environmental Science, 2015, 2(3): 565-576. doi: 10.3934/environsci.2015.3.565

    Related Papers:

    [1] M. D. König, Stefano Battiston, M. Napoletano, F. Schweitzer . On algebraic graph theory and the dynamics of innovation networks. Networks and Heterogeneous Media, 2008, 3(2): 201-219. doi: 10.3934/nhm.2008.3.201
    [2] Regino Criado, Julio Flores, Alejandro J. García del Amo, Miguel Romance . Structural properties of the line-graphs associated to directed networks. Networks and Heterogeneous Media, 2012, 7(3): 373-384. doi: 10.3934/nhm.2012.7.373
    [3] Valérie Dos Santos, Bernhard Maschke, Yann Le Gorrec . A Hamiltonian perspective to the stabilization of systems of two conservation laws. Networks and Heterogeneous Media, 2009, 4(2): 249-266. doi: 10.3934/nhm.2009.4.249
    [4] Alexandre M. Bayen, Alexander Keimer, Nils Müller . A proof of Kirchhoff's first law for hyperbolic conservation laws on networks. Networks and Heterogeneous Media, 2023, 18(4): 1799-1819. doi: 10.3934/nhm.2023078
    [5] Steinar Evje, Kenneth H. Karlsen . Hyperbolic-elliptic models for well-reservoir flow. Networks and Heterogeneous Media, 2006, 1(4): 639-673. doi: 10.3934/nhm.2006.1.639
    [6] Alessandro Gondolo, Fernando Guevara Vasquez . Characterization and synthesis of Rayleigh damped elastodynamic networks. Networks and Heterogeneous Media, 2014, 9(2): 299-314. doi: 10.3934/nhm.2014.9.299
    [7] Dilip Sarkar, Shridhar Kumar, Pratibhamoy Das, Higinio Ramos . Higher-order convergence analysis for interior and boundary layers in a semi-linear reaction-diffusion system networked by a $ k $-star graph with non-smooth source terms. Networks and Heterogeneous Media, 2024, 19(3): 1085-1115. doi: 10.3934/nhm.2024048
    [8] Nathaniel J. Merrill, Zheming An, Sean T. McQuade, Federica Garin, Karim Azer, Ruth E. Abrams, Benedetto Piccoli . Stability of metabolic networks via Linear-in-Flux-Expressions. Networks and Heterogeneous Media, 2019, 14(1): 101-130. doi: 10.3934/nhm.2019006
    [9] José Ignacio Alvarez-Hamelin, Luca Dall'Asta, Alain Barrat, Alessandro Vespignani . K-core decomposition of Internet graphs: hierarchies, self-similarity and measurement biases. Networks and Heterogeneous Media, 2008, 3(2): 371-393. doi: 10.3934/nhm.2008.3.371
    [10] Alessia Marigo, Benedetto Piccoli . A model for biological dynamic networks. Networks and Heterogeneous Media, 2011, 6(4): 647-663. doi: 10.3934/nhm.2011.6.647
  • The study aimed to examine whether urinary 8-hydroxy-2'-deoxyguanosine (8-OHdG) could serve as a biomarker for assessing sperm DNA integrity. Urine and semen samples were collected from 124 coke-oven workers, who had chronically been exposed to polycyclic aromatic hydrocarbons (PAHs), from a steel plant in Kaohsiung, Taiwan. The coke-oven workers were divided into two subgroups: topside-oven workers and side-oven workers. Sperm DNA integrity was assessed in terms of DNA fragmentation and 8-oxo-7,8-dihydro-2'-deoxyguanosine (8-oxodGuo). Urine samples were used to detect 1-hydroxypyrene (1-OHP) and urinary 8-OHdG, which served to assess exposure to PAHs and the whole body oxidative stress, respectively. Urinary 1-OHP concentrations were significantly higher in the topside-oven workers than the side-oven workers (p = 0.02). Sperm 8-oxodGuo concentrations were correlated with percentages of sperm fragmentation (p = 0.044), but urinary 8-OHdG levels failed to correlate with sperm 8-oxodGuo levels and with percentages of fragmentation. In conclusion, exposure to PAHs was linked to increased oxidative stress on sperm. However, urinary 8-OHdG may not be a suitable biomarker for examining sperm DNA damage associated with oxidative stress.


    In this paper we study systems of linear hyperbolic equations on a bounded interval, say, $ (0, 1) $, sometimes referred to as port-Hamiltonians, see for example [14], and the role is played by the boundary conditions coupling the incoming and outgoing Riemann invariants determined by the system at the endpoints $ x = 0 $ and $ x = 1 $. In particular, we consider the set of equations

    $ t(υϖ)=(C+00C)x(υϖ),0<x<1,t>0, $ (1a)
    $ Ξ(υ(0,t),ϖ(1,t),υ(1,t),ϖ(0,t))T=0,t>0, $ (1b)
    $ υ(x,0)=˚υ(x),ϖ(x,0)=˚ϖ(x)0<x<1, $ (1c)

    where $ \boldsymbol{ \upsilon} $ and $ \boldsymbol{ \varpi} $ are the Riemann invariants flowing from $ 0 $ to $ 1 $ and from $ 1 $ to $ 0 $, respectively, $ \mathcal{C}_+ $ and $ \mathcal{C}_- $ are $ m_+\times m_+ $ and $ m_-\times m_- $ diagonal matrices with positive entries and, with $ 2m = m_++m_- $, $ \boldsymbol{\Xi} $ is a $ 2m\times 4m $ matrix relating outgoing $ \boldsymbol{{ \upsilon}}(0), \boldsymbol{{ \varpi}}(1) $ and incoming $ \boldsymbol{{ \upsilon}}(1), \boldsymbol{{ \varpi}}(0) $ boundary values so that (1b) can be written as

    $ Ξout(υ(0,t),ϖ(1,t))T+Ξin(υ(1,t),ϖ(0,t))T=0,t>0. $ (2)

    An important class of such problems arises from dynamical systems on metric graphs. Let $ \Gamma $ be a graph with $ r $ vertices $ \{ \mathbf{v}_j\}_{1\leq j\leq r} = :\Upsilon $ and $ m $ edges $ \{\boldsymbol{e_j}\}_{1\leq j\leq m} $ (identified with $ (0, 1) $ through a suitable parametrization). The dynamics on each edge $ \boldsymbol{e_j} $ is described by

    $ tpj+Mjxpj=0,0<x<1,t>0,1jm, $ (3)

    where $ \boldsymbol{p}^j = (p^j_1, p^j_2)^T $ and $ \mathcal{M}^j = (M^j_{lk})_{1\leq k, l\leq 2} $ are defined on $ [0, 1] $ and $ \mathcal{M}^j(x) $ is a strictly hyperbolic real matrix for each $ x\in [0, 1] $ and $ 1\leq j\leq m $. System (3) is complemented with initial conditions and suitable transmission conditions coupling the values of $ \boldsymbol{p}^j $ at the vertices which the edges $ \boldsymbol{e_j} $ are incident to. Then, (1) can be obtained from (3) by diagonalization so that (suitably re-indexed) $ \boldsymbol{ \upsilon} $ and $ \boldsymbol{ \varpi} $ are the Riemann invariants (see [7,Section 1.1]) of $ \boldsymbol{p} = (\boldsymbol{p}^j)_{1\leq j\leq m} $.

    Such problems have been a subject of intensive research, both from the dynamics on graphs, [1,10,5,4,11,17,19], and the 1-D hyperbolic systems, [7,21,15,14], points of view. However, there is hardly any overlap, as there seems to be little interest in the network interpretation of the results in the latter, while in the former the conditions on the Riemann invariants seem to be "difficult to adapt to the case of a network", [11,Section 3].

    The main aim of this paper, as well as of the preceding one [2] is to bring together these two approaches. In [2] we have provided explicit formulae allowing for a systematic conversion of Kirchhoff's type network transmission conditions to (1b) in such a way that the resulting system(1) is well-posed. We also gave a proof of the well-posedness on any $ L_p $, $ 1\leq p<\infty $, which, in contrast to [21], is purely based on an operator semigroup approach. For notational clarity, we focused on $ 2\times 2 $ hyperbolic systems on each edge but the method works equally well for systems of arbitrary (finite) dimension. In this paper we are concerned with the reverse question, that is, to determine under what assumptions on $ \boldsymbol{\Xi} $, (1) describes a network dynamics given by $ 2\times 2 $ hyperbolic systems on each edge, coupled by Kirchhoff's transmission conditions at incident vertices.

    To briefly describe the content of the paper, we observe that if the matrix $ \boldsymbol{\Xi} = \{\xi_{ij}\}_{1\leq i\leq 2m, 1\leq j\leq 4m} $ in (1b) describes transmission conditions at the vertices of a graph, say $ \Gamma, $ on whose edges we have $ 2\times 2 $ systems of hyperbolic equations, then we should be able to group the indices $ j $ into pairs $ \{j', j"\} $ corresponding to the edges of $ \Gamma $ on which we have $ 2\times 2 $ systems for the components $ j' $ and $ j" $ of $ (\boldsymbol{ \upsilon}, \boldsymbol{ \varpi}). $ Thus, in a sense, the columns of $ \boldsymbol{\Xi} $ determine the edges of $ \Gamma. $ It follows that it is easier to split the reconstruction of $ \Gamma $ into two steps and first build a digraph $ \boldsymbol{\Gamma}, $ where each column index $ j $ of $ \boldsymbol{\Xi} $ is associated with an arc, say $ \boldsymbol{{\varepsilon}}^j, $ on which we have a first order system for either $ \upsilon_j $ or $ \varpi_j. $ Thus, the main problem is to construct vertices of $ \Gamma $ (and $ \boldsymbol{\Gamma} $) which should be somehow determined by a partition of the row indices of $ \boldsymbol{\Xi} $. To do this, we observe that the coefficients of $ \boldsymbol{\Xi} $ represent a map of connections of the edges in the sense that, roughly speaking, if $ \xi_{ij}\neq 0 $ and $ \xi_{ik}\neq 0 $, then arcs $ \boldsymbol{{\varepsilon}}^j $ and $ \boldsymbol{{\varepsilon}}^k $ are incident to the same vertex and, if they are incoming to it, then they cannot be incoming to any other vertex. A difficulty here is that while for the flow to occur from, say, $ \boldsymbol{{\varepsilon}}^j $ to $ \boldsymbol{{\varepsilon}}^k, $ these arcs must be incident to the same vertex but the converse may not hold, that is, for $ \boldsymbol{{\varepsilon}}^j $ incoming to and $ \boldsymbol{{\varepsilon}}^k $ outgoing from the same $ \mathbf{v}, $ the flow from $ \boldsymbol{{\varepsilon}}^j $ may not enter $ \boldsymbol{{\varepsilon}}^k $ but go to other outgoing arcs. To avoid such a case, in this paper we formulate conditions ensuring that the flow connectivity at each vertex is the same as the graph connectivity. This assumption yields a relatively simple criterion for the reconstruction of $ \boldsymbol{\Gamma}, $ which is that

    $ \widehat {\left(\widehat{\boldsymbol{\Xi}_{out}}\right)^T\widehat{\boldsymbol{\Xi}_{in}}} $

    is the adjacency matrix of a line graph (where for a matrix $ \mathcal{A} $, $ \widehat{\mathcal{A}} $ is obtained by replacing non-zero entries of $ \mathcal{A} $ by 1.) This, together with some technical assumptions, allows us to apply the theory of [13], see also [6,Theorem 4.5.1], to construct first $ \boldsymbol{\Gamma} $ and then $ \Gamma $ in such a way that (2) can be localized at each vertex of $ \Gamma $ in a way which is consistent with (1a).

    The main idea of this paper is similar to that of [3]. However, [3] dealt with first order problems with (2) solved with respect to the outgoing data. Here, we do not make this assumption and, while(1) technically is one-dimensional, having reconstructed $ \boldsymbol{\Gamma}, $ we still have to glue together its pairs of arcs to obtain the edges of $ \Gamma $ in such a way that the corresponding pairs of solutions of (1a) are Riemann invariants of $ 2\times 2 $ systems on $ \Gamma $. Another difficulty in the current setting is potential presence of sources and sinks in $ \Gamma. $ Their structure is not reflected in the line graph, [3], and reconstructing them in a way consistent with a system of $ 2\times 2 $ equations on $ \Gamma $ is technically involved.

    The paper is organized as follows. In Section 2 we briefly recall the notation and relevant results from [2]. Section 3 contains the main result of the paper. In Appendix we recall basic results on line graphs in the interpretation suitable for the considerations of the paper.

    We consider a network represented by a finite, connected and simple (without loops and multiple edges) metric graph $ \Gamma $ with $ r $ vertices $ \{ \mathbf{v}_j\}_{1\leq j\leq r} = :\Upsilon $ and $ m $ edges $ \{\boldsymbol{e_j}\}_{1\leq j\leq m} $. We denote by $ E_{ \mathbf{v}} $ the set of edges incident to $ \mathbf{v}, $ let $ J_{ \mathbf{v}} = \{j;\; \boldsymbol{e_j} \in E_{ \mathbf{v}}\} $ and $ |E_{ \mathbf{v}}| = |J_{ \mathbf{v}}| $ be the valency of $ \mathbf{v} $. We identify the edges with unit intervals through sufficiently smooth invertible functions $ l_j: \boldsymbol{e_j} \mapsto [0, 1] $. In particular, we call $ \mathbf{v} $ with $ l_j( \mathbf{v}) = 0 $ the tail of $ \boldsymbol{e_j} $ and the head if $ l_j( \mathbf{v}) = 1 $. On each edge $ \boldsymbol{e_j} $ we consider the system (3). Let $ \lambda^j_-< \lambda^j_+ $ be the eigenvalues of $ \mathcal{M}^j, 1\leq j\leq m $ (the strict inequality is justified by the strict hyperbolicity of $ \mathcal{M}^j $). The eigenvalues can be of the same sign as well as of different signs. In the latter case, we have $ \lambda^j_-<0< \lambda^j_+ $. By $ f^j_\pm = (f^j_{\pm, 1}, f^j_{\pm, 2})^T $ we denote the eigenvectors corresponding to $ \lambda^j_\pm $, respectively, and by

    $ \mathcal{F}^j = \left(fj+,1fj,1fj+,2fj,2\right), $

    the diagonalizing matrix on each edge. The Riemann invariants $ \boldsymbol{u}^j = (u^j_1, u^j_2)^T, 1\leq j\leq m, $ are defined by

    $ uj=(Fj)1pjandpj=(fj+,1uj1+fj,1uj2fj+,2uj1+fj,2uj2). $ (4)

    Then, we diagonalize (3) and, discarding lower order terms, we consider

    $ tuj=Ljxuj=(λj+00λj)xuj, $ (5)

    for each $ 1\leq j\leq m $.

    Remark 1. We refer an interested reader to [7,Section 1.1] for a detailed construction of the Riemann invariants for a general $ 1-D $ hyperbolic system and the explanation of the name.

    The most general linear local boundary conditions at $ \mathbf{v}\in \Upsilon $ are given by

    $ Φvp(v)=0, $ (6)

    where $ \boldsymbol{p} = ((p_1^j, p_2^j)_{1\leq j\leq m})^T $ and the real matrix $ \boldsymbol{\Phi}_{ \mathbf{v}} $ is given by

    $ Φv:=(ϕj1v,1φj1v,1ϕj|Jv|v,1φj|Jv|v,1ϕj1v,kvφj1v,kvϕj|Jv|v,kvφj|Jv|v,kv),  $ (7)

    where $ J_{ \mathbf{v}} = \{j_1, \ldots, j_{|J_{ \mathbf{v}}|}\} $ and $ k_{ \mathbf{v}} $ is a parameter determined by the problem. The difficulty with such a formulation is that it is not immediately clear what properties $ \boldsymbol{\Phi}_{ \mathbf{v}} $ should have to ensure well-posedness of the hyperbolic problem for which (6), $ \mathbf{v}\in \Upsilon $, serve as boundary conditions. There are various ways tackling this difficulty. For example, in [20,11] conditions are imposed directly on $ \boldsymbol{\Phi}_{ \mathbf{v}} $ to ensure specific properties, such as dissipativity, of the resulting initial boundary value problem. However, we follow the paradigm introduced in [7,Section 1.1.5.1] and require that at each vertex all outgoing data must be determined by the incoming data. Since for a general system (3) it is not always obvious which data are outgoing and which are incoming at a vertex, we write (6) in the equivalent form using the Riemann invariants $ \boldsymbol{u} = \mathcal{F}^{-1} \boldsymbol{p} $, as

    $ Ψvu(v):=ΦvF(v)u(v)=0. $ (8)

    For Riemann invariants, we can define their outgoing values at $ \mathbf{v} $ as follows.

    Definition 2.1. Let $ \mathbf{v}\in \Upsilon $. The following values $ u^j_k( \mathbf{v}), j\in J_{ \mathbf{v}}, k = 1, 2, $ are outgoing at $ \mathbf{v} $.

    If $ λ^j_+>λ^j_->0 $ $ λ^j_+>0>λ^j_- $ $ 0>λ^j_+>λ^j_- $
    $ l_j(\bf v) =0 $ $ u^j_1(\bf v), u^j_2(\bf v) $ $ u^j_1(\bf v) $ none
    $ l_j(\bf v) =1 $ none $ u^j_2(\bf v) $ $ u^j_1(\bf v), u^j_2(\bf v) $

     | Show Table
    DownLoad: CSV

    Denote by $ \alpha_j $ the number of positive eigenvalues on $ \boldsymbol{e_j} $. Then, we see that for a given vertex $ \mathbf{v} $ with valence $ |J_{ \mathbf{v}}| $ the number of outgoing values is given by

    $ kv:=jJv(2(1αj)lj(v)+αj). $ (9)

    Definition 2.2. We say that $ \mathbf{v} $ is

    ● a sink if either $ \alpha_j = 2 $ and $ l_j( \mathbf{v}) = 1 $ or $ \alpha_j = 0 $ and $ l_j( \mathbf{v}) = 0 $ for all $ j \in J_{ \mathbf{v}} $;

    ● a source if either $ \alpha_j = 0 $ and $ l_j( \mathbf{v}) = 1 $ or $ \alpha_j = 2 $ and $ l_j( \mathbf{v}) = 0 $ for all $ j \in J_{ \mathbf{v}} $;

    ● a transient (or internal) vertex if it is neither a source nor a sink.

    We denote the sets of sources, sinks and transient vertices by $ \Upsilon_s, \Upsilon_z $ and $ \Upsilon_t $, respectively.

    We observe that if $ \mathbf{v}\in \Upsilon_z $, then $ k_{ \mathbf{v}} = 0 $ (so that no boundary conditions are imposed at a sink), while if $ \mathbf{v}\in \Upsilon_s $, then $ k_{ \mathbf{v}} = 2|J_{ \mathbf{v}}| $.

    A typical example of (8) is Kirchhoff's law that requires that the total inflow rate into a vertex must equal the total outflow rate from it. Its precise formulation depends on the context, we refer to [8,Chapter 18] for a detailed description in the context of flows in networks. Since it provides only one equation, in general it is not sufficient to ensure the well-posedness of the problem. So, we introduce the following definition.

    Definition 2.3. We say that $ \boldsymbol{p} $ satisfies a generalized Kirchhoff conditions at $ \mathbf{v} \in \Upsilon\setminus\Upsilon_z $ if, for $ \boldsymbol{u} = \mathcal{F}^{-1} \boldsymbol{p}, $ (8) is satisfied for some matrix $ \boldsymbol{\Phi}_{ \mathbf{v}} = \boldsymbol{\Psi}_{ \mathbf{v}}\mathcal{F}^{-1} $ with $ k_{ \mathbf{v}} $ given by (9).

    To realize the requirement that the outgoing values should be determined by the incoming ones, we have to analyze the structure of $ \boldsymbol{\Psi}_{ \mathbf{v}} $. Let us introduce the partition

    $ {1,,m}=:J1J2J0, $ (10)

    where $ j \in J_1 $ if $ \alpha_j = 1, $ $ j \in J_2 $ if $ \alpha_j = 2 $ and $ j\in J_0 $ if $ \alpha_j = 0 $. This partition induces the corresponding partition of each $ J_{ \mathbf{v}} $ as

    $ J_{ \mathbf{v}} : = J_{ \mathbf{v}, 1}\cup J_{ \mathbf{v}, 2}\cup J_{ \mathbf{v}, 0}. $

    We also consider another partition $ J_{ \mathbf{v}} = J_{ \mathbf{v}}^0\cup J_{ \mathbf{v}}^1, $ where $ j\in J_{ \mathbf{v}}^0 $ if $ l_j( \mathbf{v}) = 0 $ and $ j\in J_{ \mathbf{v}}^1 $ if $ l_j( \mathbf{v}) = 1 $. Then, we can give an alternative expression for $ k_{ \mathbf{v}} $ as

    $ kv=jJ0vαj+jJ1v(2αj)=|Jv,1|+2(|J0vJv,2|+|J1vJv,0|). $ (11)

    Then, by [2,Lemma 3.6],

    $\rm (i) $ $ u^j_1(0) $ is outgoing if and only if $ j \in (J_{ \mathbf{v}, 1}\cup J_{ \mathbf{v}, 2})\cap J^0_{ \mathbf{v}}, $

    $\rm (ii) $ $ u^j_2(0) $ is outgoing if and only if $ j \in J_{ \mathbf{v}, 2}\cap J^0_{ \mathbf{v}}, $

    $ \rm(iii) $ $ u^j_1(1) $ is outgoing if and only if $ j \in J_{ \mathbf{v}, 0}\cap J^1_{ \mathbf{v}}, $

    $\rm (iv) $ $ u^j_2(1) $ is outgoing if and only if $ j \in (J_{ \mathbf{v}, 1}\cup J_{ \mathbf{v}, 0})\cap J^1_{ \mathbf{v}}. $

    We introduce the block diagonal matrix

    $ ˜Fout(v)=diag{˜Fjout(v)}jJv, $ (12)

    where

    $ \mathcal{{\tilde{F}}}_{out}^j( \mathbf{v}) = \left\{\begin{array} {ccc} \left(\begin{array}{cc}0&0\\0&0\end{array}\right)&\text{if} &j\in (J_{ \mathbf{v}, 0} \cap J_{ \mathbf{v}}^0)\cup (J_{ \mathbf{v}, 2}\cap J_{ \mathbf{v}}^1), \\ \left(fj+,1(lj(v))fj,1(lj(v))fj+,2(lj(v))fj,2(lj(v))\right)&\text{if}& j\in (J_{ \mathbf{v}, 0} \cap J_{ \mathbf{v}}^1)\cup (J_{ \mathbf{v}, 2}\cap J_{ \mathbf{v}}^0), \\ \left(fj+,1(0)0fj+,2(0)0\right)&\text{if} &j\in J_{ \mathbf{v}, 1} \cap J_{ \mathbf{v}}^0, \\ \left(0fj,1(1)0fj,2(1)\right)&\text{if}& j\in J_{ \mathbf{v}, 1} \cap J_{ \mathbf{v}}^1. \end{array} \right. $

    Further, by $ \mathcal{F}_{out}( \mathbf{v}) $ we denote the contraction of $ \mathcal{{\tilde{F}}}_{out}( \mathbf{v}) $; that is, the $ 2|J_{ \mathbf{v}}| \times k_{ \mathbf{v}} $ matrix obtained from $ \mathcal{{\tilde{F}}}_{out}( \mathbf{v}) $ by deleting $ 2|J_{ \mathbf{v}}| - k_{ \mathbf{v}} $ zero columns, and then define $ \mathcal{F}_{in}( \mathbf{v}) $ as the analogous contraction of $ \mathcal{F}( \mathbf{v})-\mathcal{{\tilde{F}}}_{out}( \mathbf{v}) $.

    In a similar way, we extract from $ \boldsymbol{u}( \mathbf{v}) $ the outgoing boundary values $ \widetilde{\boldsymbol{u}}_{out}( \mathbf{v}) = (\widetilde{\boldsymbol{u}}^j_{out}( \mathbf{v}))_{j\in J_{ \mathbf{v}}} $ by

    $ \widetilde{\boldsymbol{u}}^j_{out}( \mathbf{v}) = \left\{(0,0)Tifj(Jv,0J0v)(Jv,2J1v),(uj1(lj(v)),uj2(lj(v)))Tifj(Jv,0J1v)(Jv,2J0v),(uj1(0),0)TifjJv,1J0v,(0,uj2(1))TifjJv,1J1v, \right. $

    and $ \widetilde{\boldsymbol{u}}_{in}( \mathbf{v}) = {\boldsymbol{u}( \mathbf{v})}-\widetilde{\boldsymbol{u}}_{out}( \mathbf{v}). $ As above, we define $ \boldsymbol{u}_{out}( \mathbf{v}) $ to be the vector in $ \mathbb{R}^{k_{ \mathbf{v}}} $ obtained by discarding the zero entries in $ \widetilde{\boldsymbol{u}}_{out}( \mathbf{v}) $, as described above and, similarly, $ \boldsymbol{u}_{in}( \mathbf{v}) $ is the vector in $ \mathbb{R}^{2|J_{ \mathbf{v}}|-k_{ \mathbf{v}}} $ obtained from $ \widetilde{\boldsymbol{u}}_{in}( \mathbf{v}) $.

    Proposition 1. [2,Proposition 3.8] The boundary system (8) at $ \mathbf{v}\in \Upsilon\setminus \Upsilon_z $ is equivalent to

    $ ΦvFout(v)uout(v)+ΦvFin(v)uin(v)=0 $ (13)

    and hence it uniquely determines the outgoing values of $ \boldsymbol{u}( \mathbf{v}) $ at $ \mathbf{v} $ as defined by Definition 2.1 if and only if

    $ ΦvFout(v)is nonsingular. $ (14)

    In this case,

    $ uout(v)=(ΦvFout(v))1ΦvFin(v)uin(v). $ (15)

    To pass from (3) with Kirchhoff's boundary conditions at each vertex $ \mathbf{v} \in \Upsilon\setminus \Upsilon_z $ to(1) we have to write the former in a global form. Assuming the vertices in $ \Upsilon\setminus\Upsilon_z $ are ordered as $ \{ \mathbf{v}_1, \ldots, \mathbf{v}_{r'}\} $, we define $ \boldsymbol{\Psi}' = {\rm diag} \{\boldsymbol{\Psi}_{ \mathbf{v}} \}_{ \mathbf{v}\in \Upsilon\setminus\Upsilon_z} $ and $ \gamma \boldsymbol{u} = ((\boldsymbol{u}( \mathbf{v}))_{ \mathbf{v} \in \Upsilon\setminus\Upsilon_z})^T $ and write (8) as

    $ Ψγu=0. $ (16)

    We note that the function values that are incoming at $ \mathbf{v}\in \Upsilon_z $ do not influence any outgoing data. However, to keep the track of all vertex values, we extend $ \boldsymbol{\Psi}' $ with zero columns corresponding to edges coming to sinks and denote such an extended matrix by $ \boldsymbol{\Psi} $. Since, by the hand shake lemma, we have $ 2\sum_{ \mathbf{v}\in \Upsilon} |J_{ \mathbf{v}}| = 4m $ and by [2,Section 3.2] also $ \sum_{ \mathbf{v}\in \Upsilon\setminus\Upsilon_z} k_{ \mathbf{v}} = 2m, $ $ \boldsymbol{\Psi} $ is a $ 2m\times 4m $ matrix. In the same way, we can provide a global form of (13), splitting (16) as

    $ Ψoutγuout+Ψinγuin=0, $ (17)

    where $ \boldsymbol{\Psi}^{out} = {\rm diag} \{\boldsymbol{\Phi}_{ \mathbf{v}} \mathcal{F}_{out}( \mathbf{v}) \}_{ \mathbf{v}\in \Upsilon\setminus\Upsilon_z} $ and $ \boldsymbol{\Psi}^{in} = {\rm diag} \{\boldsymbol{\Phi}_{ \mathbf{v}} \mathcal{F}_{in}( \mathbf{v}) \}_{ \mathbf{v}\in \Upsilon\setminus\Upsilon_z}, $ extended by zero columns corresponding to the incoming functions at the sinks, $ \gamma \boldsymbol{u}_{out} : = ((\boldsymbol{u}_{out}( \mathbf{v}))_{ \mathbf{v} \in \Upsilon\setminus\Upsilon_z})^T $, and $ \gamma \boldsymbol{u}_{in} $ is $ ((\boldsymbol{u}_{in}( \mathbf{v}))_{ \mathbf{v} \in \Upsilon\setminus\Upsilon_z})^T $ extended by incoming values at the sinks.

    Using the adopted parametrization and the formalism of Definition 2.1, we only need to distinguish between functions describing the flow from $ 0 $ to $ 1 $ and from $ 1 $ to $ 0 $. Accordingly, we group the Riemann invariants $ \boldsymbol{u} $ into parts corresponding to positive and negative eigenvalues and rename them as:

    $ υ:=((uj1)jJ1J2,(uj2)jJ2)=(υj)jJ+,ϖ:=((uj1)jJ0,(uj2)jJ1J0)=(ϖj)jJ, $ (18)

    where $ J^+ $ and $ J^- $ are the sets of indices $ j $ with at least 1 positive eigenvalue, and at least 1 negative eigenvalue of $ \mathcal{M}^j $, respectively. In $ J^+ $ (respectively $ J^- $) the indices from $ J_2 $ (respectively $ J_0 $) appear twice so that we rearrange them in some consistent way to avoid confusion. For instance, we can take $ J^+ = \{1, \ldots, m^u, m^u+1, \ldots, m^+\} $ and $ J^- = \{m^++1, \ldots, m_u, m_u+1, \ldots, 2m\} $ and there are bijections between $ J_1\cup J_2 $ and $ \{1, \ldots, m^u\} $, $ J_2 $ and $ \{m^u+1, \ldots, m^+\} $, $ J_0 $ and $ \{m^++1, \ldots, m_u\} $, and $ J_1\cup J_0 $ and $ \{m_u+1, \ldots, 2m\} $, respectively. We emphasize that such a renumbering is largely arbitrary and different ways of doing it result in just re-labelling of the components of(1) without changing its structure.

    In this way, we converted $ \Gamma $ into a multi digraph $ \boldsymbol{\Gamma} $ with the same vertices $ \Upsilon $, in such a way that each edge of $ \Gamma $ was split into two arcs parametrized by $ x\in [0, 1], $ where $ x = 0 $ on each arc corresponds to the same vertex in $ \Gamma $ and the same is valid for $ x = 1 $. Conversely, if we have a multi digraph $ \boldsymbol{\Gamma}, $ where all edges appear in pairs and each two edges joining the same vertex are parametrized concurrently, then we can collapse $ \boldsymbol{{\Gamma}} $ to a graph $ \Gamma $. We note that this approach originates from [16].

    Using this construction, the second order hyperbolic problem (3), (17) was transformed into first order system (1) with (17) written in the form (2). However, it is clear that (1) can be formulated with an arbitrary matrix $ \boldsymbol{\Xi} $. Thus, we arrive at the main problem considered in this paper:

    how to characterize matrices $\boldsymbol{\Xi}$ that arise from $\boldsymbol{\Psi}$ so that (1) describes a network dynamics?

    For a graph $ \Gamma $, let us consider the multi digraph $ \boldsymbol{\Gamma} $ constructed above. The sets of vertices $ \Upsilon $ are the same for $ \Gamma $ and $ \boldsymbol{\Gamma} $. For $ \mathbf{v} \in\Upsilon $ of $ \boldsymbol{{\Gamma}}, $ we can talk about incoming and outgoing arcs which are determined by $ j\in J_{ \mathbf{v}}, $ $ l_j( \mathbf{v}) $ and the signs of $ \lambda_+^j $ and $ \lambda_-^j $, as in Definition 2.1. We denote by $ \boldsymbol{J}_{ \mathbf{v}}^+ $ and $ \boldsymbol{J}_{ \mathbf{v}}^- $ the (ordered) sets of indices of arcs $ \boldsymbol{{\varepsilon}}^j $ incoming and outgoing from $ \mathbf{v} $ in $ \boldsymbol{{\Gamma}} $, respectively. We note that $ |\boldsymbol{J}_{ \mathbf{v}}^-| = k_{ \mathbf{v}} $, the number of the outgoing conditions. With this notation, the matrix $ \boldsymbol{{\Psi_{ \mathbf{v}}}} $ can be split into two matrices

    $ \boldsymbol{{\Psi}}_{ \mathbf{v}}^{out} = (\psi_{ \mathbf{v}, i}^j)_{1\leq i\leq k_{ \mathbf{v}}, j\in \boldsymbol{J}_{ \mathbf{v}}^-}, \qquad \boldsymbol{{\Psi}}_{ \mathbf{v}}^{in} = (\psi_{ \mathbf{v}, i}^j)_{1\leq i\leq k_{ \mathbf{v}}, j\in \boldsymbol{J}_{ \mathbf{v}}^+}. $

    Since no outgoing value should be missing, we adopt the following

    Assumption 1. No column or row of $ \boldsymbol{\Psi}_{ \mathbf{v}}^{out} $ is identically zero.

    These matrices provide some insight into how the arcs are connected by the flow which is an additional feature, superimposed on the geometric structure of the incoming and outgoing arcs at the vertex. In principle, these two structures do not have to be the same, that is, it may happen that the substance flowing from $ \boldsymbol{{\varepsilon}}^j $, $ j \in \boldsymbol{J}_{ \mathbf{v}}^+, $ is only directed to some of the outgoing arcs. An extreme case of such a situation is when both $ \boldsymbol{{\Psi}}_{ \mathbf{v}}^{out} $ and $ \boldsymbol{{\Psi}}_{ \mathbf{v}}^{in} $ are completely decomposable, see [9], with blocks in both matrices having the same row indices. Then, from the flow point of view, $ \mathbf{v} $ can be regarded as several nodes of the flow network, which are not linked with each other. Such cases, where the geometric structure at a vertex is inconsistent with the flow structure, may generate problems in determining the graph underlying transport problems. Thus in this paper we adopt assumptions ensuring that the map of the flow connections given by the matrices $ \boldsymbol{{\Psi}}_{ \mathbf{v}}^{out} $ and $ \boldsymbol{{\Psi}}_{ \mathbf{v}}^{in} $ coincides with the geometry at $ \mathbf{v} $. We begin with the necessary definitions.

    Definition 3.1. Let $ \mathbf{v}\in\Upsilon_t $. We say that an arc $ \boldsymbol{{\varepsilon}}^j $, $ j\in \boldsymbol{J}^+_{ \mathbf{v}}, $ flow connects to $ \boldsymbol{{\varepsilon}}^l, l \in \boldsymbol{J}^-_{ \mathbf{v}}, $ if $ \psi_{ \mathbf{v}, i}^j \neq 0 $ and $ \psi_{ \mathbf{v}, i}^l\neq 0 $ for some $ 1\leq i\leq k_{ \mathbf{v}}. $

    Using this idea, we define a connectivity matrix $ \mathsf C_{ \mathbf{v}} = (\mathsf{c}_{ \mathbf{v}, lj})_{l \in \boldsymbol{J}^-_{ \mathbf{v}}, j \in \boldsymbol{J}^+_{ \mathbf{v}}} $ by

    $ \mathsf{c}_{ \mathbf{v}, lj} = \left\{1ifεjflow connects toεl,0otherwise. \right. $

    Remark 2. We observe that

    ● the above definition implies that for $ \boldsymbol{{\varepsilon}}^j $ and $ \boldsymbol{{\varepsilon}}^l $ to be flow connected, $ \boldsymbol{{\varepsilon}}^j $ and $ \boldsymbol{{\varepsilon}}^l $ must be incident to the same vertex;

    $ \mathsf C_{ \mathbf{v}} $ can be interpreted as the adjacency matrix of the bipartite line digraph constructed from the incoming and outgoing arcs at $ \mathbf{v}, $ where the connections between arcs are defined by flow connections, see [9,Section 3].

    For an arbitrary matrix $ A = (a_{ij})_{1\leq i\leq p, 1\leq j\leq q}, $ by $ \widehat A = (\hat a_{ij})_{1\leq i\leq p, 1\leq j\leq q} $ we denote the matrix with every nonzero entry of $ A $ replaced by 1.

    Lemma 3.2. If $ \mathbf{v} $ is a transient vertex, then

    $ Cv=^(^Ψoutv)T^Ψinv. $ (19)

    Proof. Denote $ \mathsf{B} = \widehat{\left(\widehat{\boldsymbol{{\Psi}}^{out}_{ \mathbf{v}}}\right)^T \widehat{\boldsymbol{{\Psi}}^{in}_{ \mathbf{v}}}}. $ Then, $ \mathsf{b}_{ij} = 1, i\in \boldsymbol{J}_{ \mathbf{v}}^-, j \in \boldsymbol{J}_{ \mathbf{v}}^+ $ if and only if

    $ \sum\limits_{r = 1}^{k_{ \mathbf{v}}} \hat{\psi}_{ \mathbf{v}, r}^i\hat{\psi}^j_{ \mathbf{v}, r} \neq 0. $

    This occurs if and only if there is $ r = 1, \ldots, k_{ \mathbf{v}} $ such that both $ \hat{\psi}_{ \mathbf{v}, r}^i\neq 0 $ and $ \hat{\psi}^j_{ \mathbf{v}, r} \neq 0, $ which is equivalent to $ \boldsymbol{{\varepsilon}}^j $ flow connecting with $ \boldsymbol{{\varepsilon}}^i $, that is, $ \mathsf{c}_{ \mathbf{v}, ij} = 1 $.

    Let $ \mathbf{v} $ be a source (as we do not impose boundary conditions on sinks). As above, we need to ensure that the flow from a source cannot be split into several isolated subflows. Though here we do not have inflows and outflows, we use a similar idea to that for transient vertices.

    Definition 3.3. Let $ \mathbf{v}\in\Upsilon_s $. We say that $ \boldsymbol{{\varepsilon}}^i $ and $ \boldsymbol{{\varepsilon}}^j $, $ i, j \in \boldsymbol{J}_{ \mathbf{v}}^-, $ are flow connected if there is $ l \in\{1, \ldots, k_{ \mathbf{v}}\} $ such that $ \psi_{ \mathbf{v}, l}^i\neq 0 $ and $ \psi_{ \mathbf{v}, l}^j\neq 0. $

    As before, we construct a connectivity matrix $ \mathsf C_{ \mathbf{v}} = (\mathsf{c}_{ \mathbf{v}, ij})_{i, j \in \boldsymbol{J}^-_{ \mathbf{v}}}, $ where

    $ cv,ij={1ifεjandεiare flow connected,0otherwise. $ (20)

    Note that, contrary to an internal vertex, here the connectivity matrix is symmetric. We also do not stipulate that $ i\neq j $ so that $ \boldsymbol{{\varepsilon}}^j $ is always flow connected to itself and hence, by Assumption 1, each entry of the diagonal of $ \mathsf{C}_{ \mathbf{v}} $ is 1. From this, we get a result similar to Lemma 3.2.

    Lemma 3.4. If $ \mathbf{v} $ is a source,

    $ Cv=^(^Ψoutv)T^Ψoutv. $ (21)

    Proof. As before, let $ \mathsf{B} = \widehat{\left(\widehat{\boldsymbol{{\Psi}}^{out}_{ \mathbf{v}}}\right)^T \widehat{\boldsymbol{{\Psi}}^{out}_{ \mathbf{v}}}}. $ Then, $ \mathsf{b}_{ij} = 1, i, j\in \boldsymbol{J}_{ \mathbf{v}}^-, $ if and only if

    $ \sum\limits_{r = 1}^{k_{ \mathbf{v}}} \hat{\psi}_{ \mathbf{v}, r}^i\hat{\psi}^j_{ \mathbf{v}, r} \neq 0. $

    Certainly, by Assumption 1, $ \mathsf{b}_{ii} = 1, i\in \boldsymbol{J}_{ \mathbf{v}}^-. $ For $ i\neq j $, this occurs if and only if there is $ r\in\{1, \ldots, k_{ \mathbf{v}}\} $ such that both $ \hat{\psi}_{ \mathbf{v}, r}^i\neq 0 $ and $ \hat{\psi}^j_{ \mathbf{v}, r} \neq 0 $ which is equivalent to $ \boldsymbol{{\varepsilon}}^j $ and $ \boldsymbol{{\varepsilon}}^i $ being flow connected, that is, $ \mathsf{c}_{ \mathbf{v}, ij} = 1 $.

    We adopt an assumption that the structure of flow connectivity is the same as of the geometry at the vertex. Thus, if $ \mathbf{v} $ is an internal vertex and $ j\in \boldsymbol{J}_{ \mathbf{v}}^+ $ and $ i\in \boldsymbol{J}^-_{ \mathbf{v}} $, then $ \boldsymbol{{\varepsilon}}^j $ flow connects to $ \boldsymbol{{\varepsilon}}^i $. In particular, we have

    Assumption 2. For all $ \mathbf{v}\in\Upsilon_t $,

    $ \mathsf{C}_{ \mathbf{v}} = \boldsymbol{1}_{ \mathbf{v}} = \left(111111\right). $

    We observe that the dimension of $ \mathsf{C}_{ \mathbf{v}} $ is $ |\boldsymbol{J}_{ \mathbf{v}}^-|\times |\boldsymbol{J}^+_{ \mathbf{v}}| $.

    If $ \mathbf{v} \in \Upsilon_s $, then we assume that the outflow from $ \mathbf{v} $ cannot be separated into independent subflows, that is, that the arcs outgoing from $ \mathbf{v} $ cannot be divided into groups such that no arc in any group is flow connected to an arc in any other. Equivalently, for each two arcs $ \boldsymbol{{\varepsilon}}^i $ and $ \boldsymbol{{\varepsilon}}^j $, $ i, j \in \boldsymbol{J}_{ \mathbf{v}}^-, $ there is a sequence $ j = j_0, j_1, \ldots, j_k = i $ such that $ \boldsymbol{{\varepsilon}}^{j_r} $ and $ \boldsymbol{{\varepsilon}}^{j_{r+1}} $, $ r = 0, \ldots, k-1 $ are flow connected. Indeed, if such a division was possible, then it would be impossible to find such a sequence between indices $ j $ and $ i $ in different groups as some pair would have to connect arcs from these different groups. Conversely, if for some arcs $ \boldsymbol{{\varepsilon}}^i $ and $ \boldsymbol{{\varepsilon}}^j $ there is no such a sequence, then we can build two groups of indices containing $ i $ and $ j $, respectively, by considering all indices for which such sequences can be found. Clearly, no arc in the first group is flow connected to any arc in the second as otherwise there would be a sequence connecting $ \boldsymbol{{\varepsilon}}^j $ and $ \boldsymbol{{\varepsilon}}^i $. By Lemma 3.4, $ \mathsf{C}_{ \mathbf{v}} $ can be considered as the adjacency matrix of the graph with vertices given by $ \{ \boldsymbol{{\varepsilon}}^j\}_{j \in \boldsymbol{J}^-_{ \mathbf{v}}} $ and the edges determined by the flow connectivity (20). Moreover, $ \mathsf{C}_{ \mathbf{v}} $ is symmetric, which shows that the assumption described at the beginning of this paragraph is equivalent to

    Assumption 3. For all $ \mathbf{v}\in\Upsilon_s, \mathsf{C}_{ \mathbf{v}} $ is irreducible.

    Remark 3. Assumption 3 is weaker than requiring each two arcs from $ \{ \boldsymbol{{\varepsilon}}^j\}_{j \in \boldsymbol{J}^-_{ \mathbf{v}}} $ to be flow connected. Then we would have $ \mathsf{C}_{ \mathbf{v}} = \boldsymbol{1}_{|\boldsymbol{J}_{ \mathbf{v}}^-|\times |\boldsymbol{J}^-_{ \mathbf{v}}|}. $

    Proposition 2. Let $ \mathbf{v}\in \Upsilon_t $. If the system (8), that is,

    $ Ψvu(v)=0, $ (22)

    contains a Kirchhoff's condition

    $ jJv(ψjv,ruj1(v)+ψjv,ruj2(v))=0, $ (23)

    with $ \psi^j_{ \mathbf{v}, r} \neq 0 $ for all $ j\in J_{ \mathbf{v}} $ and some $ r \in\{1, \ldots, k_{ \mathbf{v}}\} $, then Assumption 2 is satisfied.

    Proof. Condition (23) ensures that each entry of the $ r $-th row of both $ \widehat{\boldsymbol{{\Psi}}_{ \mathbf{v}}^{out}} $ and $ \widehat{\boldsymbol{{\Psi}}_{ \mathbf{v}}^{in}} $ is 1 and thus the product of each column of $ \widehat{\boldsymbol{{\Psi}}_{ \mathbf{v}}^{out}} $ with each column of $ \widehat{\boldsymbol{{\Psi}}_{ \mathbf{v}}^{in}} $ is non-zero, which yields Assumption 2.

    Example 1. Consider the model of [20], analysed in the framework of our approach in [2,Example 5.12], i.e.,

    $ tpj1+Kjxpj2=0,tpj2+Ljxpj1=0, $ (24)

    for $ t>0, 0<x<1, 0\leq j\leq m, $ where $ K^j>0, L^j>0 $ for all $ j $. For a given vertex $ \mathbf{v}, $ we define $ (\boldsymbol{p_1}( \mathbf{v}), \boldsymbol{p_2}( \mathbf{v})) = ((p^j_1( \mathbf{v}), p_2^j( \mathbf{v}))_{j\in J_{ \mathbf{v}}}, $ $ \nu^j( \mathbf{v}) = -1 $ if $ l_j( \mathbf{v}) = 0 $ and $ \nu^j( \mathbf{v}) = 1 $ if $ l_j( \mathbf{v}) = 1, $ and $ T_{ \mathbf{v}} \boldsymbol{p_2} ( \mathbf{v}) = (\nu^j( \mathbf{v})p^j_2( \mathbf{v}))_{j\in J_{ \mathbf{v}}}. $ In this case $ \alpha_j = 1 $ for any $ j $ and thus for any vertex $ \mathbf{v} $ we need $ |J_{ \mathbf{v}}| $ boundary conditions. We focus on $ \mathbf{v} $ with $ |E_{ \mathbf{v}}|>1 $. Then, we split $ \mathbb{R}^{|J_{ \mathbf{v}}|} $ into $ X_{ \mathbf{v}} $ of dimension $ n_{ \mathbf{v}} $ and its orthogonal complement $ X_{ \mathbf{v}}^\perp $ of dimension $ l_{ \mathbf{v}} = |J_{ \mathbf{v}}|-n_{ \mathbf{v}} $ and require that

    $ \boldsymbol{p_1}( \mathbf{v}) \in X_{ \mathbf{v}}, \quad T_{ \mathbf{v}}\boldsymbol{p_2}( \mathbf{v}) \in X^\perp_{ \mathbf{v}}, $

    that is, denoting $ I_1 = \{1, \ldots, n_{ \mathbf{v}}\} $ and $ I_2 = \{n_{ \mathbf{v}}+1, \ldots, |J_{ \mathbf{v}}|\} $,

    $ jJvϕjrpj1(v)=0,rI2,jJvφjrνj(v)pj2(v)=0,rI1, $ (25)

    where $ (( \varphi^j_r)_{j\in J_{ \mathbf{v}}})_{r\in I_1} $ is a basis in $ X_{ \mathbf{v}} $ and $ ((\phi^j_r)_{j\in J_{ \mathbf{v}}})_{r\in I_2} $ is a basis in $ X^\perp_{ \mathbf{v}}. $ It is clear that, in general, boundary conditions (25) do not satisfy Assumption 2. Consider $ \mathbf{v} $ such that each $ \boldsymbol{e}^j $ incident to $ \mathbf{v} $ is parameterised so as $ l_j( \mathbf{v}) = 0 $ so that each $ u^j_1(0) $ is outgoing and each $ u^j_2(0) $ is incoming. If we take $ \phi_r^j = \varphi^j_r = \delta_{rj} $ and $ L^j = K^j = 1 $ for $ j\in J_{ \mathbf{v}} $, we obtain

    $ pr1(0)=ur1(0)+ur2(0)=0,r=nv+1,,|Jv|,pr2(0)=ur1(0)ur2(0)=0,r=1,,nv. $

    Thus $ \widehat{\boldsymbol{\Psi}_{ \mathbf{v}}^{out}} $ and $ \widehat{\boldsymbol{\Psi}_{ \mathbf{v}}^{in}} $ are both the identity matrices and Assumption 2 is not satisfied.

    On the other hand, the Kirchhoff condition,

    $ jJvνj(v)pj2(v)=0, $ (26)

    see [20,Eqn (4)], satisfies the assumption of Proposition 2, as we have

    $ 0=jJvνj(v)pj2(v)=jJvνj(v)(fj+,2(v)uj1(v)+fj,2(v)uj2(v))=jJvνj(v)KjLj(uj1(v)uj2(v))=jJ0vKjLjuj1(0)jJ1vKjLjuj2(1)x+jJ1vKjLjuj1(1)+jJ0vKjLjuj2(0), $

    where we used [2,Eqn 5.2]. Hence, by Proposition 2, Assumption 2 is satisfied.

    Example 2. Let us consider the linearized Saint-Venant system,

    $ tpj1=Vjxpj1Hjxpj2,tpj2=gxpj1Vjxpj2, $ (27)

    see [2,Example 1.2], assuming that on each edge we have $ \lambda^j_{\pm} = V^j \pm \sqrt{gH^j} >0 $. Then, we have

    $ (pj1pj2)=(fj+,1uj1+fj,1uj2fj+,2uj1+fj,2uj2)=(Hjuj1+Hjuj2gHjuj1gHjuj2). $ (28)

    We use the flow structure of [11,Example 5.1], shown in Fig. 1, and focus on $ \mathbf{v}_1, $ where we need $ 2N-2 $ boundary conditions which were given as

    $ p^j_1(0) = p^1_1(1), \quad p^j_2(0) = p_2^1(1), \qquad j = 2, \ldots, N. $
    Figure 1. 

    Starlike network of channels

    .

    In terms of the Riemann invariants, they can be written as

    $ (H2H20000gH2gH200000000HNHN0000gHNgHN)(u21(0)u22(0)uN1(0)uN2(0))xxxx=(H1H1gH1gH1H1H1gH1gH1)(u11(1)u12(1)) $

    and it is clear that Assumption 2 is satisfied.

    Figure 2. 

    The reconstructed multi digraph $ \boldsymbol{\Gamma} $. It is seen that it cannot describe a flow on $ \Gamma $ as $ \varpi_5 $ and $ \varpi_6 $ must flow in the same direction

    .
    Figure 3. 

    The reconstructed multi digraph $ \boldsymbol{\Gamma} $ for (46), (47)

    .
    Figure 4. 

    A network $ \Gamma $ realizing the flow (48), (49)

    .

    For a matrix $ A = (a_{ij})_{1\leq i\leq n, 1\leq j\leq m}, $ let us denote by $ \boldsymbol{a}^c_j, 1\leq j\leq m, $ the columns of $ A $ and by $ \boldsymbol{a}^r_i, 1\leq i\leq n, $ its rows. Then, we often write

    $ A=(acj)1jm=(ari)1in, $ (29)

    that is, we represent the matrix as a row vector of its columns or a column vector of its rows. In particular, we write

    $ Ξout=(ξoutij)1i2m,1j2m=(ξout,cj)1j2m=(ξout,ri)1i2m,Ξin=(ξinij)1i2m,1j2m=(ξin,cj)1j2m=(ξin,ri)1i2m. $

    For any vector $ \boldsymbol{\mu} = (\mu_1, \ldots, \mu_k) $, we define $ \text{supp}\;\boldsymbol{\mu} = \{j\in \{1, \ldots, k\};\; \mu_j\neq 0\} $.

    Definition 3.5. We say that the problem (1) is graph realizable if there is a graph $ \Gamma = \{\{ \mathbf{v}_i\}_{1\leq i\leq r}, \{\boldsymbol{e_k}\}_{1\leq k\leq m}\} $ and a grouping of the column indices of $ \boldsymbol{\Xi} $ into pairs $ (j'_k, j_k^{"})_{1\leq k\leq m} $ such that (1a) describes a flow along the edges $ \boldsymbol{e_k} $ of $ \Gamma $, which satisfies generalized Kirchhoff's condition at each vertex of $ \Gamma $. In other words, (1) is graph realizable if there is a graph $ \Gamma $ and a matrix $ \boldsymbol{\Psi} $ such that (1a), (1c) can be written, after possibly permutating rows and columns of $ \boldsymbol{\Xi}, $ as (5), (17).

    Before we formulate the main theorem, we need to introduce some notation. Let us recall that we consider the boundary system (2), i.e.,

    $ \boldsymbol{\Xi}_{out}(( \upsilon_j(0, t))_{j\in J^+}, ( \varpi_j(1, t))_{j\in J^-}) = - \boldsymbol{\Xi}_{in}(( \upsilon_j(1, t))_{j\in J^+}, ( \varpi_j(0, t))_{j\in J^-}). $

    Let us emphasize that in this notation, the column indices on the left and right hand side correspond to the values of the same function. To shorten notation, let us renumber them as $ 1\leq j\leq 2m $. As noted in Introduction, appropriate pairs of the columns would determine the edges of the graph $ \Gamma $ that we try to reconstruct, hence the first step is to identify the possible vertices of $ \Gamma $. For this, first, we will try to construct a multi digraph $ \boldsymbol{\Gamma} $ on which (2) can be written in the form (17) for $ (\boldsymbol{ \upsilon}, \boldsymbol{ \varpi}) $. Roughly speaking, this corresponds to $ \boldsymbol{\Xi}_{out} $ and $ \boldsymbol{\Xi}_{in} $ being composed (up to permutations) of non-communicating blocks corresponding to the vertices. Here, each $ j $ should correspond to an arc $ \boldsymbol{{\varepsilon}}^j $ and the column $ j $ on the left hand side corresponds to the outflow along $ \boldsymbol{{\varepsilon}}^j $ from a unique vertex, while the column $ j $ on the right hand side corresponds to the inflow along $ \boldsymbol{{\varepsilon}}^j $ to a unique vertex. The vertices of $ \boldsymbol{\Gamma} $ should be then determined by a suitable partition of the rows of $ \boldsymbol{\Xi}_{out} $ (and of $ \boldsymbol{\Xi}_{in} $).

    In the second step we will determine additional assumptions that allow $ \boldsymbol{\Gamma} $ to be collapsed into a graph $ \Gamma $ on which (2) can be written in the form (17).

    Since we do not want (2) to be under- or over-determined, we adopt

    Assumption 4. For all $ 1\leq j\leq 2m $,

    $ \boldsymbol{\xi}^{out, c}_j\neq 0\quad\mathit{\text{and}}\quad \boldsymbol{\xi}^{out, r}_j \neq 0. $

    Our strategy is to treat $ \boldsymbol{\Xi}_{out} $ and $ \boldsymbol{\Xi}_{in} $ as the outgoing and incoming incidence matrices of a multi digraph with vertices 'smeared' over subnetworks of flow connections. Thus we have

    Assumption 5. The matrix

    $ \mathsf{A}: = \widehat {\left(\widehat{\boldsymbol{\Xi}_{out}}\right)^T\widehat{\boldsymbol{\Xi}_{in}}} $

    is the adjacency matrix of the line graph of a multi digraph.

    For $ \mathsf{A} $, let $ V^{out}_j $ and $ V^{in}_i $ be groups of row and column indices, respectively, potentially outgoing from (respectively incoming to) a vertex, see Appendix A. We introduce

    $ I : = \{i\in \{1, \ldots, 2m\};\; \boldsymbol{\xi}^{in, r}_i = 0\} $

    and adopt

    Assumption 6. For all $ i\in I $ there exists $ j\in\{1, \ldots, M'\} $ such that

    $ \mathit{\text{supp}}\; \boldsymbol{\xi}^{out, r}_i \subset V^{out}_j. $

    In the next proposition we shall show that $ V^{out}_j $ and $ V^{in}_i $ determine a partition of the row indices into sets that can be used to define vertices. The idea is that if supports of columns of $ \boldsymbol{\Xi}_{out} $ or $ \boldsymbol{\Xi}_{in} $ overlap, the arcs determined by these columns must be incident to the same vertex.

    Proposition 3. If Assumptions 4, 5 and 6 are satisfied, then the sets

    $ {{Vi}1in,VS}, $ (30)

    where

    $ VS={i{1,,2m};suppξout,riVoutM}, $ (31)
    $ Vi=sVoutisuppξout,cs,1in, $ (32)

    form a partition of the row indices of both $ \boldsymbol{\Xi}_{out} $ and $ \boldsymbol{\Xi}_{in} $ such that if for any $ j $, $ \mathrm{supp\;}\boldsymbol{\xi}^{out, c}_j\cap \mathcal{V_i}\neq \emptyset $ for some $ i = 1, \ldots, n, S $, then $ \mathrm{supp\;}\boldsymbol{\xi}^{out, c}_j\subset \mathcal{V_i}, $ and if for any $ k $, $ \mathrm{supp\;}\boldsymbol{\xi}^{in, c}_k\cap \mathcal{V_l}\neq \emptyset $ for some $ l = 1, \ldots, n, S $, then $ \mathrm{supp\;}\boldsymbol{\xi}^{in, c}_k \subset \mathcal{V_l}. $

    Since the proof is quite long, we first present its outline.

    Step 1. Reconstruct a multi digraph $ \boldsymbol{\Gamma} $ from the adjacency matrix of its line graph.

    Step 2. Identify the rows of $ \boldsymbol{\Xi}_{out} $ which are zero in $ \boldsymbol{\Xi}_{in} $ and which correspond to sources and associate them with vertices.

    Step 3. Associate other rows of $ \boldsymbol{\Xi}_{out} $ which are zero in $ \boldsymbol{\Xi}_{out} $ with appropriate vertices.

    Step 4. Associate remaining rows with vertices and construct a possible partition of the row indices.

    Step 5. Check that the constructed partition has the required properties.

    Proof. Step 1. By Assumption 5, $ \mathsf{A} $ is the adjacency matrix of $ L(\boldsymbol{\Gamma}) $ for some multi digraph $ \boldsymbol{\Gamma} $. As explained in Appendix A, we can reconstruct $ \boldsymbol{\Gamma} $ with the transient vertices defined in a unique way and admissible sources and sinks. Let us fix such a construction. Then, we have the sets $ \{V^{in}_j\}_{1\leq j\leq n} $ and $ \{V^{out}_i\}_{1\leq i\leq n} $ of incoming and outgoing arcs determining any transient vertex. Further, we have (possibly) the sets $ V^{in}_{N'} $ and $ V^{out}_{M'} $ that group the arcs incoming to sink(s) and outgoing from source(s), respectively. Since $ \mathsf{A} $ represents all arcs, the same decomposition is valid for $ \boldsymbol{\Xi}_{out} $ and $ \boldsymbol{\Xi}_{in} $, that is, we have subdivisions $ \{V^{out}_i\}_{1\leq i\leq M'} $ and $ \{V^{in}_j\}_{1\leq j\leq N'} $ of the columns of $ \boldsymbol{\Xi}_{out} $ and $ \boldsymbol{\Xi}_{in} $, respectively, and hence the correspondence of the columns with the vertices. Thus, we have to show that (30) is a partition of the rows of $ \boldsymbol{\Xi}_{out} $ and $ \boldsymbol{\Xi}_{in} $ satisfying the conditions of the proposition.

    Let us recall that the entry $ (i, j) $ of $ \mathsf{A} $ is defined by

    $ \mathsf{a}_{ij} = \widehat{\widehat{\boldsymbol{{\xi}}^{out, c}_i}\cdot \widehat{\boldsymbol{{\xi}}^{in, c}_j}} $

    and if a row $ k $ of $ \mathsf{A} $ is zero, that is, it represents a source, then there is a zero row in $ {\boldsymbol{\Xi}_{in}} $. Indeed, since, by Assumption 4, $ \text{supp}\; \boldsymbol{{\xi}}^{out, c}_k \neq \emptyset $, there is a nonzero entry, say, $ \xi^{out}_{lk} $ and thus we must have $ \xi^{in}_{lj} = 0 $ for any $ j $. So, to every zero row in $ \mathsf{A}, $ there corresponds a zero row in $ \boldsymbol{\Xi}_{in} $. However, there may be other zero rows in $ \boldsymbol{\Xi}_{in} $.

    Step 2. To determine the rows in $ \boldsymbol{\Xi}_{out} $, corresponding to sources, we consider Assumption 6. First, we note that any $ j $ of that assumption, if it exists, is determined in a unique way as the sets $ V^{out}_j $ are not overlapping. Next, we observe that if $ \text{supp}\; \boldsymbol{\xi}^{out, r}_i \subset V^{out}_j $ and $ \text{supp}\; \boldsymbol{\xi}^{out, r}_i\cap \text{supp}\; \boldsymbol{\xi}^{out, r}_k \neq \emptyset $, then $ \text{supp}\; \boldsymbol{\xi}^{out, r}_k \subset V^{out}_j $. Indeed, if $ \text{supp}\; \boldsymbol{\xi}^{out, r}_k \subset V^{out}_p $ and $ l\in \text{supp}\; \boldsymbol{\xi}^{out, r}_i\cap \text{supp}\; \boldsymbol{\xi}^{out, r}_k $, then $ l \in V^{out}_j\cap V^{out}_p $ which implies $ V^{out}_j = V^{out}_p. $ Thus we can define the set

    $ \mathcal{V_S} = \{i \in \{1, \ldots, 2m\};\; \text{supp}\; \boldsymbol{\xi}^{out, r}_i \subset V^{out}_{M'}\}. $

    For any $ i\in \mathcal{V_S} $ and $ q\in \text{supp}\; \boldsymbol{\xi}^{out, r}_i, $ we have $ \text{supp}\; \boldsymbol{\xi}^{out, c}_q\subset I, $ as otherwise there would be a nonzero product $ \boldsymbol{\xi}^{out, c}_q\cdot \boldsymbol{\xi}^{in, c}_s $ for some $ s $ as $ \boldsymbol{\xi}^{in, r}_{t} \neq 0 $ for each $ t\notin I $. Then, let there be $ k $ such that $ \text{supp}\; \boldsymbol{\xi}^{out, c}_p \cap \text{supp}\; \boldsymbol{\xi}^{out, c}_j\ni k $ for some $ j \in V^{out}_{M'} $. This means, by Assumption 6, that $ \text{supp}\;\boldsymbol{\xi}^{out, r}_k \subset V^{out}_{M'} $ and hence $ p \in V^{out}_{M'} $. Consider any nonzero element of $ \boldsymbol{\xi}^{out, c}_p, $ say, $ \xi^{out}_{lp}\neq 0 $. By the above argument, $ l \in I $. If $ \text{supp}\; \boldsymbol{\xi}^{out, r}_l \subset V^{out}_{M'} $, then $ l\in \mathcal{V_S} $. If not, $ \text{supp}\; \boldsymbol{\xi}^{out, r}_l\cap V^{out}_j\neq \emptyset $ for some $ j\neq M' $ which contradicts Assumption 6. Thus, $ \mathcal{V_S} $ satisfies the first part of the statement. The second part is void as there is no $ \boldsymbol{\xi}^{in, c}_j $ with $ \text{supp}\;\boldsymbol{\xi}^{in, c}_j\cap \mathcal{V_S}\neq \emptyset. $ Therefore, all indices $ i\in \mathcal{V_S}, $ that is, such that $ \text{supp}\; \boldsymbol{\xi}^{out, r}_i \subset V^{out}_{M'}, $ determine a source as there is no connection to any inflow.

    Step 3. Now, consider the indices $ i\in I\setminus \mathcal{V_S} $. Then, again by Assumption 6, for any $ i\in I\setminus \mathcal{V_S} $ there is a unique $ j\neq M' $ such that $ \text{supp}\; \boldsymbol{\xi}^{out, r}_i \subset V^{out}_{j} $, that is, such an $ i $ belongs to the vertex determined by $ V^{out}_{j} $. This determines a partition of $ I $ corresponding to the vertices (recall that there are no zero rows in $ \boldsymbol{\Xi}_{out} $ and so each row must belong to a vertex).

    Step 4. Next, we associate the remaining rows in $ \boldsymbol{\Xi}_{out} $ and $ \boldsymbol{\Xi}_{in} $ with the vertices. Consider $ V^{out}_i $ and $ V^{in}_j $ for some $ 1\leq i\leq n $ and $ j $ defined by (A.2). The non-zero entries $ \mathsf{a}_{pq} $ of $ \mathsf{A} $, where $ p\in V^{out}_i $ and $ q\in V^{in}_j $, occur whenever $ \text{supp}\; \boldsymbol{{\xi}}^{out, c}_p \cap \text{supp}\; \boldsymbol{{\xi}}^{in, c}_q\neq \emptyset $. Hence, the rows with indices $ k \in \text{supp}\; \boldsymbol{{\xi}}^{out, c}_p \cap \text{supp}\; \boldsymbol{{\xi}}^{in, c}_q $ must belong to a vertex through which the incoming arc $ \boldsymbol{{\varepsilon}}^q $ communicates with the outgoing arc $ \boldsymbol{{\varepsilon}}^p $. Since all nonzero entries in $ \text{supp}\; \boldsymbol{{\xi}}^{out, c}_p $ and $ \text{supp}\; \boldsymbol{{\xi}}^{in, c}_q $, respectively, reflect non-zero outflow along $ \boldsymbol{{\varepsilon}}^p, $ and inflow along $ \boldsymbol{{\varepsilon}}^q $, respectively, $ \text{supp}\; \boldsymbol{{\xi}}^{out, c}_p $ and $ \text{supp}\; \boldsymbol{{\xi}}^{in, c}_q $ must belong to the same vertex. Since the same is true for any indices from $ V^{out}_i $ and $ V^{in}_j $, plausible partitions of row indices of $ \boldsymbol{\Xi}_{out} $ and $ \boldsymbol{\Xi}_{in} $ defining vertices are

    $ \mathcal{V}^{out}_i = \bigcup\limits_{s \in V^{out}_i} \text{supp}\; \boldsymbol{{\xi}}^{out, c}_s, \qquad \mathcal{V}^{in}_j = \bigcup\limits_{q \in V^{in}_j} \text{supp}\; \boldsymbol{{\xi}}^{in, c}_q. $

    We first observe that if $ V^{out}_i $ and $ V^{in}_j $ determine the same transient vertex, then

    $ Vouti{sVouti;ξin,rs=0}=Vinj. $ (33)

    Indeed, let $ p \in \mathcal{V}^{out}_i\setminus \{s \in \mathcal{V}^{out}_i;\; \boldsymbol{\xi}^{in, r}_s = 0\} $. Then, there is $ s \in V^{out}_i $ such that $ p \in \text{supp}\; \boldsymbol{{\xi}}^{out, c}_s $. Since $ p\notin \{s \in \mathcal{V}^{out}_i;\; \boldsymbol{\xi}^{in, r}_s = 0\} $, there is $ q $ such that $ \xi^{in}_{pq}\neq 0 $ and thus $ \widehat{\boldsymbol{\xi}^{out, c}_s} \cdot \widehat{\boldsymbol{\xi}^{in, c}_q} \neq 0 $. Hence, $ q \in V^{in}_j $ and consequently $ p \in \mathcal{V}^{in}_j $. The converse can be proved in the same way by using Assumption 4 since if $ p\in \mathcal{V}^{in}_j $ then, by construction, $ p $ must belong to a support of some $ \boldsymbol{{\xi}}^{in, c}_q $ and thus cannot be in $ \{s \in \mathcal{V}^{out}_i;\; \boldsymbol{\xi}^{in, r}_s = 0\} $. As we see, if $ \mathcal{V}^{out}_i $ contains rows $ \boldsymbol{\xi}^{out, r}_k $ with $ k\in I\setminus \mathcal{V_S} $, then these rows satisfy $ \text{supp}\; \boldsymbol{\xi}^{out, r}_k \subset V^{out}_i $. If we add the indices of such rows to $ \mathcal{V}^{in}_j $ with $ V^{in}_j $ determining the same vertex as $ V^{out}_i $, then such an extended $ \mathcal{V}^{in}_j $ will be equal to $ \mathcal{V}^{out}_i $ and thus we use can use (30) to denote the partition of $ \{1, \ldots, 2m\} $ into $ \mathcal{V_1}^{out}, \ldots, \mathcal{V}^{out}_n, \mathcal{V_S} $.

    Step 5. We easily check that this partition satisfies the conditions of the proposition. We have already checked this for $ \mathcal{V_S} $. So, let $ \text{supp}\;\boldsymbol{\xi}^{out, c}_q \cap \mathcal{V}^{out}_i\neq\emptyset $ for some $ 1\leq i\leq n $, then there is $ s\in V^{out}_i $ such that $ k\in \text{supp}\;\boldsymbol{\xi}^{out, c}_q \cap \text{supp}\;\boldsymbol{\xi}^{out, c}_s $. Clearly, $ k\notin \mathcal{V_S} $ by the construction of $ \mathcal{V}^{out}_i $. If $ k \in I\setminus \mathcal{V_S} $, then $ q \in V^{out}_i $ by Assumption 6 and hence $ \text{supp}\;\boldsymbol{\xi}^{out, c}_q \subset \mathcal{V}^{out}_i $. If $ k\notin I $, then $ \widehat{ \boldsymbol{\xi}^{out, c}_s}\cdot \widehat{\boldsymbol{\xi}^{in, c}_p}\neq 0 $ for some $ p $ but then also $ \widehat{\boldsymbol{\xi}^{in, c}_p}\cdot \widehat{\boldsymbol{\xi}^{out, c}_q}\neq 0 $ and hence $ p \in V^{in}_j $, yielding $ q \in V^{out}_i $ and consequently $ \text{supp}\; \boldsymbol{\xi}^{out, c}_q \subset \mathcal{V}^{out}_i. $ Similarly, if $ \text{supp}\;\boldsymbol{\xi}^{in, c}_p \cap \mathcal{V}^{out}_i\neq\emptyset $, then there is $ k \in \text{supp}\;\boldsymbol{\xi}^{in, c}_p\cap \text{supp}\;\boldsymbol{\xi}^{out, c}_q $ for some $ q \in V^{out}_i $. But then, immediately from the definition, $ \text{supp}\;\boldsymbol{\xi}^{in, c}_p\subset \mathcal{V}^{in}_j\subset \mathcal{V}^{out}_i $ by (33).

    We note that (30) does not contain rows corresponding to sinks and they must be added following the rules described in Appendix A. With such an extension, we consider the multi digraph $ \boldsymbol{\Gamma} $, determined by

    $ {{Vi}1in,VS,VZ},{{Vouti}1in,VoutM,},{{Vinji}1in,,VinN}, $ (34)

    where the association $ i\mapsto j_i $ is defined in (A.2). By construction, if we take the triple $ \mathcal{V_i}, V^{out}_i, V^{in}_{j_i} $, $ 1\leq i\leq n $, it determines a transient vertex, the outgoing arcs given by the indices of columns in $ \boldsymbol{\Xi}_{out} $ and the incoming arcs given by the indices of columns in $ \boldsymbol{\Xi}_{in} $. Similarly, the pair $ \mathcal{V_S}, V^{out}_{M'} $ determines the sources and all outgoing arcs, while the set of incoming arcs is empty. Thus, if we denote by $ \boldsymbol{\Xi}^i_{out} $ and $ \boldsymbol{\Xi}^i_{in} $ the submatrices of $ \boldsymbol{\Xi}_{out} $ and $ \boldsymbol{\Xi}_{in} $ consisting of the rows with indices in $ \mathcal{V_i} $ and columns in $ V^{out}_i $ and $ V^{in}_{j_i} $, respectively, with an obvious modification for $ \mathcal{V_S} $, then (2) decouples into $ n $ (or $ n+1 $) independent systems

    $ Ξiout((υj(0,t))jJ+Vouti,(ϖj(1,t))jJVouti)xxx=Ξiin((υj(1,t))jJ+Vinji,(ϖj(0,t))jJVinji),1in,ΞSout((υj(0,t))jJ+VoutM,(ϖj(1,t))jJVoutM)=0. $ (35)

    This system can be seen as a Kirchhoff system on the multi digraph $ \boldsymbol{\Gamma} $ but we need to collapse $ \boldsymbol{\Gamma} $ to a graph $ \Gamma $ on which (35) can be written as (17). We observe that the question naturally splits into two problems – one is about collapsing the graph, while the other is about grouping the components of $ (\boldsymbol{ \upsilon}, \boldsymbol{ \varpi}) $ into pairs compatible with the parametrization of $ \Gamma $.

    Let $ \mathsf{A} $ be the adjacency matrix of $ L(\boldsymbol{\Gamma}) $, with Assumptions 4 and 6 satisfied. As in Appendix A, we can construct outgoing and incoming incidence matrices $ A^+ $ and $ A^- $ but these are uniquely determined only if there are no sources and sinks. However, we have an additional piece of information about sources.

    If we grouped all sources into one node, as before Proposition 5, then, by Lemma 3.4, the flow connectivity in this source would be given by

    $ \mathsf{C}_{ \mathbf{v}}: = \widehat{\left(\widehat{\boldsymbol{\Xi}_{out}^{S}}\right)^T\widehat{\boldsymbol{\Xi}_{out}^{S}}}. $

    However, such a matrix would not necessarily satisfy Assumption 3. Thus, we separate the arcs into non-communicating groups, each determining a source satisfying Assumption 3. For this, by simultaneous permutations of rows and columns, $ \mathsf{C}_{ \mathbf{v}} $ can be written as

    $ ΞS=diag{ΞSi}1ik, $ (36)

    where $ k $ may equal 1. Since the simultaneous permutation of rows and columns is given as $ P\mathsf{C}_{ \mathbf{v}}P^T, $ where $ P $ is a suitable permutation matrix, [18,p. 140], we see that $ \boldsymbol{\Xi}^S $ is a symmetric matrix, along with $ \mathsf{C}_{ \mathbf{v}} $. By [12,Sections III $ \S $ 1 and Ⅲ $ \S $ 4], $ \boldsymbol{\Xi}^S $ is irreducible if and only if it cannot be transformed by simultaneous row and column permutations to the form (36) with $ k>1 $ (since $ \boldsymbol{\Xi}^S $ is symmetric, all off-diagonal blocks must be zero). Then, $ \boldsymbol{\Xi}^S $ can be reduced to the canonical form, [12,Section Ⅲ,Eq. (68)], in which each $ \boldsymbol{\Xi}^S_i $ is irreducible. If (36) is in the canonical form, then we say that $ \boldsymbol{\Xi}^S $ allows for $ k $ sources, each satisfying Assumption 3. The indices of the columns contributing to the blocks define the $ k $ non-communicating sources $ \mathcal{V}_{S_1}, \ldots, \mathcal{V}_{S_k} $ in $ \boldsymbol{\Gamma} $, which we denote $ V^{out}_{S_1}, \ldots, V^{out}_{S_k} $, $ V^{out}_{M'} = V^{out}_{S_1}\cup \ldots\cup V^{out}_{S_k} $. Finally, define $ \boldsymbol{\xi}^{out, r} = (\xi^{out, r}_{S_i j})_{1\leq i\leq k, 1\leq j\leq 2m} $ by

    $ \xi^{out, r}_{S_i j} = \left\{1ifjVoutSi,0otherwise.\right. $

    For the sinks, it is simpler as there is no constraining information from (2). We have columns with indices in $ V^{in}_{N'} $ corresponding to sinks. These are zero columns in $ \boldsymbol{\Xi}_{in} $ but the columns with these indices in $ \boldsymbol{\Xi}_{out} $ have nonempty supports and thus we can determine from which vertices they are outgoing. Let us denote

    $ VinVi={jVinN;suppξout,cjVi},i=1,,n,S1,,Sk. $ (37)

    For each $ i $ we consider a partition

    $ VinVi=Vini,1Vini,li,i=1,,n,S1,,Sk, $ (38)

    where $ l_i\leq |V^{in}_{\mathcal{V_i}}|, $ into non-overlapping sets $ V^{in}_{i, {l}} $, $ 1\leq l\leq l_i $. Then, we define sinks $ \mathcal{V}_{i, {l}} $ as the heads of the arcs with indices from $ V^{in}_{i, {l}} $; we have $ n_z = l_1+\cdots + l_{S_k} $ sinks. Then, as above, define $ \boldsymbol{\xi}^{in, r} = (\xi^{in, r}_{\{i, {l}\}, q})_{i\in \{1, \ldots, S_k\}, l\in \{1, \ldots, l_i\}, q\in \{1, \ldots, 2m\}} $ by

    $ \xi^{in, r}_{\{i, {l}\}, q} = \left\{1ifqVini,l,0otherwise.\right. $

    Remark 4. We expect $ |V^{in}_{\mathcal{V_i}}|, i = 1, \ldots, n, S_1, \ldots, S_k, $ to be even numbers and (38) to represent a partition of $ V^{in}_{\mathcal{V_i}} $ into pairs so that $ l_i = |V^{in}_{\mathcal{V_i}}|/2. $

    Then, as in Remark 5, the incoming and outgoing incidence matrices are

    $ A^+ = \left(A+0ξin,r\right), \qquad A^- = \left(Aξout,r0\right) $

    which, by a suitable permutation of columns moving the sources and the sinks to the last positions, can be written as

    (39)

    respectively. Both matrices have $ 2m $ columns and $ \mathcal{n} : = n + k + n_z $ rows. Hence, as shown in Remark 5, the adjacency matrix of the full multi digraph $ \boldsymbol{\Gamma} $ is given by

    (40)

    where the dimensions of the blocks in the first row are, respectively, $ n\times n $, $ n\times k $ and $ n\times n_z $, in the second row, $ k\times n $, $ k\times k $ and $ k\times n_z $ and in the last one, $ n_z\times n $, $ n_z\times k $ and $ n_z\times n_z $. Thus, if $ a_{ij} $ is in the block $ (p, q), 1\leq p, q\leq 3, $ then $ a_{ji} $ will be in the block $ (q, p) $.

    Consider a nonzero pair $ (a_{ij}, a_{ji}) $ of entries of $ A(\boldsymbol{\Gamma}) $. If, say, $ a_{ij} = h $, then it means that the $ i $-th row of $ A^+ $ and $ j $-th row of $ A^- $ have entry 1 in the same $ h $ columns, that is, there are exactly $ h $ arcs coming from $ \mathbf{v}_j $ to $ \mathbf{v}_i $. Similarly, if $ a_{ji} = e $, then there are exactly $ e $ arcs coming from $ \mathbf{v}_i $ to $ \mathbf{v}_j $. Conversely, if there are $ h $ arcs from $ \mathbf{v}_j $ to $ \mathbf{v}_i $ and $ e $ arcs from $ \mathbf{v}_i $ to $ \mathbf{v}_j $, then $ (a_{ij}, a_{ji}) = (h, e). $ In particular, $ h+e = 2 $ if and only if there are two arcs between $ \mathbf{v}_j $ and $ \mathbf{v}_i $ running either concurrently or countercurrently. Since the columns of $ A^+ $ and $ A^- $ are indexed in the same way as that of $ \boldsymbol{\Xi}_{out} $ and $ \boldsymbol{\Xi}_{in} $, the pair $ (a_{ij}, a_{ji}) $ determines the rows $ \boldsymbol{a}^{+, r}_i $ and $ \boldsymbol{a}^{+, r}_j $ of $ A^+ $ and thus the indices

    $ (aij{kij1,,kijh},aji{kji1,,kjie})=(aijsuppa+,risuppa,rj,ajisuppa+,rjsuppa,ri) $ (41)

    of columns of $ \boldsymbol{\Xi}_{out} $ (and of $ \boldsymbol{\Xi}_{in} $). Finally, we are ready to formulate the main result of this paper.

    Theorem 3.6. System (2) is graph realizable with generalized Kirchhoff's conditions satisfying Assumptions 1 and 2 for $ \mathbf{v}\in \Upsilon_t $ and Assumption 3 for $ \mathbf{v}\in\Upsilon_s $ if and only if, in addition to Assumptions 4, 5 and 6, there is a partition (38) such that $ A(\boldsymbol{\Gamma}) $ defined by (40) satisfies

    1. for any $ 1\leq i, j\leq \mathcal{n} $, $ a_{ii} = 0 $ and $ (a_{ij}, a_{ji}) $ is in one of the following form:

    $ (2,0),(1,1),(0,2)or(0,0); $ (42)

    2. if $ (a_{ij}, a_{ji}) $ determines the indices $ k $ and $ l $ according to (41), then

    $ if(aij,aji)=(2,0)or(0,2),thenk,lJ+ork,lJandckcl,if(aij,aji)=(1,1),thenkJ+andlJorkJandlJ+.$ (43)

    Proof. Necessity. Let us consider the Kirchhoff system (17). By construction, both matrices $ \boldsymbol{\Psi}^{out} $ and $ \boldsymbol{\Psi}^{in} $ are in block diagonal form with equal row dimensions of the blocks. We consider the problem already transformed to $ \boldsymbol{\Gamma} $. We note that each arc's index must appear twice in $ \boldsymbol{\Psi} $ – once in $ \boldsymbol{\Psi}^{out} $ and once in $ \boldsymbol{\Psi}^{in} $ (if there are sinks, the indices of incoming arcs will correspond to the zero columns). Further, whenever column indices $ k $ and $ l $ appear in the blocks of $ \boldsymbol{\Psi}^{out} $ and $ \boldsymbol{\Psi}^{in} $, respectively, then $ \boldsymbol{{\varepsilon}}^l $ is incoming to, while $ \boldsymbol{{\varepsilon}}^k $ is outgoing from the same vertex (and not any other). Thus, by Assumption 2, the matrix

    $ \tilde A = (\tilde a_{ij})_{1\leq i, j\leq 2m} = \widehat{(\widehat{\boldsymbol{\Psi}^{out}})^T\widehat{\boldsymbol{\Psi}^{in}}} $

    is block diagonal with blocks of the form $ \mathsf{C}_{ \mathbf{v}} = \boldsymbol{1}_{ \mathbf{v}}, $ except for zero rows corresponding to the sources and zero columns corresponding to the sinks. However, in general, the column indices in $ \boldsymbol{\Psi}^{out} $ and $ \boldsymbol{\Psi}^{in} $ do not correspond to the indices of the arcs they represent. Precisely, $ \tilde a_{ij} = 1 $ if and only if there is a vertex for which the arc $ \boldsymbol{{\varepsilon}}^{j'} $ is incoming and $ \boldsymbol{{\varepsilon}}^{i'} $ is outgoing, where $ j' $ and $ i' $ are the indices of the arcs that correspond to the columns $ j $ and $ i $ of $ \boldsymbol{\Psi}^{in} $ and $ \boldsymbol{\Psi}^{out} $, respectively. To address this, we construct $ \boldsymbol{\Xi}_{out} $ and $ \boldsymbol{\Xi}_{in} $ as

    $ \boldsymbol{\Xi}_{out} = \boldsymbol{\Psi}^{out}P, \qquad \boldsymbol{\Xi}_{in} = \boldsymbol{\Psi}^{in}Q, $

    where $ P $ and $ Q $ are permutation matrices, so that in both matrices the column indices $ 1, \ldots, 2m $ correspond to $ \boldsymbol{{\varepsilon}}^1, \ldots, \boldsymbol{{\varepsilon}}^{2m} $. Hence

    $ A = (a_{ij})_{1\leq i, j\leq 2m} : = \widehat{(\widehat{\boldsymbol{\Xi}_{out}})^T\widehat{\boldsymbol{\Xi}_{in}}} = \widehat{(\widehat{\boldsymbol{\Psi}^{out}}P)^T\widehat{\boldsymbol{\Psi}^{in}}Q} = P^T\tilde AQ $

    is a matrix where the indices $ 1, \ldots, 2m $ of both the columns and the rows correspond to $ \boldsymbol{{\varepsilon}}^1, \ldots, \boldsymbol{{\varepsilon}}^{2m} $. Since $ \Gamma $ did not have loops, it is clear that $ a_{ii} = 0 $ for all $ i = 1, \ldots, 2m. $ It is also clear that any two columns (or rows) of $ \tilde A $ are either equal or orthogonal and this property is preserved by permutations of columns and of rows. Hence, by Proposition 4, $ A $ is the adjacency matrix of a line digraph and hence Assumption 5 is satisfied. Since the arcs' connections given by $ A $ and $ \tilde A $ are the same, we see that $ A $ is equal to the adjacency matrix of the line graph of $ \boldsymbol{\Gamma} $. Therefore, the transient vertices determined by $ A $ are the same as in $ \boldsymbol{\Gamma} $ (and hence in $ \Gamma $). On the other hand, as we know, $ A $ does not determine the structure of sources and sinks in $ \boldsymbol{\Gamma} $. The fact that Assumption 4 is satisfied is a consequence of Assumption 1. For Assumption 6, we recall, see Appendix A, that the sets $ V_j^{out} $ group together the indices representing arcs $ \boldsymbol{{\varepsilon}}^k $ outgoing from a single vertex, thus they correspond to the blocks $ \boldsymbol{\Psi}^{out}_{ \mathbf{v}} $ in the matrix $ \boldsymbol{\Psi}^{out} $ and therefore Assumption 6 is satisfied, even for any $ i $. Next, since $ \boldsymbol{\Gamma} $ has been constructed from $ \Gamma $, the structure of the blocks in $ \boldsymbol{\Psi}^{out} $ corresponding to sources ensures that, after permutations, their entries will coincide with $ \boldsymbol{\Xi}^S_{out} $ and thus (36) will hold with the blocks in (36) exactly corresponding to the sources in $ \boldsymbol{\Gamma}, $ on account of Assumption 3. Similarly, in $ \Gamma $ the sinks are determined and thus we have groupings of pairs of the arcs (a partition of the set of indices corresponding to arcs incoming to sinks) coming from transient vertices or sources to sinks and thus the constructions (39) and (40) are completely determined. Then, we observe that whenever we have a source $ \mathbf{v} $, then the arcs outgoing from $ \mathbf{v} $ must be coming in pairs with a single pair coming from $ \mathbf{v} $ to any other possible vertex, meaning that the respective entry in $ A(\boldsymbol{\Gamma}) $ must be either $ (2, 0) $ or $ (0, 2). $ Similar argument holds for the sinks. Since the problem comes from a graph, by construction, the orientation of the flows is consistent with the parametrization.

    Sufficiency. Given (2), we have flows $ (( \upsilon_j)_{j\in J^+}, ( \varpi_j)_{j\in J^-}) $ defined on $ (0, 1). $ Assumptions 4, 5 and 6 ensure, by Proposition 3, the existence of a multi digraph $ \boldsymbol{\Gamma} $ on which (2) can be localized to decoupled systems at vertices and written as (35). Precisely speaking, Assumption 5 associates the indices of incoming components of the solution at a vertex with incoming arcs and similarly for the indices of the outgoing components. Therefore, if an arc $ \boldsymbol{{\varepsilon}}^p $ runs from $ \mathbf{v}_j $ to $ \mathbf{v}_i $, then the flow occurs from $ \mathbf{v}_j $ to $ \mathbf{v}_i $, that is, if $ p\in J^+ $, then the flow on $ \boldsymbol{{\varepsilon}}^p $ is given by $ \upsilon_p $ with $ \upsilon_p(0) $ at $ \mathbf{v}_j $ and $ \upsilon_p(1) $ at $ \mathbf{v}_i $ and analogous statement holds for $ p\in J^- $. In other words, the index $ p $ of the arc $ \boldsymbol{{\varepsilon}}^p $ running from $ \mathbf{v}_j $ to $ \mathbf{v}_i $ determines the orientation of the parametrization: $ 0\mapsto \mathbf{v}_j $ and $ 1\mapsto \mathbf{v}_i $ if $ p\in J^+ $ and $ 0\mapsto \mathbf{v}_i $ and $ 1\mapsto \mathbf{v}_j $ if $ p\in J^-. $

    Now, (42) ensures that there are no loops at vertices and that between any two vertices there are either two arcs or none. If $ a_{ij} $ is an entry in $ \mathsf{A}^+_S\mathsf{S}^T $, $ \mathsf{Z}(\mathsf{A_Z}^-)^T $ or $ \mathsf{Z_S}(\mathsf{S_Z})^T $, then $ a_{ij} = 0 $ or $ a_{ij} = 2 $ and, by the dimensions of the blocks, $ a_{ji}\in\{0, 2\} $ and $ a_{ji} = 0 $, respectively. On the other hand, if $ a_{ij} $ is an entry in $ \mathsf{A}^+_T(\mathsf{A}^-_T)^T $, then it can take any value $ 0, 1 $ or $ 2 $ and the $ a_{ji} $ equals $ 2 $ or $ 0 $, $ 1, \, 0 $, respectively. Thus, double arcs indexed, say, by $ (k, l) $ between vertices could be combined into edges of an undirected graph (with no loops and multiple edges). However, in this way we construct a combinatorial graph which does not take into account that if $ \boldsymbol{{\varepsilon}}^k $ and $ \boldsymbol{{\varepsilon}}^l $ are combined into one edge $ \boldsymbol{e} $ of $ \Gamma $, their orientations must be the same. Thus, if $ (a_{ij}, a_{ji}) = (2, 0) $ determines the pair of indices $ k, l $ according to (43), then both arcs $ \boldsymbol{{\varepsilon}}^k $ and $ \boldsymbol{{\varepsilon}}^l $ of $ \boldsymbol{\Gamma} $ run from $ \mathbf{v}_j $ to $ \mathbf{v}_i $ and the components $ k $ and $ l $ of the solution flow concurrently along $ \boldsymbol{e} $. By assumption, $ k, l \in J^+ $ or $ k, l \in J^- $. In the first case, we associate $ \mathbf{v}_j $ with 0 and $ \mathbf{v}_i $ with 1 and we have $ ( \upsilon_{k}, \upsilon_{l}) $ on $ \boldsymbol{e} $, in agreement with the orientation. Otherwise, we associate $ \mathbf{v}_j $ with 1 and $ \mathbf{v}_i $ with 0 and we have $ ( \varpi_{k}, \varpi_{l}) $ on $ \boldsymbol{e} $. On the other hand, if $ (a_{ij}, a_{ji}) = (1, 1) $ then, by assumption, either $ k\in J^+ $ and $ l\in J^- $ or $ k\in J^- $ and $ l\in J^+ $ and the components $ k $ and $ l $ flow countercurrently. Again, in the first case, $ k\in J^+ $ and $ \boldsymbol{{\varepsilon}}^k $ running from $ \mathbf{v}_j $ to $ \mathbf{v}_i $ requires $ \mathbf{v}_j $ to be associated with 0 and $ \mathbf{v}_i $ with 1, while $ l\in J^- $ and $ \boldsymbol{{\varepsilon}}^l $ running from $ \mathbf{v}_i $ to $ \mathbf{v}_j $ also requires $ \mathbf{v}_j $ to be associated with 0 and $ \mathbf{v}_i $ with 1. Thus, we have $ ( \upsilon_{k}, \varpi_{l}) $ on $ \boldsymbol{e} $. Otherwise, we associate $ \mathbf{v}_j $ with 1 and $ \mathbf{v}_i $ with 0 and we have $ ( \varpi_{k}, \upsilon_{l}) $ on $ \boldsymbol{e} $.

    Finally, the assumption $ c_k\neq c_l $ in the first case of assumption (43) ensures that the resulting system is hyperbolic on each edge.

    Example 3. Let us consider the system

    $ tυj+cjxυj=0,1j4,tϖjcjxϖj=0,5j6, $ (44)

    where $ c_j>0, $ with boundary conditions

    $ (011000100000110000001100000010000001)(υ1(0)υ2(0)υ3(0)υ4(0)ϖ5(1)ϖ6(1))=(000000000000000000000000110001001110)(υ1(1)υ2(1)υ3(1)υ4(1)ϖ5(0)ϖ6(0)). $ (45)

    Thus,

    $ A = \widehat{(\widehat{\boldsymbol{\Xi}_{out}})^T\widehat{\boldsymbol{\Xi}_{in}}} = \left(000000000000000000000000110001001110 \right). $

    Thus, there is a multi digraph $ \boldsymbol{\Gamma} $ for which $ A $ is the adjacency matrix of $ L(\boldsymbol{\Gamma}) $. There is no sink and to determine the structure of the sources, we observe that

    $ \boldsymbol{\Xi}^S_{out} = \left(0110100011000011\right) \quad\text{and so}\quad \boldsymbol{\Xi}^S = \widehat{\left(\widehat{\boldsymbol{\Xi}_{out}^{S}}\right)^T\widehat{\boldsymbol{\Xi}_{out}^{S}}} = \left(1100111001110011\right). $

    This matrix is irreducible and thus we have one source. Therefore

    $ A^+ = \left(110001001110000000\right), \qquad A^- = \left(000010000001111100\right) $

    and consequently

    Further,

    $ suppa+,r1={1,2,6},suppa+,r2={3,4,5},suppa,r1={5},suppa,r2={6},suppa,r3={1,2,3,4}, $

    hence, by (41),

    $ (a_{12}\mapsto \{6\}, a_{21}\mapsto \{5\}), \quad (a_{13}\mapsto \{1, 2\}, a_{31}\mapsto \emptyset), \quad (a_{23}\mapsto \{3, 4\}, a_{32}\mapsto \emptyset). $

    To reconstruct $ \Gamma $, we see that $ \boldsymbol{{\varepsilon}}^5 $ and $ \boldsymbol{{\varepsilon}}^6 $ should be combined into a single edge $ \boldsymbol{e} $. However, since $ J^+ = \{1, 2, 3, 4\}, J^- = \{5, 6\}, $ the flow along $ \boldsymbol{{\varepsilon}}^5 $ runs from 1 to 0 and hence $ \mathbf{v}_1 $ should correspond to 1 in the parametrization, while $ \mathbf{v}_2 $ to 0. On the other hand, $ \boldsymbol{{\varepsilon}}^6 $ runs also from 1 to 0 but from $ \mathbf{v}_2 $ to $ \mathbf{v}_1 $ and hence $ \mathbf{v}_2 $ should correspond to 1, while $ \mathbf{v}_1 $ to 0. This contradiction is in agreement with the violation of assumption (43) as $ (a_{12}, a_{21}) = (1, 1) $ but in the corresponding $ (k, l) = (6, 5) $, both $ k $ and $ l $ belong to $ J^- $.

    Consider a small modification of (44), (45),

    $ tυj+cjxυj=0,1j5,tϖ6c6xϖ6=0, $ (46)

    $ c_j>0 $, with the last two boundary conditions of (45) accordingly changed to

    $ υ5(0)υ1(1)υ2(1)ϖ6(0)=0,ϖ6(1)υ3(1)υ4(1)υ5(1)=0. $ (47)

    The matrices $ \boldsymbol{\Xi}_{out}, \boldsymbol{\Xi}_{in}, \boldsymbol{\Xi}^S, A^+, A^- $ and $ A $ are the same as above and thus the multi digraph $ \boldsymbol{\Gamma} $ is the same as before. However, this time on $ \boldsymbol{{\varepsilon}}^5 $ we have the flow $ \upsilon_5 $, occurring from 0 to 1 and thus $ \boldsymbol{{\varepsilon}}^5 $ and $ \boldsymbol{{\varepsilon}}^6 $ can be combined with a parametrization running from $ 0 $ at $ \mathbf{v}_1 $ to $ 1 $ at $ \mathbf{v}_2 $. Assuming $ c_1>c_2 $, $ c_3>c_4 $, we identify $ u^1_1 = \upsilon_1, u^1_2 = \upsilon_2, u^3_1 = \upsilon_3, u^3_2 = \upsilon_4, u^2_1 = \upsilon_5, u^2_2 = \varpi^6 $ and write (46) as a system of $ 2\times 2 $ hyperbolic systems on a graph $ \Gamma = (\{ \mathbf{v}_1, \mathbf{v}_2, \mathbf{v}_3\}, \{\boldsymbol{e_1}, \boldsymbol{e_2}, \boldsymbol{e_3}\}) $ of the form

    $ tu11+c1xu11=0,tu12+c2xu12=0,tu21+c5xu21=0,tu22c6xu22=0,tu31+c3xu31=0,tu32+c4xu32=0, $ (48)

    with boundary conditions at

    (49)

    Consider a digraph $ G $ (possibly with multiple arcs but with no loops) and its line graph $ L(G) $. For both $ G $ and $ L(G) $ we consider their adjacency matrices $ \mathsf{A}(G) $ and $ \mathsf{A}(L(G)) $. The matrix $ \mathsf{A}(L(G)) $ is always binary, with zeroes on the diagonal. Not any binary matrix is the adjacency matrix of a line graph, see [3,6]. In fact, we have

    Proposition 4. [6,Thm. 2.4.1] A binary matrix $ \mathsf{A} $ is the adjacency matrix of a line digraph of a multi digraph if and only if all diagonal entries are 0 and any two columns (equivalently rows) of $ \mathsf{A} $ are either equal or orthogonal (as vectors).

    For our analysis, it is important to understand the reconstruction of $ G $ from a matrix $ \mathsf{A} = (\mathsf{{a}}_{ij})_{1\leq i, j\leq m} $ satisfying the above conditions. As in (29), we write

    $ \mathsf{A} = (\mathsf{{a}}_{ij})_{1\leq i, j\leq m} = (\boldsymbol{{a}}^c_j)_{1\leq j\leq m} = (\boldsymbol{{a}}^r_i)_{1\leq i\leq m}. $

    If for some $ i_1 $ we have $ \mathsf{a}_{i_1j_1} = \ldots = \mathsf{a}_{i_1j_k} = 1, $ then it means that $ \boldsymbol{e}_{j_1}, \ldots, \boldsymbol{e}_{j_k} $ join $ \boldsymbol{e}_{i_1} $ and thus they must be incident to the same vertex $ \mathbf{v} $ and all $ \boldsymbol{e}_{i_l} $ for which $ \mathsf{a}_{i_lj_1} = 1 $ (and thus all $ \mathsf{a}_{i_lj_p} = 1 $ for $ p = 1, \ldots, k $) are outgoing from $ \mathbf{v} $. We further observe that all zero rows can be identified with source(s). Similarly, zero columns correspond to sinks. If $ \boldsymbol{a}^c_j = \boldsymbol{a}^r_j = 0 $ for some $ j $, then $ \boldsymbol{e_j} $ connects a source to a sink.

    Using the adjacency matrix of a line digraph, we cannot determine how many sources or sinks the original graph could have without additional information. We can lump all potential sources and sinks into one source and one sink, we can have as many sinks and sources as there are zero columns and rows, respectively, or we can subdivide the arcs into some intermediate arrangement. We describe a construction with one source and one sink and indicate its possible variants.

    We introduce $ V^{in}_1 = \{r\in\{1, \ldots, m\};\; \boldsymbol{a}^c_r = \boldsymbol{a}^c_1\} $ and, inductively, $ V^{in}_k = \{r\in\{1, \ldots, m\};\; \boldsymbol{a}^c_r = \boldsymbol{a}^c_{j_k}, j_k = \min\{j;\; j\notin \bigcup_{1\leq p\leq k-1} V^{in}_p\}\} $ and the process terminates at $ N' $ such that $ \bigcup_{1\leq p\leq N'} V^{in}_p = \{1, \ldots, m\} $. In the same way, $ V^{out}_1 = \{l\in\{1, \ldots, m\};\; \boldsymbol{a}^r_l = \boldsymbol{a}^r_1\} $ and $ V^{out}_k = \{l\in\{1, \ldots, m\};\; \boldsymbol{a}^r_l = \boldsymbol{a}^r_{j_k}, j_k = \min\{j;\; j\notin \bigcup_{1\leq p\leq k-1} V^{out}_p\}\} $ and the process terminates at $ M' $ such that $ \bigcup_{1\leq p\leq M'} V^{out}_p = \{1, \ldots, m\} $. In other words, $ \{V^{in}_j\}_{1\leq j\leq N'} $ and $ \{V^{out}_i\}_{1\leq i\leq M'} $ represent the vertices of $ G $ through incoming and outgoing arcs, respectively. If there are any zero rows in $ G $, then we swap the corresponding set $ V^{out}_{j_0} $ with the last set $ V^{out}_{M'} $. In this way, $ V^{out}_{M'} $ represents all arcs outgoing from sources (if they exist). For this construction, we represent them as coming from a single source but other possibilities are allowed, see Remark 5. Similarly, if there are any zero columns, we swap the corresponding set $ V^{in}_{i_0} $ with $ V^{in}_{N'} $, that is, $ V^{in}_{N'} $ represents the arcs incoming to sink(s). Then, we denote

    $ M:={MifVoutM={j;arj0},M1ifVoutM={j;arj=0},N:={NifVinN={j;acj0},N1ifVinN={j;acj=0}. $ (A.1)

    Thus, we see that the number of internal (or transient) vertices, that is, which are neither sources nor sinks is $ n: = M = N $. For such vertices it is important to note that, in general, $ V^{out}_j $ and $ V^{in}_j, $ $ 1\leq j\leq n, $ do not represent the same vertex. To combine $ V^{out}_i $ and $ V^{in}_j $ into the same vertex, we have, for $ 1\leq i, j\leq n $,

    $ vj={Vinj,Vouti},aipjr=1for some/any ipVouti,jrVinj. $ (A.2)

    With this notation, we present a more algorithmic way of reconstructing $ G $ from $ \mathsf{A} $. First, we collapse equal rows of $ \mathsf{A} $ into a single row of $ \mathsf{A}^+ $ and equal columns of $ \mathsf{A} $ into a single column and then take the transpose to get $ \mathsf{A}^- $. Mathematically, let $ \mathbb{I}^+ $ be a set of indices such that $ \mathbb{I}^+\cap V^{out}_i $ consists of exactly one point for each $ 1\leq i\leq M' $ and ordered by the order of $ \{V^{out}_i\} $. Similarly, let $ \mathbb{I}^- $ be a set of indices such that $ \mathbb{I}^-\cap V^{in}_i $ consists of exactly one point for each $ 1\leq i\leq N'. $ We order $ \mathbb{I}^- $ consistently with $ \mathbb{I}^+ $, namely, if $ i_k\in \mathbb{I}^+ $ and $ j_k \in \mathbb{I}^- $ are the $ k $-th indices in $ \mathbb{I}^+ $ and $ \mathbb{I}^- $, respectively, then $ i_k \in V^{out}_k $ and $ j_k\in V^{in}_j $ with $ j $ related to $ k $ by (A.2), that is, $ \{V^{out}_k, V^{in}_j\} $ determines the same vertex $ \mathbf{v}_k $. As mentioned above, possible zero rows correspond to the highest indices. With this, we define

    $ A+=(ari)iI+,A=((acj)jI)T. $ (A.3)

    We see now that each row of $ \mathsf{A}^+ $ corresponds to a vertex and each column of $ \mathsf{A}^+ $ corresponds to an incoming arc. If there are zero rows in $ \mathsf{A} $, there is a zero row at the bottom of $ \mathsf{A}^+ $ showing the presence of a (single) source. The presence of a sink is indicated by zero columns in $ \mathsf{A}^+ $. Similarly, each row of $ \mathsf{A}^- $ corresponds to a vertex with arcs outgoing from it represented by nonzero entries in this row, in columns with indices corresponding to the indices of the arcs. If there are zero columns in $ \mathsf{A} $, they appear as a zero row in $ \mathsf{A}^- $, which represents a (single) sink. Possible sources are visible in $ \mathsf{A}^- $ as zero columns. However, what is important is that even though we lumped all sources and sinks into one single source and a single sink, the zero columns in $ \mathsf{A}^+ $ and $ \mathsf{A}^- $ keep track of the arcs going into the sink or out of the source, respectively. Unless there are no sources and sinks, $ \mathsf{A}^+ $ and $ \mathsf{A}^- $ are not the incoming and outgoing incidence matrices of a graph (for the definition of these, see e.g. [3]). Indeed, $ \mathsf{A}^+ $ does not contain sinks that, clearly, are part of the incoming incidence matrix. Similarly, $ \mathsf{A}^- $ does not include sources. If we keep our requirement that there is only one sink and one source, then we add one row to $ \mathsf{A}^+ $ and one to $ \mathsf{A}^- $ to represent the sink and the source, respectively. We use convention that, if both the sink and the source are present, the source is the last but one row and the sink is the last one. To determine the entries we use the required property of the incidence matrices, that there is exactly one non-zero entry in each column (expressing the fact that each arc has a unique tail and a unique head). Thus, we put 1 in the added rows in any column that was zero in $ \mathsf{A}^+ $ (resp. $ \mathsf{A}^- $). We denote such extended matrices by $ A^+ $ and $ A^- $. It is easy to see that the following result is true.

    Proposition 5. $ A^+ $ and $ A^- $ are, respectively, incoming and outgoing incidence matrices of a multi digraph $ G $ having $ \mathsf{A} $ as the adjacency matrix of $ L(G) $.

    Proof. Since each column of $ A^+ $ and $ A^- $ contains 1 only in one row, we can construct a multi digraph $ G $ from them using $ \mathsf{A}(G) = A^+(A^-)^T $ as its adjacency matrix. Since we allow $ G $ to be a multi digraph, the entries of $ \mathsf{A}(G) $ give the number of arcs joining the vertices. A $ (k, l) $ entry in $ \mathsf{A}(G) $ is given by $ \boldsymbol{a_k}^{+, r}\cdot \boldsymbol{a_l}^{-, r} $ and, by construction, $ \boldsymbol{a_k}^{+, r} $ is a row in $ \mathsf{A} $ belonging to $ \mathbf{v}_k $ and $ \boldsymbol{a_l}^{-, r} $ is a column in $ \mathsf{A} $ corresponding to $ \mathbf{v}_l $. Nonzero entries in $ \boldsymbol{a_k}^{+, r} $ correspond to the arcs incoming to $ \mathbf{v}_k $ and nonzero entries in $ \boldsymbol{a_l}^{-, r} $ correspond to the arcs outgoing from $ \mathbf{v}_l $ so the value of $ \boldsymbol{a_k}^{+, r}\cdot \boldsymbol{a_l}^{-, r} $ is the number of nonzero entries occurring at the same places in both vectors and thus the number of arcs from $ \mathbf{v}_l $ to $ \mathbf{v}_k $.

    The adjacency matrix $ \mathsf{A}(L(G)) $ is determined as $ (A^-)^T A^+ $. The entries of this product are given by $ \boldsymbol{a_k}^{-, c}\cdot \boldsymbol{a_l}^{+, c} $. Since each column has only one nonzero entry (equal to 1), the product will be either 0 or 1. It is 1 if and only if there is $ i $ (exactly one) such that the entry 1 appears as the $ i $-th coordinate of both $ \boldsymbol{a_k}^{-, c} $ and $ \boldsymbol{a_l}^{+, c}. $ Now, by construction, $ a_{ik}^{-} = 1 $ if and only if $ k \in V^{out}_i $ and $ a_{il}^{+} = 1 $ if and only if $ l \in V^{in}_j, $ where the correspondence between $ j $ and $ i $ is determined by (A.2). This is equivalent to $ \mathsf{a}_{kl} = 1 $.

    Remark 5. Assume that $ \mathsf{A} $ has $ k $ zero rows and $ l $ zero columns. We cannot identify the numbers of sinks and sources from $ \mathsf{A} $ without additional information. Above, we lumped all sources and all sinks into one source and one sink, respectively, but sometimes we require more flexibility. As we know, the $ k $ zero rows in $ \mathsf{A} $ become $ k $ zero columns in $ \mathsf{A}^- $ associated with the arcs outgoing from sources. In a similar way, the $ l $ zero columns in $ \mathsf{A} $ stay to be $ l $ zero columns in $ \mathsf{A}^+ $ associated with the arcs incoming to sinks. We can group these arcs in an arbitrary way, with each group corresponding to a source or a sink, respectively. Assume we wish to have $ \bar k $ sources and $ \bar l $ sinks. Then, we build the corresponding matrix $ A^+ $ by adding $ \bar k-1 $ zero rows for the sources to $ \mathsf{A}^+ $ and $ \bar l $ rows corresponding to sinks, which will consist of zeroes everywhere apart from the columns that were zero columns in $ \mathsf{A}^+ $; in these columns we put 1s in such a way that each column contains only one nonzero entry (and zeroes elsewhere). Then, columns having 1 in a particular row will represent the arcs incoming to a given sink. In exactly the same way we extend $ \mathsf{A}^- $, by creating $ \bar l-1 $ zero rows for the sinks and $ \bar k $ rows for the sources. In this way, we construct the following incoming and outgoing incidence matrices, respectively, $ A^+ $ and $ A^- $ that, by a suitable permutation of columns, can be written as

    where $ P $ is the required permutation matrix and, in both cases, the first group of columns have indices corresponding to $ V^{in}_i, 1\leq i\leq n $ (resp. $ V^{out}_j, 1\leq j\leq n $), the second group corresponds to the arcs incoming from the sources to the transient vertices, the third group combines arcs connecting sources and sinks and the last group corresponds to the sinks fed by the transient vertices. We observe that the number of columns in each group in $ A^- $ and $ A^+ $ is the same. Since

    $ \bar A^+(\bar A^-)^T = (A^+P)(A^-P)^T = A^+PP^T(A^-)^T = A^+A^-, $

    as $ P^T = P^{-1} $, see [18,p. 140], for such a digraph $ G $ we have

    (A.4)

    Example 4 Consider the networks $ G_1 $ and $ G_2 $ presented on Fig. 5. We observe that grouping of sources and sinks does not affect the line graph, $ L(G_1) = L(G_2), $ see Fig. 6. To illustrate the discussion above, we have

    $ A=(0000000000000011010000010100000000000101000000000). $ (A.5)
    Figure 5. 

    Multi digraphs $ G_1 $ with 3 sources and two sinks and $ G_2 $ with all sources and all sinks grouped into a single source and a single sink

    .
    Figure 6. 

    The line digraph for both $ G_1 $ and $ G_2 $

    .

    Then,

    $ A+=(110100000101000000000) $ (A.6)

    and we see that there are two transient (internal) vertices $ \mathbf{v}_1 $ and $ \mathbf{v}_2 $ with arcs $ \boldsymbol{e_1}, \boldsymbol{e_2} $ and $ \boldsymbol{e_4} $ incoming to $ \mathbf{v}_1 $ and arcs $ \boldsymbol{e_3} $ and $ \boldsymbol{e_5} $ are incoming to $ \mathbf{v}_2. $ The last row in $ \mathsf{A}^+ $ corresponds to source(s) with outgoing arcs $ \boldsymbol{e_1}, \boldsymbol{e_2}, \boldsymbol{e_5} $ and $ \boldsymbol{e_7} $. We also note that the zero columns in $ \mathsf{A}^+ $ correspond to arcs $ \boldsymbol{e_6} $ and $ \boldsymbol{e_7} $ that are incoming to sinks. To build $ \mathsf{A}^- $, we first collapse the identical columns of $ \mathsf{A} $ and take the transpose. We see from $ \mathsf{A} $ that the first row of the transpose corresponds to the incoming arcs $ \boldsymbol{e_1}, \boldsymbol{e_2} $ and $ \boldsymbol{e_4} $ and thus also to vertex $ \mathbf{v}_1 $ of $ \mathsf{A}^+. $ Hence, there is no need to re-order the rows and so (A.3) gives

    $ A=(001000000010100000000). $ (A.7)

    The last row corresponds to sinks and the zero columns inform us that arcs $ \boldsymbol{e_1}, \boldsymbol{e_2}, \boldsymbol{e_5} $ and $ \boldsymbol{e_7} $ emanate from sources.

    If we want to reconstruct the original graph with one source and one sink, then

    $ A^+ = \left(1101000001010000000000000011\right), \qquad A^- = \left(0010000000101011001010000000\right) $

    and

    $ A^+(A^-)^T = \left(0120101000000110\right), $

    which describes the right multi digraph in Fig. 5. On the other hand, we can consider two sinks (maximum number, as there are two zero columns in $ \mathsf{A} $) and, say, three sources. Then,

    $ A^+ = \left(1101000001010000000000000000000000000000010000010\right), \qquad A^- = \left(0010000000101010000010100000000010000000000000000\right) $

    and

    $ A^+(A^-)^T = \left(0111000100010000000000000000000000000100000100000\right) $

    which describes the left multi digraph in Fig. 5.

    It is easily seen that both digraphs have the same line digraph, shown on Fig. 6, whose adjacency matrix is $ \mathsf{A} $.

    [1] IARC Monographs Polynuclear Aromatic Compounds, Part 1, Chemical, Environmental and Experimental Data. Vol. 32. Lyon, France: International Agency for Research on Cancer (1983).
    [2] Knuckles ME, Inyang F, Ramesh A. (2001) Acute and subchronic oral toxicities of benzo[a]pyrene in F-344 rats. Toxicol Sci 61: 382-388. doi: 10.1093/toxsci/61.2.382
    [3] Inyang F, Ramesh A, Kopsombut P, et al. (2003) Disruption of testicular steroidogenesis and epididymal function by inhaled benzo[a]pyrene. Reprod Toxicol 17: 527-537. doi: 10.1016/S0890-6238(03)00071-6
    [4] Perera FP, Santella RM, Brenner D, et al. (1989) Application of biological markers to the study of lung cancer causation and prevention. IARC Sci 89: 451-459.
    [5] Chao MR, Wang CJ, Wu MT, et al. (2008) Repeated measurements of urinary methylated/oxidative DNA lesions, acute toxicity, and mutagenicity in coke-oven workers. Cancer Epidemiol Biomarkers Prev 17: 3381-3389. doi: 10.1158/1055-9965.EPI-08-0721
    [6] Liu AL, Lu WQ, Wang ZZ (2006) Elevated levels of urinary 8-hydroxy-2-deoxyguanosine, lymphocytic micronuclei, and serum glutathione S-transferase in workers exposed to coke oven emissions. Environ Health Perspect 114: 673-677.
    [7] Park JH, Troxel AB, Harvey RG, et al. (2006) Polycyclic aromatic hydrocarbon (PAH) o-quinones produced by the aldo-keto-reductases (AKRs) generate abasic sites, oxidized pyrimidines, and 8-oxo-dGuo via reactive oxygen species. Chem Res Toxicol 19: 719-728. doi: 10.1021/tx0600245
    [8] Cavalieri EL, Rogan EG (1995) Central role of radical cations in metabolic activation of polycyclic aromatic hydrocarbons. Xenobiotica 25: 677-678. doi: 10.3109/00498259509061885
    [9] Cooper CS, Grover PL, Sim P (1983) The metabolism and activation of benzo[a]pyrene. Prog Drug Metab 7: 297-396.
    [10] Grollman AP, Moriya M (1993) Mutagenesis by 8-oxoguanine: an enemy within. Trends Genet 9: 246-249. doi: 10.1016/0168-9525(93)90089-Z
    [11] Nuernberg AM, Boyce PD, Cavallari JM, et al. (2008) Urinary 8-isoprostane and 8-OHdG concentrations in boilermakers with welding exposure. J Occup Environ Med 50: 182-189. doi: 10.1097/JOM.0b013e31815cf6cc
    [12] Singh R, Sram RJ, Binkova B, et al. (2007) The relationship between biomarkers of oxidative DNA damage, polycyclic aromatic hydrocarbon DNA adducts, antioxidant status and genetic susceptibility following exposure to environmental air pollution in humans. Mutat Res 620: 83-92. doi: 10.1016/j.mrfmmm.2007.02.025
    [13] Cheng KC, Cahhill DS, Kasai H, et al. (1992) 8-hydroxyguanine, an abundant form of oxidative DNA damage, causes G-T and A-C substitutions. J Biol Chem 267: 166-172.
    [14] Singh R, Sram RJ, Binkova B, et al. (2007) The relationship between biomarkers of oxidative DNA damage, polycyclic aromatic hydrocarbon DNA adducts antioxidant status and genetic susceptibility following exposure to environmental air pollution in humans. Mutat Res 620: 83-92. doi: 10.1016/j.mrfmmm.2007.02.025
    [15] Guz J, Gachowski F M, Rozalski R, et al. (2013) Comparison of Oxidative Stress/DNA Damage in Semen and Blood of Fertile and Infertile Men. PLoS One 8: e68490. doi: 10.1371/journal.pone.0068490
    [16] Evans MD, Olinski R, Loft S, et al. (2010) Toward consensus in the analysis of urinary 8-oxo-7, 8-dihydro-29-deoxyguanosine as a noninvasive biomarker of oxidative stress. FASEB J 24: 1249-1260. doi: 10.1096/fj.09-147124
    [17] Zenzes MT, Bielechi R, Reed TE (1999). Detection of benzo(a)pyrene diol epoxides-DNA adducts in sperm of men exposed to cigarette smoke. Fertil Steril 72: 330-335. doi: 10.1016/S0015-0282(99)00230-7
    [18] Binková B, Lewtas J, Misková I, et al. (1995) DNA adducts and personal air monitoring of carcinogenic polycyclic aromatic hydrocarbons in an environmentally exposed population. Carcinogenesis 16: 1037-1046. doi: 10.1093/carcin/16.5.1037
    [19] Revel A, Raanani H, Younglai E, et al. (2001) Resveratrol, a natural aryl hydrocarbon receptor antagonist, protects sperm from DNA damage and apoptosis caused by benzo(a)pyrene. Reprod Toxicol 15: 479-486. doi: 10.1016/S0890-6238(01)00149-6
    [20] Park JH, Troxel AB, Harvey RG, et al. (2006) Polycyclic aromatic hydrocarbon (PAH) o-quinones produced by the aldo-keto-reductases (AKRs) generate abasic sites, oxidized pyrimidines, and 8-oxo-dGuo via reactive oxygen species. Chem Res Toxicol 19: 719-728. doi: 10.1021/tx0600245
    [21] Cheng KC, Cahhill DS, Kasai H, et al. (1993) 8-hydroxyguanine, an abundant form of oxidative DNA damage, causes G-T and A-C substitutions. J Biol Chem 267: 166-172.
    [22] Hu CW, Chao MR, Sie CH (2010) Urinary analysis of 8-oxo-7,8-dihydroguanine and 8-oxo-7,8-dihydro-2'-deoxyguanosine by isotope-dilution LC-MS/MS with automated solid-phase extraction: Study of 8-oxo-7,8-dihydroguanine stability. Free Radic Biol Med 48: 89-97. doi: 10.1016/j.freeradbiomed.2009.10.029
    [23] Sawyer DE, Mercer BG, Wiklendt AM, et al. (2003) Quantitative analysis of gene-specific DNA damage in human spermatozoa. Mutat Res 529: 21-34. doi: 10.1016/S0027-5107(03)00101-5
    [24] Paracchini V, Chang SS, et al. (2005) GSTM1 [corrected] deletion modifies the levels of polycyclic aromatic hydrocarbon-DNA adducts in human sperm. Mutat Res 586: 97-101. doi: 10.1016/j.mrgentox.2005.06.008
    [25] Gaspari L, Chang SS, Santella RM, et al. (2003) Polycyclic aromatic hydrocarbon-DNA adducts in human sperm as a marker of DNA damage and infertility. Mutat Res 535: 155-160. doi: 10.1016/S1383-5718(02)00297-8
    [26] Spano M, Bonde JP, Hjollund HI (2000) Sperm chromatin damage impairs human fertility. Fertil Steril 73: 43-50. doi: 10.1016/S0015-0282(99)00462-8
    [27] Klungland A, Bjelland S (2007) Oxidative damage to purines in DNA: role of mammalian Ogg1. DNA Repair (Amst) 6: 481-488. doi: 10.1016/j.dnarep.2006.10.012
    [28] Gedik CM, Collins A (2005) ESCODD (European Standards Committee on Oxidative DNA Damage) Establishing the background level of base oxidation in humanlymphocyte DNA: results of an interlaboratory validation study. FASEB J 19: 82-84.
    [29] Guetens G, De Boech, G, Highley M, et al. (2002) Oxidative DNA damage: biological significance and methods of analysis. Crit Rev Clin Lab Sci 39:3310-457.
    [30] Holmberg I, Stal P, Hamberg M (1999) Quantitative determination of 8-hydroxy-2'-deoxyguanosine in human urine by isotope dilution mass spectrometry normal levels in hemochromatosis. Free Radical Bio Med 26: 126-129.
    [31] Dizdaroglu M, Jaruga P, Birincioglu M, et al. (2002) Free radical-induced damage to DNA mechanisms and measurement. Free Radical Biol Med 32: 1102-1115. doi: 10.1016/S0891-5849(02)00826-2
    [32] Cooke MS, Singh R, Hall GK, et al. (2006) Evaluation of enzyme-linked immunosorbent assay and liquid chromatography-tandem mass spectrometry methodology for the analysis of 8-oxo-7, 8-dihydro-2'-deoxyguanosine in saliva and urine. Free Radic Biol Med 41: 1829-1836. doi: 10.1016/j.freeradbiomed.2006.09.009
    [33] Ravanat JL, Duretz B, Guiller A, et al. (1998) Isotope dilution high-performance liquid chromatography-electrospray tandem mass spectrometry assay for the measurement of 8-oxo-7,8-dihydro-2'-deoxy-guanosine in biological samples. J Chromatogr B 715: 349-356. doi: 10.1016/S0378-4347(98)00259-X
    [34] World Health Organization. (2010) WHO laboratory manual for the examination and processing of human semen. WHO, 5.
    [35] Evenson DP, Wixon R (2005) Environmental toxicants cause sperm DNA fragmentation as detected by the sperm chromatin structure assay. Toxicol Appl Pharmacol 207: S532-537. doi: 10.1016/j.taap.2005.03.021
    [36] Evenson DP, Larson KL, Jost LK (2002) Sperm chromatin structure assay: its clinical use for detecting sperm DNA fragmentation in male infertility and comparisons with other techniques. J Androl 23: 25-43.
    [37] Gedik CM, Collins A. ESCODD (European Standards Committee on Oxidative DNA Damage) (2005) Establishing the background level of base oxidation in human lymphocyte DNA: results of an interlaboratory validation study. FASEB J 19: 82-84.
    [38] Chao MR, Yen CC, Hu CW (2008) Prevention of artifact oxidation in determination of cellular 8-oxo-7,8-dihydro-2'-deoxyguanosine by isotope-dilution LC-MS/MS with automated solid-phase extraction. Free Radic Biol Med 44: 464-473. doi: 10.1016/j.freeradbiomed.2007.10.003
    [39] Rubes J, Selevan SG, Evenson D P, et al. (2005). Episodic air pollution is associated with increased DNA fragmentation in human sperm without other changes in semen quality. Hum Reprod 20: 2776-2783. doi: 10.1093/humrep/dei122
    [40] Zribi N, Chakroun NF, Elleuch H, et al. (2011) Sperm DNA fragmentation and oxidation are independent of malondialdheyde. Reprod Biol Endocrinol 9: 47. doi: 10.1186/1477-7827-9-47
    [41] Singh S, Sharma S, Jain M, et al. (2011) Importance of papanicolaou staining for sperm morphologic analysis: comparison with an automated sperm quality analyzer. Am J Clin Pathol 136: 247-251. doi: 10.1309/AJCPCLCSPP24QPHR
    [42] Sheridan J, Wang LM, Tosetto M, et al. (2009) Nuclear oxidative damage correlates with poor survival in colorectal cancer. British J Cancer 100: 381-388. doi: 10.1038/sj.bjc.6604821
    [43] Twigg J, Irvine DS, Houston P, et al. (1998) Iatrogenic DNA damage induced in human psermatozoa during sperm preparation: protective significance of seminal plasma. Mol Hum Reprod 4: 439-445. doi: 10.1093/molehr/4.5.439
    [44] Kao SH, Chao HT, Chen HW, et al. (2008) Increase of oxidative stress in human sperm with lower motility. Fertil Steril 89: 1183-1190. doi: 10.1016/j.fertnstert.2007.05.029
    [45] Aitken RJ, Baker MA, De Iuliis GN, et al. (2010) New insights into sperm physiology and pathology. Handb Exp Pharmacol 198: 99-115. doi: 10.1007/978-3-642-02062-9_7
    [46] De Vries M, Ramos L, Housein Z, et al. (2012) Chromatin remodelling initiation during human spermiogenesis. Biol Open 1: 446-457. doi: 10.1242/bio.2012844
    [47] Köberle B, Masters JRW, Hartley JA, et al. (1999) Defective repair of cisplatin-induced DNA damage caused by reduced XPA protein in testicular germ cell tumors. Current Bio 9: 273-278. doi: 10.1016/S0960-9822(99)80118-3
    [48] Smith R, Vantman D, Ponce J, et al. (1996) Total antioxidant capacity of human seminal plasma. Hum Reprod 11:1655-60. doi: 10.1093/oxfordjournals.humrep.a019465
    [49] Thiele JJ, Friesleben HJ, Fuchs J, et al. (1995) Ascorbic acid and urate in human seminal plasma: determination and interrelationships with chemiluminescence in washed semen. Hum Reprod 10: 110-115. doi: 10.1093/humrep/10.1.110
  • This article has been cited by:

    1. Jacek Banasiak, Adam Błoch, Telegraph systems on networks and port-Hamiltonians. Ⅲ. Explicit representation and long-term behaviour, 2022, 11, 2163-2480, 2165, 10.3934/eect.2022016
    2. Adam Bobrowski, Elżbieta Ratajczyk, Pairs of complementary transmission conditions for Brownian motion, 2024, 388, 0025-5831, 4317, 10.1007/s00208-023-02613-x
  • Reader Comments
  • © 2015 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(8736) PDF downloads(948) Cited by(0)

Figures and Tables

Tables(3)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog