We present an approach to the determination of causal connectivities and part of the kinetics of complex reaction systems. Our approach is based on analytical and computational methods for studying the effects of a pulse change of concentration of a chemical species in a reaction network, either at equilibrium or in a nonequilibrium stationary state. Such disturbances generally propagate through a few species, depending on the values of the kinetic coefficients, before being broadened and dissipated. This short range gives a local probe of the kinetics and connectivity of the reaction network. The range of propagation also indicates species to perturb in further experiments. From piecing together these local connectivities, the global structure of the network can be constructed. The experimental design allows deduction of both reaction orders and rate constants in many cases. An example of the usefulness of the approach is illustrated on a model of a part of glycolysis.